CMB Parameter Estimation
CMB Parameter Estimation
CMB Parameter Estimation
Effect of noncircularity of
experimental beam on CMB
parameter estimation
Santanu Das,a Sanjit Mitraa and Sonu Tabitha Paulsonb
a
b
Abstract. Measurement of Cosmic Microwave Background (CMB) anisotropies has been playing a
lead role in precision cosmology by providing some of the tightest constrains on cosmological models
and parameters. However, precision can only be meaningful when all major systematic effects are
taken into account. Non-circular beams in CMB experiments can cause large systematic deviation in
the angular power spectrum, not only by modifying the measurement at a given multipole, but also
introducing coupling between different multipoles through a deterministic bias matrix. Here we add
a mechanism for emulating the effect of a full bias matrix to the Planck likelihood code through
the parameter estimation code SCoPE. We show that if the angular power spectrum was measured
with a non-circular beam, the assumption of circular Gaussian beam or considering only the diagonal
part of the bias matrix can lead to huge error in parameter estimation. We demonstrate that, at
least for elliptical Gaussian beams, use of scalar beam window functions obtained via Monte Carlo
simulations starting from a fiducial spectrum, as implemented in Planck analyses for example, leads
to only few percent of sigma deviation of the best-fit parameters. However, we notice more significant
differences in the posterior distributions for some of the parameters, which would in turn lead to
incorrect errorbars. These differences can be reduced, so that the errorbars match within few percent,
by adding an iterative reanalysis step, where the beam window function would be recomputed using
the best-fit spectrum estimated in the first step.
Contents
1 Introduction
2
2
3
3
4
5
5
5
6
6
Introduction
Measurements of Cosmic Microwave Background (CMB) opened a whole new era in theoretical
physics. Not only it made the standard cosmological model widely acceptable, but it offered precise measurements of cosmological parameters. Several high precision ground based and space based
experiments were carried out in the past decades to measure CMB anisotropies. Recently, WMAP
[14] and Planck [5, 6] mapped the full CMB sky with few arcmin level resolutions. However, high
precision measurements demand accurate accounting of the systematic errors. Systematic errors can
arise at different stages in CMB analysis, such as, foreground cleaning [7, 8], instrument calibration [9, 10], measurement of beam response [11, 12] and estimation of the angular power spectra [13]
from the observed maps. In this paper we focus on the effect of noncircular beams [1421], one of the
most important and challenging sources of systematic error.
Observed CMB sky is a convolution of the underlying true sky with the instrumental beam
response function. Accounting for the beam is thus necessary for measuring the statistical properties
of the true sky. This turns out to be trivial if the beam is symmetric about the pointing direction, an
unbiased estimator of the true angular power spectrum can be readily obtained. However, an actual
beam response function is never so symmetric. The intrinsic optics of the instruments, aberrations due
to the placement of detectors away from the principle optic axis, non-uniform distribution of pointing
directions in a pixel, finite sampling duration in scanning and many other effects can distort a beam
to make it non-circular. The asymmetries become progressively important in different analyses as
the experiments strive to extract almost all the information embedded in the anisotropies. Beam
asymmetry may lead to spurious effects in the measured CMB data. It can bias the angular power
spectrum in a non-trivial way by introducing coupling between different multipoles through a bias
matrix [16]. It can also cause statistical isotropy violation in the CMB anisotropy maps [2225].
Here we study how the distortions in the angular power spectra caused by non-circular beams affect
cosmological parameter estimation.
If the bias matrix that couples different multipoles was available, it appears obvious that one
would use it to obtain an unbiased estimator of angular power spectra and perform parameter estimation starting from the unbiased estimator. However, this process would encounter two major hurdles.
First, attempt to unbias the estimator will introduce covariances among different multipoles. Second,
perhaps the more important one, it is highly non-trivial to obtain the bias matrix for an experiment
with non-circular beam shapes and complex scanning strategy. In literature the near-diagonal components of the bias matrix has been computed to leading order for simple scanning strategies only for
temperature anisotropies [16, 18], though efforts are on to extend the method to polarization also [26].
Our present work thus, in fact, evaluates the need for investing (enormous) effort in getting the full
bias matrix for experiments like Planck in order to perform unbiased analyses.
Here we develop a mechanism to include the effect of a bias matrix in the parameter estimation
code SCoPE [27] in conjunction with the Planck likelihood code [13, 28]. We use Planck likelihood
code in order to get realistic representation of noise. However, instead of deconvolving the observed
power spectrum to get an unbiased estimator, we apply inverse distortion to the theoretical spectra
in the likelihood code to get the desired posterior distributions. This step alleviates the need for incorporating a non-trivial covariance matrix. We use different scalar transfer functions to do parameter
estimation from the same convolved spectra and study the differences in posterior distributions with
the correct estimates.
It is very common in CMB analysis to use an effective scalar window function, which inherently
assumes the beam to be azimuthally symmetric about the pointing direction. However, if the effective
window function was derived using Monte Carlo simulations [19], including the full details of beams
and scanning strategy, as was done for Planck analyses, one may be able to deceive an analysis by
closely emulating the effect of a non-circular beam, as long as the fiducial power spectrum used for
the simulations is close to the true one. Here we use WMAP and Planck best-fit power spectra to get
two different scalar transfer functions and use them for estimating parameters from a test observed
power spectrum obtained by convolving the Planck best fit spectrum with a bias matrix. However,
since a full bias matrix for a mission like Planck does not exist, for numerical computation we must
limit ourselves to a case where it is available. We use an elliptical Gaussian beam of similar size and
ellipticity as one of the Planck high frequency detectors and non-rotating scan pattern.
This paper is organised as follows. We present a brief primer on the connection between noncircular beams and window function in section 2. The general strategy for studying the effect of
non-circular beam on cosmological parameters and the results of our study are included in section 3.
Section 4 presents the conclusion and discussions.
Different experiments observe the CMB anisotropy field by scanning the sky through an instrumental
beam of finite resolution. The observed time ordered data (TOD) is passed through a refined pixel
binning procedure which is generally kept unaware of the beam shapes. This is because deconvolution
of all the TOD samples (few trillions for Planck detectors) is computationally prohibitive. The
resultant map can then be expressed as a convolution of the true sky with an effective beam function.
If the effective beam for every direction was the same and symmetric about each direction, one could
el is trivially biased, C
el = B 2 Cl , where
show that the angular power spectrum of the measured map C
l
Cl is the angular power spectrum of the true sky and Bl is the Legendre transform of the azimuthally
symmetric beam function. It is indeed very common in CMB analysis to use an effective scalar window
function Bl2 . However, perhaps for every CMB experiments, beams are non-circular. Non-circular
el = P 0 All0 Cl0 , as we briefly review below.
beams make the effective window function tensorial, C
l
2.1
g (
can be expressed as the convolution of
The observed CMB anisotropy T
q) in a given direction q
0 ) plus additive noise n(
the true sky T (
q) with the instrumental beam function B(
q, q
q),
Z
g (
0 ) T (
q, q
q0 ) + n(
q) ,
(2.1)
T
q) =
dq 0 B(
S2
S2
S2
and assuming Statistical Isotropy (SI) of the true sky, it can be shown [16] that the expected power
spectrum of the observed map can be expressed as
el i =
hC
Z
Z
l
X
1
1
2
g (
g (
2)
h|e
alm | i =
dq 1
dq 2 hT
q1 ) T
q2 )i Pl (
q1 q
2l + 1
4 S2
S2
(2.3)
m=l
lX
max
(2.4)
l0 =0
Here, lmax is the maximum harmonic multipole, ClN is the angular power spectrum of noise and the
bias matrix,
Z
Z
2l0 + 1
2 )Wl0 (
2) ,
d
All0 :=
dq 2 Pl (
q1 q
q1 , q
(2.5)
1
q
16 2 S2
S2
where,
Z
2 ) :=
Wl (
q1 , q
S2
dq 10
S2
10 ) B(
20 ) Pl (
20 ) .
dq 20 B(
q1 , q
q2 , q
q10 , q
(2.6)
It is interesting to note that the noise-free two-point correlation function of the observed CMB
anisotropy sky can be expressed as
X
(2l + 1)
g (
g (
hT
q) T
q0 )i =
l=0
2.2
0) .
Cl Wl (
q, q
(2.7)
Circular Beams
It is fairly straightforward to show that if the beam is azimuthally symmetric about the pointing
0 ) B(
0 ), so that the beam can be expanded in terms of Legendre
direction, that is, B(
q, q
qq
polynomials,
lmax
1 X
0) ,
0 ) B(
0) =
(2l + 1) Bl Pl (
qq
(2.8)
B(
q, q
qq
4
l=0
the observed angular power spectrum is trivially biased, All0 = ll0 Bl2 ,
el = Bl2 Cl + ClN .
C
(2.9)
Thus in this case it is easy to get an unbiased estimator using the scalar window function Bl2 ,
el C N ] ,
Clscalar = Bl2 [C
l
(2.10)
assuming that the noise power spectrum, ClN , can be precisely estimated independently from instrument noise characteristics. The co-variance of the unbiased estimators is given by
Cov(Clscalar , Clscalar
) =
0
2
2ll0
Cl + Bl2 ClN .
2l + 1
(2.11)
The above equations imply that in case of circular beams, there is no coupling between power spectrum
at different multipoles.
Note that, unless otherwise mentioned, we will use a twiddle ( e ) to denote an observed quantity,
a bar ( ) on a quantity to denote an estimator and a superscript scalar to denote that the estimator
is derived using a scalar transfer function.
2.3
Non-circular Beams
location could be transformed to other directions [14]. The final expression for the observed power
spectrum can then be expressed in terms of the spherical harmonic transforms of the beam (blm )
pointing at the z-axis of spherical polar coordinates. Numerical computation then shows that, for a
trivial scanning strategy and elliptical Gaussian beams, to the leading order the bias matrix All0 has
a large number of small off diagonal components, which can imply a significantly large difference in
power spectrum if estimated assuming a circular beam [16]. Hence, from Eq. (2.4), one should define
the true unbiased estimator as
X
N
e0
Cl =
A1
(2.12)
ll0 [Cl Cl0 ] .
l0
Even though the approximate calculations provide a good estimate of the level of the effect, in
order to fully account for the effect one needs the true bias matrix computation involving the actual
non-circular beam shape, the exact scan pattern and beyond leading order computation. This is a
challenging task and has not yet been accomplished in literature. Efforts are on to compute effective
bl2 for a complex experiment, incorporating the effect of scanning strategy, in turn enabling one to
compute the true bias matrix to a reasonable accuracy using only leading order computation [29].
2.4
Above we described that circular beam window functions can lead to significant systematic error in
power spectrum estimation. But, realistic bias matrices for non-circular beams are not available. This
situation might seem discouraging ! In practice though there is a middle-ground. One can compute a
scalar beam window function through Monte Carlo simulations, incorporating beam asymmetry and
scanning strategy. A large number of maps are simulated from a fiducial angular power spectrum,
which are convolved with actual non-circular beams. The scalar window function can then be estimated by taking the average ratio of the convolved to unconvolved power spectra. However, a scalar
window function corresponds to a circular beam (blm m0 ), as non-vanishing blm lead to off-diagonal
terms in the bias matrix. So this method essentially replaces the complex effective beams, which can
also vary across the sky, by one effective circular beam. The goal of this paper is to study the accuracy
and adequacy of this scalar beam window functions in cosmological parameter estimation.
el )
The effective beam window function is defined as the average of the ratio of the convolved (C
fid
to unconvolved fiducial (Cl ) simulated maps
el /Clfid i .
Wl := hC
(2.13)
In this paper we simulate a convolved power spectrum by multiplying an unconvolved spectrum (Cl )
with a given bias matrix
lX
max
el :=
All0 Cl0 .
(2.14)
C
l0 =0
This method alleviates the need for performing the Monte Carlo simulations for estimating the scalar
transfer function, as neither All0 nor the fiducial spectrum Cl are random realizations. However, this
is possible here because we are restricting ourselves to a case for which the bias matrix is available,
which is a primary requirement for this study.
Note that, in this paper we will be using two different unconvolved spectra, one is the true
spectrum of the sky, which the analysis aims to recover, the other one is a fiducial spectrum (Clfid )
for estimating the scalar window function via Eq. (2.13). In both the cases convolution is done through
Eq. (2.14). If one includes the the full bias matrix in the unbiased estimator, as in Eq. (2.12), the
estimator is truly unbiased. This is because, combining with Eq. (2.4), one gets
X
N
e0
hCl i =
A1
(2.15)
ll0 [hCl i Cl0 ] = Cl .
l0
However, this need not be the case if one uses a scalar transfer function Wl to define the unbiased
estimator, as in Eq. (2.10), because
el i ClN ] = W 1
hClscalar i = Wl1 [hC
l
lX
max
l0 =0
All0 Cl0 .
(2.16)
Inserting the convention for true and fiducial described above one gets
Plmax
true
0 =0 All0 Cl0
scalar
fid
.
hCl
i = Cl Pllmax
fid
0
l0 =0 All Cl0
(2.17)
It is easy to see from the above equation that if Clfid = Cltrue , one would get hClscalar i = Cltrue . However,
this still does not ensure that the posterior distribution of the parameters can be correctly recovered
by using the best-fit Cl as Clfid . This is because parameter estimation codes compute likelihood for
different sets of Cltrue , obtained by sampling the parameter space, to get the posterior distributions.
Hence Clfid can match Cltrue at most for one realisation of the set of parameters. Moreover, since the
aim of an experiment is to estimate the best-fit Cl one can not guess Cltrue a priori. Hence, there is no
reason to assume that Clfid = Cltrue . Nevertheless, the difference between Clfid and Cltrue need not be
very large either, as we do have approximate knowledge of Cltrue from previous experiments and only
the nearby regions of the parameter space are sampled. In this work, as a realistic test case, we
first take Clfid to be the WMAP best-fit spectrum and Cltrue to be the Planck best-fit spectrum and
study whether this difference causes significant deviation in cosmological parameter estimation. We
then study whether the true posterior distribution can be recovered (not only the best-fit parameter
values) by taking Clfid and Cltrue to be the same, the Planck best-fit power spectrum.
The broad approach we follow here is to generate non-circular beam convolved power spectrum,
analyse it including the corresponding bias matrix in the parameter estimation codes and compare
the results with the ones obtained using different scalar transfer function.
3.1
the ellipticity of the beam used here is 1/ 1 2 = 1.4, which is on a slightly higher side compared
to the Planck detectors, leading to a more conservative estimate which will become evident later.
We chose a trivial scan pattern that keeps the alignment of the beam fixed with respect to the
local meridian over the whole sky. We could not include a realistic scan strategy as the bias matrix for
such a case has not been numerically computed in literature. In fact one of the main aim of this work
is to study the need for investing enormous effort in computing the realistic bias matrix. However, this
is not a big assumption either for this work. We are studying whether the scalar transfer functions
can mimic a typical bias matrix in parameter estimation. Had we incorporated the full scan in the
bias matrix, the matrix would be different, but we believe the broad characteristics would remain the
same, as can be seen in the bias matrix plots for two different toy scan patterns presented in [16].
Hence this study should remain valid for realistic scan patterns.
3.2
We study four different variants of scalar transfer functions (Wl ) for calculating the estimator of the
el /Wl as listed below:
power spectrum Clscalar := C
1. Wl is taken as Bl2 , where Bl is the Legendre transform of a circular Gaussian beam, as defined
in Eq. (2.8), whose radius is the geometric mean of the semi-major and semi-minor axis of the
chosen elliptical Gaussian beam, preserving the total collection solid angle of the beam.
2. Wl is taken as the diagonal components of the of the bias matrix (All0 ), i.e. Wl = All .
3. Wl is computed from a fiducial Cl using Eq.(2.13). We use two different Clfid
(a) Clfid is taken as the best fit WMAP-9 Cl ,
(b) Clfid is taken as the best fit Planck Cl .
3.3
We compute the posterior distributions of cosmological parameters from the original Cl and the
power spectrum estimators (Cl ) obtained by using different scalar and tensor transfer functions described above. Parameter estimation is done using the code SCoPE [27]. We add log-likelihoods
from commander v4.1 lm49.clik, lowlike v222.clik and CAMspec v6.2TN 2013 02 26.clik [28].
The effect of beam non-circularity is insignificantly small at the low multipoles. Since commander
and lowlike only use Cl s from low multipoles, modifications to Cl with bias matrix or either of the
above scalar transfer functions have negligible effect on these likelihoods. The only likelihood that
gets affected in this process is that from CAMspec.
CAMspec likelihood provides 2 = (Cl Clscalar )T [C]1 (Cl Clscalar ), where [C] is the covariance
matrix and Cl is a theoretical spectrum for a given set of parameters. In most situations, to perform
parameter estimation for a specific cosmological model, one runs MCMC chains over the parameter
space and for each set of parameters Cl and 2 are computed. This procedure does not question
the validity of the estimator Clscalar and its covariance, which is in direct contrast to the situation
considered in this work. Here we would like to modify Clscalar to capture the effect of P
using different
,
scalar and tensor beam window functions. However, if we modify Clscalar by say Wl1 l0 All0 Clscalar
0
the covariance matrix also gets modified to [Wl1 All0 ][C][Wl1 All0 ]T . Therefore, for running MCMC
with CAMSpec likelihood we need to change both Clscalar and the covariance matrix and feed those
into the CAMSpec likelihood code. Rather, we use an alternate approach. For each MCMC realisation
we P
multiply the Cl s obtained from CAMB[30] with Wl A1
ll0 , that is, replacing theoretical Cl by
scalar and [C] fixed. This step provides us the intended 2 at a reduced
0
Wl l0 A1
ll0 Cl keeping Cl
complication. We then estimate the posterior distributions of cosmological parameters from the exact
estimator Cl incorporating a bias matrix and the estimators Clscalar for different approximations to
the beam window described in Section 3.2.
Note that, the above procedure is justified since the aim of this paper is not a reanalysis of
Planck data, which would require a nearly exact treatment of beams and noise, the aim here is
to verify if the scalar transfer functions are adequate for representing the distortions described by
a bias matrix. So the above procedure in a way assumes that Planck measurements are correctly
represented by the supplied beam transfer functions and noise covariance matrix and the default
likelihood produces correct values. We distort the theoretical Cl by an inverse bias matrix and check
if the effect can be compensated in the likelihood codes by multiplying it with a scalar transfer function.
We use Planck likelihood code here only to get realistic noise characteristics and resolution. If one
wanted to reanalyse Planck power spectra with bias matrices (assuming that they have somehow
been made available for each channel), one would have to modify the observed spectra and the noise
covariance matrices separately for each frequency.
3.4
Results
Wl
A1
ClPlanck
0
ll0
Planck
Cl
l0
best fit Cl and Wl is estimated in three different (inexact) ways. This plot shows how much error
is introduced in the power spectrum estimator for not considering the non-circularity of the beam
properly. It can be seen that if we consider the beam to be circular Gaussian then the error in Cl
at high multipoles is very large, more than 3%. However, if we consider the diagonal components
of the All0 as the scalar transfer function then the error involved at high multipoles is lesser, 1%.
These levels of error in Cl can affect the cosmological parameters very significantly. On the other
hand Wl computed through the forward method from the fiducial Cl , as in Eq.(2.17), performs much
better. In the plot we choose the fiducial Cl as ClWMAP . Here the error at high multipoles is less
than 0.1%. Of course, if Wl was estimated using ClPlanck as the fiducial spectrum, Fl would be 1 at
every l. However, as mentioned before, ClPlanck is used as the true spectrum here, which the analysis
aims to recover, it can not be guessed a priori. We could also use some other set of parameters to
generate the true spectrum. The Planck best fit parameters are chosen here as they promise to
be the closest to reality in the history of CMB measurements.
1.005
1
0.995
0.99
Fl
Wl =
0.985
Wl =
W M AP
0
l0 All C l0
W
MAP
Cl
All0 ll0
l0
0.98
0.97
0.965
Wl
P lanck
A1
ll0 Cl0
ClPlanck
l0
500
1000
l
1500
2000
Figure 1.
Comparison of different effective window functions. A convolved CMB power spectrum is
generated by multiplying the Planck best-fit with a bias matrix for an elliptical Gaussian beam. It is then
corrected by a scalar transfer function and its ratio to the original Cl is plotted. The scalar transfer function
is obtained in different ways: by assuming a circular Gaussian beam (red), the diagonal of the bias matrix
(blue) and using a forward approach with WMAP best-fit as the fiducial Cl (green). If we used Planck
best-fit as the fiducial spectrum, the correction would be perfect and we would get Fl = 1. The good news is
that even with WMAP best-fit fiducial Cl , Fl is close to unity. The other two scalar transfer functions lead to
significant deviation of Fl from unity, implying highly incorrect estimation of the unbiased power spectrum.
We have done parameter estimation for the standard six parameter LCDM cosmology, namely
{c ,b ,h, ,ns ,As } using SCoPE [27]. In figure (2) we show the likelihood contours (in blue) for
the six parameters obtained from Planck best fit Cl , without introducing any instrumental effect.
According to the scheme we followed to introduce the bias matrix, this case is equivalent of accounting
for the full bias matrix in the analysis through the exact unbiased estimator [Eq. (2.12)]. Hence, this
is the correct distribution we aim to recover. The contours corresponding to circular Gaussian beam
approximation, the first Clscalar estimator discussed in Section 3.2, are overlaid (in red). It can be
clearly seen that circular Gaussian approximation leads to very different posteriors, in some of the
cases the contours get shifted by more than 2, implying that the mean values of the parameters
can be highly incorrect as well. Using Wl = All also leads to similar, but little better, results. The
likelihood contours are not shown here, the marginalised distributions are shown in figure 4 (in blue).
Planck data analysis uses scalar transfer functions obtained via Monte Carlo simulations starting from a fiducial Cl as mentioned in Section 3.2. We first use the best fit WMAP Cl as the fiducial
Cl and use the resulting scalar transfer function corrected power spectrum Clscalar for parameter estimation. We find that though the mean parameters are recovered to few percent of sigma accuracy, the
distributions are significantly different, which would lead to wrong errorbars. Likelihood contours are
not shown for this case, but marginalised distributions are plotted in figure 4. Next we show results
m h2
b h2
0.470
0.946
0.266
0.849
0.045
0.493
0.953
0.277
0.846
0.026
0.696
0.214
0.458
0.116
0.696
0.213
0.449
0.093
0.70
0.279
0.843
0.073
0.65
0.296
0.839
0.062
0.275
0.972
0.285
0.964
0.023
0.022
0.15
0.05
0.98
n
0.064
0.049
log(As 1010)
0.94
3.20
3.00
0.02
0.14 0.15 0.02
2
2
b h
m h
0.65
0.70
0.05
0.150.94
0.98
3.2
10
log(As 10 )
Figure 2. We plot the likelihood contours for the six cosmological parameters starting from the correct
l and a scalar transfer function corrected one (C
lscalar )
(unbiased using full bias matrix) pseudo-Cl estimator C
obtained for circular Gaussian beam. It clearly shows that circular Gaussian approximation for a non-circular
beam can lead to very large deviation in the posterior distribution. The correlation coefficients between pairs
of parameters are shown in upper triangle of the plot.
with Planck best-fit as the fiducial spectrum. In figure 3 the likelihood contours obtained from the
full bias matrix corrected Cl are plotted in blue and those from the estimator Clscalar are in red. The
contours are close for most parameters, but not for all, the 1-D marginalised distributions provide a
clearer picture.
In figure 4 we overlay the 1-D marginalized probability distributions for the four different cases
discussed in section 3, along with the distribution from the correct bias matrix corrected Cl (green).
It shows that if the WMAP best fit Cl is used as the fiducial spectrum for deriving the transfer function
(black) then also the average values of the distribution are almost the same, but the distributions are
narrower, i.e. the standard deviations are smaller than that of the original distribution sometimes
by as much as 20%. This implies that the choice of wrong fiducial Cl for obtaining scalar transfer
function may lead to wrong distribution of the parameter space, but not significantly changing the
mean. When the fiducial spectrum matches the best-fit (yellow), the differences in distributions are
reduced and the standard deviations differ by few percent, as in figure 3.
Thus, we can conclude that the method adapted by the Planck team for obtaining the parameters can provide reasonably accurate posteriors, at least if the beams are not too different from
elliptical gaussian and the fiducial spectrum is close to the actual best-fit. However, the actual best-fit
may not be known a priori. To address this issue, we suggest running the analysis first with a transfer
function derived from a fiducial spectrum constructed from a reasonable initial guess of parameters,
obtain the best-fit parameters, recompute the transfer function with the best-fit spectrum and finally
get the posterior distribution and update the best-fit by using the new transfer function.
m h2
b h2
0.441
0.947
0.266
0.837
0.021
0.493
0.953
0.277
0.846
0.026
0.666
0.169
0.378
0.070
0.696
0.213
0.449
0.093
0.275
0.821
0.048
0.296
0.839
0.062
0.273
0.965
0.285
0.964
0.023
0.022
0.70
0.65
0.14
log(As 1010)
0.06
0.96
0.042
0.94
0.049
3.15
3.00
0.02
0.14 0.15 0.02
2
2
h
m h
b
0.65
0.70
0.06
0.14 0.94
0.96
3.15
log(As 1010)
Figure 3. We plot the likelihood contours for the six cosmological parameters starting from the correct
l (in blue) and a scalar transfer function corrected one (C
lscalar ),
(bias matrix corrected) pseudo-Cl estimator C
where the transfer function was obtained using Planck best-fit as the fiducial Cl (in red). The correlation
coefficients between pairs of parameters are shown in upper triangle of the plot. The plot shows that using a
scalar transfer function derived via Monte Carlo simulations provides reasonably close best-fit values, which
is clearer in the 1-D marginalized distributions. However, the posterior distributions do show deviations from
the actual distributions, which become more significant if WMAP best-fit was used as the fiducial spectrum
[see figure 4] instead of Planck, as the later may not be known a priori. These differences can be reduced
by redoing the analysis using transfer functions derived from the best-fit spectrum obtained in the first step,
the likelihood contours would then correspond to the red contours in this figure.
We have estimated cosmological parameters from observed CMB temperature power spectrum by
correctly accounting for non-circular beams with the help of a bias matrix and also by using different
scalar beam window functions. We show that if the proper beam window function is not used, the
power spectrum estimator gets biased and the posterior distribution of parameters becomes significantly different. Namely, assuming a circular Gaussian beam creates serious departure from the true
posterior, while using a scalar transfer function estimated through Monte Carlo simulations involving
non-circular beams does well in recovering the best fit values to few percent of sigma. However,
the posterior distributions show significant differences, leading to incorrect errorbars, if the fiducial
spectrum is very different from the actual best-fit. The posterior distributions can be brought closer
to actual by introducing an iterative step, where the scalar transfer function is reestimated using the
best-fit spectrum obtained in the first step.
Correct
Circular beam
0.8
0.6
0.4
0.2
Diagonal
0.8
0.6
0.4
0.2
0.11
0.13
0.02
0.02
(0.932T)
(0.345T)
( 0.000T)
( 0.004T)
(0.010T)
0.097 0.0126
0.090 0.0115
0.087 0.0119
0.087 0.0113
0.087 0.0095
0.15
0.689 0.0108
0.680 0.0105
0.675 0.0109
0.675 0.0103
0.675 0.0086
( 1.275T)
( 0.458T)
( 0.000T)
(0.012T)
( 0.004T)
0.8
0.6
0.4
0.2
0.94
0.98
ns
( 0.814T)
( 0.224T)
( 0.000T)
(0.052T)
(0.054T)
0.70
h
0.8
0.6
0.4
0.2
0.65
bh2
0.8
0.6
0.4
0.2
0.05
WMAP
0.8
0.6
0.4
0.2
ch2
0.118 0.0024
0.119 0.0023
0.120 0.0025
0.120 0.0023
0.120 0.0019
Planck
0.968 0.0059
0.959 0.0057
0.955 0.0059
0.955 0.0056
0.955 0.0046
( 2.240T)
( 0.747T)
( 0.000T)
(0.020T)
(0.021T)
3.00
log(As 1010)
3.099 0.0242
3.087 0.0220
3.083 0.0228
3.082 0.0217
3.082 0.0182
3.20
( 0.703T)
( 0.182T)
( 0.000T)
(0.051T)
(0.055T)
Figure 4. One dimensional marginal probability distributions for all the four scalar transfer functions mentioned in Section 3.2 are plotted along with the correct/true distribution that corresponds to the case where
the bias matrix is accounted for. We show the average and the standard error below the plots. The deviation
of the mean from the correct mean is shown in the bracket in terms of T , where T is the standard deviation
of the correct distribution for the parameter in context [that is, in the plot for ns , T is the standard deviation
of the distribution for ns (green)]. Notice that both the yellow and the black marginalised distributions have
means shifted only by few percent T from the correct analysis, illustrating that the scalar transfer functions
derived through Monte Carlo simulations give accurate mean parameters, while the other approximations
badly fail. However, the posterior distributions are different, the differences are more significant when the
fiducial spectrum used for transfer function estimation is not close to the best fit (compare black and yellow
with green).
We have considered here an elliptical Gaussian beam and a non-rotating scanning strategy, as
the bias matrix is available in only in such cases. If the bias matrix was available for a more realistic
beam shape and scanning strategy, this work could be immediately repeated for that case. However,
we believe that even in such complex situations, the bias matrix would be different but still would not
show any dramatically different features. Since we have not used any special characteristics of the bias
matrix, this work should still remain valid. Testing that hypothesis of course requires the whole bias
matrix to be computed, though this work suggests that, at least for temperature and low-multipole
polarization analysis, a full bias matrix may not be necessary for cosmological parameter estimation.
Obtaining the bias matrix for polarization analysis is a even bigger challenge, though efforts are
on [26]. We emphasise that the systematic effects of this kind become important when the errorbars
are small even beyond the multipoles corresponding to the beam width. Which happens, for example,
in CMB temperature anisotropy measurement with Planck. However, for polarization measurement
the errorbars are still large, one may not need such fine corrections for current experiments.
Here we would also like to mention that for our analysis we only change ClT T , whereas a beam
10
function that affects ClT T should affect ClT E . However the effect on ClT E is not significant in the
present analysis, considering that the effects of the noncircular beam only affect the high multipoles
and presently the polarization and the cross power spectrum measurement is limited low multipoles.
It is well known that masking of foreground contaminated regions in the maps couples the lower
multipoles. The full bias matrix for masking with circular beam is routinely used in CMB analysis [31].
Masking and non-circular beam effects are strong in different multipole ranges, still there can be small
coupling between these two effects [18]. However, the effect of this on parameter estimation should be
small and may not be significant at the level of precision attainable with current CMB experiments.
The bottomline is that the need for incorporating full bias matrices in parameter estimation
may not be crucial with current sensitivities of CMB measurements. Use of scalar transfer functions
derived through Monte Carlo simulations, similar to what was done in Planck analysis, provides
sufficient accuracy in best-fit parameter estimation and can also provide somewhat accurate posterior
distributions if handled carefully. However, efforts should be invested in estimating the full bias
matrices in these complex situations for putting more precise bounds on the cosmological parameters.
Moreover, polarised bias matrices must be computed, to the leading order to start with, at least for
performing a study like this one to make a statement about the adequacy of scalar transfer functions
in high resolution polarisation analysis in Planck and post-Planck era.
Acknowledgments
We would like to thank Krzysztof G
orski and Tarun Souradeep for useful discussions. We have
used the HPC facility at IUCAA for the required computation. SD acknowledge Council for Science
and Industrial Research (CSIR), India, for the financial support as Senior Research Fellows. SM
acknowledges the support of the Science and Engineering Research Board (SERB), India through the
Fast Track grant SR/FTP/PS-030/2012.
References
[1] G. Hinshaw et al., Three-year wilkinson microwave anisotropy probe (wmap*) observations:
Temperature analysis, Astrophys. J. Suppl. 170 (2007) 288, [astro-ph/0603451].
[2] C. L. Bennett et al., Seven-year wilkinson microwave anisotropy probe (wmap) observations: Are there
cosmic microwave background anomalies?, Astrophys. J. Suppl. 192 (2011) 19.
[3] G. Hinshaw et al., Nine-year wilkinson microwave anisotropy probe (wmap) observations: cosmological
parameter results, Astrophys. J. Suppl. 208 (2013) 19, [arXiv:1212.5226].
[4] C. L. Bennett et al., Nine-year wilkinson microwave anisotropy probe (wmap) observations: final maps
and results, Astrophys. J. Suppl. 208 (2013) 20, [arXiv:1212.5225].
[5] Planck Collaboration, P. A. R. Ade, et al., Planck 2013 results. xvi. cosmological parameters,
Astronomy and Astrophysics 571 (2013), no. A16 66, [arXiv:1303.5076].
[6] Planck Collaboration, P. A. R. Ade, et al., Planck 2013 results. xxiii. isotropy and statistics of the cmb,
Astronomy and Astrophysics 571 (2013), no. A23 48, [arXiv:1303.5083].
[7] M. Tegmark, A. de Oliveira-Costa, and A. J. S. Hamilton, High resolution foreground cleaned cmb map
from wmap, Phys. Rev. D 68 (2003) 123523.
[8] Planck Collaboration, P. A. R. Ade, et al., Planck 2013 results. XII. Diffuse component separation,
Astronomy and Astrophysics 571 (Nov., 2014) A12, [arXiv:1303.5072].
[9] Planck Collaboration, N. Aghanim, et al., Planck 2013 results. V. LFI calibration, Astronomy and
Astrophysics 571 (Nov., 2014) A5, [arXiv:1303.5066].
[10] Planck Collaboration, P. A. R. Ade, et al., Planck 2013 results. VIII. HFI photometric calibration and
mapmaking, Astronomy and Astrophysics 571 (Nov., 2014) A8, [arXiv:1303.5069].
[11] K. M. Huffenberger et al., Measuring Planck beams with planets, Astronomy and Astrophysics 510
(Feb., 2010) A58+.
11
[12] G. Rocha et al., Markov chain beam randomization: a study of the impact of PLANCK beam
measurement errors on cosmological parameter estimation, Astronomy and Astrophysics 513 (Apr.,
2010) A23+, [arXiv:0907.5254].
[13] Planck Collaboration, P. A. R. Ade, et al., Planck 2013 results. XV. CMB power spectra and likelihood,
Astronomy and Astrophysics 571 (Nov., 2014) A15, [arXiv:1303.5075].
[14] T. Souradeep and B. Ratra, Window function for noncircular beam cosmic microwave background
anisotropy experiment, Astrophys. J. 560 (2001) 28, [astro-ph/0105270].
[15] P. Fosalba, O. Dore, and F. R. Bouchet, Elliptical beams in cmb temperature and polarization
anisotropy experiments: An analytic approach, Phys. Rev. D 65 (2002) 063003.
[16] S. Mitra, A. S. Sengupta, and T. Souradeep, Cmb power spectrum estimation using noncircular beams,
Phys. Rev. D 70 (2004) 103002, [astro-ph/0405406].
[17] T. Souradeep, S. Mitra, A. Sengupta, S. Ray, and R. Saha, Non-circular beam correction to the cmb
power spectrum, New Astron. Rev. 50 (Dec., 2006) 1030, [astro-ph/0608505].
[18] S. Mitra, A. S. Sengupta, S. Ray, R. Saha, and T. Souradeep, Cosmic microwave background power
spectrum estimation with non-circular beam and incomplete sky coverage, MNRAS 394 (2009) 1419,
[astro-ph/0702100].
[19] S. Mitra et al., Fast pixel space convolution for cosmic microwave background surveys with asymmetric
beams and complex scan strategies: Febecop, ApJS 193 (2011) 5, [arXiv:1005.1929].
[20] Planck Collaboration, N. Aghanim, et al., Planck 2013 results. IV. Low Frequency Instrument beams
and window functions, Astronomy and Astrophysics 571 (Nov., 2014) A4, [arXiv:1303.5065].
[21] Planck Collaboration, P. A. R. Ade, et al., Planck 2013 results. VII. HFI time response and beams,
Astronomy and Astrophysics 571 (Nov., 2014) A7, [arXiv:1303.5068].
[22] D. Hanson, A. Lewis, and A. Challinor, Asymmetric beams and cmb statistical anisotropy, Phys. Rev. D
81 (May, 2010) 103003.
[23] S. Das and T. Souradeep, Leakage of power from dipole to higher multipoles due to non-symmetric
wmap beam, arXiv:1307.0001.
[24] S. Das, S. Mitra, A. Rotti, N. Pant, and T. Souradeep, Statistical isotropy violation in WMAP CMB
maps due to non-circular beams, arXiv:1401.7757.
[25] N. Joshi, S. Das, A. Rotti, S. Mitra, and T. Souradeep, Revealing non-circular beam effect in wmap-7
cmb maps with biposh measures of statistical isotropy, arXiv:1210.7318 [astro-ph.CO] (2012)
[arXiv:1210.7318].
[26] F. A. Ramamonjisoa, S. Ray, S. Mitra, and T. Souradeep, CMB polarization TE power spectrum
estimation with non-circular beam, arXiv:1309.4784.
[27] S. Das and T. Souradeep, SCoPE: an efficient method of Cosmological Parameter Estimation, JCAP 7
(July, 2014) 18, [arXiv:1403.1271].
[28] Planck Collaboration, The explanatory supplement to the planck 2013 results.
http://wiki.cosmos.esa.int/planckpla/index.php/Main_Page, 2013.
[29] N. Pant et al., Estimation of average asymmetry of effective beams using cmb biposh measurements, .
[30] Antony Lewis, Anthony Challinor, Camb. http://camb.info/.
[31] E. Hivon et al., MASTER of the Cosmic Microwave Background Anisotropy Power Spectrum: A Fast
Method for Statistical Analysis of Large and Complex Cosmic Microwave Background Data Sets,
Astrophys. J. 567 (Mar., 2002) 217, [astro-ph/0105302].
12