# Cross-correlation of gravitational lensing from DES Science Verification data with SPT and Planck lensing

###### Abstract

We measure the cross-correlation between weak lensing of galaxy images and of the cosmic microwave background (CMB) on the same patch of sky. The effects of gravitational lensing on different sources will be correlated if the lensing is caused by the same mass fluctuations. We use galaxy shape measurements from 139 deg of the Dark Energy Survey (DES) Science Verification data and overlapping CMB lensing from the South Pole Telescope (SPT) and Planck. The DES source galaxies have a median redshift of , while the CMB lensing kernel is broad and peaks at . The resulting cross-correlation is maximally sensitive to mass fluctuations at . Assuming the Planck 2015 best-fit cosmology, the amplitude of the DESSPT cross-power is found to be and that from DESPlanck to be , where corresponds to the theoretical prediction and the errors are 68% confidence limits. These are consistent with the expected signal and correspond to significances of and respectively. We demonstrate that our results are robust to a number of important systematic effects including the shear measurement method, estimator choice, photometric redshift uncertainty and CMB lensing systematics. Significant intrinsic alignment of galaxy shapes would increase the cross-correlation signal inferred from the data; we calculate a value of for DESSPT when we correct the observations with a simple IA model. With three measurements of this cross-correlation now existing in the literature, there is not yet reliable evidence for any deviation from the expected LCDM level of cross-correlation, given the size of the statistical uncertainties and the significant impact of systematic errors, particularly IAs. Future data sets, including those from upcoming releases of DES and SPT, will cover more sky area and provide both greater depth and higher resolution, making this correlation a potentially very powerful cosmological tool. We provide forecasts for the expected signal-to-noise of the combination of the five-year DES survey and SPT-3G.

###### keywords:

gravitational lensing: weak; methods: data analysis; cosmic microwave background^{†}

^{†}pagerange: Cross-correlation of gravitational lensing from DES Science Verification data with SPT and Planck lensing–Cross-correlation of gravitational lensing from DES Science Verification data with SPT and Planck lensing

^{†}

^{†}pubyear: 2015

## 1 Introduction

Weak lensing (WL) refers to the coherent bending, by gravity, of light from distant sources as it passes through the large-scale massive structures of the Universe. WL is a powerful tool for studying the distribution and evolution of large-scale structure in the Universe because it is directly sensitive to both dark matter and baryonic matter. Here we focus on WL of two background sources: lensing of galaxies (which we will refer to as GWL) and lensing of the cosmic microwave background (referred to as CMBWL).

Observations of GWL rely on accurately measuring the shapes of a large number of small, faint galaxies, which are often at the edge of our survey magnitude detection limits. GWL surveys have steadily improved since the first measurements at the turn of the millennium (Bacon et al., 2000; Kaiser et al., 2000; Wittman et al., 2000; van Waerbeke et al., 2000), most notably with the Canada-France-Hawaii Telescope Lensing Survey (CFHTLenS, Heymans et al., 2012), and are now a viable probe of cosmology. Detection of CMBWL exploits mode couplings in the temperature and polarisation fields of the CMB, which are negligible in absence of WL (Blanchard & Schneider, 1987; Okamoto & Hu, 2003). This technique has been used in numerous studies to date (Hirata & Seljak, 2003; Smith et al., 2007; van Engelen et al., 2012; Planck Collaboration, 2014; Das et al., 2014; Planck Collaboration, 2015a; Story et al., 2015). Lewis & Challinor (2006) provide a comprehensive review of CMBWL.

GWL and CMBWL signals from the same patch of sky are expected to be correlated as both are, in part, sourced by the gravitational potentials of the same large-scale mass fluctuations. The cross-correlation of the lensing measurements from two such different sources offers a number of important applications. First, it provides a powerful check of systematics for cosmic shear measurements. For example, Vallinotto (2012) suggested that this cross-correlation can be used to mitigate shear measurement bias, whether from noise or complex galaxy morphologies (Jarvis et al., 2015), to which CMBWL is insensitive. The same is true of other observational and astrophysical systematics such as modelling the point spread function (PSF) and galaxy intrinsic alignments (IAs), though the estimation of lensing potential from CMB maps is prone to some multiplicative errors of its own. Taking GWL and CMBWL together, there is considerable scope for calibration of these bias terms through cross-correlation of the GWL and CMBWL signals. In addition, CMB lensing offers an extra high-redshift source bin that can be included in joint analyses of late-universe probes (Vallinotto, 2013) to study late-time dark energy or modifications to gravity.

Two measurements of the GWLCMBWL cross-correlation have previously been reported: Atacama Cosmology Telescope (ACT) lensing data crossed with 121 deg of galaxy lensing convergence as measured by the Canada-France-Hawaii Telescope (CFHT) Stripe 82 Survey (Hand et al., 2013), and Planck lensing crossed with 140 deg galaxy lensing data as measured by CFHTLenS (Liu & Hill, 2015). Both report low detected signals compared with expectations for the Planck best-fit cosmology, with Liu & Hill (2015) reporting a particularly low signal at roughly half the expected amplitude, a discrepancy. In this paper we aim to obtain a measurement of the GWLCMBWL cross-correlation using new data from DES Science Verification (SV) and CMBWL maps from the South Pole Telescope (SPT) and Planck. Our measurement has slightly deeper CMBWL data (SPT compared with ACT and Planck), similar sky coverage, but slightly shallower GWL measurements compared to those used in previous GWLCMBWL results. Importantly, with DES and SPT, we employ different GWL and CMBWL data to those used before for this cross-correlation and cover a different patch of sky to the previous analyses. Our results therefore serve as an independent check on the measurements made by Hand et al. (2013) and Liu & Hill (2015).

We begin by describing the relevant theory and formalism for GWL and CMBWL in Section 2. In Section 3 we describe our data from DES, SPT and Planck. We describe our methods in Section 4 and present our measurements of the GWLCMBWL cross-correlation in Section 5. In Section 6 we demonstrate that our results are robust to a variety of important systematic effects and consistency checks. In Section 7 we compare the power of the current data with that expected from the full DES survey and SPT-3G. We discuss the implications of our measurements, their relation to previous results and the future potential of this cross-correlation in Section 8. Throughout this paper we employ Planck 2015 cosmology (TT+TE+EE+lowP+lensing+ext) with , , , , .

## 2 Theory

In this paper we consider two light sources that experience weak lensing: galaxies and the CMB. Two particularly useful quantities associated with the distortion of light are the spin-0 convergence field, , and spin-2 shear field, (see for example Bartelmann &
Schneider, 2001; Munshi et al., 2008; Hoekstra &
Jain, 2008, for details and definitions of the lensing quantities). Both are derivatives of the lensing potential, which describes the strength of lensing for a given configuration of source, lens and observer.
In GWL, the main observable is shear, which is measured by the distortions of the source galaxy shapes.^{1}^{1}1Or, more correctly, the reduced shear, . For WL, and . Convergence, a measure of the magnification of the image, can be reconstructed from the shear. In CMBWL both shear and convergence can be reconstructed from the temperature map. The analytic expressions given in this section are equally applicable to shear or convergence power spectra.

Since their means vanish, it is convenient to quantify the fluctuations in both GWL and CMBWL with angular two-point functions, in particular auto- and cross-power spectra in harmonic space. Under the Limber approximation (Kaiser, 1992), these take the form of integrals over the non-linear matter power spectrum, , and a pair of appropriately chosen window functions. We are interested in the cross-correlation between GWL and CMBWL,

(1) |

where is the comoving distance to redshift , and is the distance to the horizon. Here and are the GWL and CMBWL window functions.

The GWL window function, also known as the lensing efficiency function or lensing kernel, takes the form

(2) |

where is the Hubble parameter, the speed of light, the total matter density and is the galaxy redshift distribution. We have assumed a flat universe, as we will throughout the paper.

The CMBWL window function takes a similar form but is somewhat simpler due to the single source plane,

(3) |

where is the comoving distance to the last scattering surface. Although the CMBWL weight function peaks at (Lewis & Challinor, 2006), it is sensitive to the integrated gravitational potential between the source and the observer. The window functions for CMBWL and GWL corresponding to the DES sources used in this work are shown in Figure 1. The DES source galaxies are sensitive to mass fluctuations at lower redshift than the peak of the CMB lensing kernel () but there is sufficient overlap to expect a significant cross-correlation in the gravitational lensing signal of both sources.

During the epoch of structure formation, galaxies experience tidal forces due to the gravitational potential of the surrounding mass distribution. The presence of such forces may cause the ellipticity and orientation of neighbouring galaxies to become aligned (Heavens et al., 2000; Catelan et al., 2001; Hirata & Seljak, 2004). This effect, known as intrinsic alignment (IA), produces a correlation between the intrinsic shapes of galaxies and will be present as an extra term in our GWLCMBWL cross-correlation measurement (Troxel & Ishak, 2014; Hall & Taylor, 2014; Chisari et al., 2015),

(4) |

The IACMBWL cross-correlation can be calculated, like the other power spectra, using Equation 1 with the appropriate IA weight function instead of the GWL weight function. In this paper we assume IAs are described by the widely-used non-linear alignment (NLA) model (Bridle & King, 2007; Hirata & Seljak, 2010), which means the weight function is given by

(5) |

Here is the critical density at and , a normalisation constant based on the SuperCOSMOS measurement at low redshift (Brown et al., 2002). is the linear growth function, normalised to unity at . Current measurements of IAs for different galaxy types over different redshifts still leave a significant uncertainty as to the expected level of the IA contribution for any given sample of source galaxies. The DES SV cosmology analysis (Dark Energy Survey Collaboration, 2015) was consistent with an IA signal between zero and four times within 2 confidence limits. We treat the range of possible IA contributions in Section 6.5 and discuss the significance of IAs for our measurement and future GWLCMBWL analyses in Section 8. First though, we neglect the IA contribution in our main analysis, fitting for the pure GWLCMBWL signal only. Note that the negative sign in Equation 5 means that the IA contribution subtracts from the total observed GWLCMBWL cross-correlation. This is because the galaxy shape alignment sourced by IAs is of the opposite sense to that sourced by GWL. This means any measured cross-correlation amplitude will be lower than it would be if IAs were taken into account and that the IA contribution, if significant, would therefore increase the significance of the detection if it were included in the analysis.

## 3 Data

In this section we describe the data sets used in this paper. The DES SV data products are introduced in Section 3.1, including the photometric redshift (photo-) and shear catalogues. Next, we describe the CMB data sets and their lensing maps, including those from SPT (Section 3.2.1) and Planck (Section 3.2.2). Finally, in Section 3.3, we describe the GWL and CMBWL simulations and mock catalogues used for pipeline testing and covariance estimation.

### 3.1 The Dark Energy Survey

The Dark Energy Survey (DES) is an optical survey, currently in progress, which will cover 5,000 deg in five filters () using DECam (Flaugher et al., 2015) over five years, reaching a limiting magnitude of in the i-band (Sánchez et al., 2014). The DES footprint was designed to overlap significantly with the region observed by the SPT (described in the next section), enabling many interesting cross-correlation measurements (Carlstrom et al., 2011; Giannantonio et al., 2015; Saro et al., 2015). In this paper we use the DES Science Verification (SV) data, which was taken during the period November 2012 – February 2013, before the start of the main survey in late August 2013. Specifically, we use the 139 deg contiguous area of the DES SV data that overlaps with the SPT East field; this is known as the SPT-E region (Jarvis et al., 2015) and is centred on RA 77.5 deg, DEC -51 deg. All the DES data products in this paper have been reduced from the raw survey data by the DES Data Management pipeline (Mohr et al., 2012; Desai et al., 2012).

#### 3.1.1 The DES SV photo- catalog

Large optical surveys like DES use photometry to estimate the redshift of source galaxies; this estimation technique is known as photo- (Hildebrandt et al., 2010). Fluxes in multiple broad filter bands are measured as a kind of very low resolution spectrum and a variety of methods are employed to estimate the corresponding true redshift, exploiting broad features of the spectral distribution rather than spectral lines. This is obviously less accurate than spectroscopic approaches, which produce high resolution spectra for each object, but the photo- approach is faster and cheaper for large imaging surveys. For our cross-correlation measurement it is essential that the overall redshift distribution, , is well characterised as this affects the theoretical cross-correlation power spectrum that we fit to our measured values.

The primary photo- catalogue used in this work is produced by the SkyNet2 neural network algorithm described in Graff & Feroz (2013); Bonnett (2015); Bonnett et al. (2015). Four different photo- methods were extensively tested and characterised in Bonnett et al. (2015) for the galaxy sample used in the DES SV shear catalogs. In addition to SkyNet2, these were the ANNz2 (Sadeh et al., 2015), BPZ (Benítez, 2000; Coe et al., 2006), and TPZ (Carrasco Kind & Brunner, 2013, 2014) photo- estimation pipelines, selected because they performed well in the analysis of Sánchez et al. (2014). Bonnett et al. (2015) found that the catalogues agreed to better than in the mean photo- of the distribution and we explore the effect of using these alternate photo- estimators in Section 6.6.

SkyNet2 produces a full probability density function, , for each galaxy. We use the mean of the for each galaxy, , as its point-estimate redshift and select/bin the galaxies according to . We use this point estimate to select galaxies in the redshift range , as this is the range over which the SkyNet2 algorithm gives reliable photo-z estimates according to Bonnett et al. (2015). Other point estimates were tested but the mean produced the best results for our application as it reduces the impact of unphysical features in individual galaxy . However, for the calculation of the theoretical prediction of our cross-correlation measurements, we use the full stacked from each galaxy to estimate the total redshift distribution for our forecast. Figure 1 shows the resulting redshift distribution of the DES sources used in this cross-correlation as well as the corresponding GWL lensing kernel, Equation 2, and the CMBWL kernel, Equation 3, for comparison.

#### 3.1.2 The DES SV shear catalog

The DES shear catalogs, described in detail by Jarvis et al. (2015), are based on two independent galaxy shape measurement algorithms: ngmix (Sheldon, 2014) and im3shape (Zuntz et al., 2013). Note that although the coadd catalogues were used in the initial processing of the shear measurement pipeline, the galaxy shape measurements were carried out at the single-exposure level. The shape (or shear) of each galaxy is estimated by jointly fitting a galaxy model to multiple single-exposure images of that same galaxy (Jarvis et al., 2015), thus reducing the impact of instrumental artefacts and PSF variation between exposures.

The analysis of Jarvis et al. (2015) showed that both catalogues passed all requirements on contamination by systematic effects and that the two catalogues were consistent with each other under a range of statistics. Becker et al. (2015) and Dark Energy Survey Collaboration (2015) demonstrated that the two catalogues remain consistent at the level of two-point statistics and inferred cosmological constraints. We use the ngmix catalogue for our main analysis because it has a larger effective source galaxy number density, 5.7 arcmin as opposed to 3.7 arcmin for im3shape (see Jarvis et al. (2015) for details of the weighting used to calculate this effective number density). This choice is consistent with the DES SV GWL two-point analysis (Becker et al., 2015) and the DES SV GWL cosmology analysis (Dark Energy Survey Collaboration, 2015). We repeat our analysis using the im3shape catalogue as a consistency check; these results can be found in Section 6.1. Both catalogues have demonstrated that they can provide shear measurements with systematic uncertainties (whether from astrophysical, observational, or measurement effects) subdominant to the statistical uncertainty in the SV data for different cosmological probes including two-point statistics in real and harmonic space, galaxy-galaxy lensing and mass-mapping (Jarvis et al., 2015; Chang et al., 2015; Vikram et al., 2015; Becker et al., 2015; Clampitt et al., 2015).

We produce maps of our shear catalogues, described in more detail in Section 4, using the HEALPix pixelisation scheme at (Górski et al., 2005). This corresponds to a pixel area of arcmin or a pixel scale of kpc at , our redshift of maximal sensitivity.

Shape estimates from both pipelines were ‘blinded’ during our analysis to avoid experimenter bias (Klein & Roodman, 2005). This meant that a constant scaling factor (between 0.9 and 1) was applied to all ellipticities. This would slightly alter the amplitude of the cross-correlation, preventing over-fitting to results from other papers or to any given cosmology. Our analysis procedure was finalised and fixed before de-blinding.

### 3.2 CMBWL Maps

#### 3.2.1 SPT Lensing Maps

The CMB maps are based on temperature measurements made with the SPT (Carlstrom et al., 2011), which is a 10 metre diameter telescope located at the National Science Foundation Amundsen-Scott South Pole Station in Antarctica. During 2008-2011, this telescope was used to conduct a tri-band (90, 150, 220 GHz) wide-field survey covering 2540 (Story et al., 2013). The survey area is composed of 19 subfields, all of which were scanned in a similar fashion, reaching minimum depths (maximum noise levels) of 40 -arcmin (90 GHz), 18 -arcmin (150 GHz) and 70 -arcmin (220 GHz) with roughly arcminute resolution.

The maps are produced by applying quadratic estimators (Okamoto & Hu, 2003) on the filtered 150 GHz temperature map in the range on a region extracted from the full survey area, centred on the DES SPT-E field. Modes with are also removed to minimise extra noise along the scanning direction (van Engelen et al., 2012). Positive sources brighter than (corresponding to approximately 10 mJy) and clusters detected above 6 were masked with a aperture, and the masked regions were filled using Wiener filter interpolation. The maps are produced on a HEALPix grid of .

The signal-to-noise in our cross-correlation drops significantly at high , and the contamination to the reconstruction from emissive sources and galaxy clusters that are not masked becomes increasingly important at high . For these reasons, we impose a conservative cut on the map of This cut does not lead to a significant loss of signal. We have tested this choice of and found our result robust to values between and . All of these issues are even more relevant to the Planck CMBWL data (which has higher noise per mode than SPT, particularly at high , as well as potentially higher contamination from point sources and clusters because of the larger Planck beams), and we apply the cut to the Planck data (described in the next section) as well.

#### 3.2.2 Planck

We use the Planck lensing maps from the second data release, which were made public in 2015^{2}^{2}2http://pla.esac.esa.int/pla/. These maps are produced by using filtered temperature and polarisation measurements from the Planck satellite (Planck
Collaboration, 2015a). The temperature and polarisation maps are both constructed by taking linear combinations of multi-frequency data (30 - 857 GHz for temperature and 30 - 353 GHz for polarisation) with scale-dependent coefficients using the SMICA method to produce foreground cleaned minimum variance maps (Planck
Collaboration, 2015b). Planck HFI beams range from 4 (857 GHz) to 10 (100 GHz) arcminutes in resolution, compared to SPT’s 1 arcminute. Planck covers the full sky, while SPT is focused on a smaller patch.

Similar to the SPT map, the CMBWL potential is produced using the quadratic estimators from Okamoto & Hu (2003). The main difference here is the availability of E and B-mode polarisation, which allows for additional estimators (, , , ) in addition to . These estimators are combined to form a minimum-variance estimate of the lensing potential , which is provided in the form of spherical harmonic coefficients of the CMB lensing convergence , filtered to , along with the analysis mask. The map is in HEALPix format with resolution .

### 3.3 Simulations

We test our estimators using simulated data sets, constructed specifically to mimic the noise and other statistical properties of each of our observables: GWL shape catalogues for DES and CMBWL convergence maps for SPT and Planck.

For DES, we use two sets of mock catalogues in addition to the data itself. The first are based on N-body simulations and are the same set of simulated galaxy catalogues described in Becker et al. (2015), consisting of 126 realisations of the SV SPT-E patch with lensing fields calculated by ray-tracing. We sub-sample the galaxies in the simulations to match the galaxy number density, redshift distribution and noise properties of the data, drawing a corresponding shape value from the data and adding it to the cosmological shear signal from the simulation to ensure a realistic shape noise distribution.

When estimating the noise properties of our cross-correlations we also employ a separate set of 100 DES mock catalogues produced, not from N-body simulations, but simply by applying a random rotation to the orientation of each source galaxy in our DES shape catalogue. This retains the spatial and redshift distribution of galaxies as well as the overall intrinsic ellipticity distribution across the sample, while destroying any cosmological information. Therefore these randomised catalogues act as noise-only realisations.

The SPT collaboration has produced 100 lensed and 100 unlensed simulated sky realisations. Instrumental noise is added to these sky realisations and the result is converted into simulated time streams, then processed in the same way as the real data, including the effects of masking and filtering. The 100 output reconstructions were then cross-correlated with the input lensing potential, , to obtain the lens transfer function. This transfer function was then applied to reconstructions from the unlensed CMB realisations, providing 100 realistic realisations of noise in the map. Noise in the CMB convergence maps comes primarily from fluctuations in the CMB temperature field itself, rather than instrumental noise. More discussion of these noise simulations can be found in Giannantonio et al. (2015).

Planck has produced 100 simulated lensing maps using their Full Focal Plane 8 (FFP8) Monte Carlo simulations (Planck Collaboration, 2015c). The noise level in these simulations has been tuned by the Planck collaboration to match the amplitude of the Planck CMB power spectra.

## 4 Methods

In this section we describe our main analysis pipeline - how we make our cross-correlation measurement and calculate the associated error. This main analysis is conducted using the PolSpice code (Szapudi et al., 2000; Chon et al., 2004) to estimate projected angular power spectra in harmonic space using the ngmix shape catalogues and the CMB convergence maps from SPT and Planck. The results themselves are presented in Section 5.

We produce shear maps by averaging shear estimates for individual galaxies from the DES shape catalogues into HEALPix pixels at a resolution of , applying the standard cuts, weighting and bias corrections as described in Jarvis et al. (2015). The two shear components, and , are treated separately to produce two maps. We then combine and maps from DES with the maps from SPT/Planck into a triplet. We use the PolSpice (Szapudi et al., 2000; Chon et al., 2004) code in polarisation mode to estimate projected angular power spectra for the auto- and cross-correlations of and DES . Our input are analogous to the CMB triplet. We have tested this process on the DES, SPT and Planck simulations described in Section 3.3 to confirm it is capable of returning an unbiased estimate in the presence of the relatively restrictive DES SV survey mask.

The shear maps can be decomposed into spin spherical harmonics

(6) |

where and are the spin spherical harmonics and their coefficients. Of particular interest to us are the linear combinations of these spherical harmonics which produce curl-free E-mode and divergence-free B-mode components of the shear field (Bartelmann, 2010; Schneider et al., 2002; Crittenden et al., 2002):

(7) | ||||

(8) |

To first order, WL only generates E-modes because the lensing potential is a real scalar (Schneider et al., 2002; Crittenden et al., 2002). We therefore aim to measure the cross-correlation. Any measurable correlation is unphysical and evidence of some untreated systematic effect in the DES data which is correlated with CMBWL (see Section 6.4 below for analysis). The decomposition is handled by PolSpice, which is designed to make the same split into curl- and divergence-free components in the context of CMB polarisation studies (Munshi et al., 2002). We have used our simulation catalogues to confirm that PolSpice can recover an unbiased estimate of the cosmic shear power spectrum via this method.

We apply simple binary masks to our data. For DES, HEALPix pixels are set to zero if they do not contain any DES SV source galaxies after all our quality cuts have been applied; the remaining region covers 139 deg. For SPT we have data over a wider contiguous patch of 600 deg, overlapping the DES SV SPT-E region. The mask is one inside this region and zero outside. For Planck we use the publicly available mask designed for power spectrum estimation, where point sources and galactic emission are masked out, removing roughly one third of the full sky.

We take a hybrid approach to estimating the noise in our cross-correlation, describing the noise with an analytic expression that includes estimates of noise in the constituent GWL and CMBWL parts acquired from simulations or randomised realisations of the data. Alternate noise estimators and the considerations governing our choice are described in Section 6.3. The analytic form of the noise is given by

(9) |

where , and are the theory power spectra of the WL auto-correlation, the CMB lensing auto-correlation and the WL/CMB cross-correlation respectively, assuming the Planck 2015 best-fit cosmology, and is the fraction of sky covered by the survey. These are calculated using CAMB sources (Challinor & Lewis, 2011).

The auto-correlation contributions to the noise estimate include measurement noise terms, and . The noise term for DES is the mean of the auto-correlations of the maps produced from the 100 randomised realisations of our shape catalogue. The CMBWL noise term is calculated as the mean of the auto-correlations of the 100 noise realisations provided by both the SPT and Planck collaborations. Each of these sets of noise simulations are described in Section 3.3 above. Note that we have estimated the uncertainty using the mean noise level from Planck, whereas in reality the SV patch covers a region where Planck noise is lower than average, hence we are slightly over-estimating the noise in the DESPlanck case.

We calculate using our most restrictive mask, in this case that of the DES shape catalogues. As mentioned above, the DES SV SPT-E patch which we are using covers 139 deg of sky, though the exact sky fraction varies slightly depending on the shape catalogue and redshift range considered. We take this variation into account in our calculations. In the case of DESPlanck there are a small number of additional regions inside the DES SV mask but excluded from the Planck mask. We take the product of the two masks when estimating , reducing the sky fraction by 1.5% compared to DESSPT.

We also tested a number of alternative techniques to estimate the cross-correlation signal, covariance and noise as consistency checks on our main analysis. These are detailed in Section 6 below.

## 5 Results

Figure 2 shows our results using PolSpice to correlate with in harmonic space. The measurement is averaged into 16 linearly-spaced bins over the multipole range . We then use our forecast cross-correlation power spectrum to fit a single free parameter, the cross-correlation amplitude , by a simple minimisation:

(10) |

where the error, , is calculated according to Equation 9. The fits to the cross-correlation amplitude are detailed in Table 1. If our measurement were consistent with the expectation from theory, assuming the Planck cosmology and that we have correctly modelled the DES galaxy redshift distribution, then we would expect a result consistent with .

Redshift Range | ||
---|---|---|

ngmix SPT | 0.93 | |

ngmix Planck | 1.52 |

Our measurement shows E-mode cross-correlations with best-fit amplitudes of for DESSPT and for DESPlanck, giving a significance of and respectively. We estimate the goodness-of-fit by calculating per degree of freedom (15, the number of bins minus one), finding good fits in both cases, with for DESSPT and for DESPlanck. The measurements with SPT and Planck are consistent with each other and with the theoretical expectation. The DESPlanck cross-correlation has a relatively high , with the probability to exceed such being ; a similar result was found in Giannantonio et al. (2015) for DES SV LSSPlanck CMBWL. These measurements of fix all other cosmological parameters at the Planck 2015 best-fit cosmology and ignore IAs.

We consider in detail the impact of a variety of systematic effects in Section 6. Some effects, including uncertainty in shear measurement bias and estimation of photo-z could change the measured amplitude of the cross-correlation. The result quoted here should be considered the ‘bare’ constraint on , when we assume our best estimates for shear measurement bias, photo-z. This is a reasonable approach because any deviation in these quantities will scale both the best-fit amplitude and the error bars, leaving the significance of detection unaffected.

## 6 Consistency and Systematics Tests

In this section we summarise a number of checks carried out to ensure that our analysis is accurate and robust to observational and astrophysical systematic effects. A substantial amount of work has been done quantifying the contribution of systematics to our data sets (Jarvis et al., 2015; Bonnett et al., 2015; Becker et al., 2015; Vikram et al., 2015; Leistedt et al., 2015). In particular, Giannantonio et al. (2015) dealt with a number of systematics that could potentially manifest as spurious signal in the cross-correlation of DES SV galaxy number density with CMBWL; all were found to be of negligible importance. In this paper we will concentrate on those of particular relevance to the GWLCMBWL cross-correlation.

### 6.1 Shape Measurement Pipelines

As mentioned in Section 3.1.2, the DES collaboration has produced two independent shape measurement catalogues: ngmix (Sheldon, 2014), which we use for our main analysis, and im3shape (Zuntz et al., 2013).

We have repeated our cross-correlation measurement using the im3shape shape catalogue. The results are in good agreement with those from ngmix and our forecasts but the errors are larger due to the lower effective source number density in im3shape (3.7/arcmin compared to 5.7/arcmin for ngmix). With im3shape, we measure a cross-correlation amplitude of for DESSPT and for DESPlanck. Like our main results, the im3shape measurements are consistent with the expected signal within errors.

### 6.2 Alternate Estimators

Besides the aforementioned PolSpice pipeline, we tested an additional flat-sky implementation of the same calculation, also known as the Kaiser-Squires method (KS, Kaiser & Squires, 1993). The KS method was used in both Hand et al. (2013) and Liu & Hill (2015) and shown to perform well. With an eye on the larger sky coverages in future data sets, we have performed our main calculations based on a curved sky analysis but we also checked whether our results are consistent with the those from the flat-sky KS analysis used in previous literature.

In the flat-sky KS approach and for each galaxy in the shear catalogue are projected onto a zenith equal area projection coordinate grid, then averaged over all the galaxies that fall on the same square grid. Since the shear measurements are made with respect to spherical coordinates, we apply the rotation

(11) | ||||

(12) |

where is the local angle between equal right ascension to the -axis of the image, such that the shear measurements are described with respect to the flat-sky -coordinates.

The maps are converted into by Fourier transforming using

(13) |

and are cross-correlated with CMB in Fourier space,

(14) |

where . In this process, we also apply a mask due to the finite survey area, which inevitably induces mode-coupling between the bins. To account for the mask, we follow Hivon et al. (2002) to decouple the effect of the mask from the computed spectra. Using this method, we obtain amplitudes of for DESSPT and for DESPlanck, with a of 1.18 and 1.17 respectively, which is consistent with the PolSpice pipeline.

We also tested an alternate method which uses the pseudo-Cl estimation technique for CMB polarization developed in Kogut et al. (2003). Spectra are computed from pseudo-multipole coefficients and related to the true spectra through a coupling matrix, which we bin before inversion (see also Efstathiou, 2006). Using this approach, we obtain amplitudes of for DESSPT and for DESPlanck, which are in good agreement with the other two methods.

### 6.3 Alternate covariance estimates

We have checked the hybrid noise estimates used in our main analysis by also estimating the full covariance of our cross-correlation using both our N-body simulations and jack-knife resampling. Both estimate the covariance according to

(15) |

where for the N-body case and for the jack-knife and counts over the separate cross-correlation realisations, . For the N-body method , the number of DES/SPT simulations used, and for the jack-knife , the number of equal area regions (3.5 deg) into which the data is split. is the mean of the cross-correlations in each case. For more details of both approaches to covariance estimation see the excellent review by Norberg et al. (2009).

All three approaches give consistent estimates of the error, with both the jack-knife and N-body errors agreeing with our hybrid estimate to within 10% at all . There is little off-diagonal power in the resulting covariance matrices using either method. When inverting the covariance matrices we apply the correction of Hartlap et al. (2007) to account for the effect of having a finite number of realisations. We are confident that our hybrid approach remains the most accurate available noise estimate for this cross-correlation measurement. The simulations used are not correlated between the GWL and CMBWL, therefore the N-body method, which would otherwise be preferred, misses power due to both correlated non-linear structure growth and correlated cosmological signal. In addition the assumed cosmology in the simulations produces a lower amplitude signal than the Planck best-fit assumed throughout the rest of this work. The jack-knife approach is useful as a consistency check but is known to over-estimate errors on small scales (Norberg et al., 2009). In future work, with the full DES survey data set, the production of correlated N-body simulations of GWL and CMBWL should be a main priority so that we can use the full covariance estimate and capture off-diagonal power from mode-mixing due to correlated structure formation and mask effects.

### 6.4 B-modes

In Section 4 we describe how the GWL shear field, , can be decomposed into curl-free (E-mode) and divergence-free (B-mode) contributions. WL will only produce E-mode signal to first order so, in cosmic shear measurements, detection of a significant B-mode auto-correlation is a diagnostic for systematics in the observations or the measurement, for example a poorly reconstructed PSF. While we do not necessarily expect these processes to generate B-modes that are positively correlated with CMBWL, it is possible that our estimator could introduce spurious power to our measurement (which would appear equally in the E/B-mode) or allow power to leak from the to the signal. Testing for significant cross-correlation is therefore still a useful check on the efficacy of our estimators. Note that although the CMB lensing signature may contain some B-mode fluctuations, the magnitude is substantially less than the gravitationally-induced lensing signal (Planck Collaboration, 2014) and we ignore this effect in our current study.

In Section 5 we presented the cross-correlation of our GWL E-mode signal with lensing from the CMB, . In Figure 3 and Table 2 we show the equivalent cross-correlation. We estimate the significance of the signal by fitting the expected E-mode signal with, as before, a varying constant cross-correlation coefficient, . We estimate the best-fit value of that constant by minimising in the same way as with the theoretical forecasts in our main measurement. A significant signal of this kind would indicate leakage of power from E- to B-modes. Our best-fit B-mode cross-correlation amplitudes are consistent with zero for both DESSPT, , , and DESPlanck, , . We also checked the for zero cross-correlation, , finding this to be a good fit to the data in both the DESSPT, , and DESPlanck, , cases. This gives us confidence that our estimator is not suffering from either spurious power or E to B leakage at a level that could bias our measurement significantly, given the size of our errors.

Redshift Range | |||
---|---|---|---|

ngmix SPT | 0.79 | 0.83 | |

ngmix Planck | 0.92 | 0.95 |

### 6.5 Intrinsic Alignments

As described in Section 2, GWL measurements are contaminated by the alignment of unlensed galaxy shapes determined by large-scale gravitational potentials during galaxy formation, known as intrinsic alignments (IA) (Heavens et al., 2000; Croft & Metzler, 2000; Catelan et al., 2001; Crittenden et al., 2001; Hirata & Seljak, 2004). The IAs reduce the overall observed correlation because they tend to align galaxy shapes and the matter distribution with the opposite sign to the GWL alignment.

We examine the potential contamination of our GWLCMBWL cross-correlation by IAs using the NLA model as given in Equation 5. Figure 4 shows the effect of IAs on the cross power spectrum. For our fiducial SkyNet2 source redshift distribution the decrement in expected observed power spectrum, , can be as much as 18% around but is substantially lower at other scales.

This is a relatively simplistic approach to the modelling of IAs. We have assumed that the NLA model applies to all types of galaxy even though it is only designed to explain alignment of elliptical galaxies and there is, as yet, no positive detection of IAs in spiral galaxies for large surveys in the field (Mandelbaum et al., 2011). Furthermore we have assumed that the amplitude of alignment is set exactly by the normalisation of Brown et al. (2002); Bridle & King (2007), though the DES SV cosmic shear analysis was equally consistent, at the 2 level, with there being no IAs or IAs at four times this assumed amplitude (Dark Energy Survey Collaboration, 2015).

Our aim is to determine the significance of IAs in the detection of the cross-correlation and check that our measurement is robust to IAs. We defer attempts to make precise measurements of the IA signal to later work when the volume and quality of DES and SPT data will have greatly increased, producing a significant improvement in precision; see Figure 6 and discussion in Section 7 for more details.

Assuming all galaxies are affected by the NLA model, with normalisation given by , is a conservative way to model a ‘significant’ IA effect that can be compared to our main analysis where IAs are ignored entirely. Chisari et al. (2015) have recently made more careful models of IA contamination in the context of the Hand et al. (2013) data sets. Their estimated potential levels of contamination for both red and blue galaxies are slightly lower than our assumed model ( vs. ), confirming that our implementation is a realistic, conservative example.

Including IAs in this way shifts our best-fit cross-correlation amplitude from to . The effect of IAs reduces the expected cross-correlation signal and accounting for this effect increases our measured cross-correlation amplitude. The significance of the detection remains unchanged at and the result is entirely consistent with forecast expectations and the measurement without IAs within the errors. Nevertheless the shift due to intrinsic alignment is at level, indicating that future, higher precision measurements of this cross-correlation have the potential to be a powerful probe of IAs.

### 6.6 Photometric Redshift Uncertainties

Our main analyses are based on the SkyNet2 photo- catalog but we also cross-check our main results with three other photo- estimation pipelines validated by the DES collaboration: BPZ2, ANNz2 and TPZ. See Section 3.1.1 and Bonnett (2015) for more discussion of these estimators.

We recalculate our measurement of the cross-power spectrum in exactly the same way as our main analysis, but using the different photo- estimates. The photo- estimates are used to select the galaxies which enter our maps and to create the , which is in turn used to produce the theory used to fit the cross-correlation amplitude, . See Bonnett et al. (2015) for more detail on the testing of these photo- pipelines in the context of the DES SV GWL analysis.

Figure 5 shows the difference in estimated for each choice of estimator, relative to the fiducial SkyNet2. It is clear that the variation in measured cross-correlation due to different photo- estimation codes is significantly smaller than the error on the same measurement. The scatter in best-fit is well within the confidence limits.

Bonnett et al. (2015) used their analysis of multiple photo-z pipelines to define a Gaussian prior on the mean of the photo-z distribution (of tomographic bins) of width . Any variation in the mean of the photo-z distribution would scale both our best-fit value of and our errors, leaving the significance of our detection unchanged. We have checked that the shift in best-fit cross-correlation amplitude due to a change of in the mean of the photo-z distribution is well within our one sigma errors and we quote the ‘bare’ constraint on as our headline result, where we have assumed the accuracy of the SkyNet2 derived by stacking the of individual galaxies.

### 6.7 Systematic uncertainties in the map

We have also tested for the degree of contamination in the SPT -map due to point-sources and the tSZ effect by applying a more stringent mask than the one used for lensing reconstruction. The more stringent mask removes point sources detected between (corresponding to approximately 6 mJy) and using a radius circular aperture in addition to the mask applied to sources detected above using the main mask. Clusters catalogued in Bleem et al. (2015) with detections between and are also masked with a radius disk in addition to the mask applied to clusters detected above . We obtain an amplitude of with when applying this mask, which is entirely consistent with our main result, suggesting that our kappa maps are minimally contaminated by these sources.

## 7 Forecasts

The volume of GWL surveys will greatly increase over the coming years. Dark Energy Survey^{3}^{3}3http://www.darkenergysurvey.org/ (DES) will deliver an unprecedented of lensing data by 2018, with projects including Hyper Suprime Cam (HSC)^{4}^{4}4http://www.naoj.org/Projects/HSC/, Kilo-Degree Survey (KiDS)^{5}^{5}5http://kids.strw.leidenuniv.nl/, Euclid^{6}^{6}6http://www.euclid-ec.org/ and Large Synoptic Survey Telescope (LSST)^{7}^{7}7 http://www.lsst.org/ also producing data over the next decade. These surveys will push deeper than previous efforts, increasing the overlap with the CMB lensing kernel, which is broad and peaks at (Lewis &
Challinor, 2006). On the CMB side, Planck^{8}^{8}8http://www.cosmos.esa.int/web/planck has set a new standard for all-sky CMB surveys from space but it will also be important to maximise the overlap of galaxy surveys with high-resolution CMB surveys using the upgraded cameras on the South Pole Telescope (SPT)^{9}^{9}9https://pole.uchicago.edu/ and the Atacama Cosmology Telescope (ACT)^{10}^{10}10http://www.princeton.edu/act/, as well as the next-generation PolarBear^{11}^{11}11http://bolo.berkeley.edu/polarbear/ instrument.

This paper represents an important test of the WL measurement pipelines in both the DES and SPT collaborations, and allows us to look forward with confidence to more scientifically ambitious analyses in the future when these more powerful data sets become available, particularly the full DES survey and the SPT third generation camera (SPT-3G) (Benson et al., 2014). The principal gain will be increased sky coverage, with of overlapping area expected from the full DES five year survey (Y5) and SPT-3G. This represents an -fold increase over the data used in this work. In addition, the SPT-3G upgrade will significantly decrease the noise level compared to current SPT measurements (SPT SZ). Estimates from the SPT collaboration foresee a factor of decrease in effective noise between SPT SZ and SPT-3G when temperature measurements alone are used to reconstruct the CMBWL convergence map, and a factor of between SPT SZ and SPT-3G when the SPT-3G reconstruction also uses CMB polarisation measurements.

Figure 6 shows the expected signal-to-noise (S/N) from DES Y5 and SPT-3G, compared to that from the DES SV and SPT SZ data used in this paper. We have restricted this forecast to the expected overlapping area available by DES Y5 and assumed moderate improvements in number density and GWL shape noise for DES (see figure caption for details). We can confidently expect a detection of GWLCMBWL from DES Y5 SPT-3G with a S/N of . This huge increase in measurement power over the coming years will allow us to move beyond detection of the cross-correlation and to exploit this measurement to answer a number of science questions. Note that there is a turnover in the cross-correlation power spectrum at low ell. We have excluded this turnover from these forecasts by retaining a minimum scale of . Increased coverage of this feature would further improve the power of this particular cross-correlation.

The very different observational properties of the two surveys means that the cross-correlation is an extremely useful discriminant of measurement systematics. Both CMBWL and GWL are affected by multiplicative biases in the measurement of the lensing signal. For example, uncertainties in measuring galaxy shapes leads to a shear measurement bias in DES GWL, currently marginalised over in the cosmology analysis (Dark Energy Survey Collaboration, 2015). As both probes are estimated from different types of data using very different techniques, there is considerable scope for calibration of these bias terms through cross-correlation of the GWL and CMBWL signals.

As the precision of our cross-correlation increases, the systematic effects will become more significant. In this work we estimated the order of magnitude effect of galaxy IAs, finding that the presence of IAs could shift our best-fit measurement of by a significant fraction of the errors. However, this did not alter our level of agreement with theory, given the size of our error bars. The much higher S/N measurement we can expect from future data means that the impact of IAs will be much more significant. On the one hand this means that we need to improve our modelling of IAs, paying particular attention to the impact of galaxy type and luminosity; at the same time, it means that the GWLCMBWL cross-correlation has the potential to make precision measurements of the IA signal in exactly the data sets we want to use for cosmic shear analyses, possibly allowing us to discriminate between competing IA models.

Using this cross-correlation to measure the amplitude of the IA signal for different types of source galaxy is a real possibility. This will be important for cosmic shear studies and the first robust test of the modelling done in Troxel & Ishak (2014); Hall & Taylor (2014); Chisari et al. (2015) and Larsen & Challinor (2015), who have suggested that the type of spiral galaxies that dominate WL data sets may experience an IACMBWL correlation of opposite sign to that expected in the linear alignment model.

## 8 Summary

We have found evidence for the cross-correlation of GWL measured in 139 deg of DES SV galaxy shape catalogues (effective source number density 5.7 arcmin) and CMBWL from both SPT and Planck at a significance of and 2.2 respectively. When we fit an amplitude, , as a free parameter to the DESSPT cross-correlation, we measure (68% confidence limit), using the Planck 2015 best-fit cosmology to calculate the expected power spectrum. The cross-correlation amplitude for DESPlanck is . Therefore we can conclude that our measurement is consistent with the expected level of cross-correlation.

Two previous works reported detections of the GWL-CMBWL cross-correlation. Hand et al. (2013) used 121 deg of Canada-France-Hawaii Telescope (CFHT) Stripe 82 GWL (12.3 source galaxies per arcmin) and ACT CMBWL to measure the cross-correlation with confidence and a best-fit amplitude of using best-fit 2013 Planck+lensing+WP+high-ell cosmology and using a WMAP9 best-fit cosmology i.e. a result consistent with either the WMAP9 or Planck expectation. In contrast, Liu & Hill (2015) found significantly lower best-fit amplitudes for their cross-correlation of 140 deg of CFHTLenS GWL (12.5 source galaxies per arcmin) and Planck CMBWL. They report and using the 2013 and 2015 Planck releases respectively. The authors speculate about a number of potential systematic effects, including IAs and photo- errors, which could produce the unexpectedly low cross-correlation amplitude. None seemed sufficient to account for the observed discrepancy.

Our results are consistent with expectations from the Planck 2015 cosmology and with Hand et al. (2013), though our errors are larger due to the smaller number density of sources in the DES SV catalogues. We see no evidence of the significantly low cross-correlation amplitudes reported by Liu & Hill (2015), though our result is also consistent with their measurement as the errors overlap, even with their low best-fit amplitude. This is true for both our cross-correlation of DES with SPT and with Planck.

With three measurements of this cross-correlation now existing in the literature, there is not yet reliable evidence for any deviation from the expected LCDM level of cross-correlation, given the size of the statistical uncertainties and the significant impact of systematic errors, particularly IAs. The low best-fit amplitude found by Liu & Hill (2015) seems to be an outlier. We have demonstrated that IAs can shift the expected result by significant fractions of the errors of the current experiments. This underlines the fact that accurate modelling of IAs (Troxel & Ishak, 2014; Hall & Taylor, 2014; Chisari et al., 2015; Larsen & Challinor, 2015) must be an immediate priority as new data increases the precision of this measurement in the coming years.

We tested the various tools and procedures used to produce our measurement for systematic variations that could bias the result. This included checks on our measurement estimator, noise calculation and the re-making of the measurement using a different DES shape catalogue. Results when alternate estimators for both the cross-correlation and the covariance were substituted into our analysis pipeline were found to be consistent with our primary measured result. We have also tested two independent shear measurement pipelines and four independent photometric redshift estimators produced by the DES collaboration. Results from each were entirely consistent with our fiducial analysis choices. No change in the analysis procedure had a significant impact on our measured cross-correlation strength, within the errors.

CMBWL also presents us with physical phenomena that could impact future cross-correlation measurements. One example is the bias on the estimated CMB convergence arising from galaxies and clusters. These non-Gaussian foreground objects induce mode coupling of the observed sky in a way that mimics gravitational lensing. Amblard et al. (2004) estimated that this contamination could be as high as 10% based on large-scale simulations. van Engelen et al. (2014) show that biases arise from galaxy-lensing correlation and cluster-lensing correlation. The exact size of the effect depends on the choice of mask radius and maximal temperature multipole. This level of contamination does not impact our current results greatly, given the magnitude of our uncertainties, so we left it untreated in the current analysis. It will become increasingly important as we move on to DES year one data and beyond, meaning accurate estimation of the strength of this effect should be a priority.

Thinking beyond the calibration of systematic effects, as a cosmological tool GWL-CMBWL cross-correlation will be particularly important as an ingredient in a joint analysis framework with tomographic GWL and large-scale structure (LSS) data sets. Here it will help constrain cosmological parameters and calibrate systematic nuisance parameters, including galaxy bias in LSS measurements, while acting as an additional high redshift ‘bin’. This will increase our sensitivity to phenomena at the upper end of the redshift range of late-time probes. As the power of the combined data sets increases over the coming years, novel techniques to optimise cosmological information should be considered, including optimal weighting schemes to maximise the overlap of the respective sensitivity kernels of GWL, CMBWL and LSS surveys.

## Acknowledgements

We are grateful for the extraordinary contributions of our CTIO colleagues and the DECam Construction, Commissioning and Science Verification teams in achieving the excellent instrument and telescope conditions that have made this work possible. The success of this project also relies critically on the expertise and dedication of the DES Data Management group.

Funding for the DES Projects has been provided by the U.S. Department of Energy, the U.S. National Science Foundation, the Ministry of Science and Education of Spain, the Science and Technology Facilities Council of the United Kingdom, the Higher Education Funding Council for England, the National Center for Supercomputing Applications at the University of Illinois at Urbana-Champaign, the Kavli Institute of Cosmological Physics at the University of Chicago, the Center for Cosmology and Astro-Particle Physics at the Ohio State University, the Mitchell Institute for Fundamental Physics and Astronomy at Texas A&M University, Financiadora de Estudos e Projetos, Fundação Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro, Conselho Nacional de Desenvolvimento Científico e Tecnológico and the Ministério da Ciência, Tecnologia e Inovação, the Deutsche Forschungsgemeinschaft and the Collaborating Institutions in the Dark Energy Survey.

The Collaborating Institutions are Argonne National Laboratory, the University of California at Santa Cruz, the University of Cambridge, Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas-Madrid, the University of Chicago, University College London, the DES-Brazil Consortium, the University of Edinburgh, the Eidgenössische Technische Hochschule (ETH) Zürich, Fermi National Accelerator Laboratory, the University of Illinois at Urbana-Champaign, the Institut de Ciències de l’Espai (IEEC/CSIC), the Institut de Física d’Altes Energies, Lawrence Berkeley National Laboratory, the Ludwig-Maximilians Universität München and the associated Excellence Cluster Universe, the University of Michigan, the National Optical Astronomy Observatory, the University of Nottingham, The Ohio State University, the University of Pennsylvania, the University of Portsmouth, SLAC National Accelerator Laboratory, Stanford University, the University of Sussex, and Texas A&M University.

The DES data management system is supported by the National Science Foundation under Grant Number AST-1138766. The DES participants from Spanish institutions are partially supported by MINECO under grants AYA2012-39559, ESP2013-48274, FPA2013-47986, and Centro de Excelencia Severo Ochoa SEV-2012-0234. Research leading to these results has received funding from the European Research Council under the European Union’s Seventh Framework Programme (FP7/2007-2013) including ERC grant agreements 240672, 291329, and 306478.

This research used resources of the Calcul Quebec computing consortium, part of the Compute Canada network.

The South Pole Telescope program is supported by the National Science Foundation through grant PLR-1248097. Partial support is also provided by the NSF Physics Frontier Center grant PHY-0114422 to the Kavli Institute of Cosmological Physics at the University of Chicago, the Kavli Foundation, and the Gordon and Betty Moore Foundation through Grant GBMF#947 to the University of Chicago.

OL acknowledges support from a European Research Council Advanced Grant FP7/291329, which also supported DK.

ABL thanks CNES for financial support through its post-doctoral programme.

PL is funded jointly by the Royal Society of New Zealand Rutherford Foundation Trust and the Cambridge Commonwealth Trust.

CR acknowledges support from the University of Melbourne.

ES is supported by DOE grant DE-AC02-98CH10886.

## References

- Amblard et al. (2004) Amblard A., Vale C., White M., 2004, New A, 9, 687
- Bacon et al. (2000) Bacon D. J., Refregier A. R., Ellis R. S., 2000, MNRAS, 318, 625
- Bartelmann (2010) Bartelmann M., 2010, Classical and Quantum Gravity, 27, 233001
- Bartelmann & Schneider (2001) Bartelmann M., Schneider P., 2001, Phys. Rep., 340, 291
- Becker et al. (2015) Becker M. R., et al., 2015, preprint, (arXiv:1507.05598)
- Benítez (2000) Benítez N., 2000, ApJ, 536, 571
- Benson et al. (2014) Benson B. A., et al., 2014, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series. p. 1 (arXiv:1407.2973), doi:10.1117/12.2057305
- Blanchard & Schneider (1987) Blanchard A., Schneider J., 1987, A&A, 184, 1
- Bleem et al. (2015) Bleem L. E., et al., 2015, ApJS, 216, 27
- Bonnett (2015) Bonnett C., 2015, Monthly Notices of the Royal Astronomical Society, 449, 1043
- Bonnett et al. (2015) Bonnett C., et al., 2015, preprint, (arXiv:1507.05909)
- Bridle & King (2007) Bridle S., King L., 2007, New Journal of Physics, 9, 444
- Brown et al. (2002) Brown M. L., Taylor A. N., Hambly N. C., Dye S., 2002, MNRAS, 333, 501
- Carlstrom et al. (2011) Carlstrom J. E., et al., 2011, PASP, 123, 568
- Carrasco Kind & Brunner (2013) Carrasco Kind M., Brunner R. J., 2013, MNRAS, 432, 1483
- Carrasco Kind & Brunner (2014) Carrasco Kind M., Brunner R. J., 2014, MNRAS, 442, 3380
- Catelan et al. (2001) Catelan P., Kamionkowski M., Blandford R. D., 2001, MNRAS, 320, L7
- Challinor & Lewis (2011) Challinor A., Lewis A., 2011, Phys. Rev. D, 84, 043516
- Chang et al. (2015) Chang C., et al., 2015, Physical Review Letters, 115, 051301
- Chisari et al. (2015) Chisari N. E., Dunkley J., Miller L., Allison R., 2015, preprint, (arXiv:1507.03906)
- Chon et al. (2004) Chon G., Challinor A., Prunet S., Hivon E., Szapudi I., 2004, MNRAS, 350, 914
- Clampitt et al. (2015) Clampitt J., et al., 2015, in preparation
- Coe et al. (2006) Coe D., Benítez N., Sánchez S. F., Jee M., Bouwens R., Ford H., 2006, AJ, 132, 926
- Crittenden et al. (2001) Crittenden R. G., Natarajan P., Pen U.-L., Theuns T., 2001, ApJ, 559, 552
- Crittenden et al. (2002) Crittenden R. G., Natarajan P., Pen U.-L., Theuns T., 2002, ApJ, 568, 20
- Croft & Metzler (2000) Croft R. A. C., Metzler C. A., 2000, ApJ, 545, 561
- Dark Energy Survey Collaboration (2015) Dark Energy Survey Collaboration 2015, preprint, (arXiv:1507.05552)
- Das et al. (2014) Das S., et al., 2014, J. Cosmology Astropart. Phys, 4, 14
- Desai et al. (2012) Desai S., et al., 2012, ApJ, 757, 83
- Efstathiou (2006) Efstathiou G., 2006, MNRAS, 370, 343
- Flaugher et al. (2015) Flaugher B., et al., 2015, AJ, 150, 150
- Giannantonio et al. (2015) Giannantonio T., et al., 2015, preprint, (arXiv:1507.05551)
- Górski et al. (2005) Górski K. M., Hivon E., Banday A. J., Wandelt B. D., Hansen F. K., Reinecke M., Bartelmann M., 2005, ApJ, 622, 759
- Graff & Feroz (2013) Graff P., Feroz F., 2013, SkyNet: Neural network training tool for machine learning in astronomy, Astrophysics Source Code Library (ascl:1312.007)
- Hall & Taylor (2014) Hall A., Taylor A., 2014, MNRAS, 443, L119
- Hand et al. (2013) Hand N., et al., 2013, preprint, (arXiv:1311.6200)
- Hartlap et al. (2007) Hartlap J., Simon P., Schneider P., 2007, A&A, 464, 399
- Heavens et al. (2000) Heavens A., Refregier A., Heymans C., 2000, MNRAS, 319, 649
- Heymans et al. (2012) Heymans C., et al., 2012, MNRAS, 427, 146
- Hildebrandt et al. (2010) Hildebrandt H., et al., 2010, A&A, 523, A31
- Hirata & Seljak (2003) Hirata C. M., Seljak U., 2003, Phys. Rev. D, 68, 083002
- Hirata & Seljak (2004) Hirata C. M., Seljak U., 2004, Phys. Rev. D, 70, 063526
- Hirata & Seljak (2010) Hirata C. M., Seljak U., 2010, ArXiv Astrophysics e-prints,
- Hivon et al. (2002) Hivon E., Górski K. M., Netterfield C. B., Crill B. P., Prunet S., Hansen F., 2002, ApJ, 567, 2
- Hoekstra & Jain (2008) Hoekstra H., Jain B., 2008, Annual Review of Nuclear and Particle Science, 58, 99
- Jarvis et al. (2015) Jarvis M., et al., 2015, preprint, (arXiv:1507.05603)
- Kaiser (1992) Kaiser N., 1992, ApJ, 388, 272
- Kaiser & Squires (1993) Kaiser N., Squires G., 1993, ApJ, 404, 441
- Kaiser et al. (2000) Kaiser N., Wilson G., Luppino G. A., 2000, arXiv:astro-ph/0003338,
- Klein & Roodman (2005) Klein J. R., Roodman A., 2005, Annual Review of Nuclear and Particle Science, 55, 141
- Kogut et al. (2003) Kogut A., et al., 2003, ApJS, 148, 161
- Larsen & Challinor (2015) Larsen P., Challinor A., 2015, preprint, (arXiv:1510.02617)
- Leistedt et al. (2015) Leistedt B., et al., 2015, preprint, (arXiv:1507.05647)
- Lewis & Challinor (2006) Lewis A., Challinor A., 2006, Phys. Rep., 429, 1
- Liu & Hill (2015) Liu J., Hill J. C., 2015, preprint, (arXiv:1504.05598)
- Mandelbaum et al. (2011) Mandelbaum R., et al., 2011, MNRAS, 410, 844
- Mohr et al. (2012) Mohr J. J., et al., 2012, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series. (arXiv:1207.3189), doi:10.1117/12.926785
- Munshi et al. (2002) Munshi D., Challinor A. D., Doré O., Mortlock D. J., Efstathiou G. P., 2002, ArXiv Astrophysics e-prints,
- Munshi et al. (2008) Munshi D., Valageas P., van Waerbeke L., Heavens A., 2008, Phys. Rep., 462, 67
- Norberg et al. (2009) Norberg P., Baugh C. M., Gaztañaga E., Croton D. J., 2009, MNRAS, 396, 19
- Okamoto & Hu (2003) Okamoto T., Hu W., 2003, Phys. Rev. D, 67, 083002
- Planck Collaboration (2014) Planck Collaboration 2014, A&A, 571, A17
- Planck Collaboration (2015b) Planck Collaboration 2015b, preprint, (arXiv:1502.05956)
- Planck Collaboration (2015c) Planck Collaboration 2015c, preprint, (arXiv:1509.06348)
- Planck Collaboration (2015a) Planck Collaboration 2015a, preprint, (arXiv:1502.01591)
- Sadeh et al. (2015) Sadeh I., Abdalla F. B., Lahav O., 2015, preprint, (arXiv:1507.00490)
- Sánchez et al. (2014) Sánchez C., et al., 2014, MNRAS, 445, 1482
- Saro et al. (2015) Saro A., et al., 2015, preprint, (arXiv:1506.07814)
- Schneider et al. (2002) Schneider P., van Waerbeke L., Mellier Y., 2002, A&A, 389, 729
- Sheldon (2014) Sheldon E. S., 2014, MNRAS, 444, L25
- Smith et al. (2007) Smith K. M., Zahn O., Doré O., 2007, Phys. Rev. D, 76, 043510
- Story et al. (2013) Story K. T., et al., 2013, ApJ, 779, 86
- Story et al. (2015) Story K. T., et al., 2015, ApJ, 810, 50
- Szapudi et al. (2000) Szapudi I., Prunet S., Pogosyan D., Szalay A. S., Bond J. R., 2000, ArXiv Astrophysics e-prints,
- Troxel & Ishak (2014) Troxel M. A., Ishak M., 2014, Phys. Rev. D, 89, 063528
- Vallinotto (2012) Vallinotto A., 2012, ApJ, 759, 32
- Vallinotto (2013) Vallinotto A., 2013, ApJ, 778, 108
- Vikram et al. (2015) Vikram V., et al., 2015, preprint, (arXiv:1504.03002)
- Wittman et al. (2000) Wittman D. M., Tyson J. A., Kirkman D., Dell’Antonio I., Bernstein G., 2000, Nature, 405, 143
- Zuntz et al. (2013) Zuntz J., Kacprzak T., Voigt L., Hirsch M., Rowe B., Bridle S., 2013, MNRAS, 434, 1604
- van Engelen et al. (2012) van Engelen A., et al., 2012, ApJ, 756, 142
- van Engelen et al. (2014) van Engelen A., Bhattacharya S., Sehgal N., Holder G. P., Zahn O., Nagai D., 2014, ApJ, 786, 13
- van Waerbeke et al. (2000) van Waerbeke L., et al., 2000, Astron. Astrophys., 358, 30

Astrophysics Group, Department of Physics and Astronomy, University College London, 132 Hampstead Road, London, NW1 2PS, United Kingdom

Department of Physics, McGill University, 3600 rue University, Montreal, QC, H3A 2T8, Canada

Sorbonne Universités, UPMC Univ Paris 6 et CNRS, UMR 7095, Institut d’Astrophysique de Paris, 98 bis bd Arago, 75014 Paris, France

Fermi National Accelerator Laboratory, Batavia, IL 60510-0500, USA

Kavli Institute for Cosmological Physics, University of Chicago, 933 East 56th Street, Chicago, IL 60637, USA

Department of Astronomy and Astrophysics, University of Chicago, Chicago, IL 60637, USA

Department of Physics, ETH Zurich, Wolfgang-Pauli- Strasse 16, CH-8093 Zurich, Switzerland

Kavli Institute for Cosmology, University of Cambridge, Madingley Road, Cambridge CB3 0HA, United Kingdom

Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, United Kingdom

Institute of Cosmology & Gravitation, University of Portsmouth, Portsmouth, PO1 3FX, UK

Department of Physics, University of Chicago, 5640 South Ellis Avenue, Chicago, IL, USA 60637

Institut de Ciències de lâEspai, IEEC-CSIC, Campus UAB, Carrer de Can Magrans, s/n, 08193 Bellaterra, Barcelona, Spain

Centre for Theoretical Cosmology, DAMTP, University of Cambridge, CB3 0WA, United Kingdom

Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 19104, USA

Jodrell Bank Center for Astrophysics, School of Physics and Astronomy, University of Manchester, Oxford Road, Manchester, M13 9PL, UK

Department of Physics, Stanford University, 382 Via Pueblo Mall, Stanford, CA 94305, USA

Department of Astronomy and Department of Physics, University of Illinois, 1002 West Green St., Urbana, IL 61801

Argonne National Laboratory, High-Energy Physics Division, 9700 S. Cass Avenue, Argonne, IL, USA 60439

Cerro Tololo Inter-American Observatory, National Optical Astronomy Observatory, Casilla 603, La Serena, Chile

Department of Physics and Electronics, Rhodes University, PO Box 94, Grahamstown, 6140, South Africa

Kavli Institute for Particle Astrophysics & Cosmology, P. O. Box 2450, Stanford University, Stanford, CA 94305, USA

Carnegie Observatories, 813 Santa Barbara St., Pasadena, CA 91101, USA

Institut de Física d’Altes Energies (IFAE), The Barcelona Institute of Science and Technology, Campus UAB, 08193 Bellaterra (Barcelona) Spain

SLAC National Accelerator Laboratory, Menlo Park, CA 94025, USA

Laboratório Interinstitucional de e-Astronomia - LIneA, Rua Gal. José Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil

Observatório Nacional, Rua Gal. José Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil

National Center for Supercomputing Applications, 1205 West Clark St., Urbana, IL 61801, USA

School of Physics and Astronomy, University of Southampton, Southampton, SO17 1BJ, UK

Excellence Cluster Universe, Boltzmannstr. 2, D-85748 Garching bei München, Germany

Universitäts-Sternwarte, Fakultät für Physik, Ludwig-Maximilians Universität München, Scheinerstr. 1, 81679 München, Germany

Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Dr., Pasadena, CA 91109, USA

Department of Astronomy, University of Michigan, Ann Arbor, MI 48109, USA

Department of Physics, University of Michigan, Ann Arbor, MI 48109, USA

Department of Astronomy, University of California, Berkeley, 501 Campbell Hall, Berkeley, CA 94720, USA

Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA

Max Planck Institute for Extraterrestrial Physics, Giessenbachstrasse, 85748 Garching, Germany

Center for Cosmology and Astro-Particle Physics, The Ohio State University, Columbus, OH 43210, USA

Department of Physics, The Ohio State University, Columbus, OH 43210, USA

Australian Astronomical Observatory, North Ryde, NSW 2113, Australia

Departamento de Física Matemática, Instituto de Física, Universidade de São Paulo, CP 66318, CEP 05314-970, São Paulo, SP, Brazil

Department of Astronomy, The Ohio State University, Columbus, OH 43210, USA

Department of Astrophysical Sciences, Princeton University, Peyton Hall, Princeton, NJ 08544, USA

Institució Catalana de Recerca i Estudis Avançats, E-08010 Barcelona, Spain

School of Physics, University of Melbourne, Parkville, VIC 3010, Australia

Department of Physics, University of Arizona, Tucson, AZ 85721, USA

Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas (CIEMAT), Madrid, Spain

Instituto de Física Teórica, Universidade Estadual Paulista, Rua Dr. Bento T. Ferraz 271, Sao Paulo, SP 01140-070, Brazil