[astro-ph/0209489] Analyzing weak lensing of the cosmic microwave background using the likelihood function

Authors:  Christopher M. Hirata, Uros Seljak
Abstract:  Future experiments will produce high-resolution temperature maps of the cosmic microwave background (CMB) and are expected to reveal the signature of gravitational lensing by intervening large-scale structures. We construct all-sky maximum-likelihood estimators that use the lensing effect to estimate the projected density (convergence) of these structures, its power spectrum, and cross-correlation with other observables. This contrasts with earlier quadratic-estimator approaches that Taylor-expanded the observed CMB temperature to linear order in the lensing deflection angle; these approaches gave estimators for the temperature-convergence correlation in terms of the CMB three-point correlation function and for the convergence power spectrum in terms of the CMB four-point correlation function, which can be biased and non-optimal due to terms beyond the linear order. We show that for sufficiently weak lensing, the maximum-likelihood estimator reduces to the computationally less demanding quadratic estimator. The maximum likelihood and quadratic approaches are compared by evaluating the root-mean-square (RMS) error and bias in the reconstructed convergence map in a numerical simulation; it is found that both the RMS errors and bias are of order 1 percent for the case of Planck and of order 10--20 percent for a 1 arcminute beam experiment. We conclude that for recovering lensing information from temperature data acquired by these experiments, the quadratic estimator is close to optimal, but further work will be required to determine whether this is also the case for lensing of the CMB polarization field.
[PDF]  [PS]  [BibTex]  [Bookmark]

Discussion related to specific recent arXiv papers
Post Reply
Antony Lewis
Posts: 1941
Joined: September 23 2004
Affiliation: University of Sussex
Contact:

[astro-ph/0209489] Analyzing weak lensing of the cosmic micr

Post by Antony Lewis » September 24 2004

This seems like a nice paper, however some things are rather confusing:

1. p6, "unit determinant of Lambda". Surely the lensing matrix (in any given realisation) does not have unit determinant? Also the claim that the noise and the CMB power are rarely of the same order seems spurious because most of the information for lensing reconstruction will come from the smallest scales that are not noise dominated (i.e. where the noise and cmb are about the same)

2. Section D. The authors use the first order series expansion in this section. Since it is not at all true that this is a good approximation on most scales this may be a very bad approximation. Also setting H[Theta] to H[Theta_hat] is also extremely bad for l>300.

3. p19, estimator map of lensing potential: How do steps (i) and (ii) accomplish a good approximation to [tex]C^{-1}\hat{\Theta}[/tex]?

More generally, can anyone explain why the first order expansion in the lensing potential actually gives a quadratic estimator that is quite accurate? (given that the expansion itself is rather bad). Seems like the expansion is encapsulating most of the non-Gaussianity even though the Gaussian signal is quite wrong on most scales.

Post Reply