Issue 
A&A
Volume 571, November 2014
Planck 2013 results



Article Number  A17  
Number of page(s)  39  
Section  Cosmology (including clusters of galaxies)  
DOI  https://doi.org/10.1051/00046361/201321543  
Published online  29 October 2014 
Planck 2013 results. XVII. Gravitational lensing by largescale structure
^{1}
APC, AstroParticule et Cosmologie, Université Paris Diderot,
CNRS/IN2P3, CEA/lrfu, Observatoire de Paris, Sorbonne Paris Cité, 10 rue Alice Domon et Léonie
Duquet, 75205
Paris Cedex 13,
France
^{2}
Aalto University Metsähovi Radio Observatory, ,
Metsähovintie 114, 02540
Kylmälä,
Finland
^{3}
African Institute for Mathematical Sciences, ,
68 Melrose Road, Muizenberg,
Cape Town, South
Africa
^{4}
Agenzia Spaziale Italiana Science Data Center, ,
via del Politecnico snc,
00133
Roma,
Italy
^{5}
Agenzia Spaziale Italiana, Viale Liegi 26, Roma, Italy
^{6}
Astrophysics Group, Cavendish Laboratory, University of
Cambridge, J J Thomson
Avenue, Cambridge
CB3 0HE,
UK
^{7}
Astrophysics & Cosmology Research Unit, School of Mathematics,
Statistics & Computer Science, University of KwaZuluNatal, ,
Westville Campus, Private Bag
X54001, 4000
Durban, South
Africa
^{8}
Atacama Large Millimeter/submillimeter Array, ALMA Santiago Central
Offices, Alonso de Cordova 3107,
Vitacura, Casilla
763 0355, Santiago, Chile
^{9}
CITA, University ofToronto, 60 St. George St., Toronto, ON
M5S 3H8,
Canada
^{10}
CNRS, IRAP, 9 Av.
colonel Roche, BP
44346, 31028
Toulouse Cedex 4,
France
^{11}
California Institute of Technology, Pasadena, California, USA
^{12}
Centre for Theoretical Cosmology, DAMTP, University of
Cambridge, Wilberforce
Road, Cambridge
CB3 0WA,
UK
^{13}
Centro de Estudios de Física del Cosmos de Aragón
(CEFCA), Plaza San Juan, 1, planta
2, 44001
Teruel,
Spain
^{14}
Computational Cosmology Center, Lawrence Berkeley National
Laboratory, Berkeley,
California,
USA
^{15}
Consejo Superior de Investigaciones Científicas
(CSIC), Madrid,
Spain
^{16}
DSM/Irfu/SPP, CEASaclay, 91191
GifsurYvette Cedex,
France
^{17}
DTU Space, National Space Institute, Technical University of
Denmark, Elektrovej
327, 2800
Kgs. Lyngby,
Denmark
^{18}
Département de Physique Théorique, Université de
Genève, 24 quai E.
Ansermet, 1211
Genève 4,
Switzerland
^{19}
Departamento de Física Fundamental, Facultad de Ciencias,
Universidad de Salamanca, 37008
Salamanca,
Spain
^{20}
Departamento de Física, Universidad de Oviedo, ,
Avda. Calvo Sotelo s/n,
Oviedo,
Spain
^{21}
Department of Astronomy and Astrophysics, University of
Toronto, 50 Saint George Street,
Toronto, Ontario,
Canada
^{22}
Department of Astrophysics/IMAPP, Radboud University
Nijmegen, PO Box
9010, 6500 GL
Nijmegen, The
Netherlands
^{23}
Department of Electrical Engineering and Computer Sciences,
University of California, Berkeley, California, USA
^{24}
Department of Physics & Astronomy, University of British
Columbia, 6224 Agricultural Road,
Vancouver, British
Columbia, Canada
^{25}
Department of Physics and Astronomy, Dana and David Dornsife College
of Letter, Arts and Sciences, University of Southern California, ,
Los Angeles, CA
90089,
USA
^{26}
Department of Physics and Astronomy, University College
London, London
WC1E 6BT,
UK
^{27}
Department of Physics and Astronomy, University of
Sussex, Brighton
BN1 9QH,
UK
^{28}
Department of Physics, Carnegie Mellon University, ,
5000 Forbes Ave, Pittsburgh, PA
15213,
USA
^{29}
Department of Physics, Florida State University, ,
Keen Physics Building, 77 Chieftan
Way, Tallahassee,
Florida,
USA
^{30}
Department of Physics, Gustaf Hällströmin katu 2a, University of
Helsinki, Helsinki,
Finland
^{31}
Department of Physics, Princeton University, ,
Princeton, New Jersey, USA
^{32}
Department of Physics, University of California, ,
Berkeley, California, USA
^{33}
Department of Physics, University of California, ,
One Shields Avenue, Davis, California, USA
^{34}
Department of Physics, University of California, ,
Santa Barbara, California, USA
^{35}
Department of Physics, University of Illinois at
UrbanaChampaign, 1110 West Green
Street, Urbana,
Illinois,
USA
^{36}
Dipartimento di Fisica e Astronomia G. Galilei, Università degli
Studi di Padova, via Marzolo
8, 35131
Padova,
Italy
^{37}
Dipartimento di Fisica e Scienze della Terra, Università di
Ferrara, via Saragat
1, 44122
Ferrara,
Italy
^{38}
Dipartimento di Fisica, Università La Sapienza, ,
P.le A. Moro 2, 00185
Roma,
Italy
^{39}
Dipartimento di Fisica, Università degli Studi di
Milano, via Celoria,
16, 20133
Milano,
Italy
^{40}
Dipartimento di Fisica, Università degli Studi di
Trieste, via A. Valerio
2, 34127
Trieste,
Italy
^{41}
Dipartimento di Fisica, Università di Roma Tor
Vergata, via della Ricerca Scientifica,
1, 00133
Roma,
Italy
^{42}
Discovery Center, Niels Bohr Institute, Blegdamsvej 17, 2100
Copenhagen,
Denmark
^{43}
Dpto. Astrofísica, Universidad de La Laguna (ULL), ,
38206, La Laguna, Tenerife,
Spain
^{44}
European Southern Observatory, ESO Vitacura, Alonso de Cordova 3107, Vitacura,
Casilla
19001
Santiago,
Chile
^{45}
European Space Agency, ESAC, Planck Science Office, Camino bajo del
Castillo, s/n, Urbanización Villafranca del Castillo, Villanueva de la
Cañada, Madrid,
Spain
^{46}
European Space Agency, ESTEC, Keplerlaan 1, 2201 AZ
Noordwijk, The
Netherlands
^{47}
Finnish Centre for Astronomy with ESO (FINCA), University of
Turku, Väisäläntie
20, 21500
Piikkiö,
Finland
^{48}
Haverford College Astronomy Department, 370 Lancaster Avenue, Haverford, Pennsylvania, USA
^{49}
Helsinki Institute of Physics, Gustaf Hällströmin katu 2, University
of Helsinki, Helsinki,
Finland
^{50}
INAF − Osservatorio Astrofisico di Catania, ,
via S. Sofia 78, Catania, Italy
^{51}
INAF − Osservatorio Astronomico di Padova, ,
Vicolo dell’Osservatorio 5,
Padova,
Italy
^{52}
INAF − Osservatorio Astronomico di Roma, ,
via di Frascati 33, Monte Porzio Catone,
Italy
^{53}
INAF − Osservatorio Astronomico di Trieste, ,
via G.B. Tiepolo 11, Trieste, Italy
^{54}
INAF Istituto di Radioastronomia, via P. Gobetti 101, 40129
Bologna,
Italy
^{55}
INAF/IASF Bologna, via Gobetti 101, Bologna, Italy
^{56}
INAF/IASF Milano, via E. Bassini 15, Milano, Italy
^{57}
INFN, Sezione di Bologna, via Irnerio 46, 40126
Bologna,
Italy
^{58}
INFN, Sezione di Roma 1, Università di Roma Sapienza, ,
P.le Aldo Moro 2, 00185
Roma,
Italy
^{59}
IPAG:Institut de Planétologie et d’Astrophysique de Grenoble,
Université Joseph Fourier, Grenoble 1/CNRSINSU, UMR 5274, 38041
Grenoble,
France
^{60}
IUCAA, Post Bag 4, Ganeshkhind, Pune University
Campus, 411 007
Pune,
India
^{61}
Imperial College London, Astrophysics group, Blackett
Laboratory, Prince Consort
Road, London,
SW7 2AZ,
UK
^{62}
Infrared Processing and Analysis Center, California Institute of
Technology, Pasadena,
CA
91125,
USA
^{63}
Institut Néel, CNRS, Université Joseph Fourier Grenoble
I, 25 rue des
Martyrs, Grenoble,
France
^{64}
Institut Universitaire de France, 103 bd SaintMichel, 75005
Paris,
France
^{65}
Institut d’Astrophysique Spatiale, CNRS (UMR 8617) Université
ParisSud 11, Bâtiment
121, Orsay,
France
^{66}
Institut d’Astrophysique de Paris, CNRS (UMR 7095), ,
98bis boulevard Arago, 75014,
Paris,
France
^{67}
Institute for Space Sciences, BucharestMagurale,
Romania
^{68}
Institute of Astronomy and Astrophysics, Academia
Sinica, Taipei,
Taiwan
^{69}
Institute of Astronomy, University of Cambridge, ,
Madingley Road, Cambridge
CB3 0HA,
UK
^{70}
Institute of Theoretical Astrophysics, University of
Oslo, Blindern,
Oslo,
Norway
^{71}
Instituto de Astrofísica de Canarias, C/Vía Láctea s/n, La Laguna, Tenerife, Spain
^{72}
Instituto de Física de Cantabria (CSICUniversidad de
Cantabria), Avda. de los Castros
s/n, Santander,
Spain
^{73}
Jet Propulsion Laboratory, California Institute of
Technology, 4800 Oak Grove
Drive, Pasadena,
California,
USA
^{74}
Jodrell Bank Centre for Astrophysics, Alan Turing Building, School
of Physics and Astronomy, The University of Manchester, Oxford Road, Manchester, M13
9PL, UK
^{75}
Kavli Institute for Cosmology Cambridge, ,
Madingley Road, Cambridge, CB3 0HA, UK
^{76}
LAL, Université ParisSud, CNRS/IN2P3, Orsay, France
^{77}
LERMA, CNRS, Observatoire de Paris, 61 avenue de
l’Observatoire, Paris,
France
^{78}
Laboratoire AIM, IRFU/Service d’Astrophysique − CEA/DSM − CNRS −
Université Paris Diderot, Bât. 709,
CEASaclay, 91191
GifsurYvette Cedex,
France
^{79}
Laboratoire Traitement et Communication de l’Information, CNRS (UMR
5141) and Télécom ParisTech, 46 rue
Barrault, 75634
Paris Cedex 13,
France
^{80}
Laboratoire de Physique Subatomique et de Cosmologie, Université
Joseph Fourier Grenoble I, CNRS/IN2P3, Institut National Polytechnique de
Grenoble, 53 rue des
Martyrs, 38026
Grenoble Cedex,
France
^{81}
Laboratoire de Physique Théorique, Université ParisSud 11 &
CNRS, Bâtiment 210,
91405
Orsay,
France
^{82}
Lawrence Berkeley National Laboratory, Berkeley, California, USA
^{83}
MaxPlanckInstitut für Astrophysik, KarlSchwarzschildStr. 1, 85741
Garching,
Germany
^{84}
McGill Physics, Ernest Rutherford Physics Building, McGill
University, 3600 rue
University, Montréal,
QC, H3A 2T8, Canada
^{85}
MilliLab, VTT Technical Research Centre of Finland, ,
Tietotie 3, Espoo, Finland
^{86}
National University of Ireland, Department of Experimental
Physics, Maynooth,
Co. Kildare,
Ireland
^{87}
Niels Bohr Institute, Blegdamsvej 17, Copenhagen, Denmark
^{88}
Observational Cosmology, Mail Stop 36717, California Institute of
Technology, Pasadena,
CA
91125,
USA
^{89}
Optical Science Laboratory, University College London, ,
Gower Street, London, UK
^{90}
SBITPLPPC, EPFL, 1015
Lausanne,
Switzerland
^{91}
SISSA, Astrophysics Sector, via Bonomea 265, 34136
Trieste,
Italy
^{92}
School of Physics and Astronomy, Cardiff University, ,
Queens Buildings, The Parade,
Cardiff, CF24 3AA, UK
^{93}
School of Physics and Astronomy, University of
Nottingham, Nottingham
NG7 2RD,
UK
^{94}
Space Research Institute (IKI), Russian Academy of
Sciences, Profsoyuznaya Str,
84/32, 117997
Moscow,
Russia
^{95}
Space Sciences Laboratory, University of California, ,
Berkeley, California, USA
^{96}
Special Astrophysical Observatory, Russian Academy of
Sciences, Nizhnij Arkhyz,
Zelenchukskiy region, 369167
KarachaiCherkessian Republic,
Russia
^{97}
Stanford University, Dept of Physics, Varian Physics Bldg, 382 via Pueblo
Mall, Stanford,
California,
USA
^{98}
SubDepartment of Astrophysics, University of Oxford, ,
Keble Road, Oxford
OX1 3RH,
UK
^{99}
Theory Division, PHTH, CERN, 1211
Geneva 23,
Switzerland
^{100}
UPMC Univ Paris 06, UMR7095, 98bis boulevard Arago, 75014
Paris,
France
^{101}
Université de Toulouse, UPSOMP, IRAP, 31028
Toulouse Cedex 4,
France
^{102}
University of Granada, Departamento de Física Teórica y del Cosmos,
Facultad de Ciencias, Granada, Spain
^{103}
Warsaw University Observatory, Aleje Ujazdowskie 4, 00478
Warszawa,
Poland
Received: 21 March 2013
Accepted: 16 March 2014
On the arcminute angular scales probed by Planck, the cosmic microwave background (CMB) anisotropies are gently perturbed by gravitational lensing. Here we present a detailed study of this effect, detecting lensing independently in the 100, 143, and 217 GHz frequency bands with an overall significance of greater than 25σ. We use thetemperaturegradient correlations induced by lensing to reconstruct a (noisy) map of the CMB lensing potential, which provides an integrated measure of the mass distribution back to the CMB lastscattering surface. Our lensing potential map is significantly correlated with other tracers of mass, a fact which we demonstrate using several representative tracers of largescale structure. We estimate the power spectrum of the lensing potential, finding generally good agreement with expectations from the bestfitting ΛCDM model for the Planck temperature power spectrum, showing that this measurement at z = 1100 correctly predicts the properties of the lowerredshift, latertime structures which source the lensing potential. When combined with the temperature power spectrum, our measurement provides degeneracybreaking power for parameter constraints; it improves CMBalone constraints on curvature by a factor of two and also partly breaks the degeneracy between the amplitude of the primordial perturbation power spectrum and the optical depth to reionization, allowing a measurement of the optical depth to reionization which is independent of largescale polarization data. Discarding scale information, our measurement corresponds to a 4% constraint on the amplitude of the lensing potential power spectrum, or a 2% constraint on the rootmeansquared amplitude of matter fluctuations at z ~ 2.
Key words: gravitational lensing: weak / methods: data analysis / cosmic background radiation / largescale structure of Universe
© ESO, 2014
1. Introduction
This paper, one of a set of papers associated with the 2013 release of data from the Planck^{1} mission (Planck Collaboration I 2014), describes our reconstruction of the cosmic microwave background (CMB) lensing potential based on 15 months of data, estimation of the lensing potential power spectrum, and a first set of associated science results.
When Blanchard and Schneider first considered the effect of gravitational lensing on the CMB anisotropies in 1987, they wrote with guarded optimism that although “such an observation is far from present possibilities [...] such an effect will not be impossible to find and to identify in the future” (Blanchard & Schneider 1987). In the proceeding years, and with the emergence of the concordance ΛCDM cosmology, a standard theoretical picture has emerged, in which the largescale, linear structures of the Universe which intercede between ourselves and the CMB lastscattering surface induce small but coherent (Cole & Efstathiou 1989) deflections of the observed CMB temperature and polarisation anisotropies, with a typical magnitude of 2′. These deflections blur the acoustic peaks (Seljak 1996), generate smallscale power (Linder 1990; Metcalf & Silk 1997), nonGaussianity (Bernardeau 1997), and convert a portion of the dominant Emode polarisation to Bmode (Zaldarriaga & Seljak 1998). Gravitational lensing of the CMB is both a nuisance, in that it obscures the primordial fluctuations (Knox & Song 2002), as well as a potentially useful source of information; the characteristic signatures of lensing provide a measure of the distribution of mass in the Universe at intermediate redshifts (typically 0.1 <z< 5). In the ΛCDM framework, there exist accurate methods to calculate the effects of lensing on the CMB power spectra (Challinor & Lewis 2005), as well as optimal estimators for the distinct statistical signatures of lensing (Hu & Okamoto 2002; Hirata & Seljak 2003a).
In recent years there have been a number of increasingly sensitive experimental measurements of CMB lensing. Lensing has been measured in the data of the WMAP satellite both in crosscorrelation with largescalestructure probed by galaxy surveys (Hirata et al. 2004, 2008; Smith et al. 2007; Feng et al. 2012a), as well as internally at lower signaltonoise (Smidt et al. 2011; Feng et al. 2012b). The current generation of lownoise, highresolution groundbased experiments has done even better; the Atacama Cosmology Telescope (ACT) has provided an internal detection of lensing at 4.6σ (Das et al. 2011, 2013), and the South Pole Telescope detects lensing at 8.1σ in the temperature power spectrum, and 6.3σ from a direct reconstruction of the lensing potential (Keisler et al. 2011; van Engelen et al. 2012; Story et al. 2013). Significant measurements of the correlation between the reconstructed lensing potential and other tracers of largescale structure have also been observed (Bleem et al. 2012; Sherwin et al. 2012).
Planck enters this field with unique fullsky, multifrequency coverage. Nominal map noise levels for the first data release (approximately 105, 45, and 60 μK arcmin for the three CMB channels at 100, 143, and 217 GHz respectively) are approximately five times lower than those of WMAP (or twenty five times lower in power), and the Planck beams (approximately 10′, 7′ and 5′ at 100, 143, and 217 GHz), are small enough to probe the deflections typical of lensing. Full sky coverage is particularly beneficial for the statistical analysis of lensing effects, as much of the “noise” in temperature lens reconstruction comes from CMB fluctuations themselves, which can only be beaten down by averaging over many modes.
Lensing performs a remapping of the CMB fluctuations, such that the observed temperature anisotropy in direction is given in terms of the unlensed, “primordial” temperature anisotropy as (e.g. Lewis & Challinor 2006) (1)where is the CMB lensing potential, defined by (2)Here χ is conformal distance (with χ_{∗} ≈ 14 000 Mpc denoting the distance to the CMB lastscattering surface) and is the (Weyl) gravitational potential at conformal distance χ along the direction at conformal time η (the conformal time today is denoted as η_{0}). The angulardiameter distance f_{K}(χ) depends on the curvature of the Universe, and is given by (3)The lensing potential is a measure of the integrated mass distribution back to the lastscattering surface. To first order, its effect on the CMB is to introduce a correlation between the lensed temperature and the gradient of the unlensed temperature, a property which can be exploited to make a (noisy) reconstruction of the lensing potential itself.
In Fig. 1 we plot the noise power spectrum for reconstruction of the lensing potential using the three Planck frequencies which are most sensitive to the CMB anisotropies on the arcminute angular scales at which lensing effects become apparent.
Fig. 1
Skyaveraged lens reconstruction noise levels for the 100, 143, and 217 GHz Planck channels (red, green, and blue solid, respectively), as well as for experiments that are cosmicvariance limited to a maximum multipole ℓ_{max} = 1000, 1500, and 1750 (upper to lower solid grey lines). A fiducial ΛCDM lensing potential using bestfit parameters to the temperature power spectrum from Planck Collaboration XVI (2014) is shown in dashed black. The noise level for a minimumvariance (“MV”) combination of 143 + 217 GHz is shown in black (the gain from adding 100 GHz is negligible). 
Fig. 2
Overview of forecasted contributions to the detection significance as a function of lensing multipole L for the power spectrum (solid black), as well as for several other mass tracers, at the noise levels of our MV lens reconstruction. Our measurement of the power spectrum is presented in Sect. 6, The ISWφ correlation believed to be induced by dark energy is studied in Sect. 6.2. The NVSSφ correlation is studied (along with other galaxy correlations) in Sect. 6.3. The CIBφ prediction (dashed cyan) uses the linear SSED model of Hall et al. (2010), assuming no noise or foreground contamination. A full analysis and interpretation of the CIBφ correlation is performed in Planck Collaboration XVIII (2014). 
The angular size of the Planck beams (5′ FWHM and greater) does not allow a high signaltonoise ratio (S/N) reconstruction of the lensing potential for any individual mode (our highest S/N on an individual mode is approximately 2/3 for the 143 and 217 GHz channels, or 3/4 for a minimumvariance combination of both channels), however with fullsky coverage the large number of modes that are probed provides considerable statistical power. To provide a feeling for the statistical weight of different regions of the lensing measurement, in Fig. 2 we plot (forecasted) contributions to the total detection significance for the potential power spectrum as a function of lensing multipole L. In addition to the power spectrum of the lensing potential, there is tremendous statistical power in crosscorrelation of the Planck lensing potential with other tracers of the matter distribution. In Fig. 2 we also plot forecasted S/N contributions for several representative tracers.
This paper describes the production, characterization, and first science results for two Planckderived lensing products:

1.
A map of the CMB lensing potential over a large fraction of the sky (approximately 70%). This represents an integrated measure of mass in the entire visible Universe, with a peak sensitivity to redshifts of z ~ 2. At the resolution of Planck, this map provides an estimate of the lensing potential down to angular scales of 5′ at L = 2048, corresponding to structures on the order of 3 Mpc in size at z = 2.

2.
An estimate of the lensing potential power spectrum and an associated likelihood, which is used in the cosmological parameter analysis of Planck Collaboration XVI (2014). Our likelihood is based on the lensing multipole range 40 ≤ L ≤ 400. This multipole range (highlighted as a dark grey band in Fig. 2), was chosen as the range in which Planck has the greatest sensitivity to lensing power, encapsulating over 90% of the anticipated signaltonoise, while conservatively avoiding the lowL multipoles where meanfield corrections due to survey anisotropy (discussed in Appendix C) are large, and the highL multipoles where there are large corrections to the power spectra from Gaussian (disconnected) noise bias. Distilled to a single amplitude, our likelihood corresponds to a 4% measurement of the amplitude of the fiducial ΛCDM lensing power spectrum, or a 2% measurement of the amplitude of the matter fluctuations (neglecting parameter degeneracies).
Our efforts to validate these products are aided by the frequency coverage of the three Planck channels that we employ, which span a wide range of foreground, beam, and noise properties. For the mask levels that we use, the rootmeansquared (RMS) foreground contamination predicted by the Planck sky model (Delabrouille et al. 2013) has an amplitude of 14, 22, and 70 μK at 100, 143, and 217 GHz, which can be compared to a CMB RMS for the Planck bestfitting ΛCDM power spectrum of approximately 110 μK. The dominant foreground component at all three CMB frequencies is dust emission, both from our Galaxy as well as the cosmic infrared background (CIB), although at 100 GHz freefree emission is thought to constitute approximately 15% of the foreground RMS. Contamination from the thermal SunyaevZeldovich (tSZ) effect is a potential worry at 100 and 143 GHz, but negligible at 217 GHz (Sunyaev & Zeldovich 1980). On the instrumental side, these frequency channels also span a wide range of beam asymmetry, with typical ellipticities of 19%, 4%, and 18% at 100, 143, and 217 GHz. The magnitude of correlated noise on small scales (due to deconvolution of the bolometer time response) also varies significantly. The ratio of the noise power (before beam deconvolution) at ℓ = 1500 to that at ℓ = 500 is a factor of 1.5, 1.1, and 1.0 at 100, 143, and 217 GHz. The agreement of lens reconstructions based on combinations of these three channels allows a powerful suite of consistency tests for both foreground and instrumental biases. We will further validate the robustness of our result to foreground contamination using the componentseparated maps from the Planck consortium (Planck Collaboration XII 2014).
At face value, the 4% measurement of in our fiducial likelihood corresponds to a 25σ detection of gravitational lensing effects. In fact, a significant fraction (approximately 25% of our error bar) is due to sample variance of the lenses themselves, and so the actual “detection” of lensing effects (under the null hypothesis of no lensing) is significantly higher. We have also been conservative in terms of mask and multipole range in the construction of our fiducial lensing likelihood. As we will show in Sect. 7.1, we obtain consistent results on sky fractions larger than our fiducial 70% sky mask.
The Planck lensing potential is part of a significant shift for CMB lensing science from the detection regime to that of precision cosmological probe. The NVSS quasar catalogue, for example, has been a focus of previous lensing crosscorrelation studies with WMAP (Hirata et al. 2004; Smith et al. 2007; Hirata et al. 2008), where evidence for crosscorrelation was found at approximately 3.5σ. As we will see in Sect. 6.3, the significance for this correlation with Planck is now 20σ. Notably, this is less than the significance with which lensing may be detected internally with Planck. The lensing potential measured by Planck now has sufficient signaltonoise that shot noise of the NVSS quasar catalogue is the limiting source of noise in the crosscorrelation.
The majority of this paper is dedicated to the production and testing of the Planck lensing map and power spectrum estimate. Our focus here is on extracting the nonGaussian signatures of lensing, although we note that lensing effects are also apparent at high significance (10σ) as a smoothing effect in the Planck temperature power spectra (Planck Collaboration XV 2014). We begin in Sect. 2, where we describe and motivate our methodology for producing unbiased estimates of the lensing potential and its power spectrum. The Planck maps and data cuts that are used for this purpose are described in Sect. 3, and the simulations that we use to characterize our reconstruction and its uncertainties are described in Sect. 4. In Sect. 5 we give an overview of our error budget, and discuss the various sources of systematic and statistical uncertainty for our lensing estimates. In Sect. 6 we present our main results: the first Planck lensing map and a corresponding estimate of the lensing potential power spectrum. The likelihood based on this power spectrum is combined with the likelihood for the temperature anisotropy power spectrum (Planck Collaboration XV 2014) to derive parameter constraints in Planck Collaboration XVI (2014). In Sect. 6.1 we highlight a subset of parameter results where the information provided by the lensing likelihood has proven particularly useful. In the concordance ΛCDM cosmology, there is believed to be a correlation between the CMB lensing potential and the lowℓ temperature anisotropies, driven by the effects of dark energy. We also present a measurement of this correlation in Sect. 6.2. Finally, we connect our lensing potential map to other tracers of largescale structure with several illustrative crosscorrelations using galaxy, quasar, cluster and infrared source catalogues in Sect. 6.3. These main results are followed in Sect. 7 by a large suite of systematic and consistency tests, where we perform null tests against a variety of different data cuts and processing. We conclude in Sect. 8. A series of appendices provide further details on some technical aspects of our methodology and lensing potential estimates.
Throughout this paper, when we refer to the concordance or fiducial ΛCDM cosmology we are referring to a model with baryon density ω_{b} = Ω_{b}h^{2} = 0.0221, cold dark matter density ω_{c} = Ω_{c}h^{2} = 0.1199, density parameter for the cosmological constant Ω_{Λ} = 0.6910, Hubble parameter H_{0} = 100 h km s^{1} Mpc^{1} with h = 0.6778, spectral index of the power spectrum of the primordial curvature perturbation n_{s} = 0.96, amplitude of the primordial power spectrum (at k = 0.05 Mpc^{1}) A_{s} = 2.21 × 10^{9}, and Thomson optical depth through reionization τ = 0.093. These values were determined from a prepublication analysis of the Planck temperature power spectrum, but are consistent with the bestfit values quoted in Planck Collaboration XVI (2014). Throughout this work we shall frequently quote lensing bandpower amplitudes relative to this fiducial model for ease of comparison, although as discussed in Sect. 2 our lensing likelihood itself is designed to be insensitive to the choice of fiducial model.
2. Methodology
In this section, we detail our methodology for reconstructing the lensing potential and estimating its angular power spectrum. These are both obtained by exploiting the distinctive statistical properties of the lensed CMB.

1.
If we consider a fixed lensing potential applied to multiple realizations of the CMB temperature anisotropies, then lensing introduces statistical anisotropy into the observed CMB; the fluctuations are still Gaussian, however the covariance varies as a function of position and orientation on the sky. We use this idea to obtain a (noisy) estimate of . The noise of this map is a combination of instrumental noise and statistical noise due to the fact that we only have a single realization of the CMB to observe, analogous to shape noise in galaxy lensing.

2.
If we consider averaging over realizations of both the lensing potential and the CMB fluctuations, then lensing introduces nonGaussianity into the observed CMB. This appears at lowest order in the connected part of the CMB 4point function, or trispectrum^{2}. We use this to measure the lensing power spectrum .
The estimators that we use are derived from maximizing the likelihood function of the lensed CMB under the assumption that the instrumental noise is Gaussian and the lensed CMB is perturbatively Gaussian and statistically isotropic. These estimators are optimal (in the minimumvariance sense). In cases where we have made suboptimal analysis choices, we provide estimates of the loss of signaltonoise.
2.1. Lens reconstruction
To gain intuition for the process of lens reconstruction, it is useful to consider the effect of lensing on a small patch of the sky. Lensing remaps the temperature fluctuations by a deflection field . The part of that is constant over our patch is not an observable effect; it describes only a recentering of the map. The variation of the deflection field across the patch is observable, however. This can be usefully decomposed into convergence (κ) and shear modes (γ_{+}, γ_{−}) as (4)If we observe a patch that is small enough that these quantities can be taken as constant, then the observational consequences are simple. The convergence mode causes a local change of scale, either magnifying or demagnifying the fluctuations. Taking the local power spectrum of our small patch, we would find that the CMB peaks would shift to larger or smaller scales, relative to the fullsky average. The shear modes also describe changes of scale, however they are now orientation dependent. On a small patch, convergence and shear estimators can be constructed from local estimates of the (orientationdependent) power spectrum and then stitched together to recover the lensing potential φ (Zaldarriaga & Seljak 1998; Bucher et al. 2012). This procedure describes a quadratic estimator for the local convergence and shear.
From the description above, it is not immediately clear how to go about stitching together estimates of convergence and shear in different regions of the sky, or what weight to give the local power spectrum estimates as a function of scale. These questions can be resolved by considering a generic form for the quadratic estimator, and optimizing its weight function for sensitivity to lensing (Okamoto & Hu 2003). To first order in the lensing potential, the statistical anisotropy introduced by lensing appears as an offdiagonal contribution to the covariance matrix of the CMB: (5)where the average ⟨⟩ is taken over CMB realizations with a fixed lensing potential. Here the bracketed term is a Wigner 3j symbol, is the harmonic transform of the lensing potential, and the weight function is given by (6)Here is the ensembleaverage power spectrum of the lensed CMB. In our analysis, we will use the fiducial model described at the end of Sect. 1 to determine , however our lensing likelihood can be renormalized to account for uncertainties in this model. In this approach, errors in the fiducial model do not bias our lensing bandpower estimates, they only result in slight suboptimality. Note that we use the lensed power spectrum here, rather than the unlensed spectrum that is sometimes used in the literature, as this is accurate to higher order in φ (Lewis et al. 2011), an improvement which is necessary at Planck sensitivity (Hanson et al. 2011). Use of the unlensed spectrum would lead to biases on the order of 15% at L< 200.
Now we construct a quadratic estimator to search for the covariance which is introduced by lensing. We will use several different estimators for the lensing potential, as well as to probe possible pointsource contamination, and so it will be useful to keep this discussion as general as possible. A completely generic quadratic estimator for the lensing potential can be written as (7)where ℛ^{xφ} is a normalization matrix, and is a quadratic “building block” which takes in a pair of filtered sky maps and , and sums over their empirical covariance matrix with a weight function : (8)The “meanfield” term accounts for all known sources of statistical anisotropy in the map, which could otherwise bias the lensing estimate. It is given by (9)where the ensemble average here is taken over realizations of the CMB and noise.
We may now optimize the generic quadratic estimator above. If the primordial CMB fluctuations and instrumental noise are Gaussian and the lensing potential is fixed, then the likelihood for the observed CMB fluctuations is still a Gaussian, which may be maximized with respect to the lensing potential modes φ_{LM} (Hirata & Seljak 2003a). The optimal quadratic estimator is the first step of an iterative maximization of this likelihood, and it has been shown that additional iterations of the estimator are not necessary for temperature lens reconstruction (Hirata & Seljak 2003a; Okamoto & Hu 2003). The optimal quadratic estimator has the following choices for the weight function and filtering.

1.
The weight function W^{x} should be a matched filter for the covariance induced by lensing (i.e., one should use , with weight function given by Eq. (6)). We shall use this weight function for all of our fiducial results, although for consistency tests we will also use “biashardened” estimators, which have weight functions constructed to be orthogonal to certain systematic effects (Namikawa et al. 2013). This is discussed further in Sect. 7.4.

2.
The filtered temperature multipoles should be given by , where T is a beamdeconvolved sky map and C is its total signal+noise covariance matrix. We describe our approximate implementation of this filtering in Appendix B. When combining multiple frequencies for our minimumvariance estimator, all of the available data are combined into a single map which is then filtered and used for both input multipoles of the quadratic estimator. It can be desirable to use different pairs of maps however, and we use this for several consistency tests. For example, we feed maps with independent noise realizations into the quadratic estimator to avoid possible noise biases in Sect. 7.3.
In the quadratic maximumlikelihood estimator, the meanfield correction emerges from the determinant term in the likelihood function, and it can be seen that the normalization matrix ℛ is the Fisher matrix for the φ_{LM}; this means that the normalization is the same as the covariance matrix of the lens reconstruction, and so the unnormalized lensing estimate is equivalent to an inversevarianceweighted lens reconstruction, which is precisely the quantity needed for most statistical analysis. This is why we have denoted it with an overbar, in analogy to .
We choose to treat the map noise as if it were homogeneous when constructing the filtered , and do not account for variation with hit count across the sky. This is a slightly suboptimal filtering choice; in Appendix. B we estimate that it leads to a 5% loss of total signaltonoise when constraining the power spectrum of the lensing potential. The advantage of this approach, however, is that far from the mask boundaries our filtering asymptotes to a simple form, given by (10)where is the temperature power spectrum and is the power spectrum of the homogeneous noise level that we use in our filtering. For the purposes of compact notation, in the following equations we combine both of these elements in the “filter function” F_{ℓ}. The asymptotic form of our filtering, Eq. (10), will prove useful, as it means that the normalization of our estimator, as well as its variance and response to various systematic effects, may be accurately modelled analytically. It allows us to propagate uncertainties in the beam transfer function and CMB power spectrum, for example, directly to our lens reconstruction. This filtering choice also means that the normalization does not vary as a function of position on the sky, which simplifies the analysis of crosscorrelations between the lensing potential map and external tracers. Under the approximation of Eq. (10), the estimator normalization is given by (11)where the response function for filtered maps and is (12)This can be read as “the response of estimator x to lensing on scale L”. The filter functions F_{ℓ} are those used for and respectively. In cases where the filter functions are obvious, we will drop the indices above.
Putting all of the above together, for a chosen quadratic estimator we obtain normalized, meanfielddebiased estimates of the lensing potential φ as (13)Note that our normalization function is only approximate, but we will verify its accuracy in Sect. 4. For the standard lensing estimator of Okamoto & Hu (2003; which uses the weight function of Eq. (6)), we use x = φ. This estimator is denoted simply as .
2.2. Lensing power spectrum estimation
We form estimates for the power spectrum of the lensing potential by taking spectra of the lensing estimates from Sect. 2.1, using a simple pseudoC_{ℓ} estimator. In order to reduce mode coupling, as well as to downweight regions near the analysis boundary where the meanfield due to masking can be large, we take the power spectrum from an apodized version of our lensing estimate, given by (14)where is an apodized version of the analysis mask used in our filtering and P_{L} ≡ L(L + 1) is an approximate prewhitening operation. The construction of is described in Sect. 3. Our fiducial apodization occurs over a band of approximately 5deg, and effectively reduces the sky fraction by 9%.
The power spectrum of probes the 4point function of the observed CMB, which contains both disconnected and connected parts. We model it as being due to a combination of Gaussian CMB fluctuations, lensing effects and unresolved pointsource shot noise, and estimate the power spectrum of the lensing potential with (15)where is the average value of the square of the apodizing mask. The first line of Eq. (15)isolates the connected part of the CMB 4point function, or trispectrum, which would be zero for Gaussian fluctuations. The second line contains corrections which isolate the part of the trispectrum which is directly proportional to the nonGaussianity induced by . In the following paragraphs, we explain these terms in more detail.
The first correction term subtracts the (large) disconnected contribution to the power spectrum of . To determine this term, we use the datadependent subtraction which emerges for maximumlikelihood estimators of the CMB trispectrum (Regan et al. 2010; see also Appendix D). For lensing, this procedure has the additional advantage of reducing the correlation between different multipoles L ≠ L′ of the lens reconstruction to a level which is negligible at Planck resolution and noise levels (Schmittfull et al. 2013), as well as reducing sensitivity to uncertainties in our model of the CMB and noise covariance matrices (Namikawa et al. 2013). Writing the power spectrum of explicitly as a function of the four inversevariance filtered temperature maps (16)the disconnected contribution reads (17)where indicates a MonteCarlo simulation of the corresponding map. The ensemble average is taken over two sets of independent realizations mc and mc′. Note that because of the way we have used pairs of MonteCarlo simulations and data with independent CMB and noise realizations, the meanfield correction is zero for all of the terms above.
The term corrects for the “N^{(1)}” bias due to secondary contractions of the lensing trispectrum (Hu 2001; Kesden et al. 2003). It is only a large effect at L> 100, and so we calculate it using a flatsky expression as (18)where l_{1} + l_{2} = l_{3} + l_{4} = L and is a fiducial model for the lensing potential power spectrum. The W(l,l′) are flatsky analogues of the fullsky weight functions. The flatsky lensing weight function, for example, is (19)The N^{(1)} term is proportional to the lensing potential power spectrum, and so in principle it should be used to improve our constraints on rather than subtracted as an additive bias. However the statistical power of this term is relatively small at Planck noise levels. From a Fisher matrix calculation, the trispectrum contractions which source the N^{(1)} term are only detectable in the Planck data at 4σ significance, compared to approximately 25σ for the primary contractions. We choose simply to subtract the N^{(1)} term from our power spectrum estimates. There is a small cosmological uncertainty in the N^{(1)} correction due to uncertainty in the power spectrum, which we discuss in Sect. 5.3.
The term is a correction for the bias induced by the shotnoise of unresolved point sources (including SZ clusters), and will be discussed in more detail in Sect. 2.4.
Finally, the term is a small correction that we obtain by estimating following the procedure above on a number of lensed CMB realizations, and then subtracting the input power spectrum. This term can be nonzero due to pseudoC_{ℓ} leakage effects from masking, which we have not accounted for other than apodization, errors in our calculation of the N^{(1)} term, or errors in the normalization at the power spectrum level. We will find that is sufficiently small that in practice it does not matter whether we account for it as a renormalization or an additive offset, and we choose to treat it as an offset for simplicity.
We will ultimately characterize the uncertainty of by MonteCarlo, however the following analytical expression is a useful approximation (20)where is the reconstruction noise level. We take it to be (21)where is the ensembleaverage crossspectrum between and , given by (22)where f_{sky,(i)(j)} is the fraction of sky common to both and . We will use this analytical estimate of the variance for weighting our estimates in statistical analyses.
2.3. Lensing power spectrum likelihood
Based on our measurements of the lensing potential power spectrum, we construct a Gaussian likelihood based on bins in . Our likelihood has the form (23)where ℬ represents a binning function, Σ is the covariance matrix between bins, and sums are performed over paired upper/lower indices. As the shape and amplitude of are constrained strongly in the concordance ΛCDM model, with which our results are broadly consistent, we choose a binning function designed to maximize our sensitivity to small departures from the fiducial ΛCDM expectation. This is given by (24)where V_{L} is defined in Eq. (20)and and define the multipole range of the bin. Our binned results correspond to an estimate of the amplitude of a fiducial power spectrum in a given multipole range, normalized to unity for the fiducial model. We denote these amplitude estimates explicitly for a given lens reconstruction as (25)In principle there are several reasons why a likelihood approach such as the one we have described above could fail; the usual issues with a pseudoC_{ℓ} likelihood, such as the nonGaussianity of the are compounded by the fact that the φ estimates themselves are nonGaussian, and derived from the temperature data itself, which means that the measurement uncertainties on the lens reconstruction are potentially correlated with those of the temperature power spectrum. In Appendix D, we validate the above approach to the likelihood by considering these issues in more detail.
2.4. Unresolved pointsource correction
At the high multipoles (ℓ > 1000) that provide most of the modes which are useful for lens reconstruction, the contribution from unresolved (and therefore undetected and unmasked) extragalactic foregrounds becomes apparent in the Planck power spectra. For measurements of the lensing potential power spectrum, we should therefore be concerned with possible biases from the trispectrum nonGaussianity of these unresolved sources. The work of Osborne et al. (2014) indicates that at Planck resolution and sensitivity, the contribution of unresolved foregrounds to lens reconstruction is small but potentially nonnegligible, with the largest contribution coming from the shot noise trispectrum of unresolved radio and tSZ sources. Given uncertainties in modelling the nonGaussianity of unresolved point sources, our approach is heavily data dependent. In addition to tests for the consistency of lens reconstructions at 143 and 217 GHz, we use two additional methods to measure and correct for pointsource contamination in our analysis.
In our fiducial lensing power spectrum analysis, we measure the amplitude of the shotnoise contribution to the data trispectrum and correct the measured lensing spectrum accordingly. We measure the shotnoise amplitude using the power spectrum of a quadratic estimator designed to detect the “noise” due to sources, similar to the approach advocated in Munshi et al. (2011a). The trispectrum (or connected fourpoint function) for pointsource shot noise is defined in position space as (26)or in harmonic space: (27)The quadratic estimator that we use to measure the amplitude S^{4} is denoted as , and is defined by the weight function (28)We measure the amplitude of the shotnoise trispectrum S^{4} as (29)where is defined analogously to Eq. (15)and the sums are taken over . The power spectrum weighting S_{L} is given by (30)In the limit that this estimator is equivalent to measuring the shotnoise trispectrum using the 1point kurtosis of the map. The advantage of the trispectrumrelated power spectrum approach is that it allows us to separate out regions of the trispectrum measurement that are contaminated by lensing, as well as to look for evidence of clustering, which would appear as a deviation of the measured spectrum from the shape expected for unclustered sources. We then calculate and remove an estimated bias to the measured lensing potential given by (31)The shotnoise correction described above does not take into account the correlation of sources with the dark matter distribution. And hence the lensing potential. In Sect. 7.4 we therefore additionally perform tests using a “point source biashardened estimator”, constructed using the weight function of Eq. (28). This biashardened estimator has zero response to both the pointsource shotnoise trispectrum (S^{4}), as well as to the primary trispectrum contractions involving the correlation between point sources and the lensing potential (S^{2}φ).
3. Data and cuts
Planck sky maps:
The majority of the results in this paper are based on the Planck nominalmission frequency maps at 100, 143, and 217 GHz, built from the first 15.5 months of data. These are in HEALPix^{3} format, with resolution parameter N_{side} = 2048, corresponding to pixels with a typical width of . These have effective noise levels of approximately 105 μK arcmin at 100 GHz, 45 μK arcmin at 143 GHz, and 60 μK arcmin at 216 GHz. The beam widths shrink with frequency, and are 10′ at 100 GHz, 7′ at 143 GHz, and 5′ at 217 GHz. Our primary products – a map of the lensing potential and an associated power spectrum likelihood – are based on a minimumvariance (hereafter MV) combination of the 143 and 217 GHz maps. Although lensing may be detected at a significance of approximately 10σ at 100 GHz, the CMB modes used in this reconstruction are already samplevariance dominated at 143 and 217 GHz, and so adding it to the MV combination would lead to negligible improvement. In addition to these maps, we also use the 857 GHz Planck map as a dust template, which is projected out of both maps independently in our filtering procedure (described further in Appendix B). This projection is primarily intended to remove diffuse Galactic dust contamination, although it also removes a portion of the CIB fluctuations, which have a similar spectral index to that of Galactic dust over these frequency bands (Planck Collaboration XXIV 2011; Planck Collaboration XVIII 2011). As a simple approach, this template projection overlooks several potential difficulties, including variation of Galactic dust spectral indices across the sky as well as mismatch between the beams at 100/143/217 and 857 GHz, although we find it is adequate for our purposes. As lens reconstruction is most sensitive to smallscale modes, the coupling to largescale Galactic foregrounds is relatively weak. In Sect. 7 we will also perform lens reconstruction using the more rigorously componentseparated maps of Planck Collaboration XII (2014), finding good agreement with our baseline results.
In analyzing these maps, we use the fiducial beam transfer functions described in Planck Collaboration VII (2014) and Planck Collaboration (2013). There are uncertainties associated with these transfer functions, which we propagate to an uncertainty in the lensing estimator normalization in Sect. 5.2.
Galaxy mask:
We avoid the majority of Galactic foreground power using the temperature analysis masks described in Planck Collaboration (2013). These are constructed using a combination of the 30 GHz and 353 GHz maps, corrected for an estimate of the CMB contribution, smoothed to 5′ and thresholded until a desired sky fraction is obtained (Planck Collaboration XII 2014). For our baseline results, we use the 70% masks (which remove 30% of the sky), although in Sect. 7 we will show that we obtain consistent results with both larger and smaller masks. When computing the power spectrum, we additionally multiply our φ estimates by an apodized version of the Galaxy mask. Each pixel outside of the masked region is multiplied by an apodized weight varying between zero at the mask boundary and unity at a distance greater than 5deg from the closest masked pixel. We use a sinusoid weight function, similar to the one used in Namikawa et al. (2013) and BenoitLevy et al. (2013). Note that due to this apodization, the effective sky fraction used in our power spectrum analysis is approximately 9% lower than the sky fraction of our reconstruction.
CO and extendedobject masks:
We mask regions believed to be contaminated by carbonmonoxide (CO) emission lines at 100 and 217 GHz using the type 2 CO map described in Planck Collaboration XIII (2014). We mask all pixels above 60 μK in the map after degrading to N_{side} = 256, and then restore by hand all isolated pixels that are removed by this cut but do not appear to be extended CO regions. This lowresolution mask is prograded back to N_{side} = 2048 and smoothed at 20′, slightly larger than the map resolution of 15′. We additionally exclude extended nearby objects (the two Magellanic Clouds, M31, M33, and M81) by cutting out disks centred on their locations. The radii of these disks range from 250′ for the Large Magellanic Cloud to 30′ for M81.
Pointsource masks:
We remove detected point sources (or otherwise compact objects) using a mask constructed from a combination of sources identified in the Planck Early Release Compact Source Catalogue (ERCSC; Planck Collaboration Int. VII 2013), the SZ clusters from the Planck Cluster Catalogue (PCC; Planck Collaboration XXIX 2014), and the Planck Catalogue of Compact Sources (PCCS; Planck Collaboration XXVIII 2014). We produce individual masks for 100, 143, and 217 GHz, cutting out detected sources with disks having radii of either 3σ or 5σ depending on their flux level, where σ describes the Gaussian beamwidth of the given channel and is taken to be 4.1′ at 100 GHz, 3.1′ at 143 GHz, and 2.1′ at 217 GHz. SZ clusters are not masked at 217 GHz. For the PCCS we mask all sources with S/N ≥ 5 in a given band, or S/N ≥ 10 in either adjacent frequency band, using a 3σ disk, extending to a 5σ disk for any sources with S/N ≥ 10 in the target band. For the ERCSC, we make a 3σ cut for every source detected at the target frequency, as well as for every source that is detected with a flux greater than 1 Jy at 70 GHz, 0.7 Jy at 100 GHz, 0.4 Jy at 143 GHz, 0.5 Jy at 217 GHz, or 0.9 Jy at 353 GHz. For sources that satisfy the brightness criterion at the target frequency, we also enlarge the cut to 5σ. For SZ clusters we make a cut between 3σ and 5σ depending on the cluster angular size θ_{500} for all sources in the MMF1 catalogue with S/N ≥ 5. We further mask cool cores (CC) using the Early Cold Cores catalogue (Planck Collaboration VII 2011) in a similar fashion, using a maximum 15′ excluding radius.
In Sect. 7.2 we will test our sensitivity to this SNR cut, varying it from S/N ≥ 5 to S/N ≥ 4.2 and S/N ≥ 10.
4. Simulations
We require simulations of the data described in Sect. 3 for several aspects of our analysis: to determine the lensingmean field φ^{MF}, to determine the and correction terms of Eq. (15), and to determine our measurement error bars, as well as to validate our reconstruction methodology. Here we outline the basic algorithm we use to generate these simulations.

1.
Simulate the unlensed CMB and lensing potential. Using the fiducial ΛCDM theoretical unlensed power spectra for the cosmology given at end of Sect. 1, we generate realizations and φ_{LM} of the unlensed temperature and lensing potential with appropriate correlations. By default we simulate these fields up to ℓ_{max} = 2560. In addition to these stochastic components, we also include a dipolar lensing mode in our simulations, to account for the Doppler aberration due to our motion with respect to the CMB (Kamionkowski & Knox 2003; Challinor & van Leeuwen 2002): (32)where β is in the direction of Galactic coordinates (l,b) = (264.4,48.4), with amplitude  β  = 0.00123, corresponding to a velocity of 369 km s^{1} (Hinshaw et al. 2009).

2.
Lens the temperature field. We use a fast spherical harmonic transform to compute the temperature and deflection fields (33)This computation is performed on an equicylindrical pixelization (ECP) grid with (N_{θ},N_{φ}) = (16 384,32 768) equallyspaced points in (θ,φ). We determine the lensed temperature on this ECP grid as (34)This equation is evaluated following the same algorithm as the LENSPix code (Lewis 2005), with 2D cubic Lagrange interpolation on the ECP grid, which has conveniently uniform pixel spacing. We then obtain the lensed temperature multipoles by a reverse harmonic transform. The average power spectrum of these lensed simulations agrees with the implementation in CAMB (Challinor & Lewis 2005) to better than 0.1% at l< 2048. We have also used the algorithm of Basak et al. (2009), which gives consistent results. In addition to the lensing dipole, our motion with respect to the CMB also causes a small modulation of the observed fluctuations, which we determine from the lensed temperature field as (35)where f_{ν} is a frequency dependent boost factor that we approximate as 1.5, 2, and 3 at 100, 143, and 217 GHz, respectively. There is an extended discussion of this factor in Planck Collaboration XXVII (2014).

3.
Generate a Gaussian realization of extragalactic foreground power . We simulate Gaussian foreground power at 100, 143, and 217 GHz such that the auto and crossspectra of the foreground power at these frequencies has the form (36)where A and B label frequency bands. The template power spectra are based on the modelling of Planck Collaboration XV (2014); Planck Collaboration XVI (2014), and are given by (37)We use the parameters The parameters above were obtained by fitting the auto and crossspectra of our inversevariancefiltered temperature maps to the templates above. Our point source masks differ from those used in Planck Collaboration XVI (2014), and so these numbers cannot be directly compared to the values there. The purpose of including a Gaussian foreground component in our simulations is primarily so that our measurement error bars will include the additional scatter due to the foreground power (we will see in Sect. 5 that foreground power constitutes a few percent of our total error budget). Our lensing power spectrum results are otherwise unaffected by this component of our simulations. Our lensing estimates do include an overall correction for pointsource shotnoise, although we do not include the nonGaussian covariance of the shotnoise in our error budget, given that point sources are a small fraction of the total power at the multipoles probed by Planck. We test this pointsource treatment further in Sect. 7.2.

4.
Convolve T^{φ} + T^{fg} + T^{mod} with the (asymmetric) instrumental beam using Eq. (C.7). This results in a HEALPix map of the beamconvolved CMB , from which we extract beamconvolved multipoles .

5.
Make a final map of the observed sky with the HEALPix N_{side} = 2048 pixelization, with (38)where is the HEALPix N_{side} = 2048 pixel window function, which accounts for the effective smearing of the beam due to the binning of hits into pixels. Although this smearing is an inadequate description of the pixelizationsmearing effect for any individual pixel, it is a good approximation when averaged over all pixels on the sky. We will discuss pixelization effects further in Sect. C.4. The quantity is a simulated realization of the instrumental noise. We use noise realizations from the FFP6 simulation set (Planck Collaboration I 2014). In order to make the power spectrum of our simulations more closely match the power spectrum of the data, we find that it is necessary to include an additional white noise component with an amplitude of 20, 10, and 10 μK arcmin at 100, 143, and 217 GHz respectively. Indeed, the FFP6 simulation are produced on the basis of the halfring estimates of the noise which are known to underestimate the map noise (Planck Collaboration VI 2014).
The simulations described above are used throughout this work both to debias our lensing estimates, as well as to characterize their uncertainty. Here we briefly use them to validate the normalization of our lensing estimates. There are two relevant normalizations: the maplevel normalization which describes how our reconstruction traces the input φ realization, as well as the analogous normalization for the lensing potential power spectrum. Throughout this work, we use the analytical normalization estimates of Eq. (12)to normalize our map estimates. In Fig. 3 we plot the average crosscorrelation between the input and reconstructed φ_{LM} for our minimumvariance reconstruction, divided by a simple f_{sky} factor to account for missing power in the mask. We can see that our analytical maplevel normalization is accurate at the 1% level. This is small enough that we do not further renormalize our φ map. At the power spectrum level, we find similar accuracy of our normalization in the 40 <L< 400 region which forms the basis of our fiducial lensing likelihood. The power spectrum estimates in Fig. 3 are obtained by averaging the first line of Eq. (15)over a set of lensed but otherwise Gaussian simulations. This 1% error in the power spectrum normalization is accounted for in our likelihood by the term. Given the good agreement of the power spectrum used in our simulations with that measured from data, it does not matter whether we incorporate the MonteCarlo correction as a multiplicative or additive term. We have chosen to incorporate it as additive for simplicity. At higher multipoles, there is additional scatter in our power spectrum normalization test (up to 10% for the plotted bins). It is difficult to know from this test whether it is due to actual scatter in our power spectrum normalization (which seems unlikely, given that the maplevel normalization is quite accurate), a small misestimation of the N^{(1)} term (which begins to dominate over above L = 500), or simply pseudoC_{ℓ} leakage issues. In any case, it is still small, and is absorbed into our ΔC_{L}_{MC} term when estimating power spectra at these multipoles. Although we have only presented results for our minimumvariance reconstruction here, those for individual frequency results are similar.
As a visual illustration, and a preview of our data results in Sect. 6, in Fig. 4, we show a simulated lens reconstruction as well as the input φ map. This gives a visual impression of the signaltonoise in our lens reconstruction.
Fig. 3
Validation of our estimator normalization for simulations of the MV reconstruction at the map and power spectrum levels. The map normalization (plotted as ) is tested by taking the crossspectrum of the input φ with the reconstruction averaged over MonteCarlo simulations, divided by an f_{sky} factor to account for missing power in the mask. The power spectrum normalization (plotted as ) is obtained by averaging the first line of Eq. (15)over simulations, and then comparing it to the expected value, which is because our simulations do not contain pointsource nonGaussianity. 
Fig. 4
Simulation of the Wienerfiltered lensing potential estimate for the MV reconstruction (left), and the input φ realization (right; filtered by to be directly comparable to the Wiener estimate). Both maps show the southern Galactic sky in orthographic projection. The lensing reconstruction is noise dominated on all scales, however correlations between the two maps can still be seen visually. 
5. Error budget
In this section, we describe the measurement and systematic error budget for our estimation of the lensing potential power spectrum. This is broken down into three sections; in Sect. 5.1 we describe our measurement (or “statistical”) error bars, which are due to the fact that we have only a single noisy sky with a finite number of modes to observe. In Sect. 5.2 we consider uncertainty in the instrumental beam transfer function, which we will see propagates to a normalization uncertainty for our lensing estimates. Finally, in Sect. 5.3 we discuss the effect of cosmological uncertainty; possible errors in the fiducial model for result in a normalization uncertainty for our lensing estimates, and uncertainties in the fiducial power spectrum lead to uncertainties in the correction. As a guide to the relative size and scale dependence of these terms, in Fig. 5 we summarize the error budget for our fiducial minimumvariance lens reconstruction, based on 143 and 217 GHz. Individual frequency bands, as well as 100 GHz are qualitatively similar.
Fig. 5
Relative contribution of statistical measurement uncertainty (grey band, representing the 1σ uncertainty on for the plotted bins) and the systematic errors we assign due to uncertainty in the N^{(1)} bias correction, pointsource shotnoise correction, and the beam transfer function. The contributions represent error eigenmodes in , and so are completely correlated between L. Note that these contributions are much smaller than our measurement uncertainty, and they have been multiplied by a factor of 10 for clarity. 
5.1. Measurement
Although our measurement uncertainties are ultimately assigned by Monte Carlo, we can use the analytical expression of Eq. (21)to gain intuition for how they are sourced by various components. Our simple model of the sky after masking and dust cleaning is that it consists of three uncorrelated signals: CMB, instrumental noise, and unresolved isotropic foreground power. The noise variance of the lens reconstruction in Eq. (21)involves two power spectra, and so we can think of the noise contribution as the sum of six possible terms involving pairs of the CMB, noise, and foreground power spectra. In Fig. 6 we combine these contributions into three representative contributions to the reconstruction noise: “pure CMB” in which both spectra are due to CMB fluctuations; the “noise” contribution in which either both spectra are those for noise power, or one is noise and one is CMB; and, finally, the “foreground” contribution in which either one or both of the spectra are due to unresolved foreground power. We can see that for most reconstruction multipoles, the pure CMB contribution constitutes the largest part of the reconstruction noise, followed by noise. The unresolved foreground power is a fairly small contribution to our measurement error. Note that the dominant terms for both the “noise” and “foreground” contributions are the ones in which one of the spectra is a CMB fluctuation. For this reason, we will focus less on the use of crossspectra to avoid noise biases than is done for the usual CMB power spectra (Planck Collaboration XV 2014), although we will perform consistency tests using crossspectra of data to avoid noise biases. Note that our realizationdependent method for removing the disconnected noise bias (Eq. (17)) means that the majority of this contribution is estimated directly from the data itself, reducing our sensitivity to uncertainty in the noise and foreground power.
Fig. 6
Relative contributions of CMB, instrumental noise, and foreground power terms discussed in Sect 5.1 to the approximate lens reconstruction variance of Eq. (21), previously plotted in Fig. 1. As discussed in Sect. 5.1, the finite number of CMB modes observed by Planck is the dominant source of variance for the lens reconstruction. 
5.2. Beam transfer function
Errors in the effective beam transfer function appear as an error in the normalization of our lensing estimates. For simplicity here we will describe the case for a single standard quadratic lensing estimator that uses the same map for both of its inputs, although when dealing with combinations of channels for our actual results we account for differences in the beam transfer function and errors between “legs” of the estimator (or pair of estimators, in the case of the lensing potential power spectrum).
We model the effect on from a discrepancy between the “measured” and “true” effective beam transfer function as (39)where the average is over CMB fluctuations for fixed lenses. The statistics of B^{true}/B^{meas} are discussed in Planck Collaboration VII (2014), where it is found that the beam uncertainty may be described by a small number N_{eig} of eigenmodes (typically five), (40)Here, are the uncertainty eigenmodes (indexed by i) for the beam determination, and g_{i} are independent Gaussian random variables with unit variance. The frequencyband eigenmodes are small for 100−217 GHz (no more than 1% at l< 2048), and so we can equate (41)We can also ignore terms of order E^{2} and propagate the independent beam uncertainties directly to a corresponding set of independent eigenmodes for φ. These eigenmodes are given as (42)with the model that (43)
Fig. 7
Propagation of eigenmodes for beam transfer function uncertainty (lefthand panels) to an uncertainty in the lensing normalization (righthand panels). The normalization eigenmodes are flat, with amplitudes given by approximately twice the beam uncertainty at the beam and noise cutoff scale (at ℓ ≈ 1000 for 100 GHz and ℓ ≈ 1500 for 143 and 217 GHz). 
In Fig. 7 we plot the beam uncertainty eigenmodes as well as the corresponding normalization eigenmodes. The normalization eigenmodes are flatter than the beam uncertainty eigenmodes, and do not converge to zero at lowL, reflecting the fact that at all scales the lensing estimate takes most of its weight from modes of ℓ_{1} and ℓ_{2} close to the noise/beam cutoff (at ℓ ≈ 1000 for 100 GHz and ℓ ≈ 1500 for 143 and 217 GHz). As can be seen in Fig. 7, the beam normalization eigenmode is essentially twice the value of the beam uncertainty eigenmode on these scales.
Owing to the small size of the eigenmodes, the propagation to a fractional uncertainty on is accomplished essentially by applying a factor of two to the eigenmodes of Eq. (42). The normalization uncertainty due to beams for the lensing potential power spectrum is therefore approximately 1% at 100 GHz and 0.5% at 143 and 217 GHz. When computing the beam uncertainty for our minimumvariance result, we take the 143 and 217 GHz beam eigenmodes to be uncorrelated.
5.3. Cosmological
Similarly to the case with the beam transfer function of the previous subsection, uncertainty in the cosmological model may also be propagated to an uncertainty in the estimator normalization. There are two aspects in which cosmological uncertainty can enter our analysis. First, when estimating the contribution of Eq. (15)we assume a fiducial , obtained from the parameters of the fiducial ΛCDM cosmology described in Sect. 1. We take a random subset of samples from Planck MCMC chains for the ΛCDM cosmology and calculate the corresponding bias. Compared to our fiducial correction, we see smallscale fluctuations with an amplitude of 0.2% and a typical scale of ΔL = 150, superimposed on a larger overall amplitude scatter with a standard deviation of approximately 6%. This result is reasonable: in the ΛCDM model the PlanckTT measurements have the power to detect the smoothing effects of gravitational lensing on the acoustic peaks at a significance level of 10σ, corresponding to a 10% constraint on the overall amplitude of , consistent with the result from analysis of individual samples above. To account for cosmological uncertainty in the correction, we therefore assign it a 10% overall amplitude uncertainty. This correction is small for the L range considered in our fiducial likelihood, and so this uncertainty is essentially negligible compared to our statistical error bars (although we include it in our error budget nevertheless).
The second point at which cosmological uncertainty enters our analysis is in the normalization of our estimates. In Eq. (12), as well as in our tests of the normalization in Sect. 4, we have assumed that the statisticallyanisotropic covariance induced by lensing is equal to that assumed when weighting the estimator, given by Eq. (6). If the true cosmological power spectrum differs from the fiducial one used in Eq. (6)then we will misestimate φ as (44)where is given by Eq. (6), with the modification that it is evaluated with the “true” ensembleaverage lensed CMB power spectrum rather than the fiducial one . We have chosen to account for this normalization uncertainty coherently with the temperature likelihood. Given a proposed , pair, we renormalize in our likelihood (Eq. (23)) so that it is directly comparable to our measurement.
6. Results
In this section, we present the primary results of this paper; a measurement of the CMB lensing potential over approximately 70% of the sky and its associated power spectrum estimate . Our lensing potential map and power spectrum have several science implications, which we will discuss in proceeding subsections.

1.
The lensing power spectrum probes the matter power spectrum integrated back to the lastscattering surface, and provides additional lowredshift leverage for CMBalone parameter constraints with a broad maximum at z ~ 2. Highlights of the additional information that the Planck lensing likelihood contributes to the parameter constraints of Planck Collaboration XVI (2014) are presented in Sect. 6.1.

2.
Our lensing potential map correlates with the CMB temperature via the latetime integrated SachsWolfe (ISW) effect sourced by dark energy. We measure the crossspectrum between our lens reconstruction and the lowℓ temperature multipoles in Sect. 6.2, which we find is present at the level predicted in ΛCDM, and discrepant with the null hypothesis of no correlation at approximately 2.5σ.

3.
Our lensing potential map correlates with the fluctuations of the cosmic infrared background, which also trace largescale structure and have good redshift overlap with the CMB lensing potential. We do not discuss this correlation here, but instead refer to the detailed analysis of Planck Collaboration XVIII (2014).

4.
Finally, our lensing potential map correlates with external tracers of largescale structure. In Sect. 6.3 we present crosscorrelations with several representative galaxy catalogues.
Our primary results are based on a minimumvariance (MV) combination of the 143 and 217 GHz Planck maps, although as an initial consistency check we will also present most of our results in this section for the individual frequency maps as well. In Sect. 7 we will present a more extensive suite of consistency tests on our power spectrum estimates and lensing map. In particular in Sect. 7.1, we will compare with results obtained on different masks and on the foregroundcleaned Planck maps described in Planck Collaboration XII (2014).
Our MV lensing map is shown in Fig. 8, where we plot the Wienerfiltered reconstruction of the lensing potential. At the noise levels of our reconstruction, approximately half of the modes in this map are noise. At a visual level, the lensing map of Fig. 8 has comparable signaltonoise to the COBEDMR maps of the CMB temperature anisotropy (Bennett et al. 1996). Future measurements, including the Planck full mission release which contains approximately twice the amount of data used here, will improve on this first fullsky map of the CMB lensing potential. As is illustrated in the simulated reconstruction of Fig. 4, there will be clear visual correlations between this map and future measurements. For comparison with the full MV result, in Fig. 9 we also plot the individual 143 and 217 GHz reconstructions for the southern hemisphere.
Fig. 8
Wienerfiltered lensing potential estimate for our MV reconstruction, in Galactic coordinates using orthographic projection. The reconstruction is bandpass filtered to L ∈ [ 10,2048 ]. The Planck lens reconstruction has S/N ≤ 1 for individual modes on all scales, so this map is noise dominated. Comparison between simulations of reconstructed and input φ in Fig. 4 show the expected level of visible correlation between our reconstruction and the true lensing potential. 
Fig. 9
Wienerfiltered lensing potential estimates, as in Fig. 8, for the individual 143 and 217 GHz maps. Note that the “noise” due to CMB fluctuations is correlated between these two estimates. 
In Fig. 10 we plot the power spectra of our individual 100, 143, and 217 GHz reconstructions as well as the minimumvariance reconstruction. The agreement of all four spectra is striking.
Fig. 10
Lensing potential power spectrum estimates based on the individual 100, 143, and 217 GHz sky maps, as well our fiducial minimumvariance (MV) reconstruction which forms the basis for the Planck lensing likelihood. The black line is for the bestfit ΛCDM model of Planck Collaboration XVI (2014). 
Overall, our power spectrum measurement is reasonably consistent with the ΛCDM prediction, given our measurement error bars. Dividing the L ∈ [ 1,2048 ] multipole range into bins of ΔL = 64 and binning uniformly in , we obtain a reduced χ^{2} for the difference between our power spectrum estimate and the model of 40.7 with 32 degrees of freedom. The associated probability to exceed is 14%. On a detailed level, there are some discrepancies between the shape and amplitude of our power spectrum and the fiducial model however. Our likelihood is based on the multipole range 40 ≤ L ≤ 400, which captures 90% of the forecasted ensembleaverage signaltonoise for an amplitude constraint on in our fiducial cosmological model (see Fig. 2). This range was chosen as the region of our spectrum least likely to be contaminated by systematic effects (primarily uncertainties in the meanfield corrections at lowL, and uncertainties in the Gaussian and pointsource bias corrections at highL). Estimating an average amplitude for the fiducial lensing power spectrum for a single bin over this multipole range using Eq. (25)we find an amplitude of Â_{40 → 400} = 0.94 ± 0.04 relative to the fiducial model (which has A = 1). The power in this region is consistent with the fiducial model, although 1.5σ low (the corresponding probabilitytoexceed for the χ^{2} of this difference is 15%). The low and highL extent of our likelihood were deliberately chosen to have enough expected lensing signal to enable a 10σ detection of lensing on either side, bookending our likelihood with two additional consistency tests. On the lowL side, we have a good agreement with the expected power. As will be discussed in Sect. 7.4, our measurement at L< 10 fails some consistency tests at a level comparable to the expected signal. The L< 10 modes, which we suspect are somewhat contaminated by errors in the meanfield subtraction, are nevertheless consistent with the fiducial expectation, as can be seen in Fig. 10; we measure Â_{1 → 10} = 0.44 ± 0.54. Extending to the lower limit of our likelihood, with a single bin from 10 ≤ L ≤ 40 we measure Â_{10 → 40} = 1.02 ± 0.12. On the highL side of our fiducial likelihood, there is tension however. Extending from the final likelihood multipole at L = 400 to the maximum multipole of our reconstruction, we find Â_{400 → 2048} = 0.68 ± 0.13, which is in tension with A = 1 at a level of just over 2.4σ. The relatively low power in our reconstruction is driven by a dip relative to the ΛCDM model spectrum between 500 <L< 750, as can be seen in Fig. 10. We show this feature more clearly in the residual plot of Fig. 11. This deficit of power is in turn driven by the 143 GHz data. For an estimate of the power spectrum using only 143 GHz, we measure . The 217 GHz reconstruction is more consistent with the model, having . These two measurements are in tension; we have , which is a 2.5σ discrepancy. The error bar on this difference accounts for the expected correlation between the two channels due to the fact that they see the same CMB sky. A larger set of consistency tests will be presented in Sect. 7. We note for now that the bins from 40 <L< 400 used in our likelihood pass all consistency tests, and show better agreement between 143 and 217 GHz. Although L< 40 and L> 400 are not included in our nominal likelihood, when discussing the use of the lensing likelihood for cosmological parameter constraints in the following section we will perform additional crosschecks using these bins to ascertain whether they would have any significant implications for cosmology.
In addition to the Planck power spectrum measurements, in Fig. 11 we have overplotted the ACT and SPT measurements of the lensing potential power spectrum (Das et al. 2013; van Engelen et al. 2012). It is clear that all are very consistent. The Planck data provides the largest signaltonoise of these measurements; as we have already discussed the 40 <L< 400 lensing likelihood provides a 4% constraint on the amplitude of the lensing potential power spectrum, while the constraint from current ACT and SPT measurements are 32% and 16% respectively. These measurements are nevertheless quite complementary. As a function of angular scale, the fullsky Planck power spectrum estimate has the smallest uncertainty per multipole of all three experiments at L< 500, at which point the additional smallscale modes up to ℓ_{max} = 3000 used in the SPT lensing analysis lead to smaller error bars. The good agreement in these estimates of is reassuring; in addition to the fact that the experiments and analyses are completely independent, these measurements are sourced from fairly independent angular scales in the temperature map, with ℓ ≲ 1600 in the case of Planck, ℓ < 2300 in the case of ACT, and ℓ < 3000 in the case of SPT. Crosscorrelation of the Planck lensing map with these independent measures of the lensing potential will provide an additional crosscheck on their consistency, however at the power spectrum level they are already in good agreement.
Fig. 11
Replotting of Fig. 10, removing 100 GHz for easier comparison of 143 and 217 GHz. Also plotted are the SPT bandpowers from van Engelen et al. (2012), and the ACT bandpowers from Das et al. (2013). All three experiments are very consistent. The lower panel shows the difference between the measured bandpowers and the fiducial bestfit ΛCDM model. 
6.1. Parameters
Weak gravitational lensing of the CMB provides sensitivity to cosmological parameters affecting the latetime growth of structure that are otherwise degenerate in the primary CMB anisotropies imprinted around recombination. Examples include the dark energy density in models with spatial curvature and the mass of neutrinos that are light enough (m_{ν}< 0.5 eV) still to have been relativistic at recombination.
To connect our measurement of the lensing power spectrum to parameters, we construct a lensing likelihood nominally based on the multipole range 40 ≤ L ≤ 400, cut into eight equalwidth bins with ΔL = 45 to maintain parameter leverage from shape information in addition to our overall amplitude constraint. In Table 1 we present bandpowers for these eight bins using the individual 100, 143, and 217 GHz reconstructions as well as the MV reconstruction that is the basis for our nominal likelihood. The bandpower estimates and their uncertainties are broken down into constituent parts as discussed in Sect. 2. Based on these bandpowers, we form a likelihood following Eq. (23). The measurement errors on each bin are measured by MonteCarlo using 1000 simulations, and the bins are sufficiently wide that we can neglect any small covariance between them (this is discussed further in Appendix D). We analytically marginalize over uncertainties that are correlated between bins, including them in the measurement covariance matrix. This includes beam transfer function uncertainties (as described in Sect. 5.2), uncertainties in the point source correction (Sect. 7.2) and uncertainty in the N^{(1)} correction.
As the lensing likelihood is always used in conjunction with the PlanckTT power spectrum likelihood, we coherently account for uncertainty in by renormalizing our lensing potential measurement for each sample, as described in Sect. 5.3.
The lensing likelihood is combined with the main PlanckTT likelihood (Planck Collaboration XV 2014) – constructed from the temperature (pseudo) crossspectra between detector sets at intermediate and high multipoles, and an exact approach for Gaussian temperature anisotropies at low multipoles – in Planck Collaboration XVI (2014) to derive parameter constraints for the sixparameter ΛCDM model and wellmotivated extensions. Lensing also affects the power spectrum, or 2point function, of the CMB anisotropies, and this effect is accounted for routinely in all Planck results. On the angular scales relevant for Planck, the main effect is a smoothing of the acoustic peaks and this is detected at around 10σ in the Planck temperature power spectrum (Planck Collaboration XVI 2014). The information about that is contained in the lensed temperature power spectrum for multipoles ℓ ≲ 3000 is limited to the amplitude of a single eigenmode (Smith et al. 2006). In extensions of ΛCDM with a single additional latetime parameter, lensing of the power spectrum itself can therefore break the geometric degeneracy (Stompor & Efstathiou 1999; Sherwin et al. 2011; van Engelen et al. 2012; Planck Collaboration XVI 2014). As discussed in Appendix D and Schmittfull et al. (2013), cosmic variance of the lenses produces weak correlations between the CMB 2point function and our estimates of , but they are small enough that ignoring the correlations in combining the two likelihoods should produce only subpercent underestimates of the errors in physical cosmological parameters.
In the following, we illustrate the additional constraining power of our measurements in ΛCDM models and oneparameter extensions, highlighting those results from Planck Collaboration XVI (2014) where the lensing likelihood is influential.
6.1.1. Sixparameter ΛCDM model
In the sixparameter ΛCDM model, the matter densities, Hubble constant and spectral index of the primordial curvature perturbations are tightly constrained by the Planck temperature power spectrum alone. However, in the absence of lensing the amplitude A_{s} of the primordial power spectrum and the reionization optical depth τ are degenerate, with only the combination A_{s}e^{− 2τ}, which directly controls the amplitude of the anisotropy power spectrum on intermediate and small scales, being well determined. This degeneracy is broken by largeangle polarization since the power from scattering at reionization depends on the combination A_{s}τ^{2}. In this first release of Planck data, we use the WMAP nineyear polarization maps (Bennett et al. 2013) in combination with Planck temperature data. With this data combination, is rather tightly constrained in the ΛCDM model (see Fig. 12) and the direct measurements reported here provide a nontrivial consistency test of the model.
Fig. 12
Upper left: Planck measurements of the lensing power spectrum compared to the ΛCDM mean prediction and 68% confidence interval (dashed lines) for models fit to Planck+WP+highL (see text). The eight bandpowers are those used in the Planck lensing likelihood; they are renormalized, along with their errors, to account for the small differences between the lensed in the bestfit model and the fiducial model used throughout this paper. The error bars are the ± 1σ errors from the diagonal of the covariance matrix. The colour coding shows how varies with the optical depth τ across samples from the ΛCDM posterior distribution. Upper right: as upperleft but using only the temperature power spectrum from Planck. Lower left: as upperleft panel but in models with spatial curvature. The colour coding is for Ω_{K}. Lower right: as upperleft but in models with three massive neutrinos (of equal mass). The colour coding is for the summed neutrino mass ∑ m_{ν}. 
The eight bandpowers used in the lensing likelihood are compared to the expected spectrum in Fig. 12 (upperleft panel). For the latter, we have used parameter values determined from the main Planck likelihood in combination with WMAP polarization (hereafter denoted WP) and smallscale power spectrum measurements (hereafter highL) from ACT (Das et al. 2013) and SPT (Reichardt et al. 2012)^{4}. In this plot, we have renormalized the measurements and their error bars (rather than the theory) using the bestfit model with a variant of the procedure described in Sect. 5.3. Since the lensed temperature power spectrum in the bestfit model is very close to that in the fiducial model used to normalise the power spectrum estimates throughout this paper, the power spectrum renormalisation factor of Eq. (44) differs from unity by less than 0.5%. The predicted in the bestfit model differs from the fiducial model by less than 2.5% for L< 1000. The bestfit model is a good fit to the measurements, with χ^{2} = 10.9 and the corresponding probability to exceed equal to 21%^{5}. Significantly, we see that the ΛCDM model, calibrated with the CMB fluctuations imprinted around z = 1100, correctly predicts the evolution of structure and geometry at much lower redshifts. The 68% uncertainty in the ΛCDM prediction of is shown by the dashed lines in the upperleft panel of Fig. 12. We can assess consistency with the direct measurements, properly accounting for this uncertainty, by introducing an additional parameter that scales the theory in the lensing likelihood. (Note that we choose not to alter the lensing effect in .) As reported in Planck Collaboration XVI (2014) we find
in excellent agreement with . This analysis differs in detail from the simple oneparameter analysis with the fiducial template, reported in Table 1 as A = 0.943 ± 0.04 for the MV reconstruction. The latter does not account for uncertainty in the amplitude or shape of the template between bandpowers due to variations in the six ΛCDM parameters, nor small variations in the lensed temperature power spectrum (which affect the normalisation of the reconstruction power spectrum).
Fig. 13
Marginalised constraints on the optical depth in ΛCDM models from the Planck temperature power spectrum (Planck; solid black), and additionally including the lensing likelihood (Planck+lensing; dashed red) or WMAP polarization (Planck+WP; dasheddotted blue). We use a prior τ> 0.01 in all cases. 
An alternative route to breaking the A_{s}τ degeneracy is possible for the first time with Planck. Since is directly proportional to A_{s}, the lensing power spectrum measurements and the smoothing effect of lensing in (which at leading order varies as ) can separately constrain A_{s} and τwithout largeangle polarization data. The variation of with τ in ΛCDM models constrained only by the Planck temperature power spectrum is illustrated in the upperright panel of Fig. 12, and suggests that the direct measurements may be able to improve constraints on τ further. This is indeed the case, as shown in Fig. 13 where we compare the posterior distribution of τ for the Planck temperature likelihood alone with that including the lensing likelihood. We find At 95% confidence, we can place a lower limit on the optical depth of 0.04 (Planck+lensing). This is very close to the optical depth for instantaneous reionization at z = 6, providing further support for reionization being an extended process.
The τ constraints via the lensing route are consistent with, though weaker, than those from WMAP polarization. However, since the latter measurement requires very aggressive cleaning of Galactic emission (see e.g., Fig. 17 of Page et al. 2007), the lensing constraints are an important crosscheck.
6.1.2. Effect of the large and small scales on the sixparameter ΛCDM model
Before exploring the further parameters that can be constrained with the lensing likelihood, we test the effect on the ΛCDM model of adding the largescale (10 ≤ L ≤ 40) and smallscale (400 ≤ L ≤ 2048) lensing data to our likelihood. Adding additional data will produce random shifts in the posterior distributions of parameters, but these should be small here since the multipole range 40 ≤ L ≤ 400 is designed to capture over 90% of the signaltonoise (on an amplitude measurement). If the additional data are expected to have little statistical power, i.e., the error bars on parameters do not change greatly, but their addition produces large shifts in the posteriors, this would be symptomatic either of internal tensions between the data or an incorrect model.
Fig. 14
Marginalised posteriors for the sixparameter ΛCDM model, shown as box plots, for Planck+WP+highL with various lensing likelihoods. The red and blue lines are the median and mean, respectively. The box and bar correspond to 68% and 95% of the probability density, both centered on the median. The leftmost column is without the lensing likelihood and the median of these constraints is shown by the grey line. The remaining columns show the effect of adding in the fiducial lensing likelihood (second column), and further adding a lowL bin (third column), highL bins (fourth column), or both (final column). 
In Fig. 14, we compare the posterior distributions of the ΛCDM parameters for Planck+WP+highL alone with those after combining with various lensing likelihoods. Adding our fiducial lensing likelihood (second column) reduces the errors on parameters by a small amount and the median values shift by rather less than 1σ for all parameters. The largest gain is for Ω_{c}h^{2} (and H_{0}) where the errors improve by 20%. Adding further large and smallscale data produces no significant reduction in error bars, as expected. Parameter medians also change very little except for Ω_{c}h^{2}, which is dragged low by a further 0.3σ on adding the smallscale lensing information. (The shift in H_{0} is due to the anticorrelation between H_{0} and Ω_{c}h^{2} caused by the acousticscale degeneracy in the temperature power spectrum; see Planck Collaboration XVI 2014.) These findings are consistent with the power spectrum amplitude measurements discussed in Sect. 6: we can lower the lensing power by reducing the matter density, and this is favoured by the lower amplitudes measured from the smallscale lensing power spectrum.
The tension between the smallscale power and the power over the L = 40–400 range included in our fiducial likelihood, coupled with our lower confidence in the accuracy of the bias removal on small scales, is the reason that we do not include these smaller scales at this stage in the Planck lensing likelihood.
Fig. 15
Two views of the geometric degeneracy in curved ΛCDM models which is partially broken by lensing. Left: the degeneracy in the Ω_{m}Ω_{Λ} plane, with samples from Planck+WP+highL colour coded by the value of H_{0}. The contours delimit the 68% and 95% confidence regions, showing the further improvement from including the lensing likelihood. Right: the degeneracy in the Ω_{K}H_{0} plane, with samples colour coded by Ω_{Λ}. Spatiallyflat models lie along the grey dashed lines. 
6.1.3. Spatial curvature and dark energy
Inflation models with sufficient number of efolds of expansion naturally predict that the Universe should be very close to spatially flat. Constraining any departures from flatness is therefore a critical test of inflationary cosmology. However, the primary CMB anisotropies alone suffer from a geometric degeneracy, whereby models with identical primordial power spectra, physical matter densities, and angulardiameter distance to lastscattering have almost identical power spectra (Efstathiou & Bond 1999). The degeneracy is partly broken by lensing (Stompor & Efstathiou 1999), with small additional contributions from the lateISW effect (on large scales) and by projection effects in curved models (Howlett et al. 2012). In ΛCDM models with curvature, the geometric degeneracy is twodimensional, involving the curvature and dark energy density, and this limits the precision with which either can be determined from the CMB alone.
With the highsignificance detection of lensing by Planck in the temperature power spectrum (Planck Collaboration XVI 2014), and via the lens reconstruction reported here, the geometric degeneracy is partially broken, as shown in Fig. 15. The long tail of closed models with low dark energy density (and expansion rate at low redshift) allowed by the geometric degeneracy have too much lensing power to be consistent with Planck’s measured temperature and lensing power spectra (see also Fig. 12). We find marginalised constraints on the curvature parameter of so that lensing reconstruction reduces the uncertainty on Ω_{K} by more than a factor of two over limits driven by the smoothing effect on the acoustic peaks of . This improvement is consistent with the spread in in curved models constrained by the temperature power spectrum, relative to the errors on the reconstruction power spectrum; see Fig. 12. Note that the mean value of Ω_{K} also moves towards zero with the inclusion of the measurements. Adding the highL and lowL data to the likelihood brings no more than a percentlevel improvement on the constraint. We see that the CMB alone now constrains the geometry to be flat at the percent level. Previous constraints on curvature via CMB lensing have been reported by SPT in combination with the WMAP7 data: (68%; Story et al. 2013). This constraint is consistent, though almost a factor of two weaker, than that from Planck. Tighter constraints on curvature result from combining the Planck data with other astrophysical data, such as baryon acoustic oscillations, as discussed in Planck Collaboration XVI (2014).
Lensing effects provide evidence for dark energy from the CMB alone, independent of other astrophysical data (Sherwin et al. 2011; van Engelen et al. 2012). In curved ΛCDM models, we find marginalised constraints on Ω_{Λ} of Again, lensing reconstruction improves the errors by more than a factor of two over those from the temperature power spectrum alone. Note that part of the preference for closed models from the Planck + WP + highL data combination (i.e., without the 4point function) is related to the tendency of the temperature power spectrum, considered over the full multipole range, to favour models with more lensing power than in the bestfit ΛCDM model. This effect is discussed in some detail in Planck Collaboration XVI (2014).
6.1.4. Neutrino masses
The unique effect in the unlensed temperature power spectrum of massive neutrinos that are still relativistic at recombination is small. With the angular scale of the acoustic peaks fixed from measurements of the temperature power spectrum, neutrino masses increase the expansion rate at z> 1 where dark energy is negligible^{6} and so suppress clustering on scales smaller than the horizon size at the nonrelativistic transition (Kaplinghat et al. 2003). This effect reduces for L> 10 (see Fig. 12) and gives less smoothing of the acoustic peaks in . As discussed in Planck Collaboration XVI (2014), the constraint on ∑ m_{ν} from the Planck temperature power spectrum (and WMAP lowℓ polarization) is driven by the smoothing effect of lensing: ∑ m_{ν}< 0.66 eV (95%; Planck+WP+highL). Curiously, this constraint is weakened by additionally including the lensing likelihood to (45)reflecting mild tensions between the measured lensing and temperature power spectra, with the former preferring larger neutrino masses than the latter. Possible origins of this tension are explored further in Planck Collaboration XVI (2014) and are thought to involve both the measurements and features in the measured on large scales (ℓ < 40) and small scales ℓ > 2000 that are not fit well by the ΛCDM+foreground model. In regard to , Fisher estimates show that the bandpowers in the range 130 <L< 309 carry most of the statistical weight in determining the marginal error on ∑ m_{ν}, and Fig. 12 reveals a preference for high ∑ m_{ν} from this part of the spectrum. We have checked that removing the first bandpower from the lensing likelihood, which is the least stable to data cuts and the details of foreground cleaning as discussed in Sect. 7, has little impact on our neutrino mass constraints. We also note that a similar trend for lower lensing power than the ΛCDM expectation on intermediate scales is seen in the ACT and SPT measurements (Fig. 11). Adding the highL information to the likelihood weakens the constraint further, pushing the 95% limit to 1.07 eV. This is consistent with our smallscale measurement having a significantly lower amplitude. At this stage it is unclear what to make of this mild tension between neutrino mass constraints from the 4point function and those from the 2point function, and we caution overinterpreting the results. We expect to be able to say more on this issue with the further data that will be made available in future Planck data releases.
6.2. Correlation with the ISW Effect
As CMB photons travel to us from the lastscattering surface, the gravitational potentials they traverse may undergo a nonnegligible amount of evolution. This produces a net redshift or blueshift of the photons concerned, as they fall into and then escape from the evolving potentials. The overall result is a contribution to the CMB temperature anisotropy known as the latetime ISW effect, or the ReesSciama (RS) effect depending on whether the evolution of the potentials concerned is in the linear (ISW) or nonlinear (RS) regime of structure formation (Sachs & Wolfe 1967; Rees & Sciama 1968). In the epoch of dark energy domination, which occurs after z ~ 0.5 for the concordance ΛCDM cosmology, largescale potentials tend to decay over time as space expands, resulting in a net blueshifting of the CMB photons which traverse these potentials.
In the concordance ΛCDM model, there is significant overlap between the largescale structure which sources the CMB lensing potential φ and the ISW effect (greater than 90% at L< 100), although it should be kept in mind that we cannot observe the ISW component by itself, and so the effective correlation with the total CMB temperature is much smaller, on the order of 20%.
The correlation between the lensing potential and the ISW effect results in a nonzero bispectrum or threepoint function for the observed CMB fluctuations. This bispectrum is peaked for “squeezed” configurations, in which one short leg at lowℓ supported by the ISW contribution is matched to the lensinginduced correlation between two smallscale modes at highℓ. Constraints on the amplitude of the lensingISW bispectrum using several different estimators are presented in Planck Collaboration XXIV (2014). Here we will present an additional constraint, in which the bispectrum measurement is recast as an estimate for the amplitude of the crossspectrum , using the filtering and frequency map combinations of our baseline lensing reconstruction. Our measurements are in good agreement with those made in Planck Collaboration XXIV (2014); a detailed comparison of several lensingISW bispectrum estimators, including the one used here, is presented in Planck Collaboration XIX (2014).
Following Lewis et al. (2011), we begin with an estimator for the crossspectrum of the lensing potential and the ISW effect as (46)where is the Wienerfiltered temperature map and is given in Eq. (13). In Fig. 16 we plot the measured crossspectra for our individual frequency reconstructions at 100, 143, and 217 GHz as well as the MV reconstruction. We also plot the mean and scatter expected in the fiducial ΛCDM model.
Fig. 16
LensingISW bispectrumrelated cross spectra computed from Eq. (46). Black dashed lines indicate the average value for simulations, while dark/light grey filled regions indicate the expected one/two standard deviation scatter, also measured from simulations. The thin magenta line gives the expected crossspectrum for our fiducial model. The agreement of this curve with the simulation average illustrates that our estimator is accurately normalized. In all the quantitative analysis of this section we ignore L< 10, although we have plotted the crossspectra at these multipoles for interest. 
To compare quantitatively the overall level of the measured correlation to the value in ΛCDM, we estimate an overall amplitude for the crossspectrum as (47)The overall normalization is determined from MonteCarlo simulations. For our processing of the data, we find that it is well approximated (at the 5% level) by the analytical approximation (48)The estimator above is equivalent to the KSW and skewC_{ℓ} estimators of Komatsu et al. (2005); Munshi et al. (2011b) for the lensingISW bispectrum that are used in Planck Collaboration XXIV (2014; up to implementation details such as filtering). The meanfield subtraction performed when computing can be identified with the linear term of Creminelli et al. (2006), which is necessary to minimize the estimator variance. The contribution to the total S/N of this estimator as a function of the short leg L is plotted in Fig. 2, where it can be seen that the constraining power for the fiducial correlation is almost entirely at L< 100.
In Table 2 we present measured values for the amplitude of the lensingISW bispectrum using Eq. (46). The uncertainties on Â^{Tφ} are determined by MonteCarlo. We use the multipole range 10 <L< 100, given some of the potential systematic issues with L< 10 identified in Sect. 7.4, although as can be seen from Fig. 16, the inclusion of lower multipoles does not significantly affect our results. Note that for the ISWlensing measurements, inaccuracies in the meanfield subtraction do not bias the estimator although they may degrade the statistical errors on large scales. The differences between the different amplitude fits are well within the expected scatter, as we show in Table 3.
As a point of interest, we have also split our amplitude constraint into the contribution from even and odd multipoles. There are well known odd/evenmultipole power asymmetries in the temperature anisotropies on large angular scales, the study of which is somewhat limited by the small number of available modes (Land & Magueijo 2005; Kim & Naselsky 2010; Gruppuso et al. 2011; Bennett et al. 2011). The lensing potential gives a potentially new window on these power asymmetries, as a third somewhat independent measurement of power on large angular scales. As we can see in Table 2, there is a large difference between the odd/evenL contributions in our lensingISW bispectrum estimate, related to the odd/even temperature asymmetry.
In Table 3 we present differences between the amplitude measurements of Table 2, as well as the expected scatter, accounting for correlations between the estimates due to common CMB and (in the case of MV vs. 143 or 217 GHz) noise fluctuations. We see that our estimates are all very consistent. The difference in measured odd/even multipole power, while striking, is not statistically significant at greater than 2σ.
To conclude, we see a correlation between our lensing potential measurement and the largescale temperature anisotropies that is consistent with the level expected due to the presentday dark energy domination in our fiducial ΛCDM cosmology. Our amplitude measurements are in agreement with those presented in Planck Collaboration XXIV (2014); these independent measurements are compared in detail in Planck Collaboration XIX (2014).
6.3. Correlation with galaxy catalogues
Fig. 17
Crossspectra of the Planck MV lensing potential with several galaxy catalogues, scaled by the signaltonoise weighting factor defined in Eq. (53). Crosscorrelations are detected at approximately 20σ significance for the NVSS quasar catalogue, 10σ for SDSS LRGs, and 7σ for both MaxBCG and WISE. 
The CMB lensing potential is an integrated measure of all matter in the Universe back to the lastscattering surface, with a broad kernel which peaks at z ~ 2, but has significant contributions from both lower and higher redshifts. There are therefore significant correlations between the CMB lensing potential and other tracers of largescale structure. Such correlations have already been observed between the lensing potential reconstructed from WMAP data and the NVSS quasar catalogue (Smith et al. 2007; Hirata et al. 2008, on the order of 3σ in significance), as well as between the lensing potential from SPT with optical and infrared catalogues from BCS, WISE, and Spitzer/IRAC (Bleem et al. 2012, between 4 and 5σ), and the ACT lensing potential with quasars from SDSS (Sherwin et al. 2012, at 3.8σ). Here we will show several representative examples of correlations with the Planck lensing potential, between 7σ and 20σ in significance. Our goal is not to perform an in depth study, but rather to demonstrate the power of our public lensing map. In addition to the correlations presented here, there is a powerful correlation with, the cosmic infrared background (CIB). Correlation of the Planck lensing potential with the CIB fluctuations as probed by the highest frequency Planck channels is observed at greater than 40σ in significance, and has been subjected to a more detailed analysis and modelling, which is presented in an accompanying paper, Planck Collaboration XVIII (2014).
To predict expected levels of correlation for a given galaxy catalogue, we use the Limber approximation (Limber 1954) with a simple linear bias model, in which the crossspectrum between two mass tracers is given by (49)where χ is conformal distance and P(k,χ) is the 3D matter power spectrum for wavenumber k at conformal lookback time χ. K^{φ} and K^{g} are kernels associated with lensing and with the galaxy catalogue of interest respectively, and are given by (50)Here b(z) is a redshiftdependent linear bias parameter, and dN/dz describes the redshift distribution of the galaxy population. As with the lensingISW bispectrum, we can construct a simple pseudoC_{L} estimator for the cross correlation using (51)where f_{sky,j} is the sky fraction common to both the lens reconstruction and the catalogue, and g_{LM} is the harmonic transform of the galaxy fractional overdensity. Denoting the positions of the objects p in the catalogue as , the transform is given by (52)where N^{gg} is the surface density of objects in gal/steradian. With a fiducial model for the crossspectrum obtained using Eq. (49), the minimumvariance estimator for its overall amplitude is (53)where is the signal power spectrum of the catalogue, which can be estimated using Eq. (49)with K^{g} for both weight functions. Here we have defined the spectrum as a scaling at each multipole L of the crosscorrelation power spectrum. The normalization is given by (54)In Fig. 17 we plot the contributions to Â^{gφ} as a function of L for several surveys that have significant correlations with the Planck MV lensing potential: the NVSS quasar catalogue, the MaxBCG cluster catalogue, an SDSS LRG catalogue, and an infrared catalogue from the WISE satellite. The error bars for each correlation are measured from the scatter of simulated lens reconstructions correlated with each catalogue map, and are in generally good agreement (at the 20% level) with analytical expectations. These catalogues are discussed in more detail below.

1.
NVSS Quasars: the NRAO VLA Sky Survey (NVSS; Condonet al. 1998) is a catalogue ofapproximately two million sources north ofδ = −40deg which is 50% complete at 2.5 mJy. Most of the bright sources are AGNpowered radio galaxies and quasars. We process this catalogue following Smith et al. (2007), pixelizing the catalogue at HEALPix N_{side} = 256 and projecting out the azimuthallysymmetric modes of the galaxy distribution in ecliptic coordinates to avoid systematic striping effects in the NVSS dataset. We model the expected crosscorrelation for this catalogue using a constant b(z) = 1.7 and a redshift distribution centered at z_{0} = 1.1 given by (55)For this model, in the correlation with the MV lens reconstruction we measure an amplitude of .

2.
SDSS LRGs: We use the LRG catalogue of Ross et al. (2011); Ho et al. (2012) based on Sloan Digital Sky Survey Data Release 8 (SDSS DR8), which covers 25% of the sky. After cutting to select all sources with photometric redshift 0.4 ≤ z ≤ 0.8, and p_{gal}> 0.2, we are left with approximately 1.4 × 10^{6} objects with a mean redshift of z = 0.55 and a scatter of ± 0.07. Apart from the cut above, we do not perform any additional weighting on p_{gal}. We model this catalogue using dN/dz taken from the histogram of photometric redshifts, and take b(z) = 2. We measure , very consistent with expectation.

3.
MaxBCG Clusters: The MaxBCG cluster catalogue (Koester et al. 2007) is a collection of 13,823 clusters over approximately 20% of the sky selected from the SDSS photometric data, covering a redshift range 0.1 ≤ z ≤ 0.3. It is believed to be 90% pure and more than 85% complete for clusters with M ≥ 1 × 10^{14}M_{⊙}. To simplify the sky coverage, we have discarded the three southern SDSS stripes in the catalogue, which reduces the overall sky coverage to approximately 17%. There are accurate photometric redshifts (Δ_{z} ~ 0.01) for all objects in the catalogue, and so we can construct dN/dz directly from the histogram of the redshift distribution. Although these clusters are at very low redshift compared to the typical structures which source the CMB lensing potential, they are strong tracers of dark matter, with an effective bias parameter of b(z) = 3 (Hütsi 2010). We obtain a similar average bias parameter ⟨ b(M,z) ⟩ for the MaxBCG clusters if we combine the massrichness relation of Bauer et al. (2012) and the halo bias prescription of Tinker et al. (2010). Here we measure a correlation with the Planck lensing potential of . This is significantly larger than expected given the simple model above, although as can be seen in Fig. 17 the shape of the correlation is in reasonable agreement.

4.
WISE Catalogue: The Wide Field Survey Infrared Explorer (WISE) satellite (Wright et al. 2010) has mapped the full sky in four frequency bands W1 − W4 at 3.4, 4.6, 12, and 22 μm respectively. We start from the full mission catalogue, which contains over 560 × 10^{6} objects. To obtain a catalogue with roughly uniform sensitivity over the full sky and to eliminate stellar contamination we follow Kovacs et al. (2013), selecting all sources with W1 magnitudes less than 15.2 at galactic latitudes greater than 10°, and require W1 − W2 > 0.2 and W2 − W3 > 2.9. We cut all sources that are flagged as potentially spurious, or for which more than 30% of frames observing a given source are marked as contaminated by moonglow. These cuts leave us with 2 308 751 sources. We additionally remove all sources that lie outside an f_{sky} = 0.6 Galaxy mask. Following Goto et al. (2012), we model the redshift distribution of these sources using (56)with z_{0} = 0.1, σ_{z} = 0.1, and where the proportionality constant is chosen such that ^{∫}dzdN/ dz = 1. The bias is taken to be b = 1. For this model we measure an amplitude of . Note that the cuts we have used above are fairly conservative – in particular, relaxing the magnitude cut above we can obtain much larger crosscorrelation signals (on the order of 40σ in significance).
It can be seen from Fig. 17 that we have significant correlations with the objects in all four of these catalogues, with amplitudes and shapes roughly at the expected level. The correlations here are only a demonstrative sample of what is possible with the Planck lensing potential map. Note that although the lensing potential has a peak sensitivity to highredshift tracers (such as NVSS), significant correlations may also be observed with relatively low redshift catalogues. From the perspective of crosscorrelation, the Planck lensing potential represents a mass survey with a well understood redshift kernel over the fullsky, and we anticipate that correlations with it will prove useful as a direct probe of the correlation between luminous and dark matter.
7. Consistency and systematic tests
In this section we perform a number of consistency and systematic tests. Most of our tests will focus on the robustness of the MV reconstruction that forms the basis for our fiducial lensing likelihood. In Fig. 18 we summarize several of our main consistency checks: foreground tests obtained by comparing individual frequencies, adjusting the Galactic mask level, or comparing to results obtained analyzing componentseparated maps (discussed in Sect. 7.1), tests of our correction at the power spectrum level for contamination by the shotnoise from unresolved pointsource (in Sect. 7.2), tests for possible noise bias obtained by taking appropriate correlations of maps with independent noise realizations (discussed further in Sect. 7.3), and tests using alternative estimators that are bias hardened against various systematics (Sect. 7.4). In all cases, we find that our MV power spectrum is generally consistent with these alternative reconstructions.
Our baseline reconstruction method uses the methodology outlined in Sect. 2, as well as a specific choice of filtering to produce the that are the inputs to our lensing estimators. We test our robustness to these choices in Sect. 7.5, as well as the implementation of our main reconstruction pipeline, by comparing to three alternative and largely independent pipelines.
In addition to these targeted tests, we have also performed several generic null tests for unexpected signals in our data. These are deferred to Appendix A, where we study surveytosurvey discrepancies, lensing curl modes, and the maplevel statistics of our lens reconstruction.
Fig. 18
Summary of internal consistency tests between our fiducial minimumvariance (MV) reconstruction and a set of alternatives designed to test sensitivity to potential issues. The top panel shows estimates, with measurement error bars. The bottom panels show the residual with respect to the MV reconstruction in units of the MV measurement uncertainty. The grey band marks the 1σ deviation uncertainty of the MV reconstruction. The error bar on each data point in the lower panels gives the standard deviation of the scatter between each result and the MV, determined from MonteCarlo simulations that account for the correlated CMB, noise and foreground power between estimators. Comparison of the uncertainty on the scatter points and the grey band gives an indication of how constraining each test is. The various tests are described in more detail in subsections of Sect. 7. 
7.1. Foreground contamination
In this section, we test the consistency of our results with respect to the treatment of Galactic foreground contamination. We use several different tests to gain perspective on possible issues:

1.
We compare individual frequency reconstructions. The agreement of the lensing power spectra estimated at 100, 143, and 217 GHz is apparent within the error bars plotted in Fig. 10. However, because a large fraction of the uncertainty in these estimates comes from CMB fluctuations themselves they are significantly correlated. In Fig. 18 we plot the difference between the 143 and 217 GHz reconstructions (the scatter for 100 GHz is too large to provide a useful test of the MV reconstruction so we have not plotted it, although it is also consistent). Both are in reasonable agreement, although we can see that the somewhat large amplitude for our first likelihood bin, just above L = 40 is being driven by 217 GHz.

2.
We perform lens reconstruction using both more aggressive and more conservative Galaxy masks than our fiducial f_{sky} = 0.7 analysis mask, constructed following the description in Sec. 3. The results of this test are also shown in Fig. 18, for masks with f_{sky} = 0.6 and f_{sky} = 0.8. We can see a highly significant discrepancy at very lowL between our fiducial f_{sky} = 0.7 mask and the more aggressive f_{sky} = 0.8 mask, however other bins (in particular in the 40 <L< 400 range) are untouched. Our fiducial results are completely consistent with the more conservative f_{sky} = 0.6 results.

3.
We perform lens reconstruction on the componentseparated CMB maps from Planck Collaboration XII (2014). These are the product of four different foregroundremoval algorithms: Ruler, Nilc, Sevem, and Smica. Each of these methods combines the full set of nine Planck frequency maps from 30 → 857 GHz to obtain a bestestimate CMB map. This makes a detailed characterization of these maps more difficult. However analysis of these maps does provide a valuable consistency test. In Fig. 18 we compare the reconstructions obtained on these componentseparated maps to our MV reconstruction, analyzed using the same Galaxy mask. The foregroundcleaned maps are generally in good agreement with our fiducial reconstruction, although the first bin above L = 40 is consistently about one standarddeviation lower than our MV result, which as we have already pointed out is driven high by 217 GHz. The confidence masks produced in Planck Collaboration XII (2014) cover significantly more sky than the f_{sky} = 0.7 mask used in our baseline analysis.
The remainder of this subsection describes our analysis of the componentseparated maps in more detail. For the power spectrum tests of Fig. 18, we analyze these maps following the same methodology as for our fiducial results. As a mask, we use the union of our fiducial f_{sky} = 0.7 Galaxy mask, point source masks at 100, 143, and 217 GHz, as well as the “method” masks provided for each component separation algorithm. We characterize the bias terms of Eq. (15)as well as the statistical uncertainties of the lens reconstruction on each map using the same lensed CMB and FFP6 noise realizations that are discussed in Sect. 4. These have been run through each of the componentseparation algorithms above, resulting in a set of simulated foregroundcleaned maps (note that no foreground signal was included in the simulations which were run through the separation algorithms, and so these simulations only reproduce the mixing employed by each algorithm of the individual frequency maps). To match the power spectrum of these simulations to the power spectrum of the data maps, we find it is necessary to add extragalactic foreground power following the model in Sect. 4, with A_{cib} = 18 μK^{2} and A_{src} = 28 μK^{2}. The resulting simulations have a power spectrum that agrees with that of the CMB map estimate based on the data to better than 2% at l< 2048. This could be improved slightly by tailoring a specific correction for each map. We also add homogeneous pixel noise with a level of 12 μK arcmin. If we neglected this power, the agreement would be only at the 8% level, primarily due to the noise term (the A_{cib} and A_{src} contributions are each at the level of 1−2%). Due to the procedure that we use to subtract the disconnected noise bias (Eq. (17)) from our lensing power spectrum estimates, the inclusion of these components does not significantly affect our results, but comparison with the values used for our singlefrequency simulations in Sect. 4 are a useful indicator of the extent to which the componentseparation algorithms are able to remove extragalactic foreground power in the highℓ regime.
As already discussed, our results on the componentseparated CMB maps are presented in Fig. 18. Because the CMB and FFP6 noise components of the foregroundcleaned map simulations are the same as those used to characterize our fiducial lens reconstruction, we can measure the expected scatter between the componentseparated maps and our fiducial reconstruction. This scatter will be slightly overestimated because we have not attempted to model coherently the contribution to the reconstruction noise from residual diffuse extragalactic foreground power. For the eight bins in 40 ≤ L ≤ 400 on which our fiducial likelihood is based, we measure a χ^{2} for the difference between our fiducial reconstruction and the corresponding foregroundcleaned reconstruction of χ^{2} = (3.14,4.3,2.5,14.7) for Nilc, Smica, Sevem, and Ruler respectively. These χ^{2} values associated have probabilitytoexceed (PTE) values of (79%,64%,86%,2%). At the level that we are able to test, the Nilc, Smica, and Sevem foregroundcleaned maps give results that are quantitatively consistent with our fiducial reconstruction. There is more scatter between our fiducial reconstruction and the Ruler map than expected from simulations, as evidenced by a very high χ^{2} for the difference, however as can be seen in Fig. 18, there are not any clear systematic differences. Indeed, the discrepancy for the bins plotted in Fig. 18 (which differ somewhat from the linear bins used in our likelihood) is much less significant than for the bins of our fiducial likelihood.
When using the componentseparated maps above, we have used the same f_{sky} = 0.7 Galactic mask as for our MV result, although the confidence regions associated with each foregroundcleaned map allow more sky, ranging up to f_{sky} = 0.94 for the Nilc method. We have used the metis pipeline (described later in Sect. 7.5) to test whether this improved sky coverage could benefit our lens reconstruction. The same method has been used in Planck Collaboration XII (2014) to evaluate possible biases to lens reconstruction induced by these methods using the FFP6 simulated CMB realization, described in Planck Collaboration I (2014), indicating that the different componentseparation algorithms do not alter significantly the lensing signal (at the level which can be tested on a single simulation). Analyzing the Nilc map, which has the largest confidence region, we find that we can increase the usable sky area up to f_{sky} = 0.87 without encountering significant Galactic contamination. In Fig. 19 we show the striking improvement in sky coverage on the Nilc map. Smica and Sevem are very similar; we have not considered Ruler because of its larger noise level.
Power spectrum estimates at this mask level show consistency with the MV reconstruction within two standard deviations of the measurement uncertainty. The increased sky coverage does not bring significant improvements in the errorbars of the power spectrum, however. Using Eq. (20) as an estimate of the power spectrum variance, the larger sky coverage yields only a 3.5% improvement at L< 40 over the MV result, decreasing down to 0 at L = 400. This could be due to the different weighting used in the component separation compared to the one of the MV map, which results in slightly noisier maps for our purpose. While the componentseparated maps allow for a reduced mask, they lead to a marginal improvement of the power spectrum uncertainties. Nevertheless, their agreement with the MV result is reassuring.
Fig. 19
Wienerfiltered potential maps in Galactic coordinates, as in Fig. 8, plotted here in Mollweide projection. Top: MV reconstruction; Bottom: extended reconstruction on the Nilc componentseparated map. 
7.2. Point source correction
As can be seen in Table 1, the unresolved pointsource shotnoise correction in any individual band for our MV likelihood is on the order of a few percent, reaching up to 6% for the highest multipole bands. Averaged over the 40 ≤ L ≤ 400 band, the shot noise correction amounts to a 2% shift in the amplitude of , which is small but nonnegligible compared to our statistical uncertainty of 4%. Physically, the amplitudes of our source corrections are reasonable; at 143 GHz we measure . From the radio pointsource model of De Zotti et al. (2010), this corresponds to an effective flux cut of approximately 150 mJy at this frequency, roughly comparable to that expected for the S/N> 5 cut we make when masking sources in our fiducial analysis (Planck Collaboration XXVIII 2014). The shot noise measured at 217 GHz is lower, as expected given the smaller contribution from radio sources, with . The shotnoise level that we measure for the MV reconstruction lies between these two, with .
As a consistency test of our S/N> 5 source masks, we construct two new sets of source masks that only exclude S/N> 10 sources (to check the stability of our reconstruction under significantly higher source contamination), or that exclude all detected sources with S/N> 4.2 (to test the conservatism of our fiducial flux cut). The results of this test are shown in the masks panel of Fig. 18. Although not shown explicitly in the figure, the source correction term with the S/N> 10 mask is considerably larger than our fiducial correction. For this mask we measure K^{4}, now a 9σ detection of shot noise in the trispectrum. The importance of the correction grows correspondingly, becoming a 9% overall correction to the amplitude in 40 ≤ L ≤ 400. With the source correction, we measure , consistent with the measurement for our fiducial source mask. Without the source correction we would have measured , significantly discrepant. We can see that even for this extremely permissive source mask, our correction for the pointsource shotnoise trispectrum does a good job of rendering insensitive to source contamination. Turning to the S/N> 4.2 mask, we measure K^{4}, reduced from a 2σ measurement of shot noise to one very consistent with zero. For this reconstruction we measure Â_{40 → 400} = 0.93 ± 0.04, in good agreement with the MV result, although there does appear to be an overall trend for lower power above L = 400 in the S/N> 4.2 mask. At a detailed level, both of these results do fail χ^{2} tests on the difference with the MV reconstruction. As the change in sky fraction between these three reconstructions is small (the difference between the fiducial and S/N> 4.2 mask is 1%, and even smaller for the S/N> 10 mask), even the small differences quoted above are significantly larger than they should be if the source correction were perfect. This can be seen in Fig. 18, where the error bars on the difference between each of these reconstructions and the MV result is too small to be visible for the plotted points. For the 40 ≤ L ≤ 400 range on which our likelihood is based, we believe that the discrepancy between the power spectra estimated with these source masks is acceptable. Particularly given the results of the S/N> 10 cut, the source correction that we perform appears to behave reasonably, and reduces the sensitivity of our estimates to the shotnoise trispectrum of unmasked point sources. In Sect. 7.4, we will also use a modified lensing estimator that has zero response to the shotnoise trispectrum and therefore does not require a source correction. As with the tests above, we will find good overall agreement with the standard estimator used for the MV reconstruction.
7.3. Instrumental noise bias
Errors in the instrumental noise model are a source of concern for our lensing power spectrum estimates, in which all four legs of the trispectrum estimate come from a single map. In this situation there are three Wick contractions of the noise which contribute to the power spectrum estimate in two distinct ways, depending on whether they couple between the two quadratic estimators in the power spectrum estimate or not:

1.
The Wick contraction that couples noise contributions inside the same estimator leads to a meanfield term as in Eq. (9). This term is sourced by the statistical anisotropy of the noise (due to the uneven hit distribution across the sky, or the correlations along the scan direction induced by deconvolution of the bolometer time constant). The shape of this meanfield is dictated by the Planck scan strategy, which is coherent over large scales, and so this meanfield is primarily at lowL, below the L_{min} = 40 cutoff of our likelihood. It can be avoided by forming a quadratic estimator in which the two input maps have different noise realizations. We perform a null test against this contribution in Fig. 18, where the (143 × 217) × (143 × 217) panel shows a comparison of our MV power spectrum estimate to that obtained using a quadratic estimator in which one leg comes from 143 GHz and one from 217 GHz. We see generally good agreement with the MV result. For the eight bins in 40 ≤ L ≤ 400 of our fiducial likelihood, we find a χ^{2} for the difference (accounting for correlations between the two estimates) of 4.2, which has a PTE of 83%. The noise meanfield contribution can also be reduced by using an appropriate biashardened estimator. We will perform additional tests with such estimators in Sect. 7.4.

2.
The two Wick contractions that couple noise terms acrossestimators contribute to the disconnected bias . This term describes the noise contribution to the reconstruction variance. Our use of the data covariance matrix in Eq. (17)makes our estimate of the disconnected noise bias insensitive to small errors in the noise model. The noise contribution can also be avoided by taking the crossspectrum of two quadratic estimators with independent noise. We show the result of such a correlation compared to our MV reconstruction in Fig. 18, where the (143 × 143) × (217 × 217) panel shows the crossspectrum between the individual 143 and 217 GHz reconstructions. Again, we see generally good agreement with the MV result. For the eight bins in 40 ≤ L ≤ 400 of our fiducial likelihood, we find a χ^{2} for the difference (accounting for correlations between the two estimates) of 6.5, which has a PTE of 59%.
Based on these tests, we believe that noise bias cannot be a significant contaminant for our power spectrum estimates.
7.4. Biashardened estimators
Most of the results in this paper use the standard weight function to form quadratic estimates for the lensing potential. This weight function is a matched filter for lensing, and results in estimators with minimal variance. As we discuss in Appendix C, however, the standard estimator has large meanfield contributions at lowL, primary due to masking and noise inhomogeneity. As pointed out in Namikawa et al. (2013), these meanfields can be mitigated with the use of “biashardened” lensing estimators, which use weight functions specially constructed to project out these effects. Following the notation of Appendix C, we consider an uncertain bias field z_{LM} with associated weight function that gives a contribution to the offdiagonal elements of the CMB covariance matrix (and therefore the estimator meanfield) of (57)We can then form an estimator which is insensitive to z_{LM} simply as (58)The quadratic estimator (φ − z) that uses the weight function above has a significantly reduced meanfield contribution from z_{LM}. This biashardening procedure may be repeated iteratively to produce weight functions that are insensitive to several sources of meanfield simultaneously. In Fig. 20 we plot a comparison at lowL between estimates with the standard estimator and estimators that are bias hardened against noise and mask meanfields. We see generally good agreement between these sets of estimates for L ≥ 10. Below L = 10, the pseudospectrum of the difference between the standard and biashardened estimators is unexpectedly large. For this reason, we bandlimit our lensing potential map to L ≥ 10, and do not consider L< 10 for any of the quantitative results in this paper.
Fig. 20
Biashardened estimator consistency tests for our MV reconstruction. Following the discussion in Sect. 7.4 we form estimators that are biashardened against the meanfields due to inhomgeneous noise levels (φN, cyan), masking (φM, magenta), and both of these effects simultaneously (φMN, orange). The top panel shows the power spectra of the meanfield for each of these estimators. The middle panel shows the raw power spectra of the reconstructions themselves, compared to the average expected from signal+noise (dashed black). The lower panel shows the power spectrum of the difference between each biashardened estimator and the standard result. Dashed black lines give the expected average of this difference (measured on simulations). The grey line is the fiducial power spectrum. We can see that the differences are consistent with the expected scatter for all estimators at L ≥ 10. 
In addition to lowL consistency tests, we may also use these biashardened estimators to estimate . The resulting power spectra are plotted in Fig. 18. Taking the difference with the eight bins in 40 ≤ L ≤ 400 of our fiducial MV reconstruction, with scatter estimated from simulations, we find χ^{2} values of 6.6 and 4.3 for the mask/noisehardened estimators respectively. These have corresponding PTE values of 58% and 83%. The error bars on obtained with the bias hardened estimators used here are generally between 10 and 20% larger than the error bars of the standard estimator. We also construct an estimator that is bias hardened against the point source weight function of Eq. (28). This estimator has the distinction of having zero response to pointsource shotnoise, as well as any correlation between the pointsource shot noise power and the lensing potential. Again, we find consistent results with the standard estimator, with χ^{2} = 5 and a PTE of 76%.
7.5. Alternative methods
All of the primary results in this paper use a lens reconstruction pipeline based on the methodology outlined in Sect. 2. As a robustness test, both of this methodology and of its implementation, we have implemented three independent reconstruction pipelines. These independent pipelines make significantly different choices than our baseline approach, primarily in the choice of the inversevariance filter function and the calculation of the correction terms in Eq. (15). Our main result of this section is Fig. 21, where we compare the lens reconstruction power spectra for these alternative pipelines to our baseline results at 143 and 217 GHz, as well as the MV reconstruction. The agreement is excellent. The only aspect of implementation common to these four results (apart from the data that was analyzed) is the set of simulations, described in Sect. 4, which are used to estimate the various bias terms and characterize the estimator scatter.
Fig. 21
Comparison of alternative lensing pipelines. The Baseline results use the methodology of Sect. 2. Boxes are for the MV reconstruction, circles show the 143 GHz results, and triangles show the 217 GHz results. The two bottom panels show differences relative to the MV result for 143 and 217 GHz. 
We now proceed to describe these alternative pipelines in somewhat more detail. We note first two common aspects of all three alternative pipelines, which differ from our baseline results. First, rather than accounting for the Galaxy and pointsource mask in the filter function as is done for our baseline method, these methods use a completely diagonal filter, which is the same as the form (given in Eq. (10)) that our baseline filtering asymptotes to far from the mask boundaries. Instead, as we will discuss, the alternative pipelines deal with the mask using different choices of data preparation and selection. A second difference common to all three pipelines with our baseline approach is that when computing the disconnected noise bias , these methods do not use the twopoint expression of Eq. (17), but rather the approximation to it based on Eq. (21), evaluated using an estimate of the data power spectrum. Because this correction is quite large, the agreement of results that use this alternative method to calculate it is reassuring. We now proceed to describe the individual methods in more detail; the common motivation in the development of each has been the reduction of the sharp gradients induced when masking, which can induce a meanfield several orders of magnitude larger than the lensing signal at low multipoles (as discussed in Appendix C). Each method takes a different approach to mitigating this mask effect, as discussed below:

1.
The method iso consists of applying the standard quadratic lensing estimator to the sky map after multiplying by an apodized Galactic mask, and filling point source holes using local constrained Gaussian realizations of the CMB signal+noise. The mask meanfield is proportional to the power spectrum of the mask, and so as apodization smooths the mask boundary (suppressing its power spectrum on small scales), it correspondingly reduces the meanfield significantly. The combination of apodization and source filling makes this estimator very fast to apply to simulations but does require an involved set of correction terms and f_{sky} factors. Our implementation and calculation of this method is described in detail in BenoitLevy et al. (2013).

2.
The metis method consists of inpainting the Galactic mask as well as the pointsource holes, using the sparseinpainting algorithm described in Abrial et al. (2007, 2008). The resulting map resembles a fullsky CMB map, and therefore has no mask meanfield contribution. The inhomogeneous noise and beaminduced meanfields do still have to be corrected however. Our implementation is based on that described in Perotto et al. (2010), with several improvements. In Perotto et al. (2010) lens reconstruction was performed on the inpainted map and then analyzed on the full sky. However, further inspection has revealed that there are some spurious features in the lens reconstruction, localized to the inpainted region inside the Galactic mask. This is likely due to the inhomogeneous noise in Planck that was ignored in previous work and cannot be reproduced by the inpainter. We therefore remask the fullsky lens reconstruction with an apodized Galactic mask (as in Eq. (14)) to remove these regions from our analysis. We follow the same procedure when evaluating the analytical expression for the bias, prewhitening and then applying an apodized Galactic mask to the inpainted temperature multipoles to estimate their power spectrum. Small residual biases are corrected using the same procedure used in the main method.

3.
The patches method avoids the Galactic mask completely by cutting the sky into a collection of 410 small overlapping 10deg × 10deg patches centered on the locations of N_{side} = 8 HEALPix pixels, which are then analyzed under the flatsky approximation. Our implementation of this method is described in Plaszczynski et al. (2012). As with the iso method, pointsource holes are filled using constrained Gaussian realizations. The patches are extracted from a prewhitened CMB map, and apodized with a KaiserBessel window function. The Fourier modes in each patch are fitted in real space using a fast FourierToeplitz algorithm. No meanfield correction is applied. Residual biases due to noise inhomogeneity are removed using a correction that is found to be small. The patches method has been particularly useful in the early stages of our analysis, to identify outliers caused by unmasked point sources.
As can be seen in Fig. 21, all three of these methods are in good agreement with the results of our baseline method, providing reassurance that our results are insensitive to the precise details of our data filtering and reconstruction methodology.
8. Conclusions
The Planck maps have unprecedented sensitivity to gravitational lensing effects. We see significant and consistent measurements of lensing for each of the highresolution CMB channels at 100, 143 and 217 GHz. Even the noisiest channel that we have considered, 100 GHz, provides a 10σ detection of lensing, which is greater than all previous detections. Our fiducial lens reconstruction, based on a minimumvariance combination of the 143 and 217 GHz channels does even better, with a detection of lensing (relative to the null hypothesis of no lensing) at a significance of greater than 25σ. Notably, the noise on our reconstruction is low enough that it is no longer the limiting source of noise for many correlations with largescale structure catalogues (several examples of which we have given in Sect. 6.3). This marks a shift for CMB lensing, from the detection regime into that of standard cosmological probe. Our lensing potential map is publicly available, and we look forward to the uses which may be found fo it.
The percentlevel Planck lensing potential measurement pushes into the realm of precision cosmology, and requires careful validation tests that we have performed in Sect. 7. Our fiducial likelihood, based on the 40 ≤ L ≤ 400 range that is most sensitive to lensing, passes all of the tests that we have performed at an acceptable level. Most importantly, it is consistent with individual frequency reconstructions, with more aggressive or more conservative masking, and with results obtained on more rigorously componentseparated maps. Our measurement of the lensing potential power spectrum is broadly consistent with the ΛCDM expectations of Planck Collaboration XVI (2014), although there are some shape and amplitude tensions that lead to smallerthanforecasted improvements on parameters such as the sum of neutrino masses ∑ m_{ν}. As we saw in Sect. 6.1, our estimate of the potential power spectrum has direct cosmological implications, where it significantly improves CMB powerspectrumonly constraints on curvature (or alternatively, dark energy). Our lensing likelihood also breaks the A_{s} − τ degeneracy, allowing a measurement of the optical depth to reionization which is a useful independent crosscheck on the values determined with large angularscale polarization measurements.
CMB lensing science has undergone a remarkable development in the previous two years (Das et al. 2011; Sherwin et al. 2011; van Engelen et al. 2012). The next year is expected to bring even greater improvements. The 2500 deg^{2} lensing survey of the South Pole Telescope is expected have comparable statistical power to the one reported here, with a highly complementary weighting in Lspace that is more sensitive to smaller angular scale lenses. There is also additional Planck data, which will improve on the measurement reported here considerably. The results of this paper are based on the nominalmission Planck maps, collected during the first fifteen months of science operations. The cryogens used to maintain the 100 mK cooling stage for the Planck HFI instrument lasted considerably longer than the nominalmission lifetime, to a fullmission duration of approximately thirty months. This additional data, set to be released in 2014, results in map noise levels that are roughly a factor two lower in power than those used here. As instrumental noise constitutes approximately half of our error budget, we expect a corresponding decrease of approximately 25% in the variance of our lensing potential map as well as the uncertainties on . The Planck polarization data, also set for release in 2014, will provide an additional powerful probe of lensing.
Planck (http://www.esa.int/Planck) is a project of the European Space Agency (ESA) with instruments provided by two scientific consortia funded by ESA member states (in particular the lead countries France and Italy), with contributions from NASA (USA) and telescope reflectors provided by a collaboration between ESA and a scientific consortium led and funded by Denmark.
The ISWlensing correlation also introduces a nonzero bispectrum. When correlating the reconstructed with the largeangle temperature anisotropies in Sect. 6.2, we are probing this bispectrum.
As discussed in detail in Planck Collaboration XVI (2014), the primary role of the ACT and SPT data in these parameter fits is to constrain more accurately the contribution of extragalactic foregrounds, which must be carefully modelled to interpret the Planck power spectra on small scales. For ΛCDM, the foreground parameters are sufficiently decoupled from the cosmological parameters that the inclusion of the ACT and SPT data has very little effect on the cosmological constraints.
With all other physical densities held fixed to preserve the acoustic physics before recombination, increasing the neutrino mass increases the expansion rate at all times after the neutrinos become nonrelativistic (which is after recombination for the light masses considered here). However, to preserve the angular scale of the sound horizon and hence the location of the acoustic peaks and troughs, one must reduce the expansion rate at low redshift. This is done by reducing the physical density in dark energy.
Acknowledgments
The development of Planck has been supported by: ESA; CNES and CNRS/INSUIN2P3INP (France); ASI, CNR, and INAF (Italy); NASA and DoE (USA); STFC and UKSA (UK); CSIC, MICINN, JA and RES (Spain); Tekes, AoF and CSC (Finland); DLR and MPG (Germany); CSA (Canada); DTU Space (Denmark); SER/SSO (Switzerland); RCN (Norway); SFI (Ireland); FCT/MCTES (Portugal); and PRACE (EU). A description of the Planck Collaboration and a list of its members, including the technical or scientific activities in which they have been involved, can be found at http://www.sciops.esa.int/index.php?project=planck&page=Planck_Collaboration. Some of the results in this paper have been derived using the HEALPix package (Górski et al. 2005). This research used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the US Department of Energy under Contract No. DEAC0205CH11231, as well as of the IN2P3 Computer Center (http://cc.in2p3.fr) and the PlanckHFI data processing center infrastructures hosted at the Institut d’Astrophysique de Paris (France) and financially supported by CNES. We acknowledge support from the Science and Technology Facilities Council [grant number ST/I000976/1].
References
 Abrial, P., Moudden, Y., Starck, J., et al. 2007, J. Fourier Analysis and Applications, 13, 729 [CrossRef] [Google Scholar]
 Abrial, P., Moudden, Y., Starck, J.L., et al. 2008, Statistical Methodology, 5, 289 [NASA ADS] [CrossRef] [Google Scholar]
 Basak, S., Prunet, S., & Benabed, K. 2009, A&A, 508, 2213 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bauer, A. H., Baltay, C., Ellman, N., et al. 2012, ApJ, 749, 56 [NASA ADS] [CrossRef] [Google Scholar]
 Bennett, C., Banday, A., Gorski, K., et al. 1996, ApJ, 464, L1 [NASA ADS] [CrossRef] [Google Scholar]
 Bennett, C., Hill, R., Hinshaw, G., et al. 2011, ApJS, 192, 17 [NASA ADS] [CrossRef] [Google Scholar]
 Bennett, C., Larson, D., Weilard, J. L., et al. 2013, ApJS, 208, 20 [NASA ADS] [CrossRef] [Google Scholar]
 BenoitLevy, A., Dechelette, T., Benabed, K., et al. 2013, A&A, 555, A37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bernardeau, F. 1997, A&A, 324, 15 [NASA ADS] [Google Scholar]
 Blanchard, A., & Schneider, J. 1987, A&A, 184, 1 [NASA ADS] [Google Scholar]
 Bleem, L., van Engelen, A., Holder, G., et al. 2012, ApJ, 753, L9 [NASA ADS] [CrossRef] [Google Scholar]
 Bucher, M., Carvalho, C. S., Moodley, K., & Remazeilles, M. 2012, Phys. Rev. D, 85, 043016 [NASA ADS] [CrossRef] [Google Scholar]
 Challinor, A., & Lewis, A. 2005, Phys. Rev. D, 71, 103010 [NASA ADS] [CrossRef] [Google Scholar]
 Challinor, A., & van Leeuwen, F. 2002, Phys. Rev. D, 65, 103001 [NASA ADS] [CrossRef] [Google Scholar]
 Cole, S., & Efstathiou, G. 1989, MNRAS, 239, 195 [NASA ADS] [CrossRef] [Google Scholar]
 Condon, J. J., Cotton, W. D., Greisen, E. W., et al. 1998, AJ, 115, 1693 [NASA ADS] [CrossRef] [Google Scholar]
 Creminelli, P., Nicolis, A., Senatore, L., Tegmark, M., & Zaldarriaga, M. 2006, JCAP, 0605, 004 [NASA ADS] [CrossRef] [Google Scholar]
 Das, S., Sherwin, B. D., Aguirre, P., et al. 2011, Phys. Rev. Lett., 107, 021301 [NASA ADS] [CrossRef] [Google Scholar]
 Das, S., Louis, T., Nolta, M. R., et al. 2013, JCAP, submitted [arXiv:1301.1037] [Google Scholar]
 Delabrouille, J., Betoule, M., Melin, J.B., et al. 2013, A&A, 553, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 De Zotti, G., Massardi, M., Negrello, M., & Wall, J. 2010, A&ARv, 18, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Efstathiou, G., & Bond, J. R. 1999, MNRAS, 304, 75 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Feng, C., Aslanyan, G., Manohar, A. V., et al. 2012a, Phys. Rev. D, 86, 063519 [NASA ADS] [CrossRef] [Google Scholar]
 Feng, C., Keating, B., Paar, H. P., & Zahn, O. 2012b, Phys. Rev. D, 85, 043513 [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]
 Goto, T., Szapudi, I., & Granett, B. R. 2012, MNRAS, 422, L77 [NASA ADS] [CrossRef] [Google Scholar]
 Gruppuso, A., Finelli, F., Natoli, P., et al. 2011, MNRAS, 411, 1445 [NASA ADS] [CrossRef] [Google Scholar]
 Hall, N., Knox, L., Reichardt, C., et al. 2010, ApJ, 718, 632 [NASA ADS] [CrossRef] [Google Scholar]
 Hanson, D., Rocha, G., & Gorski, K. 2009, MNRAS, 400, 2169 [NASA ADS] [CrossRef] [Google Scholar]
 Hanson, D., Lewis, A., & Challinor, A. 2010, Phys. Rev. D, 81, 103003 [NASA ADS] [CrossRef] [Google Scholar]
 Hanson, D., Challinor, A., Efstathiou, G., & Bielewicz, P. 2011, Phys. Rev. D, 83, 043005 [NASA ADS] [CrossRef] [Google Scholar]
 Hinshaw, G. et al. 2009, ApJS, 180, 225 [NASA ADS] [CrossRef] [Google Scholar]
 Hirata, C. M., & Seljak, U. 2003a, Phys. Rev. D, 67, 043001 [NASA ADS] [CrossRef] [Google Scholar]
 Hirata, C. M., & Seljak, U. 2003b, Phys. Rev. D, 68, 083002 [NASA ADS] [CrossRef] [Google Scholar]
 Hirata, C. M., Padmanabhan, N., Seljak, U., Schlegel, D., & Brinkmann, J. 2004, Phys. Rev. D, 70, 103501 [NASA ADS] [CrossRef] [Google Scholar]
 Hirata, C. M., Ho, S., Padmanabhan, N., Seljak, U., & Bahcall, N. A. 2008, Phys. Rev. D, 78, 043520 [NASA ADS] [CrossRef] [Google Scholar]
 Ho, S., Cuesta, A., Seo, H.J., et al. 2012, ApJ, 761, 14 [NASA ADS] [CrossRef] [Google Scholar]
 Howlett, C., Lewis, A., Hall, A., & Challinor, A. 2012, JCAP, 4, 27 [NASA ADS] [CrossRef] [Google Scholar]
 Hu, W. 2001, Phys. Rev. D, 64, 083005 [NASA ADS] [CrossRef] [Google Scholar]
 Hu, W., & Okamoto, T. 2002, ApJ, 574, 566 [NASA ADS] [CrossRef] [Google Scholar]
 Hütsi, G. 2010, MNRAS, 401, 2477 [NASA ADS] [CrossRef] [Google Scholar]
 Kamionkowski, M., & Knox, L. 2003, Phys. Rev. D, 67, 063001 [NASA ADS] [CrossRef] [Google Scholar]
 Kaplinghat, M., Knox, L., & Song, Y.S. 2003, Phys. Rev. Lett., 91, 241301 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Keisler, R., Reichardt, C., Aird, K., et al. 2011, ApJ, 743, 28 [NASA ADS] [CrossRef] [Google Scholar]
 Kesden, M. H., Cooray, A., & Kamionkowski, M. 2003, Phys. Rev. D, 67, 123507 [NASA ADS] [CrossRef] [Google Scholar]
 Kim, J., & Naselsky, P. 2010, ApJ, 714, L265 [NASA ADS] [CrossRef] [Google Scholar]
 Knox, L., & Song, Y.S. 2002, Phys. Rev. Lett., 89, 011303 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Koester, B. et al. 2007, ApJ, 660, 239 [NASA ADS] [CrossRef] [Google Scholar]
 Komatsu, E., Spergel, D. N., & Wandelt, B. D. 2005, ApJ, 634, 14 [NASA ADS] [CrossRef] [Google Scholar]
 Kovacs, A., Szapudi, I., Granett, B. R., & Frei, Z. 2013, MNRAS Letters, 431 [arXiv:1301.0475] [Google Scholar]
 Land, K., & Magueijo, J. 2005, Phys. Rev. D, 72, 101302 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Lewis, A. 2005, Phys. Rev. D, 71, 083008 [NASA ADS] [CrossRef] [Google Scholar]
 Lewis, A., & Challinor, A. 2006, Phys.Rept., 429, 1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lewis, A., Challinor, A., & Hanson, D. 2011, JCAP, 1103, 018 [NASA ADS] [CrossRef] [Google Scholar]
 Limber, D. N. 1954, ApJ, 119, 655 [NASA ADS] [CrossRef] [Google Scholar]
 Linder, E. V. 1990, MNRAS, 243, 353 [NASA ADS] [Google Scholar]
 Metcalf, R. B., & Silk, J. 1997, ApJ, 489, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Munshi, D., Heavens, A., Cooray, A., et al. 2011a, MNRAS, 412, 1993 [NASA ADS] [CrossRef] [Google Scholar]
 Munshi, D., Valageas, P., Cooray, A., & Heavens, A. 2011b, MNRAS, 414, 3173 [NASA ADS] [CrossRef] [Google Scholar]
 Namikawa, T., Yamauchi, D., & Taruya, A. 2012, JCAP, 1201, 007 [NASA ADS] [CrossRef] [Google Scholar]
 Namikawa, T., Hanson, D., & Takahashi, R. 2013, MNRAS, 431, 609 [NASA ADS] [CrossRef] [Google Scholar]
 Okamoto, T., & Hu, W. 2003, Phys. Rev. D, 67, 083002 [NASA ADS] [CrossRef] [Google Scholar]
 Osborne, S. J., Hanson, D., & Doré, O. 2014, JCAP, 03, 024 [NASA ADS] [CrossRef] [Google Scholar]
 Page, L., Hinshaw, G., Komatsu, E., et al. 2007, ApJS, 170, 335 [NASA ADS] [CrossRef] [Google Scholar]
 Perotto, L., Lesgourgues, J., Hannestad, S., Tu, H., & Wong, Y. Y. Y. 2006, JCAP, 0610, 013 [NASA ADS] [CrossRef] [Google Scholar]
 Perotto, L., Bobin, J., Plaszczynski, S., Starck, J. L., & Lavabre, A. 2010, A&A, 519, A14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration VII. 2011, A&A, 536, A7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XVIII. 2011, A&A, 536, A18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXIV. 2011, A&A, 536, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration Int. VII. 2013, A&A, 550, A133 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration 2013, The Explanatory Supplement to the Planck 2013 results, http://www.sciops.esa.int/wikiSI/planckpla/index.php?title=Main_Page (ESA) [Google Scholar]
 Planck Collaboration I. 2014, A&A, 571, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration II. 2014, A&A, 571, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration III. 2014, A&A, 571, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration IV. 2014, A&A, 571, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration V. 2014, A&A, 571, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration VI. 2014, A&A, 571, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration VII. 2014, A&A, 571, A7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration VIII. 2014, A&A, 571, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration IX. 2014, A&A, 571, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration X. 2014, A&A, 571, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XI. 2014, A&A, 571, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XII. 2014, A&A, 571, A12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XIII. 2014, A&A, 571, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XIV. 2014, A&A, 571, A14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XV. 2014,A&A, 571, A15 [Google Scholar]
 Planck Collaboration XVI. 2014, A&A, 571, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XVII. 2014, A&A, 571, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XVIII. 2014, A&A, 571, A18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XIX. 2014, A&A, 571, A19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XX. 2014, A&A, 571, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXI. 2014, A&A, 571, A21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXII. 2014, A&A, 571, A22 [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]
 Planck Collaboration XXV. 2014, A&A, 571, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXVI. 2014, A&A, 571, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXVII. 2014, A&A, 571, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXVIII. 2014, A&A, 571, A28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXIX. 2014, A&A, 571, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXX. 2014, A&A, 571, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXXI. 2014, A&A, 571, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Plaszczynski, S., Lavabre, A., Perotto, L., & Starck, J.L. 2012, A&A, 544, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Rees, M. J., & Sciama, D. W. 1968, Nature, 217, 511 [NASA ADS] [CrossRef] [Google Scholar]
 Regan, D., Shellard, E., & Fergusson, J. 2010, Phys. Rev. D, 82, 023520 [NASA ADS] [CrossRef] [Google Scholar]
 Reichardt, C. L., Shaw, L., Zahn, O., et al. 2012, ApJ, 755, 70 [NASA ADS] [CrossRef] [Google Scholar]
 Ross, A. J., Ho, S., Cuesta, A. J., et al. 2011, MNRAS, 417, 1350 [NASA ADS] [CrossRef] [Google Scholar]
 Sachs, R. K., & Wolfe, A. M. 1967, ApJ, 147, 73 [NASA ADS] [CrossRef] [Google Scholar]
 Schmittfull, M. M., Challinor, A., Hanson, D., & Lewis, A. 2013, Phys. Rev. D, 88 [Google Scholar]
 Seljak, U. 1996, ApJ, 463, 1 [NASA ADS] [CrossRef] [Google Scholar]
 Sherwin, B. D., Das, S., Hajian, A., et al. 2012, Phys. Rev. D, 86, 083006 [NASA ADS] [CrossRef] [Google Scholar]
 Sherwin, B. D., Dunkley, J., Das, S., et al. 2011, Phys. Rev. Lett., 107, 021302 [NASA ADS] [CrossRef] [Google Scholar]
 Smidt, J., Cooray, A., Amblard, A., et al. 2011, ApJ, 728, L1 [NASA ADS] [CrossRef] [Google Scholar]
 Smith, K. M., Hu, W., & Kaplinghat, M. 2006, Phys. Rev. D, 74, 123002 [NASA ADS] [CrossRef] [Google Scholar]
 Smith, K. M., Zahn, O., Dore, O., & Nolta, M. R. 2007, Phys. Rev. D, 76, 043510 [NASA ADS] [CrossRef] [Google Scholar]
 Stompor, R., & Efstathiou, G. 1999, MNRAS, 302, 735 [NASA ADS] [CrossRef] [Google Scholar]
 Story, K., Reichardt, C., Hou, Z., et al. 2013, ApJ, 779, 86 [NASA ADS] [CrossRef] [Google Scholar]
 Sunyaev, R. A., & Zeldovich, I. B. 1980, ARA&A, 18, 537 [NASA ADS] [CrossRef] [Google Scholar]
 Tinker, J. L., Robertson, B. E., Kravtsov, A. V., et al. 2010, ApJ, 724, 878 [NASA ADS] [CrossRef] [Google Scholar]
 van Engelen, A., Keisler, R., Zahn, O., et al. 2012, ApJ, 756, 142 [NASA ADS] [CrossRef] [Google Scholar]
 Wright, E. L., Eisenhardt, P. R., Mainzer, A., et al. 2010, Astron.J., 140, 1868 [NASA ADS] [CrossRef] [Google Scholar]
 Zaldarriaga, M., & Seljak, U. 1998, Phys. Rev. D, 58, 023003 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Null tests
In this appendix we present several generic null tests of our lens reconstruction. In Appendix A.1 we perform tests using lens reconstructions on the two individual sixmonth fullsky surveys contained in the first Planck data release. In Appendix. A.2 we perform the curlmode null test which is standard practice in Galaxy lensing. Finally, in Appendix A.3 we will perform several tests on the isotropy and Gaussianity of our MV lensing potential map.
Appendix A.1: Survey consistency
The Planck nominal mission consists of approximately two fullsky sixmonth surveys. Our fiducial analysis is based on maps made from the full data range, however we may also analyze the first two surveys separately to increase our sensitivity to effects such as beam asymmetry, correlated noise along the scan direction, and zodiacal light. In Fig. A.1 we plot several consistency tests exploiting the breakdown of the data into separate surveys, for both 143 GHz and 217 GHz. We test for gross contamination of either survey by estimating the lensing power spectrum on each individually and then difference the resulting , finding results consistent with the scatter expected from simple signal+noise simulations. To test the noise model more thoroughly, we also run lens reconstruction on a map constructed by taking the halfdifference of the two surveys, filtering it in the same way as the full nominalmission map. As can be seen in Fig. A.1, the power spectrum of this reconstruction null test is small compared to the measurement errors on . Finally, to isolate the effect of systematics which couple to the CMB and which flip sign between surveys (such as odd moments of the beam asymmetry), we take the crossspectrum of lens reconstruction estimates obtained from the halfdifference of the two survey maps with those obtained from the halfsum. Again, this null test shows residuals that are negligibly small compared to our measurement errors.
Fig. A.1
Survey consistency null tests at 143 GHz and 217 GHz. Solid black points show difference of estimated separately on the two sixmonth surveys of the nominal mission. Magenta points show the autospectrum of obtained from maps made by taking the halfdifference of the two surveys (S1−S2) × (S1−S2). Yellow points show the crossspectrum of obtained on halfdifference and halfsum maps (S1−S2) × (S1 + S2). The scatter on the tests that contain difference maps is very small because (as can be seen in Fig. 6) most of the lens reconstruction “noise” is due to primary CMB fluctuations, which is nulled in these tests. 
Fig. A.2
Lensing curl power spectrum estimates for the MV reconstruction (black hatched), as well as the individual 143 (green filled) and 217 GHz (blue filled) frequency reconstructions. 
Appendix A.2: Curlmode reconstruction
In addition to the standard gradient lensing potential φ sourced by largescale structure, we may also use the Planck maps to reconstruct a curlmode lensing field ψ. The weight function associated with the curl field is given by (Namikawa et al. 2012) (A.1)Using this weight function, we may construct an estimator for the curlmode power spectrum, completely analogous to that for φ. The response function of the curl estimator to point source shot noise is zero, and so we neglect this correction. Consistency of the reconstructed curl mode with zero is a standard systematic test in weak lensing, as the curl field due to gravitational lensing should be negligible (Hirata & Seljak 2003b). One caveat for this test in the CMB lensing context is that the variance of the curl mode reconstruction is still affected by the presence of gravitational lenses, by means of an “N^{(1)}”type (Kesden et al. 2003) contraction in the trispectrum (van Engelen et al. 2012). Our curlmode power spectrum estimates are plotted in Fig. A.2, after correction for the “N^{(1)}”type bias using Eq. (18). This bias is also given explicitly for the curl mode in BenoitLevy et al. (2013). We do not see evidence for anomalous power in the curl reconstruction, although we do find in Sect. C.4 that the curl estimate is contaminated by pixelization effects at a level too small to be seen in this plot.
Appendix A.3: Statistical tests at the map level
In this section we present several sanity tests of our lensing estimates at the map level.
Fig. A.3
Histogram of pixel values for the unnormalized lensing estimate for the MV reconstruction, in units of the standard deviation of the map. A Gaussian with σ = 1 is overplotted in black. Note that the yaxis is logarithmic. The left column shows the histogram for the data itself, while the right column shows a simulated reconstruction. The reconstructed map is slightly nonGaussian, but in a way that is expected from simulations given the nonGaussian nature of the quadratic lensing estimator. 
1point PDF of :
in the left panel of Fig. A.3 we plot the histogram, or 1point function, of our unnormalized potential estimate for the MV reconstruction, bandlimited to 10 ≤ L ≤ 2048. As can be seen from the histogram, the map is slightly nonGaussian, as might be expected for an estimator that is based on the product of two Gaussian random fields. In the right panel of Fig. A.3 we plot the same histogram for a simulated MV reconstruction, which shows a comparable level of nonGaussian structure. In Appendix D we show that this small nonGaussianity is negligible for our lensing likelihood.
Power spectra on patches:
to test for gross contamination of our reconstruction from any localized region of the sky, we use the patches method (described in Sect 7.5) to build an estimate of the lensing power spectra on a collection of approximately 400 10deg × 10deg tiles. We fill pointsource holes with local constrained Gaussian realizations of the CMB signal. The patches were chosen to have no overlap with the Galaxy mask.
Fig. A.4
Bidimensional deflection power spectrum obtained by stacking the individual Fourier amplitude obtained on patches, after subtraction and smoothing for the 143 and 217 GHz channels. 
Fig. A.5
Density plot of the individual deflection power spectra constructed on approximately 400 10deg × 10deg patches for the 143 GHz and 217 GHz channels. The lines are the theoretical 1σ and 2σ error for this patch size. No strong outliers are seen. 
The distribution of the power spectra estimated for each patch can be used as a diagnostic for deficiency in the mask and dust correction. Significantly contaminated but localized regions would appear as outliers. This procedure successfully detected several point sources that were overlooked in early versions of our source masks. In Fig. A.4 we present results of this analysis using our final Galaxy and point source masks for 143 GHz and 217 GHz.
Isotropy on patches:
we can also use the same data to build a coarse test of the isotropy of the reconstructed lensing map. We perform lens reconstruction on patches which are oriented along the local spherical basis (e_{θ},e_{φ}), and then stack the corresponding 2D power spectrum estimates to build a bidimensional power spectrum for the 143 GHz and 217 GHz maps. This is shown in Fig. A.5, where we have additionally applied smoothing based on an undecimated wavelet transform to reduce the visual noise. The overplotted contour lines are placed at levels given by the azimuthally averaged power spectrum for L = 100, 200, 300 and 400.
Appendix B: Filtering
In this section, we discuss our filtering methodology, that takes us from a set of maps indexed by frequency ν and pixel p to the multipoles that are fed into the quadratic estimators used to measure lensing and possible contaminants. For an observed CMB that consists of Gaussian signal and noise, the optimal filtering (in the minimumvariance sense) for the construction of any quadratic anisotropy estimator is “C^{1}” filtering, in which the observed map is beamdeconvolved and then multiplied by the inverse of a signal+noise covariance matrix. For the purposes of lens reconstruction, the Gaussian assumption above is quite accurate for Planck. The Planck map noise is well approximated as Gaussian (Planck Collaboration VI 2014). The observed CMB signal is lessobviously well approximated as Gaussian, with the nonGaussianity due to lensing detected at high significance, however for lens reconstruction at Planck noise levels the difference in signaltonoise using a filter constructed in the Gaussian limit versus a more complicated “delensing” filter that properly accounts for lensing nonGaussianity has been shown to be negligible (Hirata & Seljak 2003a). Other sources of nonGaussianity, such as unresolved point sources, are completely subdominant to either primary CMB or noise at all of the multipoles probed by our lens reconstruction, and so we neglect them in the construction of our filters.
For a single map, with a symmetric beam transfer function, fullsky coverage and homogeneous map noise levels, the C^{1} filter is given by (B.1)where gives the value of the spherical harmonic at the center of pixel p, B_{l} gives the beam transfer function, and F_{l} is given by (B.2)where σ is the map noise level in units of μK arcmin. The maps we use must always be masked to remove bright Galactic foregrounds and point sources, and so we do not ever use the fullsky filter above, although as we will discuss the filter that we do use asymptotes to this form in regions far from the mask, and the simple form F_{l} of the filter function given in Eq. (B.2)will prove useful in many analytical calculations.
For the more realistic case of inhomogeneous map noise levels, with pixelspace noise correlations and a sky cut, the construction of is more involved. To account for the sky cut, we can take the noise covariance to be infinite for masked pixels. Reformulating the C^{1} filter using the inverse of the noise covariance matrix in pixel space (defined to be zero if either p or p′ is a masked pixel), and allowing for the jointfiltering of multiple input maps (with no noise correlations between them) we have (B.3)where ν indexes the maps which are being combined and the pointing matrix gives the value of the spherical harmonic at the center of pixel p, convolved with the appropriate beam+pixel transfer function . The matrix is given by (B.4)The filter of Eq. (B.3)is the one that we will use, with a slightly suboptimal choice for the noise covariance matrix N_{pp′} that is described below.
The Planck map noise has several sources of pixelpixel correlations for p ≠ p′. On small scales, there are correlations along the scan direction induced by the deconvolution of the bolometer time response (Planck Collaboration VII 2014). On large scales, there are correlations between scan rings, due to residual 1 /f noise fluctuations not cancelled by the baseline subtraction procedure (Planck Collaboration VI 2014). Accounting for these pixel correlations would be computationally intractable, as the full covariance matrix N_{pp′} (an sized object) is too large to work with on a modern computer, and so we neglect them in our noise filter. The remaining structure in the noise matrix is due to the variable coverage in different regions of the sky. Noise levels are lowest near the Ecliptic poles, which are visited more frequently than the region around the Ecliptic plane due to the Planck scan strategy. Although the variable noise level is tractable to include in our filtering procedure (and is indeed included in the C^{1} filtering of e.g., Planck Collaboration XXIV 2014), we have chosen to neglect this aspect of the noise, using instead a fixed effective noise level. We use an inversenoise map in pixel space given by (B.5)where M_{p} is a sky mask and σ_{ν} is an average map noise level in units of μK arcmin. In addition to this noise weighting, we also project out the five modes corresponding to the monopole, dipole, and 857 GHz Planck map (as a simple dust template). The noise level for each of these modes is taken to infinity using the Woodbury formula; putting each of these modes into a map indexed by i, our total noise covariance matrix is given in terms of the inverse noise map above by (B.6)where is the overlap matrix between templates given by (B.7)Our motivation for taking a fixed noise level is that with this approach, in regions sufficiently far from the mask boundary, our filter asymptotes to the diagonal form of Eq. (B.1). This means that the normalization of our lensing estimates can be wellapproximated analytically, which is very useful for the propagation of systematic effects, and also that the normalization of our lensing estimates does not vary across the sky with noise level, which simplifies crosscorrelation analysis. Our C^{1} filter is therefore optimally accounting for masking effects, but not for noise correlations and inhomogeneity. We estimate the suboptimality of neglecting these noise properties by calculating the quantity (B.8)where is the optimal filter and is the suboptimal filter that we have actually used. This equation gives the S/N loss as a function of lens multipole L, however in practice we find that the L dependence is small enough that it suffices to quote a single average loss. To estimate the degradation due to ignoring noise correlations we set (B.9)where is the power spectrum of the map noise. We find that the degradation due to neglect of noise correlations is small; less than 2% for all L ≤ 2048 at 100 GHz, and less than 0.1% at 143 and 217 GHz. To calculate the degradation due to ignoring noise inhomogeneity, we determine the map noise level in the 3072 regions corresponding to N_{side} = 16 HEALPix pixels, take using Eq. (B.2)with the local noise level, and estimate a resulting S/N degradation using Eq. (B.8). The neglect of noise inhomogeneity is the dominant suboptimality of our filtering, although it is still small. We find an average S/N loss (averaged over the entire sky) of approximately 4% at 100, 143, and 217 GHz, consistent with the simulationbased results of Hanson et al. (2009). We take this loss as justified, given the simpler normalization properties of our lensing estimates when neglecting variations in the map noise level.
Appendix C: Meanfields
Fig. C.1
Analytical estimates for the power spectra of the largest lowL meanfields at 143 GHz. The various components are discussed in more detail in Sect. C.1 (mask), Sect. C.2 (noise), and Sect. C.3 (beams). The meanfields all couple most strongly to even modes of the lens reconstruction, due to the approximate north/south symmetry of the scan strategy and Galactic mask. 
As discussed in Sect. 2, the quadratic lensing estimators that we use are designed to detect statistical anisotropy induced by lensing. There are a number of nonlensing sources of statistical anisotropy that can mimic the lensing signal to some extent. In our analysis, the effects we consider are

1.
The application of a sky mask, which introduces sharp gradients that may be misinterpreted as lensing.

2.
Noise inhomogeneity, which causes the overall power to fluctuate across the sky and can resemble the convergence component of lensing.

3.
Beam asymmetry, which smears the fluctuations more along one direction than another and can mimic the shear component of lensing.

4.
Pixelization, in which detector samples are accumulated into pixels, introduces a spurious deflection field on the pixel scale because the centroid of the hit distribution in each pixel does not necessarily lie at the pixel center.
In our analysis, we account for most of these effects with a corrective meanfield term, given by Eq. (9), which is determined using MonteCarlo simulations. In this appendix, we will break this meanfield down into its constituent parts and discuss each in more detail. As an overview of the results in this section, in Fig. C.1 we plot estimates for the three largest meanfields, due to masking, noise inhomogeneity, and beam asymmetry at 143 GHz (100 and 217 GHz are qualitatively similar). These meanfields all have most of their contributions on very large scales, dictated by the coherency of the scan strategy in the case of beam asymmetry and noise inhomogeneity, and of the largescale nature of the Galactic foregrounds in the case of the sky mask. We note also that the mask meanfield is concentrated near the edges of the Galaxy cut, and so for our power spectrum analysis (where we apodize the estimated φ map, removing the contributions near the edges of the mask), the effective meanfield is significantly reduced.
Our discussion will focus on constructing simple models for each source of meanfield. Following Hanson et al. (2010), we will identify each of the individual contributions to the meanfield with a tracer z_{LM} that sources a contribution to the CMB covariance matrix given by (C.1)where is a weight function describing how z_{LM} couples multipoles. Such a contaminant leads to a bias for the standard lensing estimator given by (C.2)where the response function is defined in Eq. (12). The analytical forms for the meanfields that we present here are used in Sect. 7.4 to construct “bias hardened” estimators that have smaller meanfield corrections, as proposed by Namikawa et al. (2013).
Appendix C.1: Mask
The application of a mask to the beamconvolved sky, followed by subsequent beam deconvolution, can be modelled as a statistical anisotropy “M”, with and weight function (Namikawa et al. 2013) (C.3)The mask meanfield using our fiducial Galaxy and pointsource analysis mask (described in Sect. 3) is plotted in Fig. C.1. At very lowL, the mean field induced by this mask can be much larger than the reconstruction noise level. At L = 2 for example, it can be up to 10^{4} times greater than the reconstruction noise level, although the amplitude of the meanfield is a very steep function of L.
The mask meanfield is proportional to , and so with percentlevel agreement between the power spectrum of our simulations we use to determine the meanfield and the data power spectrum (discussed in Sect. 4), we can expect errors on the power spectra of the mask meanfield at the level of one part in 10^{4}. At L = 2 for 143 and 217 GHz this corresponds to a level of , comparable to the lensing potential power spectrum. There are several good ways to mitigate the effect of the mask meanfield. By apodizing the mask before lens reconstruction, the gradients at the mask boundary can be greatly reduced, which significantly suppresses the meanfield (van Engelen et al. 2012; BenoitLevy et al. 2013). At the map level, the only issue with such an apodization for our purposes is that it results in a positiondependent normalization which can complicate crosscorrelation analyses. A related approach, which we have taken in our power spectrum analysis, is to apodize the lens reconstruction after performing the reconstruction. This amounts to removing the regions of near the mask boundary where the meanfield has all of its features. This “avoidance” approach can also be used easily in crosscorrelation analyses, as in Hirata et al. (2004). A third approach is that of Namikawa et al. (2013), who show that it is possible to construct “biashardened” lensing estimators which have smaller meanfield contributions. We have already shown consistency tests based on this approach in Sect. 7.4. In the lowL biashardened comparisons of Fig. 20, we saw discrepancies with the standard estimator of and 7 × 10^{7}, precisely the level expected above.
For our publicly released lensing map, we have chosen simply to highpass filter our map to L> 10, above which we believe that meanfield errors are subdominant to the reconstruction noise.
Appendix C.2: Instrumental noise
The instrumental noise of Planck is rich in statisticallyanisotropic structure, as discussed in Planck Collaboration VI (2014), which is sourced by the following

1.
Striping effects, due to residual noise correlations on long time scales after destriping.

2.
Correlations along the scan direction induced by the deconvolution of the detector time response and the lowpass rolloff filter.

3.
Large variations in the pixel noise level due to the larger density of pixel hits near the Ecliptic poles.
All of these effects can be avoided by forming estimators using pairs of data with independent noise realizations, as we have done in Sect. 7.3, although our baseline results do not use such crosscorrelations. Here we discuss the size of the noise meanfield. Items (1) and (2) are difficult to study analytically. For white noise with a variance in each pixel given by , item (3) can be identified with a statistical anisotropy “N” defined by Hanson et al. (2009)(C.4)Note that the shape of this bias is dictated by the hit count distribution, which is known quite well, rather than the noise level, which is more uncertain. The shape of the noise meanfield is plotted in Fig. C.1. As with the mask, it is large compared to our signal and noise but falls off quickly with multipole. For L< 10, given that there are known issues with our noise modelling at the 5% level in power (as discussed in Sect. 4), it can be important to construct quadratic estimators from correlations of maps with independent noise realizations (apart from our biashardened estimator tests in Sect. 7.4, this is also done for L = 1 in Planck Collaboration XXVII 2014), however at the L> 10 multipoles of our lensing map and power spectrum, uncertainties in the noise meanfield should be negligible.
Appendix C.3: Beam asymmetry
Planck observes the sky after convolution with a “scanning beam”, which captures its effective response to the sky as a function of displacement from the nominal pointing direction. Asymmetry in the shape of this beam (for example ellipticity) means that the fluctuations in the map are smeared depending on the beam orientation, which can mimic the shearing effects of gravitational lensing. Here we present semianalytical estimates of this beam asymmetry meanfield, and discuss how simulations of beam asymmetry are included in our simulations to correct for it. As with the noise inhomogeneity, we will see that the largescale coherency of the Planck scan strategy means that beam asymmetries only mimic largescale modes of φ_{LM}.
Decomposing the scanning beam into harmonic coefficients B_{ℓm}, each timeordered data (TOD) sample can be modelled as (neglecting the contribution from instrumental noise, which is independent of beam asymmetry) (C.5)where the TOD samples are indexed by i, and is the underlying sky signal. The spin spherical harmonic rotates the scanning beam to the pointing location (θ,φ), while the e^{− isαi} factor gives it the correct orientation. On the small scales that lensing is sensitive to, it is a good approximation to assume that the procedure of mapmaking from TOD samples is just a process of binning: (C.6)where is the centroid of pixel p, and the sum is taken over all hits assigned to pixel p. H_{p} is the number of hits in pixel p. Here we will assume that the pointing for each observation of pixel p is identical to the pixel center . The purpose of this assumption is simply to isolate the effect of beam asymmetry from the (small) effect of pixelization. We will discuss the pixelization effect in Sect. C.4. Placing all hits at pixel centers, combining Eq. (C.5)and Eq. (C.6)we have (C.7)where the “scan strategy maps” are given by (C.8)The scan strategy maps are spins objects, with corresponding harmonic decomposition (C.9)We use Eq. (C.7)to include the effect of asymmetric beams in the simulations that we use to debias our lensing estimates. We can also propagate the effect of beam asymmetry to the lensing meanfield directly. Each mode s of the scan strategy can be identified with a statistical anisotropy “b^{s}” of the type defined in Eq. (C.1), with and (C.10)In Fig. C.1 we have plotted the expected fullsky bias for the bestfit Planck beams, which are measured from planet observations in Planck Collaboration VII (2014). There are considerable differences in the meanfield between frequencies, as well as the steep scale dependence, which is dictated by the largescale coherency of the scan strategy.
Appendix C.4: Pixelization
In the process of mapmaking, TOD samples are accumulated into bins defined by HEALPix pixels. The Planck HFI maps are generated at HEALPix resolution 11, with a typical size of , which is sufficiently close to the scale of deflections expected from gravitational lensing effects that this remapping of the observed CMB is a source of concern. To model this effect, we begin from Eq. (C.5), but treat the beam as symmetric to isolate the effect of pixelization from that of beam asymmetry. Binning hits into pixels, we then have (C.11)where is the center of pixel p and in the second line we have made a gradient approximation to about the pixel center, and used the effective deflection vector (C.12)where is the pointing for observation i, which is assigned to pixel p (with center at and number of hits H_{p}). Implicit in this model is an assumption that on the arcminute pixelization scale the beamconvolved sky fluctuations are quite smooth, and can be modelled as a gradient. Pixelization therefore introduces an effective deflection field, given by d^{pix}, which biases our reconstruction of the lensing potential. We could in principle correct for the spurious lensing field which this introduces by including the remapping of Eq. (C.11)in the simulations used to estimate the meanfield correction, however as we will see it is quite small. We can decompose the spin1 deflection field as (C.13)where d_{θ} and d_{φ} are the components of d^{pix} along the θ and φ directions respectively. The associated weight function, following Eq. (C.1), is given by
Fig. C.2
Power spectra of the expected bias to our lensing estimators due to pixelization (thick solid lines), as well as crosscorrelation of the pixelization field with the estimators themselves (binned points). The theoretical lensing potential power spectrum is plotted for gradient modes as solid black, while the lower dashed black curves give respectively. The pixelization bias is always small compared to both signal and reconstruction noise. 
(C.14)The pixelization deflection field can also be usefully split into gradient and curltype parts, as (C.15)where the gradient and curl parts both satisfy the reality condition (C.16)Only the gradient mode couples to the lensing estimator, while the curl part only couples to the curlmode lensing estimator (discussed in Sect. A.2). The statistical anisotropy introduced by is almost completely degenerate with that of CMB lensing, differing only by the fact that it operates on the beamconvolved CMB sky rather than the primordial one. In Fig. C.2 we plot the expected bias to the lensing power spectrum estimator, which is given by (C.17)as well as the analogous bias for the curl nulltest estimator. To test for the presence of the pixelization lensing field in our reconstructions, we crosscorrelate the deflection fields with our lensing estimates , . In Fig. C.2 we plot (C.18)and similarly for ψ. We detect both and at greater than 15σ significance for all of our three main frequency bands, at the expected level. This detection is only possible because we know d_{LM} precisely from the satellite pointing. Except at the very highest multipoles, the gradient mode of the deflection lensing field is always significantly smaller than the CMB lensing potential. At L< 250, which provides most of the weight for the Planck lensing likelihood, the bias to the lensing potential is significantly less than 0.1%, while our statistical measurement uncertainty is at the 4% level.
Appendix D: Likelihood validation
In the construction of our lensing likelihood in Sect. 2.3 we have used several simplifying approximations to the statistics of the reconstructed lensing potential, which we discuss briefly here.
Uncorrelated binned power spectrum estimates: in Fig. D.1 we plot MonteCarlo estimates of the bandpower correlation. These plots do not include the effect of beam uncertainties, errors in the pointsource correction or cosmological uncertainty in the subtraction of the N^{(1)} bias, which we include in our likelihood covariance matrix analytically. Even for adjacent bandpowers, the correlation coefficient is constrained to be less than 10% and we ignore such correlations in constructing the likelihood. As discussed in Sect. 5, however, we do include bandpower correlations from beam uncertainties and the pointsource and N^{(1)} correction in the likelihood.
Fig. D.1
Measured correlation coefficients for the bins of Table 1, estimated from 1000 simulated lens reconstructions. There is some evidence for correlations along the diagonal at the 10% level, which is sufficiently small that we have chosen to neglect binbin correlations in our likelihood. 
Gaussianity of the binned power spectrum estimates: our estimator for the lensing potential power spectrum is a sum over the trispectrum of the observed CMB. There is no simple, rigorous argument to expect our bandpower estimates to be Gaussian distributed, however keeping in mind that we are collapsing the information in a multimillion pixel map down to eight eigenmodes, it is quite reasonable to expect that the central limit theorem will drive their distribution in this direction. With a set of simulated lens reconstructions, it is straightforward to assess the Gaussianity of our amplitude estimates. In Fig. D.2 we plot histograms of the bandpower coefficients from our set of fiducial simulations showing that the Gaussian approximation is good for the marginal distributions. Further tests of the assumption of a Gaussian likelihood are made in Schmittfull et al. (2013) for isotropic surveys. There it is shown that parameter constraints for a simple twoparameter model of scatter across simulations in a manner consistent with the likelihood width in any realization.
Fig. D.2
Histograms of amplitude coefficients for the bandpowers of Table 1, taken from 1000 simulated lens reconstructions. Each coefficient has been rescaled by its standard deviation so that it can be directly compared to a unitvariance Gaussian distribution (solid lines). 
No correlations between the lensing potential and temperature power spectra: lensing may be detected at a significance of 10σ as a smearing of the acoustic peaks in the temperature power spectrum; see Planck Collaboration XVI (2014). A potential worry when simply combining the temperature likelihood with the lensing power spectrum likelihood is that we are double counting the same information (Lewis 2005; Perotto et al. 2006; Smith et al. 2006) since we do not account for the correlated effect of sample variance of the lenses in the estimated and spectra. This effect is explored in detail in Schmittfull et al. (2013) where it is shown to derive from the connected part of the 6point function for the lensed temperature anisotropies. They estimate for Planck that the maximum correlation between and for unbinned spectra is around 0.05%, and that ignoring this leads to a misestimation of the errors on the lensing amplitude A_{L} (see Sect. 6.1) in a joint analysis of only 2%. A more pressing concern is that the statistical errors on the reconstruction power are correlated with those for the temperature power spectrum due to sample variance of the CMB anisotropies. This effect is also discussed in detail in Schmittfull et al. (2013). The correlation derives from the fully disconnected part of the CMB 6point function and leads to correlations of at most 0.2% between the unbinned spectra. However, these correlations are removed exactly by the realizationdependent aspect of the N^{(0)} debiasing procedure discussed in Sect. 2.2.
To see that the disconnected contribution to vanishes, we note that with our datadependent N^{(0)} subtraction the data enters the in the form^{7}(D.1)Here, is the inversevariance filtered temperature, with covariance , and summation over repeated indices is assumed. Equation (D.1) arises naturally for optimal trispectrum estimation of weakly nonGaussian fields (e.g., Regan et al. 2010). The expectation value of the data combination on the righthand side is simply the connected 4point function, ⟨ T_{i}T_{j}T_{k}T_{l} ⟩ _{c}, but the datadependent quadratic terms improve the variance of the estimator. The righthand side of Eq. (D.1) is explicitly symmetric under regroupings of pairs of indices so that only the primary coupling of the trispectrum need be included. For lensing, the primary coupling factors such that (D.2)see Hu (2001); Hanson et al. (2011) for explicit expressions. The first set of terms in Eq. (D.1) combine to give (D.3)where the terms in square brackets are exactly of the form of meanfieldsubtracted, unnormalized estimates of φ_{LM} and φ_{L − M}. The remaining terms are equivalent to the N^{(0)} subtraction in Eq. (17). Any quadratic estimate of the temperature power spectrum will involve a linear combination of terms of the form , where and may be different noisy estimates of the temperature anisotropies but are correlated with themselves and T_{i} through (at least) the common temperature anisotropies. If we form the covariance between and the data combination in Eq. (D.1) that enters the estimate, the disconnected terms give the contribution to the covariance between the estimated lensing and temperature power spectra due to sample variance of the temperature anisotropies. However, it is straightforward to show that by including the additional terms in Eq. (D.1) that are quadratic in the disconnected contribution vanishes. We see that, as well as mitigating against errors in the CMB and noise covariance matrices (see Sect. 2), and reducing the covariance between estimates (Hanson et al. 2011), the datadependent N^{(0)} subtraction that we adopt removes much of the covariance with the temperature power spectrum.
The final effect that correlates the lensing and temperature power spectrum estimates is the correlation between T and φ due to the latetime ISW effect in the former (Sect. 6.2). This produces only local correlations between the spectra and they fall rapidly with L being less than 0.5% at L = 40. Their impact on the errors on the lensing amplitude A_{L} in a joint estimate from the temperature and lensing likelihoods is negligible.
All Tables
All Figures
Fig. 1
Skyaveraged lens reconstruction noise levels for the 100, 143, and 217 GHz Planck channels (red, green, and blue solid, respectively), as well as for experiments that are cosmicvariance limited to a maximum multipole ℓ_{max} = 1000, 1500, and 1750 (upper to lower solid grey lines). A fiducial ΛCDM lensing potential using bestfit parameters to the temperature power spectrum from Planck Collaboration XVI (2014) is shown in dashed black. The noise level for a minimumvariance (“MV”) combination of 143 + 217 GHz is shown in black (the gain from adding 100 GHz is negligible). 

In the text 
Fig. 2
Overview of forecasted contributions to the detection significance as a function of lensing multipole L for the power spectrum (solid black), as well as for several other mass tracers, at the noise levels of our MV lens reconstruction. Our measurement of the power spectrum is presented in Sect. 6, The ISWφ correlation believed to be induced by dark energy is studied in Sect. 6.2. The NVSSφ correlation is studied (along with other galaxy correlations) in Sect. 6.3. The CIBφ prediction (dashed cyan) uses the linear SSED model of Hall et al. (2010), assuming no noise or foreground contamination. A full analysis and interpretation of the CIBφ correlation is performed in Planck Collaboration XVIII (2014). 

In the text 
Fig. 3
Validation of our estimator normalization for simulations of the MV reconstruction at the map and power spectrum levels. The map normalization (plotted as ) is tested by taking the crossspectrum of the input φ with the reconstruction averaged over MonteCarlo simulations, divided by an f_{sky} factor to account for missing power in the mask. The power spectrum normalization (plotted as ) is obtained by averaging the first line of Eq. (15)over simulations, and then comparing it to the expected value, which is because our simulations do not contain pointsource nonGaussianity. 

In the text 
Fig. 4
Simulation of the Wienerfiltered lensing potential estimate for the MV reconstruction (left), and the input φ realization (right; filtered by to be directly comparable to the Wiener estimate). Both maps show the southern Galactic sky in orthographic projection. The lensing reconstruction is noise dominated on all scales, however correlations between the two maps can still be seen visually. 

In the text 
Fig. 5
Relative contribution of statistical measurement uncertainty (grey band, representing the 1σ uncertainty on for the plotted bins) and the systematic errors we assign due to uncertainty in the N^{(1)} bias correction, pointsource shotnoise correction, and the beam transfer function. The contributions represent error eigenmodes in , and so are completely correlated between L. Note that these contributions are much smaller than our measurement uncertainty, and they have been multiplied by a factor of 10 for clarity. 

In the text 
Fig. 6
Relative contributions of CMB, instrumental noise, and foreground power terms discussed in Sect 5.1 to the approximate lens reconstruction variance of Eq. (21), previously plotted in Fig. 1. As discussed in Sect. 5.1, the finite number of CMB modes observed by Planck is the dominant source of variance for the lens reconstruction. 

In the text 
Fig. 7
Propagation of eigenmodes for beam transfer function uncertainty (lefthand panels) to an uncertainty in the lensing normalization (righthand panels). The normalization eigenmodes are flat, with amplitudes given by approximately twice the beam uncertainty at the beam and noise cutoff scale (at ℓ ≈ 1000 for 100 GHz and ℓ ≈ 1500 for 143 and 217 GHz). 

In the text 
Fig. 8
Wienerfiltered lensing potential estimate for our MV reconstruction, in Galactic coordinates using orthographic projection. The reconstruction is bandpass filtered to L ∈ [ 10,2048 ]. The Planck lens reconstruction has S/N ≤ 1 for individual modes on all scales, so this map is noise dominated. Comparison between simulations of reconstructed and input φ in Fig. 4 show the expected level of visible correlation between our reconstruction and the true lensing potential. 

In the text 
Fig. 9
Wienerfiltered lensing potential estimates, as in Fig. 8, for the individual 143 and 217 GHz maps. Note that the “noise” due to CMB fluctuations is correlated between these two estimates. 

In the text 
Fig. 10
Lensing potential power spectrum estimates based on the individual 100, 143, and 217 GHz sky maps, as well our fiducial minimumvariance (MV) reconstruction which forms the basis for the Planck lensing likelihood. The black line is for the bestfit ΛCDM model of Planck Collaboration XVI (2014). 

In the text 
Fig. 11
Replotting of Fig. 10, removing 100 GHz for easier comparison of 143 and 217 GHz. Also plotted are the SPT bandpowers from van Engelen et al. (2012), and the ACT bandpowers from Das et al. (2013). All three experiments are very consistent. The lower panel shows the difference between the measured bandpowers and the fiducial bestfit ΛCDM model. 

In the text 
Fig. 12
Upper left: Planck measurements of the lensing power spectrum compared to the ΛCDM mean prediction and 68% confidence interval (dashed lines) for models fit to Planck+WP+highL (see text). The eight bandpowers are those used in the Planck lensing likelihood; they are renormalized, along with their errors, to account for the small differences between the lensed in the bestfit model and the fiducial model used throughout this paper. The error bars are the ± 1σ errors from the diagonal of the covariance matrix. The colour coding shows how varies with the optical depth τ across samples from the ΛCDM posterior distribution. Upper right: as upperleft but using only the temperature power spectrum from Planck. Lower left: as upperleft panel but in models with spatial curvature. The colour coding is for Ω_{K}. Lower right: as upperleft but in models with three massive neutrinos (of equal mass). The colour coding is for the summed neutrino mass ∑ m_{ν}. 

In the text 
Fig. 13
Marginalised constraints on the optical depth in ΛCDM models from the Planck temperature power spectrum (Planck; solid black), and additionally including the lensing likelihood (Planck+lensing; dashed red) or WMAP polarization (Planck+WP; dasheddotted blue). We use a prior τ> 0.01 in all cases. 

In the text 
Fig. 14
Marginalised posteriors for the sixparameter ΛCDM model, shown as box plots, for Planck+WP+highL with various lensing likelihoods. The red and blue lines are the median and mean, respectively. The box and bar correspond to 68% and 95% of the probability density, both centered on the median. The leftmost column is without the lensing likelihood and the median of these constraints is shown by the grey line. The remaining columns show the effect of adding in the fiducial lensing likelihood (second column), and further adding a lowL bin (third column), highL bins (fourth column), or both (final column). 

In the text 
Fig. 15
Two views of the geometric degeneracy in curved ΛCDM models which is partially broken by lensing. Left: the degeneracy in the Ω_{m}Ω_{Λ} plane, with samples from Planck+WP+highL colour coded by the value of H_{0}. The contours delimit the 68% and 95% confidence regions, showing the further improvement from including the lensing likelihood. Right: the degeneracy in the Ω_{K}H_{0} plane, with samples colour coded by Ω_{Λ}. Spatiallyflat models lie along the grey dashed lines. 

In the text 
Fig. 16
LensingISW bispectrumrelated cross spectra computed from Eq. (46). Black dashed lines indicate the average value for simulations, while dark/light grey filled regions indicate the expected one/two standard deviation scatter, also measured from simulations. The thin magenta line gives the expected crossspectrum for our fiducial model. The agreement of this curve with the simulation average illustrates that our estimator is accurately normalized. In all the quantitative analysis of this section we ignore L< 10, although we have plotted the crossspectra at these multipoles for interest. 

In the text 
Fig. 17
Crossspectra of the Planck MV lensing potential with several galaxy catalogues, scaled by the signaltonoise weighting factor defined in Eq. (53). Crosscorrelations are detected at approximately 20σ significance for the NVSS quasar catalogue, 10σ for SDSS LRGs, and 7σ for both MaxBCG and WISE. 

In the text 
Fig. 18
Summary of internal consistency tests between our fiducial minimumvariance (MV) reconstruction and a set of alternatives designed to test sensitivity to potential issues. The top panel shows estimates, with measurement error bars. The bottom panels show the residual with respect to the MV reconstruction in units of the MV measurement uncertainty. The grey band marks the 1σ deviation uncertainty of the MV reconstruction. The error bar on each data point in the lower panels gives the standard deviation of the scatter between each result and the MV, determined from MonteCarlo simulations that account for the correlated CMB, noise and foreground power between estimators. Comparison of the uncertainty on the scatter points and the grey band gives an indication of how constraining each test is. The various tests are described in more detail in subsections of Sect. 7. 

In the text 
Fig. 19
Wienerfiltered potential maps in Galactic coordinates, as in Fig. 8, plotted here in Mollweide projection. Top: MV reconstruction; Bottom: extended reconstruction on the Nilc componentseparated map. 

In the text 
Fig. 20
Biashardened estimator consistency tests for our MV reconstruction. Following the discussion in Sect. 7.4 we form estimators that are biashardened against the meanfields due to inhomgeneous noise levels (φN, cyan), masking (φM, magenta), and both of these effects simultaneously (φMN, orange). The top panel shows the power spectra of the meanfield for each of these estimators. The middle panel shows the raw power spectra of the reconstructions themselves, compared to the average expected from signal+noise (dashed black). The lower panel shows the power spectrum of the difference between each biashardened estimator and the standard result. Dashed black lines give the expected average of this difference (measured on simulations). The grey line is the fiducial power spectrum. We can see that the differences are consistent with the expected scatter for all estimators at L ≥ 10. 

In the text 
Fig. 21
Comparison of alternative lensing pipelines. The Baseline results use the methodology of Sect. 2. Boxes are for the MV reconstruction, circles show the 143 GHz results, and triangles show the 217 GHz results. The two bottom panels show differences relative to the MV result for 143 and 217 GHz. 

In the text 
Fig. A.1
Survey consistency null tests at 143 GHz and 217 GHz. Solid black points show difference of estimated separately on the two sixmonth surveys of the nominal mission. Magenta points show the autospectrum of obtained from maps made by taking the halfdifference of the two surveys (S1−S2) × (S1−S2). Yellow points show the crossspectrum of obtained on halfdifference and halfsum maps (S1−S2) × (S1 + S2). The scatter on the tests that contain difference maps is very small because (as can be seen in Fig. 6) most of the lens reconstruction “noise” is due to primary CMB fluctuations, which is nulled in these tests. 

In the text 
Fig. A.2
Lensing curl power spectrum estimates for the MV reconstruction (black hatched), as well as the individual 143 (green filled) and 217 GHz (blue filled) frequency reconstructions. 

In the text 
Fig. A.3
Histogram of pixel values for the unnormalized lensing estimate for the MV reconstruction, in units of the standard deviation of the map. A Gaussian with σ = 1 is overplotted in black. Note that the yaxis is logarithmic. The left column shows the histogram for the data itself, while the right column shows a simulated reconstruction. The reconstructed map is slightly nonGaussian, but in a way that is expected from simulations given the nonGaussian nature of the quadratic lensing estimator. 

In the text 
Fig. A.4
Bidimensional deflection power spectrum obtained by stacking the individual Fourier amplitude obtained on patches, after subtraction and smoothing for the 143 and 217 GHz channels. 

In the text 
Fig. A.5
Density plot of the individual deflection power spectra constructed on approximately 400 10deg × 10deg patches for the 143 GHz and 217 GHz channels. The lines are the theoretical 1σ and 2σ error for this patch size. No strong outliers are seen. 

In the text 
Fig. C.1
Analytical estimates for the power spectra of the largest lowL meanfields at 143 GHz. The various components are discussed in more detail in Sect. C.1 (mask), Sect. C.2 (noise), and Sect. C.3 (beams). The meanfields all couple most strongly to even modes of the lens reconstruction, due to the approximate north/south symmetry of the scan strategy and Galactic mask. 

In the text 
Fig. C.2
Power spectra of the expected bias to our lensing estimators due to pixelization (thick solid lines), as well as crosscorrelation of the pixelization field with the estimators themselves (binned points). The theoretical lensing potential power spectrum is plotted for gradient modes as solid black, while the lower dashed black curves give respectively. The pixelization bias is always small compared to both signal and reconstruction noise. 

In the text 
Fig. D.1
Measured correlation coefficients for the bins of Table 1, estimated from 1000 simulated lens reconstructions. There is some evidence for correlations along the diagonal at the 10% level, which is sufficiently small that we have chosen to neglect binbin correlations in our likelihood. 

In the text 
Fig. D.2
Histograms of amplitude coefficients for the bandpowers of Table 1, taken from 1000 simulated lens reconstructions. Each coefficient has been rescaled by its standard deviation so that it can be directly compared to a unitvariance Gaussian distribution (solid lines). 

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.