Issue 
A&A
Volume 591, July 2016



Article Number  A97  
Number of page(s)  4  
Section  Cosmology (including clusters of galaxies)  
DOI  https://doi.org/10.1051/00046361/201424164  
Published online  21 June 2016 
Statistical isotropy violation in WMAP CMB maps resulting from noncircular beams
IUCAA, PO Bag 4, Ganeshkhind, 411007 Pune, India
email: sanjit@iucaa.in
Received: 8 May 2014
Accepted: 21 April 2016
Statistical isotropy (SI) of cosmic microwave background (CMB) fluctuations is a key observational test to validate the cosmological principle underlying the standard model of cosmology. While a detection of SI violation would have immense cosmological ramification, it is important to recognise their possible origin in systematic effects of observations. The WMAP seven year (WMAP7) release claimed significant deviation from SI in the bipolar spherical harmonic (BipoSH) coefficients A_{ll}^{20} and A_{l2l}^{20}>. Here we present the first explicit reproduction of the measurements reported in WMAP7, confirming that beam systematics alone can completely account for the measured SI violation. The possibility of such a systematic origin was alluded to in WMAP7 paper itself and other authors but not as explicitly so as to account for it accurately. We simulate CMB maps using the actual WMAP noncircular beams and scanning strategy. Our estimated BipoSH spectra from these maps match the WMAP7 results very well. It is also evident that only a very careful and adequately detailed modelling, as carried out here, can conclusively establish that the entire signal arises from noncircular beam effect. This is important since cosmic SI violation signals are expected to be subtle and dismissing a large SI violation signal as observational artefact based on simplistic plausibility arguments run the serious risk of “throwing the baby out with the bathwater”.
Key words: cosmic background radiation / cosmology: observations / methods: data analysis
© ESO, 2016
1. Introduction
Cosmic microwave background (CMB) measurements have paved the way for the transition of cosmology to a precision science. The recent observations by WMAP and Planck favour a minimal six parameter ΛCDM cosmology (Hinshaw et al. 2013; Planck Collaboration XVI 2014). An assumption of statistical isotropy (SI) of CMB is a fundamental tenet of the standard model of cosmology, which increasingly precise measurements now enable to be observationally tested. Indeed several SI violation anomalies have been detected in recent measurements of CMB temperature anisotropies (Bennett et al. 2013, 2011; Planck Collaboration XXIII 2014). Although there is no dearth of proposed theoretical models (Ackerman et al. 2007; Ramazanov & Rubtsov 2012; Rotti et al. 2011) to explain these SI violations, it is necessary to model all the known systematics and account for the biases introduced by them before probing any deeper for subtle cosmological effects.
It is well known that noncircular instrumental beams can induce artefacts in CMB measurements. The resultant effect on the CMB angular power spectrum (C_{l}) has been extensively studied in literature (Tegmark et al. 2003; Ashdown et al. 2007; Fosalba et al. 2002; Hinshaw et al. 2007; Mitra et al. 2009, 2004; Souradeep et al. 2006; Souradeep & Ratra 2001; Das & Souradeep 2015). It has been recognised that noncircular beams can lead to SI violation in the observed maps that are not captured in the angular power spectrum (Hajian & Souradeep 2003; Mitra et al. 2004). Bipolar spherical harmonic (BipoSH) expansion of the CMB two point correlation function provides a novel tool to probe SI violation in CMB maps (Hajian & Souradeep 2003). Detection of nonzero BipoSH spectra, (L> 0), is the basic probe of SI violation.
In this work, we address the high statistical significance detection of SI violation by the WMAP7 team in BipoSH coefficients and from the WMAP V and Wband maps (Bennett et al. 2011). Their paper mentions that “it seems very likely that the observed quadrupolar effect is the result of incomplete handling of beam asymmetries”, however it was never explicitly shown that the measured BipoSH spectra could be emulated by performing detailed numerical simulations. In a recent paper, we provided a formalism to study noncircular beams in BipoSH representation and show that mild noncircularity at levels comparable to WMAP beams do generate significant BipoSH spectra in observed CMB maps (Pant et al. 2016). In another approach, Hanson et al. (2010) invoked the anisotropic primordial power spectrum (aPPS) model (Pullen & Kamionkowski 2007; Ma et al. 2011) and an approximate scan pattern to estimate this effect. However, as shown in Kumar et al. (2015), it is unlikely that a surrogate estimator can be extrapolated to explain the distinct features observed in BipoSH spectra caused by noncircular beams, namely the observed shift in the location of the primary peak or the zero crossings that are clearly detected in WMAP7 analysis of in the two bands. These features arise owing to the combined effect of structures in CMB power spectra and nontrivial beam characteristics (Joshi et al. 2012; Pant et al. 2016), which are discussed in Sect. 3. Given that the detected SI violation signal caused by noncircular beams is very strong, it is important to account accurately for the effect before probing deeper into the residual signature, if any, for cosmological anomalies. Use of an appropriate estimator is thus necessary for the analysis of data from the current and future highresolution CMB experiments to avoid incorrect or inaccurate prediction of the systematic effects.
In this paper, to estimate the full systematic effect accurately, we simulate WMAP observations by incorporating published instrumental beam maps along with the real nontrivial scanning strategy for the V and W bands. BipoSH spectra obtained from these simulated maps match those recovered by an identical analysis of WMAP7 observed maps, providing the first direct demonstration that the systematic effects that are due to noncircular beams is sufficient to account for the WMAP7 measurements of quadrupolar anomaly.
2. Review of BipoSH
The standard model of cosmology assumes the Universe to be SI. The CMB temperature anisotropies can be expressed in terms of the spherical harmonics as , where a_{lm} are the coefficients of expansion on this basis. If the CMB temperature anisotropies are described as a Gaussian random field, warranted by recent precise measurements made by Planck (Planck Collaboration XXIII 2014; Planck Collaboration XXIV 2014), specifying the twopoint correlation for this field completely characterises its statistical properties. Further, if the basic assumption of SI is valid, then the twopoint correlation function is fully described by the equivalent harmonic space quantity, the angular power spectrum C_{ℓ} defined as , where the angular brackets ⟨···⟩ denote an average over an ensemble of CMB realizations.
However, in a nonSI Universe, the BipoSH spectra capture the full description of the twopoint correlation function (Hajian & Souradeep 2003). The BipoSH spectra are defined by the expression (Bennett et al. 2011), (1)where Π_{L}: = (2L + 1)^{1/2} and are the ClebschGordan coefficients. While the BipoSH spectrum is the standard angular power spectrum, detecting nonvanishing power in the remaining BipoSH spectra forms the basic criteria for probing SI violation.
3. Simulation and analysis
The WMAP satellite scanned the sky temperature in five different frequency bands at 22, 30, 40, 60, and 90 GHz. The WMAP7 team searched for SI violation signals in the 60 GHz (V) and 90 GHz (W) bands, since these channels are the least foregroundcontaminated and detected significant nonvanishing BipoSH spectra and . We simulate observed CMB maps, by closely replicating the WMAP instrument characteristics and scanning strategy for the same bands.
The V and W band are comprised of two and four detectors respectively. Since each detector has its unique beam pattern and scan path, it is important to carry out the simulations for each differencing assembly (DA) independently. We generate 30 realisations of SI CMB skies from the bestfit WMAP7 angular power spectrum using HEALPix (Górski et al. 2005). The procedure that was followed to generate the time order data (TOD) from these realisations and mapmaking is described below.
3.1. Scan pattern
Each DA of the WMAP satellite involves a pair of radiometers. The WMAP satellite measures the temperature difference between the DA radiometer pairs. Each of these antennae are inclined at an angle ~ 70.5° to the symmetry axis of the satellite. The satellite spins around this axis with 2.2 min period and slowly precesses about the SunEarth axis at an inclination angle of 22.5°, with a precession period ~ 1 h. This particular method of scanning is adapted by the WMAP satellite to reduce noise in the data. This entire system is located at L2 and is moving around the Sun with a period of one year.
We have used an analytical approximation (Das & Souradeep 2015; Moss et al. 2011) to calculate the pointing of the satellite to obtain the beam locations and their orientations on the sky as a function of time.
3.2. Beam functions
The observed temperature along a particular direction is related to the underlying sky temperature by (2)where is the beam response function. The shape and size of the intrinsic instrumental beams, elliptical Gaussian (EG) profiles fits for FWHM of W (V) band beam are ~ 13.2′ (~ 21.0′) with eccentricity ~ 0.4 (0.46) (Jarosik et al. 2007; Mitra et al. 2004). However, EG approximation of the beam profile does not provide adequately accurate model for the problem. Specifically, such an approximation cannot explain subtle zero crossing features seen in the BipoSH spectra caused by noncircular beams (Pant et al. 2016). Moreover, the effect of scanning creates varying effective beam profile across the sky (Mitra et al. 2011). Therefore, for an accurate estimation of the effect of noncircular beams on the BipoSH spectra, a full convolution of SI CMB sky with the actual beam, along with scanning strategy, is critical.
WMAP beam maps are available as 600 × 600 arrays on the LAMBDA website^{1} extending 24° on each side. Thus each of the pixels covers 2.2′ of the beam profile. Apart from the central peak, the beam consists of many small structures around the peak, as seen in Fig. 1, which span a much wider solid angle. This includes subdominant side peaks and a noisy annular region with negative pixel values, which remain significant up to a radial distance of ~ 4° from the central peak. The contribution from the outer region of the beam is critical for reproducing the amplitude and features of the observed BipoSH spectra, most importantly when explaining the flipping of signs of BipoSH spectra seen in WMAP7 year observations. We note that the negative pixels of the beam sum up to ~ 10% of the total integrated beam profile. For completeness, in our analysis we have taken a ± 6° cutoff for the beam around the pointing direction, while numerically convolving it with the sky map using Eq. (3). We also tested that ignoring the part beyond 6° does not affect our results.
3.3. Generation of timeordered data and mapmaking
The observed temperature anisotropy along a given pixel direction can be expressed as a discretised convolution as (3)where is a pointing direction, is the central direction of the ith pixel and is the set of all the pixels for which is nonzero. If the beam function were considered to be elliptical Gaussian, pixels within 3−4σ would have been sufficient for the numerical integration. However, for the actual WMAP beams, as shown in Fig. 1, where we plotted the beam map of the A side of W3 DA, this area cutoff is not valid and requires integration over the entire beam.
To compute , we include all the pixels inside the radius ± 6°, from that pointing direction and then interpolate the beam function at the centre of the pixels for convolution using Eq. (3). The WMAP team provides the temperature sky maps at N_{side} = 512. At this resolution, the distance between two adjacent pixel centres is 3.8′. As the sizes of many of the noisy peaks in the skirt of the beams are smaller than 3.8′, these features will be missed in numerical integration, hence the convolution is not precise. To overcome this problem, all the simulations in the paper are generated at a resolution N_{side} = 1024, for an accurate evaluation of the numerical beam convolution of the map, following which the pointing direction is converted to the pixel numbers that correspond to a N_{side} = 512 resolution map.
To evaluate the beam function at the centre of each pixel, it is necessary to use a fast and accurate interpolation scheme. Linear interpolation at the centre of each pixel by considering the value of the beam function at the four pixels around it only works for the central peak. However, the features in the annular region that surround the central peak of the beams being of a rapidly varying nature, the error due to linear interpolation is found to be unacceptably large. In this region, a nonlinear interpolation scheme, such as the cubic or the spline interpolation, may provide reasonable results. But these methods are computationally expensive. To overcome this issue, we upgraded the beam function from 600 × 600 grid to 2400 × 2400 grid using MATLAB’s spline2 interpolation scheme. This highresolution grid enables us to use a linear interpolation of the beam in the convolution process without compromising the precision and is computationally fast. We also checked that a further increment in the interpolated grid size is not required to improve accuracy.
We follow this procedure to generate noiseless TOD, where each sample is obtained by Eq. (3). The TOD vector for each DA can then be written as d = AT, where T is the scanned sky temperature and A is the pointing matrix, where each row consists of two nonzero components, one + 1 and one −1. We use the Jacobi iteration method (Das & Souradeep 2015) for mapmaking, although other possible methods are discussed in Hamilton (2003).
Fig. 1
Enlarged view of the Aside beam of W3 DA. The central part of the beam is truncated to show the subdominant structures in the beam. Apart from the central peak the beam consists of many features extending up to ~ 4 degree in an annular ring. The beam also takes unphysical negative values in the annular region, which perhaps arise from noisy measurements. We have separated the negative (orange) and positive (blue) sensitivity parts for visual clarity. A small portion of the beam (with both positives and negatives) is enlarged to show the presence of subtle detailed structures. Motivation for considering all these details in the shape of the beam are described in the text. 
3.4. BipoSH analysis in the presence of masks
The observed CMB maps are contaminated by foregrounds. Even though the W and V band maps are the least foreground contaminated, the region close to the galactic equator has to be masked to extract information on the genuine CMB signal. The application of masks in analyses as well as the presence of anisotropic noise in the maps, correlates different harmonic modes of CMB anisotropies, resulting in nonvanishing BipoSH spectra (L ≠ 0). To remove this bias from observed WMAP maps, masked with KQ75, firstly we evaluate BipoSH coefficients from 1000 masked realisations generated from a synfast subroutine of HEALPix for the bestfit CMB power spectrum and average beam transfer function. This estimate is used to debias and arrive at the final WMAP BipoSH spectra. The results obtained closely match the WMAP findings.
We now apply our WMAP analysis identically to our simulated maps with a noncircular beam and scan strategy. Even though we do not include foregrounds in our simulations, we use the mask in our analysis to compare them with the results derived from the observed maps.
4. Results
The WMAP7 team detected an SI violation signal in BipoSH spectra in both the V and the W bands. The primary goal of this paper is to make a detailed computation of the BipoSH spectra arising from the actual WMAP beam, coupled to the scan strategy, and establish conclusively that the measured nonzero L = 2 BipoSH spectra in WMAP7 can be entirely attributed to this systematic effect. We generate 30 simulations each for four DAs in the W band (W1, W2, W3, W4) and two DAs in the V band (V1, V2) by convolving CMB sky realisations with the actual WMAP beam and scanning strategy and performed BipoSH analysis on each set separately. The BipoSH spectra obtained from the true and simulated datasets are shown in Fig. 2.
In Fig. 2, we plot the WMAP data points (red error bars). The average s from the four DAs are plotted in black dashed lines. The location of peak of the spectra is set by both the first CMB acoustic peak, and the parameter b_{l2} characterising beam noncircularity (Pant et al. 2016). A closer inspection of the spectrum reveals a change in sign, for both bands, which can only be explained by taking cognisance of a large noisy region around the main peak in the WMAP beams shown in Fig. 1. Thus the spectra from the simulated maps match the observed spectra very well, replicating the nontrivial features, leaving very little scope for accommodating additional systematic or cosmological effects. The error bars are estimated from 1000 simulations used for estimating the mask bias after adding anisotropic noise (Rotti et al. 2011) and scaled for different cases based on the total number of simulations. The error bars are not plotted for each individual DA to avoid clutter in the plot, rather we plotted the error bars of the average spectra. We plotted two error bars, the regions shaded in orange denote standard errors estimated from the DA simulations alone, while in the grey region we also include the standard errors on the estimated mask bias from 1000 simulations in quadrature. The large errorbars seen at highl in the WMAP analysis is mostly due to noise in the data. Since we plot errorbars on spectra averaged over many simulations, they are smaller (at least by a factor of few). The nonzero measurements of spectrum reported in WMAP7 (not plotted here) can also be convincingly reproduced by our simulations. For the W band, it can be seen in the figure that BipoSH spectra for W1 and W2 are very close to each other and similarly for W3 and W4. This trend is essentially due to similarities in their respective beam transfer functions.
To highlight the reliability of our method and to emphasise the use of a proper estimator for SI violation, we plot the for V and W bands derived from the aPPS modelbased estimates (Hanson et al. 2010) in Fig. 2. This model predicts , where g_{20} lacks a multipole dependence that is crucial for characterising the beam’s noncircularity, hence fails in recovering the shift in the location of the peaks and the zero crossings of the BipoSH spectra. The match at the first peak in for the V band is due to a mere coincidence that both b_{l2} and C_{l} peak around l ≈ 220, which is not the case for the W band, where b_{l2} peaks near l ≈ 350 (see Fig. 5 Pant et al. 2016).
We also plot the ratios of observed and predicted for V and W bands for our method as well as the same derived from the aPPS modelbased estimates in Fig. 3. Our method clearly produces a close match with observation, while it is clear that the aPPS estimator is not designed to capture the observed trend.
Fig. 2
Comparison between the beam window function corrected BipoSH spectra obtained from WMAP maps (red errorbars), observed in V (top) and W (bottom) bands, and our detailed simulations of the respective channels (solid lines). The average BipoSH spectra across the DAs are denoted by the dashed lines. The saffron band shows the standard error derived from all the simulations for each band only (total 60 simulations for V and 120 for W), while the grey band also combines the error on the mask bias correction. Spectra from the simulated maps match the observed spectra very well, correctly reproducing the location of the peak and the zero crossings. Also plotted are the spectra (dashdot lines) evaluated from the aPPS model parameterised by g_{20} provided in Table 1 of Hanson et al. (2010) to emphasise that a surrogate estimator is unlikely to capture these subtle effects. 
5. Discussion and conclusion
In this paper, we have demonstrated for the first time that the quadrupole anomaly measured in WMAP7 maps can be closely reproduced by incorporating the effect of noncircular beams and scan strategy of the satellite. This clearly indicates, as was also suspected by the WMAP team, that this signal does not have any cosmological origin. The absence of similar detections in BipoSH spectra and in the more recent measurements by the Planck satellite (Planck Collaboration XXIII 2014) reinforces this claim, where an analysis strategy similar to the one described in this paper was employed to account for the noncircular beams. Prior to that, the WMAP nine year analysis team released a set of beamsymmetrised (deconvoled) maps, which also did not show this particular SI violation signature (Ramazanov & Rubtsov 2014). However, as cautioned by the WMAP team, the deconvolution procedure renders the resultant maps unsuitable for cosmology analyses. Also, new SI violation signals are present in those maps at highl (Bennett et al. 2013). However, in our socalled forward, approach we evade the complications of deconvolution by the ability to estimate a correctable bias.
The close match we obtain clearly shows that finer details like foreground residuals, coupling between beam and mask, weak lensing and anisotropic noise, etc., which, in principle, could have been important for these BipoSH spectra, cause substantially weaker SI violation in WMAP7 relative to the noncircular beam. This provides a reliable route to probing SI violation from the standard model of cosmology, which will be targeted by the current and upcoming CMB missions.
Fig. 3
Ratio of the BipoSH spectra for the V and the W bands. The trend recovered by our detailed analysis clearly matches the WMAP7 observations. The very low and high multipole regions which have poor S/N are indicated with pink shading. We also plot the ratio corresponding to Hanson et al. (2010) estimates, which is consistent with observations at lowℓs, but, not surprisingly, far off at highℓs. 
Acknowledgments
We would like to thank Saurabh Kumar, Gary Hinshaw, and Eiichiro Komatsu for useful comments and help. S.D. and A.R. acknowledge the CSIR, India for financial support through Senior Research Fellowships. S.M. acknowledges the Fast Track grant SR/FTP/PS030/2012 of SERB, India. T.S. acknowledges Swarnajayanti fellowship grant of DST India. Computations were carried out at the HPC facilities at IUCAA.
References
 Ackerman, L., Carroll, S. M., & Wise, M. B. 2007, Phys. Rev. D, 75, 083502 [NASA ADS] [CrossRef] [Google Scholar]
 Ashdown, M. A. J., Baccigalupi, C., Balbi, A., et al. 2007, A&A, 467, 761 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bennett, C. L., Hill, R. S., Hinshaw, G., et al. 2011, ApJS, 192, 17 [NASA ADS] [CrossRef] [Google Scholar]
 Bennett, C. L., Larson, D., Weiland, J. L., et al. 2013, ApJS, 208, 20 [Google Scholar]
 Das, S., & Souradeep, T. 2015, J. Cosmol. Astropart. Phys., 5, 012 [Google Scholar]
 Fosalba, P., Doré, O., & Bouchet, F. R. 2002, Phys. Rev. D, 65, 063003 [NASA ADS] [CrossRef] [Google Scholar]
 Górski, K. M., Hivon, E., Banday, A. J., et al. 2005, ApJ, 622, 759 [NASA ADS] [CrossRef] [Google Scholar]
 Hajian, A., & Souradeep, T. 2003, ApJ, 597, L5 [NASA ADS] [CrossRef] [Google Scholar]
 Hamilton, J.C. 2003, ArXiv Astrophysics eprints [arXiv:astroph/0310787] [Google Scholar]
 Hanson, D., Lewis, A., & Challinor, A. 2010, Phys. Rev. D, 81, 103003 [NASA ADS] [CrossRef] [Google Scholar]
 Hinshaw, G., Nolta, M. R., Bennett, C. L., et al. 2007, ApJS, 170, 288 [NASA ADS] [CrossRef] [Google Scholar]
 Hinshaw, G., Larson, D., Komatsu, E., et al. 2013, ApJS, 208, 19 [NASA ADS] [CrossRef] [Google Scholar]
 Jarosik, N., Barnes, C., Greason, M. R., et al. 2007, ApJS, 170, 263 [NASA ADS] [CrossRef] [Google Scholar]
 Joshi, N., Das, S., Rotti, A., Mitra, S., & Souradeep, T. 2012, ArXiv eprints [arXiv:1210.7318] [Google Scholar]
 Kumar, S., Rotti, A., Aich, M., et al. 2015, Phys. Rev. D, 91, 043501 [NASA ADS] [CrossRef] [Google Scholar]
 Ma, Y.Z., Efstathiou, G., & Challinor, A. 2011, Phys. Rev. D, 83, 083005 [NASA ADS] [CrossRef] [Google Scholar]
 Mitra, S., Sengupta, A. S., & Souradeep, T. 2004, Phys. Rev. D, 70, 103002 [NASA ADS] [CrossRef] [Google Scholar]
 Mitra, S., Sengupta, A. S., Ray, S., Saha, R., & Souradeep, T. 2009, MNRAS, 394, 1419 [NASA ADS] [CrossRef] [Google Scholar]
 Mitra, S., Rocha, G., Górski, K. M., et al. 2011, ApJS, 193, 5 [NASA ADS] [CrossRef] [Google Scholar]
 Moss, A., Scott, D., & Sigurdson, K. 2011, J. Cosmol. Astropart. Phys., 1, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Pant, N., Das, S., Rotti, A., Mitra, S., & Souradeep, T. 2016, J. Cosmol. Astropart. Phys., 3, 035 [NASA ADS] [CrossRef] [Google Scholar]
 Planck Collaboration XVI. 2014, A&A, 571, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXIII. 2014, A&A, 571, A23 [NASA ADS] [CrossRef] [EDP Sciences] [MathSciNet] [PubMed] [Google Scholar]
 Planck Collaboration XXIV. 2014, A&A, 571, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pullen, A. R., & Kamionkowski, M. 2007, Phys. Rev. D, 76, 103529 [NASA ADS] [CrossRef] [Google Scholar]
 Ramazanov, S. R., & Rubtsov, G. I. 2014, Phys. Rev. D, 89, 043517 [NASA ADS] [CrossRef] [Google Scholar]
 Ramazanov, S., & Rubtsov, G. 2012, J. Cosmol. Astropart. Phys., 2012, 033 [Google Scholar]
 Rotti, A., Aich, M., & Souradeep, T. 2011, ArXiv eprints [arXiv:1111.3357] [Google Scholar]
 Souradeep, T., & Ratra, B. 2001, ApJ, 560, 28 [NASA ADS] [CrossRef] [Google Scholar]
 Souradeep, T., Mitra, S., Sengupta, A., Ray, S., & Saha, R. 2006, New A Rev., 50, 1030 [NASA ADS] [CrossRef] [Google Scholar]
 Tegmark, M., de OliveiraCosta, A., & Hamilton, A. J. 2003, Phys. Rev. D, 68, 123523 [NASA ADS] [CrossRef] [Google Scholar]
All Figures
Fig. 1
Enlarged view of the Aside beam of W3 DA. The central part of the beam is truncated to show the subdominant structures in the beam. Apart from the central peak the beam consists of many features extending up to ~ 4 degree in an annular ring. The beam also takes unphysical negative values in the annular region, which perhaps arise from noisy measurements. We have separated the negative (orange) and positive (blue) sensitivity parts for visual clarity. A small portion of the beam (with both positives and negatives) is enlarged to show the presence of subtle detailed structures. Motivation for considering all these details in the shape of the beam are described in the text. 

In the text 
Fig. 2
Comparison between the beam window function corrected BipoSH spectra obtained from WMAP maps (red errorbars), observed in V (top) and W (bottom) bands, and our detailed simulations of the respective channels (solid lines). The average BipoSH spectra across the DAs are denoted by the dashed lines. The saffron band shows the standard error derived from all the simulations for each band only (total 60 simulations for V and 120 for W), while the grey band also combines the error on the mask bias correction. Spectra from the simulated maps match the observed spectra very well, correctly reproducing the location of the peak and the zero crossings. Also plotted are the spectra (dashdot lines) evaluated from the aPPS model parameterised by g_{20} provided in Table 1 of Hanson et al. (2010) to emphasise that a surrogate estimator is unlikely to capture these subtle effects. 

In the text 
Fig. 3
Ratio of the BipoSH spectra for the V and the W bands. The trend recovered by our detailed analysis clearly matches the WMAP7 observations. The very low and high multipole regions which have poor S/N are indicated with pink shading. We also plot the ratio corresponding to Hanson et al. (2010) estimates, which is consistent with observations at lowℓs, but, not surprisingly, far off at highℓs. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.