# Power asymmetry in CMB polarization maps from PLANCK : a local variance analysis

## Abstract

A persistent signal of power asymmetry on opposite hemispheres of CMB sky was seen in full-sky temperature measurements made so far. This asymmetry was seen in microwave sky from WMAP as well as PLANCK satellites, and calls for attention the larger question of statistical isotropy, one of the foundational principles of modern cosmology. In this work we present an analysis of polarized CMB maps from PLANCK 2015 full mission data. We apply the local variance estimator on low resolution mode maps from PLANCK 2015 polarization Commander solution. We find a significant hemispherical power asymmetry in polarization data on large angular scales, at the level of depending on the galactic mask, and the circular disc radius used for computing local variance maps. However the direction is found to be pointing broadly towards CMB kinetic dipole direction. Precise measurements of CMB polarization in future will shed light on this apparent discrepancy in the anisotropy axis seen in temperature and polarized CMB sky, and likely influence of systematics on our findings.

Korea Institute for Advanced Study (KIAS), Quantum Universe Center (QUC),

85 Hoegi-ro, Dongdaemun-gu, Seoul - 02455, Republic of Korea

Korea Astronomy and Space Science Institute (KASI),

Yuseong-gu, 776 Daedeok daero, Daejeon 34055, Republic of Korea

University of Science and Technology (UST),

Yuseong-gu 217 Gajeong-ro, Daejeon 34113, Republic of Korea

## 1 Introduction

Hemispherical power asymmetry (HPA), initially observed in WMAP first year temperature data [1] has come to be one of the outstanding anomalies that indicated violation of statistical isotropy on large angular scales of CMB sky. This anisotropic signal persisted between completely different full-sky missions viz., WMAP and PLANCK probes, that involved different analysis procedures, systematics and science teams [2, 3, 4]. HPA attracted significant attention from the cosmology community, which was estimated using a variety of methods, see for example Ref. [5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25]. Many explanations were also put forth to explain this anomaly with no final conclusion.

This hemispherical power asymmetry has come to be modelled as a dipole modulation of otherwise statistically isotropic CMB sky as [26, 27],

(1) |

where and are the modulated/observed and statistically isotropic CMB anisotropy fields respectively in the direction . and are the amplitude and direction of modulation. Here the amplitude can be scale dependent [3, 4, 18, 20, 24, 25]. This phenomenological modelling of HPA as dipole modulation of temperature field has important implications for the statistics of observed CMB signal. It induces off-diagonal correlations in the two-point function that couples a multipole to . These correlations were exploited to estimate the anisotropic signal itself [3, 4, 7, 14, 18, 24, 23, 25].

In this paper, we probe the presence of hemispherical power asymmetry in the polarized CMB signal from PLANCK full mission data. The PLANCK science team didn’t report any analysis of PLANCK 2015 CMB polarization maps in the context of isotropy studies. As the PLANCK papers inform [28], this is because of the uncertainty in the level of recovered polarization signal on large angular scales, that warrants such analysis. Accordingly the low multipoles of polarized CMB signal in the various component separated CMB polarization maps have been excluded in the full mission PLANCK data that is made public. Further the complimentary simulations to polarization data that are also made public, have a significant level of mismatch in the noise level compared to data. However in this paper, we work in the lower to intermediate range of multipoles, and suitably modifying the noise simulations, that are already made available through the second public release to compute significances.

## 2 CMB polarization maps

The CMB temperature anisotropies are conventionally expanded in terms of spherical harmonics as

(2) |

where are temperature anisotropies after subtracting the leading monopole and dipole components in the sky direction , are the coefficients of expansion, and are a suitable basis to expand the temperature anisotropies on a sphere.

Unlike the intensity field, CMB polarization stokes parameters are not scalars under rotation and are frame dependent. The combination transforms as , under rotation by an angle of the co-ordinate plane perpendicular to the direction in which are defined. Hence transform as a spin field under coordinate transformations. Thus it is expanded in terms of spin-weighted spherical harmonics, , as

(3) |

where are the spin-2 spherical harmonics, and are the coefficients of expansion in that basis [29, 30, 31, 32]. are the usual spherical harmonics that are used to decompose the CMB temperature field which is a spin-0 field.

Further, a linear combination of these spin-2 spherical harmonic coefficients, , can be defined as [30, 32]

(4) | |||||

(5) |

where and are respectively called electric and magnetic modes of polarization that are invariant under rotations. Similar to an electric and magnetic field, and respectively have even and odd parity properties under inversion ie., for , and . Thus one can obtain two spin-0 fields as

(6) |

As mentioned above, owing to their parity properties under inversion, is a scalar field while is a pseudo-scalar field. Physically these modes originate from distinct sources. Whereas, scalar perturbations give rise to purely modes, only tensor perturbations can generate modes (on large angular scales). Lensing due to large scale structure can induce modes from modes at small angular scales.

In this paper, we test for hemispherical power asymmetry in CMB mode maps derived from PLANCK full mission data.

## 3 Methodology and Data sets

### 3.1 Local variance estimator

We use the recently proposed local variance estimator (LVE) that was applied to probe hemispherical power asymmetry in temperature data [20]. This statistical approach is based on the Crossing statistic that has been proposed to address some different problems in astronomy [33, 34, 35, 36].

Hereafter we use LV or LVE to refer to local variance analysis or maps thus estimated. In this section, we briefly describe the LV method to apply on polarized CMB maps.

One can define a variance map of the observed CMB temperature sky as,

(7) |

where is the variance computed locally from all the pixels , falling inside a circular disc of chosen radius defined in the sky direction , and is the mean temperature corresponding to the same circular region.

A variance map of dipole modulated CMB sky following Eq. [1] is given by

(8) |

upto first order in the isotropy violating dipole field, where and are maps of locally computed variances of modulated and unmodulated CMB anisotropies respectively. Further, a normalized variance map is constructed as

(9) |

where is the expected mean field bias in variance estimated from CMB simulations (including inhomogeneous noise), that is subtracted from the local variance map of the observed sky. From this normalized variance map we compute the dipole amplitude (and direction) to check for consistency with random expectations for any given realization. Any significant dipole amplitude thus found in data would indicate a significant hemispherical power asymmetry, hemisphere defined by the corresponding dipole direction, and hence a breakdown of statistical isotropy. We should note that one can define different metrics to test the statistical isotropy using the normalized variance map but in this work we focus on the dipole amplitude and its direction.

In practice, one implements the local variance estimator as follows. Any cleaning procedure to clean the raw satellite data, is limited by the level of foreground emission in a given sky region. Correspondingly, owing to high emission levels in the galactic plane, the recovered CMB signal is less reliable in those regions. Bright point and extended sources of extragalactic origin also lead to the same situation. Thus these regions are masked in any CMB analysis to avoid biases due to residual foreground contamination. Hence, in computing the local variances at different sky locations, the circular disc mask defined locally is taken in union with the galactic mask to exclude potentially contaminated regions of the recovered CMB sky. Further to make the analysis robust, we can define a criterion for using these locations only if a chosen fraction of the sky pixels, say for example 50% or 80%, in comparison to the number of pixels available using the full disc mask, survive in the effective disc mask (ie., the union of galactic mask and the locally defined circular disc mask of a chosen radius). Thus, the fraction of the sky over which local variances are computed is less than the sky fraction of the galactic mask used to omit regions with potential residual contamination. Operationally, we set these invalid pixels to the HEALPix bad value . And a dipole is fit to the (normalized) local variance map thus obtained using only the valid pixels.

To analyze the polarization maps, we apply the LV estimator to the mode map (defined by Eq. [6]). The variance is computed as , where the angular brackets denote the average of the squared mode field computed using all pixels, , falling inside an effective circular disc of chosen radius, , defined at on the CMB sky.

### 3.2 Data sets used

#### mode CMB maps

We do our local variance analysis on a CMB mode map. To this extent
we use the full-sky Stokes and CMB maps from PLANCK full mission data, estimated
using the Commander algorithm. These maps are made available through PLANCK’s
second public release^{3}

The Commander estimated maps are available at HEALPix with smoothing given by a Gaussian beam of (arcmin). Since the polarized CMB maps are dominated by noise are high, in this paper we work with a low resolution mode map at with a beam resolution of (arcmin) Gaussian beam. First the HEALPix facilities - anafast and alteralm - are applied in succession to obtain ’s from Stokes maps, at the desired lower beam/pixel resolution.

Now, as mentioned earlier, various component separated maps have their low modes removed in the maps that are made publicly available. This is so because of the unsatisfactory level of systematic artefacts still present on large angular scales of PLANCK polarization maps [28]. A window function in multipole space defined as

(10) |

where and , is applied to these publicly released CMB maps. Further as discussed in Appendix A, to reduce the effect of noise dominated modes in the CMB mode map thus obtained on our analysis, we additionally apply a similar cosine filter at high as

(11) |

where and .

The hemispherical power asymmetry, as found in CMB temperature data is present significantly on the large angular scales viz., [3, 4, 14, 18, 24, 23, 25]. In order to test for similar behaviour in the polarized CMB map, we do additional high pass filtering to get two mode CMB maps using the window function defined by Eq. [10] with / and / at low. The various space windows used to produce the filtered mode CMB maps are shown in Fig. [1].

These filtered mode ’s are then converted to an mode map as

(12) |

using HEALPix facility synfast. Thus, in the current study, we use three different high pass filtered CMB mode maps as derived from PLANCK provided Stokes maps estimated using Commander algorithm. They have a beam resolution given by a Gaussian beam of (arcmin) and at , comprising different angular modes.

#### Masks

Common galactic masks for polarization analysis were provided by the PLANCK team through
their second data release^{4}

Since we work at low resolution of HEALPix , the UPB77 mask is
downgraded as follows.
The high resolution mask is available
at HEALPix . We first invert this mask and fill all island like
regions inside using the HEALPix facility `process_mask`

.
To be more conservative, these island like regions in the galactic region are excluded to begin with.
This inverted/filled high resolution mask is downgraded
to . Then a Gaussian beam of (arcmin) is applied, which is
same as that of the CMB map we analyse. Finally a cut-off of 0.05 is applied on this inverted,
filled, downgraded, and smoothed mask where those pixels with pixel-value less than this cut-off
are set to ‘1’, and ‘0’ otherwise.
We do so, so as to be able
to handle leakage due to smearing of extended/point source contamination due to smoothing.
The mask thus obtained has a sky fraction of .
This mask is referred to as mask074, and is shown in red in Fig. [2].

Now the confidence masks for individual component separated CMB polarization maps have an available sky fraction of , , and respectively for Commander, NILC, SEVEM, and SMICA derived solutions. Their mask contours are different, and hence they excise different regions of the sky. So, we further devise two more masks with larger sky fraction for use in this analysis as follows. Here we adopt the procedure of Tegmark et al. (2003) [37], used to define sky partitions for iterative cleaning of the CMB sky. To create these additional masks, we use the raw polarization maps of PLANCK 030, 070 and 353 GHz channels. These maps with different beam/pixel smoothing and map resolution are first downgraded to have a common beam, pixel and map resolution. The raw maps corresponding to each of these three frequencies are deconvolved with their respective circular beam transfer functions, and simultaneously smoothed with a Gaussian beam of (degrees), so as to have largely contiguous regions in the masks being created, while repixelizing them at a lower resolution of . Among the PLANCK frequency channels, the 070 GHz band has visibly lowest foreground levels in the raw satellite data. The raw map corresponding to this frequency band is subtracted from the lowest and highest frequency raw maps from PLANCK that are sensitive to polarization viz., the 030 and 353 GHz channels. These two difference maps, 030-070 and 353-070, in basis containing only foregrounds, are dominated by synchrotron and thermal dust respectively. The individual difference maps in are further combined to form maps, where . These two maps, are then used to create what is called a junk map which retains only the pixel-value that is largest between the two maps at any sky location. This junk map taken as representative of the maximum foreground signal in any sky direction due to any polarized non-CMB emission is thresholded at different values to obtain galactic masks with different sky fractions. Thus a suitable cut-off is applied to this junk map to obtain two more galactic masks that have a sky fraction of and , respectively.

In order to account for the point source contamination, we extend the point source masks
in polarization as provided by PLANCK team at different frequencies^{5}

Thus we have three masks labelled mask074, mask080, and mask085 which have respectively , , and unmasked sky fractions. The LVE maps are obtained using these three foreground exclusion masks. We additionally created two more masks in basis, and the whole LV analysis is repeated using these masks as test of robustness of our analysis using masks.

#### Simulations

Here we describe the generation of simulations used to supplement the clean polarized CMB maps from satellite observations for computing significances.

PLANCK team provided a set of 1000 CMB maps at each frequency band in which it made observations of polarized microwave sky, that have same instrument properties such as beam smoothing similar to the observed maps. Also provided are a set of 1000 inhomogeneous noise maps corresponding to each of the seven frequency bands (030 - 353 GHz). These are called Full Focal Plane (FFP) simulations. However the PLANCK collaboration didn’t make available, a similar set of processed CMB maps with noise, that complement polarized CMB sky recovered from data using various component separation methods. Hence we process the high resolution CMB and noise FFP simulations provided, to suitably obtain mock maps that augment the observed mode CMB maps used in the present study.

The frequency specific Stokes CMB and noise maps are available at a HEALPix
resolution of or , and have different beam smoothing, depending on the
frequency channel. Hence all these maps are brought to a lower resolution of ,
and to have beam resolution given by a Gaussian beam with (arcmin) following the procedure
outlined in Sec. [3.2.1] (ie., applying HEALPix facilities anafast
and alteralm in sequence).
In this process, we used the polarization circular beam transfer functions () for deconvolution, that are also provided with PLANCK’s second public release^{6}

The polarization (inhomogeneous) noise maps are also
subjected to the same procedure to obtain low resolution mode noise maps. Now, the
smoothed/downgraded mode noise maps are combined to obtain noise r.m.s
maps corresponding to each channel at .
These noise r.m.s maps are then used to obtain a set of 1000 mock maps
following inverse noise variance addition of mode CMB and noise maps corresponding to different
frequency channels with same simulation/seed number^{7}

However the effective noise maps thus obtained are not yet useful due to an underestimation of noise in FFP simulations compared to data. As described in Appendix A, we do a rescaling of the effective noise maps by a constant multiplicative factor . Thus a set of 1000 mock CMB mode maps with noise are generated complementing the observational data. These co-added low resolution mode CMB maps with (scaled) noise are further filtered using the high and low pass cosine filters defined in Eq. [10] and [11]. Thus from each co-added realistic mode CMB realization, we obtain three filtered maps similar to the observed data maps described in Sec. [3.2.1].

## 4 Results

In the present analysis, we use a CMB mode map derived from PLANCK 2015 Stokes CMB maps, estimated using Commander algorithm. Different window functions in multipole space are applied to this mode CMB map to get three filtered maps, as described in Sec. [3.2.1]. In order to probe any scale dependent behaviour of the amplitude of hemispherical power asymmetry, if present, that is modelled as dipole modulation of otherwise statistically isotropic CMB anisotropies, we generated these three filtered maps. All the three maps have a low-pass cosine filter applied with following Eq. [11]. A high-pass cosine filter defined by Eq. [10] with , and are applied to obtain the said three filtered mode CMB maps. Note that the filter is already applied by PLANCK team to the publicly released polarization CMB maps. We only generated the later two filtered CMB maps. As a consistency check of our results, we also performed the LV analysis on large angular scales using foreground masks derived in space directly.

### 4.1 Power asymmetry on large angular scales

Here we present results from analyzing the CMB maps containing large angular scales ie., using the filtered mode CMB map comprising the multipoles . We used four different circular disc of radii , , and (arcmin), to compute local variances. The largest angular mode that is fully available in the polarized CMB map from the second data release is , and corresponds to an angular size of (degrees). So we limit to using a maximum disc radius of (degrees) ie., a maximum circular disc of diameter (degrees). The circular discs of different radii are defined at the pixel centers of a HEALPix grid of , to compute the local variances. An map has pixels of size (arcmin) assuming them to be a square of side . Later on we also use two more filtered mode CMB maps, and given the angular scales they contain, we compute all the local variance maps using different disc radii at the same HEAPix resolution of . In practice the local variances are computed from effective circular discs meaning the circular discs defined locally are taken in conjugation with a galactic mask to avoid potential foreground bias on our analysis. We compute local variances only if the effective discs have atleast of the pixels in comparison to those in full disc so that our statistics are robust. This results in local variance maps with less sky fraction that the galactic mask itself used to exclude regions with potential galactic contamination in the data map. The three masks shown in Fig. [2] are used in the present study.

The filtered Commander 2015 mode CMB map comprising the large angular scales as available, is used to compute four local variance maps at , following Eq. [7], corresponding to the choice of four circular disc of radii , , and (arcmin) described earlier. These LV maps are further corrected for and normalized using mean field bias computed from simulations. Thus we obtain normalized local variance maps as defined by Eq. [9]. A dipole is then fit to these normalized local variance maps from data, assuming that the power asymmetry, if any, is dipolar () in nature. We eventually check for presence of any significant power in higher modes () of local variance maps, that could be inducing power asymmetry.

In order to compute the mean field bias to correct the local variance maps from data, we proceed as follows. The generation of simulations with same beam smoothing and consistent noise levels as data, are described in Section [3.2.3] and Appendix [A]. The simulations thus obtained are filtered to contain the same multipole range as data by applying the high and low pass cosine filters following Eq. [10], [11] and [12].

We then apply the local variance estimator to simulations similar to data ie., computing local variance maps at HEALPix with the same choice of circular disc radii and effective disc fraction of described above, from each simulation. Thus a set of 1000 local variance maps are derived, corresponding to each of the circular disc radii chosen and the galactic mask used. A mean field is thus estimated by taking the mean of the 1000 LV maps from simulations for each choice of disc radius and mask. This is subtracted from and also used to normalize the data LV maps (Eq. [12]). The simulations themselves are also normalized using the mean field LV estimate.

The amplitudes of dipole power asymmetry in data in comparison to simulations are shown in
Fig. [3]. The four plots in that figure correspond to the four choices
of disc radii made to compute the local variances across the sky.
The red, green and blue histograms in each plot of the figure correspond
to the three masks described in the previous section (Fig. [2]).
The amplitude of dipole in LV maps from data are denoted by triangle point types in
each plot. We used the HEALPix
subroutine `remove_dipole`

to obtain the dipole amplitudes and directions from data
and simulations. An estimate of the variance of the normalized local variance maps from
isotropic simulations is used for inverse noise variance weighing of pixels to compute
the dipole amplitudes and directions from normalized LV maps while using the subroutine
`remove_dipole`

.
The dipole amplitudes as well as probability-to-exceed the observed values (values)
corresponding to different choices of circular disc radii for mapping the local variances
in Commander 2015 CMB mode map comprising large angular scales ()
are tabulated in the top block of Table 1. The three rows in that
block correspond to the three masks used. The dispersion in recovered dipole amplitudes
from corresponding simulations are quoted as errors in braces in the table.

range | ds= | ds= | ds= | ds= |
---|---|---|---|---|

, value | , value | , value | , value | |

0.034(10), 0.000 | 0.035(10), 0.000 | 0.035(10), 0.000 | 0.036(10), 0.000 | |

0.037(10), 0.000 | 0.037(10), 0.000 | 0.038(10), 0.000 | 0.039(10), 0.000 | |

0.035(10), 0.000 | 0.035(10), 0.000 | 0.037(9), 0.000 | 0.037(9), 0.000 | |

- | 0.028(11), 0.007 | 0.028(11), 0.006 | 0.028(10), 0.003 | |

- | 0.028(10), 0.005 | 0.029(10), 0.001 | 0.029(10), 0.000 | |

- | 0.026(10), 0.006 | 0.028(10), 0.002 | 0.028(10), 0.001 | |

- | 0.019(11), 0.173 | 0.018(11), 0.166 | 0.018(11), 0.165 | |

- | 0.017(11), 0.219 | 0.018(11), 0.168 | 0.019(11), 0.134 | |

- | 0.018(10), 0.148 | 0.019(10), 0.106 | 0.020(10), 0.086 |

As can be readily noticed from Fig. [3], the amplitude of dipole in normalized LV maps from data denoted by triangle point types are significantly anomalous with a probability, . So none of the 1000 normalized local variance maps from simulations have a dipole of amplitude more than that seen in data. Further a careful observation also reveals that when using smaller disc radii, the distance between data points and the footsteps of the histogram on the right increases, in comparison to using disc radius of (arcmin). From Table 1, we also see that the amplitudes are consistent with each other, rather independent of the galactic mask used which have different masking fractions, as well as the disc radius chosen to compute local variances.

The direction of the dipoles as inferred from these normalized local variance maps are shown in top panel of Fig. [4]. They broadly point towards the CMB kinetic dipole. The dipole directions inferred from data when using galactic mask with larger sky fraction are aligned more closely towards CMB dipole. This suggests that the normalized LV maps’ dipole axes are sensitive to the available sky fraction in the mask used, and consequently to the residual contamination. still present in the CMB map. Various LV maps’ dipole directions () are listed in top block of Table [2] along with the axes proximity to the CMB kinetic dipole ( [38]) direction measured as .

range | ds= | ds= | ds= | ds= |
---|---|---|---|---|

(253.0,20.0) 29.4 | (258.9,19.7) 28.6 | (265.1,22.6) 25.4 | (267.4,22.8) 25.4 | |

(266.5,28.3) 19.8 | (268.8,28.8) 19.5 | (274.5,30.0) 19.8 | (274.6,30.7) 19.1 | |

(271.9,33.7) 15.5 | (273.9,34.9) 15.1 | (283.6,35.1) 19.4 | (284.5,36.0) 19.2 | |

- | (259.3,25.4) 22.8 | (266.5,28.5) 19.6 | (269.2,28.4) 20.0 | |

- | (265.6,37.9) 10.2 | (271.3,38.9) 10.5 | (271.1,39.2) 10.2 | |

- | (272.4,45.1) 6.5 | (284.1,44.9) 14.1 | (285.9,45.4) 15.2 | |

- | (264.6,40.2) 7.8 | (269.3,42.6) 6.5 | (270.9,42.9) 7.0 | |

- | (266.2,52.6) 4.8 | (273.7,53.8) 8.4 | (274.7,52.5) 8.2 | |

- | (279.4,54.6) 11.6 | (288.4,54.1) 16.4 | (289.4,54.1) 17.0 |

Now we check whether there is any significant power in modes other than dipole in the local variance maps. This would inform us whether the anomalous power asymmetry we found is only dipolar in nature or other modes are also contributing to the observed power asymmetry. We find that there is significant power only in the dipole component of the local variance maps from polarization data. As an example, in Fig. [5], we show the angular power spectrum plot of local variance map derived using the range from data, and estimated using disc radius of (arcmin) and the galactic mask mask074. We applied the first order correction to estimate the angular power spectrum, , where is the multipole index of LV maps.

Thus we find a signal of anomalous dipole power asymmetry in CMB mode polarization maps from PLANCK full mission data on large angular scales, broadly pointing in the CMB kinetic dipole direction. The dipole axis inferred is however found to be sensitive to the galactic mask, used to estimate the LV maps from the current release of polarization data.

### 4.2 Scale dependence of amplitude of dipolar power asymmetry

In this section we probe the dipole power asymmetry signal at small angular scales. This will allow us to evaluate any scale dependence in the amplitude of the dipole modulation signal we are studying in CMB mode polarization maps. Given the limited number of circular disc radii we could use with the current data, we exclude the large angular scales in progression to study any such scale dependence with local variance maps.

To this end, we use two filtered CMB mode maps containing small angular scales ie., comprising the multipoles and . A high and low pass filters defined in Eq. [10] and [11] are applied over a width of 20 multipoles at both ends of the multipole ranges mentioned for the two filtered maps. Since these maps contain smaller angular scales, we use only three circular discs of radii , and (arcmin) to get LV maps. The three masks shown in Fig. [2] are used as in the previous section. We again demand availability of atleast of pixels in the effective circular discs, defined as union of the galactic masks and the full circular discs defined locally on HEALPix grid to map local variances.

The amplitude of dipole of normalized local variances maps from the two filtered CMB mode maps are shown in Fig. [6]. The plots in left-hand column of the figure corresponds to using the filtered map containing , and the right-hand column corresponds to those comprising . Each plot has again three histograms in red, green and blue corresponding to use of mask074, mask080 and mask085 galactic masks. The data values are denoted by triangle point types in each plot. The three rows in Fig. [6] correspond to results obtained using the three disc radii , and (arcmin).

One can see that, the amplitude of dipole in LV maps is anomalous only when obtained from the range . The dipole amplitudes in the LV maps from the filtered map containing modes are consistent with amplitude excursions expected in isotropic maps. The dipole power asymmetry amplitudes and their values corresponding to the two filtered maps are listed in the second and third block of Table [1]. The corresponding dipole directions are shown in the middle and bottom panel of Fig. [4]. These dipole directions are still broadly pointing towards CMB Doppler boost direction, nearly independent of the choice of disc radius and galactic mask used. But the orientation of the dipole axes move towards higher latitudes around the kinetic dipole axis by excluding low multipoles. The dipole directions corresponding to LV maps from these two filtered CMB mode maps are listed in the second and third block of Table [2], along with their angular separation from the CMB kinetic dipole direction.

### 4.3 LV analysis using mode mask

Here we redo the exercise of local variance analysis described in the preceding sections, using mode masks. So far we used masks derived in basis ie. masks, where . Now we repeat the LV analysis using mode foreground masks created as described below, as a test of consistency of our results obtained using maps.

Once we obtain the difference maps 030-070 and 353-070 from raw PLANCK stokes satellite maps, as described in Sec [3.2.2], we convert these difference maps into mode maps which contain only foregrounds. The two mode difference maps are combined to form a junk map, where we retain only the largest absolute pixel-value, among the two, in the output junk map, following the procedure of Tegmark et al. (2003) [37]. This map is then thresholded such that the resultant masks contains and ie., unmasked sky fractions. The masks thus obtained, denoted here after as mask080E and mask085E, are shown in Fig. [7].

We now go on to estimate the local variance maps using these mode masks. The analysis procedure is same as described in Sec. [4.1], except for using the newly derived masks. Here we only analyze the filtered mode CMB map that includes the large angular scales ie., the data maps containing the full range , and also use a single disc radius of (arcmin) to compute local variances.

The dipole amplitudes derived by fitting a dipole to normalized local variance maps in this case are shown in left-hand panel of Fig. [8] in comparison to simulations. The two histograms in that panel correspond to the two masks used, and the data points are denoted by triangle point types. The direction of the dipole as inferred from these normalized data LVE maps are shown in the same figure in the right-hand panel. The observed dipole amplitudes and their significances using mask080E and mask085E are tabulated in Table 3.

The dipole amplitudes now obtained using space masks are found to be lower in amplitude than those obtained when using masks in Sec. [4.1]. Their significances have become and , respectively, when mask080E and mask085E are used. The dipole directions now seem to be more closely aligned with the CMB kinetic dipole direction. For reference, the dipole directions from using the masks viz., mask080 and mask085 with the same circular radius of (arcmin) to obtain LV maps are also indicated in the right panel of Fig. [8].

mask | value | |||
---|---|---|---|---|

mask080E | 0.026(10) | 0.003 | (253.7,40.2) | 10.7 |

mask085E | 0.027(9) | 0.001 | (266.8,40.8) | 7.5 |

## 5 Conclusions

In this paper, we tried to probe the presence of hemispherical power asymmetry in CMB mode polarization data, that was so far studied in CMB temperature data. To this end we used Stokes CMB maps from PLANCK full mission data estimated using Commander component separation procedure. From these full-sky Stokes maps available at HEALPix and having a Gaussian beam smoothing of (arcmin), we derived mode polarization map of CMB. This mode map is obtained at a low resolution of and with a smoothing given by a Gaussian beam of (arcmin). The low resolution mode map is further low-pass and high-pass filtered to get three maps comprising different ranges of multipoles viz., , and to test for any scale dependence of the amplitude of hemispherical power asymmetry modelled as dipole modulation.

We find a tentative evidence for the presence of dipole modulation in polarization data as well, with an amplitude of depending on the mask, and disc radius used on large angular scales (ie., using the multipole range ). The lower values are found when space masks are used for computing local variance maps, instead of masks where . This sheds light on some important analysis aspect that the mask used for analyzing CMB mode maps have to be constructed carefully, and the masks that are provided with second PLANCK data release may not be readily applicable in space. LV maps from intermediate and small angular scales of the CMB polarization map considered in the present work are found to have lower dipole amplitude. The amplitude of the dipoles on large angular scales were found to be very significant with a probability of , which goes down to on small angular scales. The direction of the dipoles inferred from fitting a dipole to the local variance maps are found to be roughly pointing towards the CMB kinetic dipole direction, and broadly independent of the angular scales involved, and masks and circular disc radii used to obtain LV maps. The dipole axis inferred when using space masks are found to be more closely aligned with the CMB dipole direction.

However, with the know issue of systematic artefacts on large angular scales of PLANCK polarization data, only precise future measurements of polarization or a future PLANCK plarization data release with carefully processed large angular scales will validate these findings, beyond any systematics still present in the multipole range studied here.

## Acknowledgements

We are very thankful to Anthony Banday for a careful reading of our draft and giving his valuable comments. We also thank Changbom Park for helpful discussions.
P.K.A. and A.S. would like to acknowledge the support of the National Research Foundation of Korea (NRF-2016R1C1B2016478).
The present work has extensively made use of the publicly available HEALPix^{8}^{9}

## Appendix

## Appendix A Treatment of noise simulations

As described in Sec. [3.2.3], effective noise realizations were obtained using
co-addition of smoothed, downgraded frequency specific FFP8 noise maps. However the
noise levels in polarized FFP realizations were under estimated compared to data.
One can use the half-mission half-difference (HMHD) map from the two Commander
polarized CMB maps derived using half mission data as proxy
to the true noise level in PLANCK Commnader mode CMB map.
These half-mission maps are also provided through PLANCK second public release^{10}

Here we compare the noise levels in simulations and data, in an attempt to be able to use the FFP polarization noise simulations with reasonable modifications. To reliably derive the noise angular power spectrum from Commnader polarized HMHD maps, and from the FFP8 polarized noise simulations, we have to first modify the masks shown in Fig. [2]. The half-mission (HM) data itself comes with additional mask owing to missing pixels in respective data sets, that is further processed and downgraded to . The original half-mission mask is made available at . To this we also add the missing pixels mask from single year maps to make it more conservative. This combined mask is downgraded to and smoothed with a Gaussian beam of (arcmin). Then a pixel cut-off of 0.95 is imposed whereby pixels with pixel-value greater than this are set to ’’ and rest to ’’. The mask thus obtained is shown in Fig. [9].

The three masks shown in Fig. [2] are then taken in union with the suitably extended, downgraded half-mission mask of Fig. [9]. However in order to compute the noise power spectrum from partial sky, it would be useful to apodize the masks to avoid ringing in the recovered noise power spectra due to sharp mask boundaries (the so-called Gibbs phenomenon). So we use a cosine profile to apodize the three union (galactic+HM) masks over a width of (arcmin) at the mask boundaries. Recall that these masks are at , which have pixels of side (arcmin), assuming them to be square shaped. So, choosing to apodize over a width of 10 pixels strips along the mask boundary equals an apodization width of .

These apodized union masks are used to mask the Commander HMHD mode map. Then the pseudo full-sky noise power spectra were obtained following the MASTER algorithm [40]. We repeat the same exercise on 1000 FFP8 co-added noise polarization maps by applying the union masks with differing sky fractions. A comparison of respective power spectra is shown in top panel of Fig. [10]. For brevity only the pseudo full-sky noise power spectrum from using the union mask mask074+HM on data and simulations is shown in that panel.

In the bottom panel of Fig. [10], we show the ratio of pseudo full-sky noise power spectrum from Commander HMHD mode map () to the corresponding mean pseudo full-sky noise spectrum from 1000 effective FFP noise maps (), obtained using the three union masks. Also shown is the mean of the three noise power spectra ratios (as dark green solid line), smoothed using a running window of (ie., average over at each multipole, ). We find that this mean of the ratio of noise power spectra is approximately constant upto . By taking the average of all three noise power spectra ratios in the range we get a factor , that is plotted as a black dashed curve in the same figure (bottom panel). Thus we scale the effective noise maps by a factor of to alleviate their noise levels in comparison to true noise in data maps.

As explained in Sec. 4.3, we devised additional masks directly in space. Since the masks we used have varying mask fractions, we expect the scaling factor found with them to still hold when estimated using masks. Nevertheless, we analyzed the HMHD noise proxy using mask085E (see Fig. [7]) following the same procedure as with masks. We find the ratio of the Commander HMHD to mean simulation noise power spectra obtained by employing mask085E to have a similar level of mismatch, and also result in approximately same factor of for scaling the simulations.

### Footnotes

- aluri@kasi.re.kr
- shafieloo@kasi.re.kr
- http://irsa.ipac.caltech.edu/data/Planck/release_2/all-sky-maps/matrix_cmb.html
- http://irsa.ipac.caltech.edu/data/Planck/release_2/ancillary-data/
- http://irsa.ipac.caltech.edu/data/Planck/release_2/ancillary-data/
- http://irsa.ipac.caltech.edu/data/Planck/release_2/ancillary-data/
- In the absence of foregrounds, an ILC procedure will reduce to an inverse noise variance weighing if the pixels are uncorrelated i.e., , where denotes the frequency channels, are the weights which are given by the inverse of r.m.s squared of smoothed/downgraded polarized noise maps at at each frequency, and denoted a pixel/sky direction in the sky. Hence ignoring the correlations between pixels due to beam, we used the diagonal approximation to the noise covariance matrix to obtain suitable simulations that complement the observed maps. This is justified following our interest in probing large angular scales of the polarized CMB sky. However the effective polarization noise maps thus obtained from FFP simulations are still not usable due to an additional complication. See text for more details.
- https://healpix.jpl.nasa.gov/
- http://www.esa.int/Planck
- http://irsa.ipac.caltech.edu/data/Planck/release_2/all-sky-maps/matrix_cmb.html

### References

- Eriksen H. K., Hansen F. K., Banday A. J., Gorski K. M., & Lilje P. B., 2004, ApJ, 605, 14
- Bennett C. L. et al., 2011, ApJS, 192, 17
- Planck Collaboration : Ade P. A. R. et al., 2014, A & A, 571, A23
- Planck Collaboration : Ade P. A. R. et al., 2016, A & A, 594, A16
- Hansen F. K., Banday A. J., & Gorski K. M., 2004, MNRAS, 354, 641
- Hajian A., Souradeep T., & Cornish N., 2005, ApJL, 618, L63
- Prunet S., Uzan J. P., Bernardeau F., & Brunier T., 2005, Phys. Rev. D, 71, 083508
- Armendariz-Picon C., 2006, JCAP, 03, 002
- Eriksen H. K., Banday A. J., Gorski K. M., Hansen F. K., & Lilje P. B., 2007, ApJL, 660, L81
- Bernui A., Mota B., Reboucas M. J., & Tavakol R., 2007, IJMPD, 16, 411
- Lew B., 2008, JCAP, 09, 023
- Bernui A., 2008, Phys. Rev. D, 78, 063531
- Hansen F. K.,, Banday A. J., Gorski K. M., Eriksen H. K., & Lilje P. B., 2009, ApJ, 704, 1448
- Hanson D., & Lewis A., 2009, Phys. Rev. D, 80, 063004
- Paci F. et al., 2010, MNRAS, 407, 399
- Paci F. et al., 2013, MNRAS, 434, 3071
- Flender S., & Hotchkiss S., 2013, JCAP, 09 033
- Rath P. K., & Jain P., 2013, JCAP, 12, 014
- Santos L. et al., 2014, A & A, 569, A75
- Akrami Y. et al., 2014, ApJL, 784 L42
- Bernui A., Oliveira A. F., & Pereira T. S., 2014, JCAP, 10, 041
- Adhikari S., 2015, MNRAS, 446, 4232
- Quartin M., & Notari A., 2015, JCAP, 01, 008
- Aiola S., Wang B., Kosowsky A., Kahniashvili T., & Firouzjahi H., 2015, Phys. Rev. D, 92, 063008
- Ghosh S., Kothari R., Jain P., & Rath P. K., 2016, JCAP, 01, 046
- Gordon C., Hu W., Huterer D., & Crawford T., 2005, Phys. Rev. D, 72, 103002
- Gordon C., 2007, ApJ, 656, 636
- Planck Collaboration: Adam R. et al., 2016, A & A, 594, A9
- Goldberg J. N., et al., 1967, Jou. Math. Phys., 8, 2155
- Zaldarriaga M., Seljak U., 1997, Phys.Rev. D, 55, 1830
- Kamionkowski M., Kosowsky A., Stebbins A., 1997, Phys.Rev. D, 55, 7368
- Zaldarriaga M., 1998, ApJ, 503, 1
- Shafieloo, A., Clifton, T. & Ferreira, P. 2011, JCAP, 9, 017
- Colin, J., Mohayaee, R., Sarkar, S., & Shafieloo, A. 2011, MNRAS, 414, 264
- Shafieloo, A. 2012, JCAP, 5, 024
- Shafieloo, A. 2012, JCAP, 8, 002
- Tegmark M., de Oliveira-Costa A., Hamilton A. J., 2003, Phys. Rev. D, 68, 123523
- Kogut A. et al., 1993, ApJ, 419, 1
- Gorski K. M. et al., 2005, ApJ, 622, 759
- Hivon E. et al., 2002, ApJ, 567, 2