ʻ ‘
The DES and SPT Collaborations
Joint analysis of DES Year 3 data and CMB lensing from SPT and Planck III:
Combined cosmological constraints
Abstract
We present cosmological constraints from the analysis of two-point correlation functions between galaxy positions and galaxy lensing measured in Dark Energy Survey (DES) Year 3 data and measurements of cosmic microwave background (CMB) lensing from the South Pole Telescope (SPT) and Planck. When jointly analyzing the DES-only two-point functions and the DES cross-correlations with SPT+Planck CMB lensing, we find and , assuming CDM. When additionally combining with measurements of the CMB lensing autospectrum, we find and . The high signal-to-noise of the CMB lensing cross-correlations enables several powerful consistency tests of these results, including comparisons with constraints derived from cross-correlations only, and comparisons designed to test the robustness of the galaxy lensing and clustering measurements from DES. Applying these tests to our measurements, we find no evidence of significant biases in the baseline cosmological constraints from the DES-only analyses or from the joint analyses with CMB lensing cross-correlations. However, the CMB lensing cross-correlations suggest possible problems with the correlation function measurements using alternative lens galaxy samples, in particular the redMaGiC galaxies and high-redshift MagLim galaxies, consistent with the findings of previous studies. We use the CMB lensing cross-correlations to identify directions for further investigating these problems.
I Introduction
The late-time large scale structure (LSS) of the Universe is sensitive to a variety of cosmological signals, ranging from the properties of dark energy and dark matter, to the masses of the neutrinos. Galaxy imaging surveys probe this structure using observations of both the positions of galaxies (which trace the LSS) and the shapes of galaxies (which are distorted by the gravitational lensing effects of the LSS). Several galaxy imaging surveys have used two-point correlations between these measurements to place constraints on cosmological models, including the Kilo Degree Survey (KiDS), the Hyper Suprime Cam Subaru Strategic Program (HSC-SSP), and the Dark Energy Survey (DES) (e.g. Heymans et al., 2021; Hamana et al., 2020; DES Collaboration et al., 2018). DES has recently presented cosmological constraints from the joint analysis of the three two-point correlation functions () between measurements of these probes from the first three years (Y3) of DES data DES Collaboration et al. (2022).
Surveys of the cosmic microwave background (CMB) are also able to probe the late-time LSS through the effects of gravitational lensing. Although CMB photons originate from the last scattering surface at redshift , their paths are perturbed by structure at late times, including the same LSS measured by galaxy surveys. CMB lensing provides a highly complementary probe of structure to galaxy surveys, and cross-correlations between the two have several appealing features. For one, current galaxy imaging surveys (like DES) identify galaxies out to , but the galaxy lensing measurements with these surveys do not have significant sensitivity beyond . Without the high-redshift lensing information, cosmological constraints from galaxy surveys at are therefore significantly degraded. CMB lensing, however, reaches peak sensitivity at . Therefore, by cross-correlating galaxy surveys with CMB lensing measurements, it is possible to obtain high-precision measurements of the evolution of the matter distribution over a broader range of redshifts than by using galaxy surveys alone. Cross-correlations of galaxy surveys with CMB lensing measurements are also expected to be robust to certain types of systematic biases. Because the galaxy survey measurements are so different from the CMB lensing measurements (e.g., they use data measured by different telescopes at different wavelengths, and use different estimators for the lensing signal), biases in the galaxy surveys are unlikely to correlate with biases in CMB lensing, making two-point functions between the two especially robust. Finally, cross-survey correlations often have different parameter dependencies than correlations within a survey, offering the possibility of improved parameter constraints via degeneracy breaking in joint analyses.
The prospect of obtaining tighter and more robust cosmological constraints from the late-time matter distribution via cross-correlations is particularly timely given recent hints of tensions between some cosmological probes. In particular, recent observations of late-time structure from galaxy surveys tend to prefer lower values of than CMB surveys Battye et al. (2015); MacCrann et al. (2015); Raveri (2016); Raveri and Hu (2019); Aghanim et al. (2020). This tension could result from physics beyond the standard cosmological constant and cold dark matter model (CDM), or it could result from systematic biases in the analyses. By cross-correlating galaxy surveys with CMB lensing, we obtain an independent handle on the late-time large scale structure measurements that can be used to investigate the origins of this possible tension Krolewski et al. (2021); Robertson et al. (2021); Chang et al. (2022); White et al. (2022). Recent analyses have also suggested the possibility of systematic biases in galaxy survey measurements Chang et al. (2019). Because cross-correlations between galaxy surveys and CMB lensing are robust to many important sources of systematic error, they provide a powerful way to ensure that late-time measurements of structure are unbiased.
This work presents the joint cosmological analysis of two-point correlations between galaxy positions and galaxy lensing measured in DES data, and CMB lensing measurements from the South Pole Telescope (SPT, Carlstrom et al., 2011) and the Planck satellite (Planck Collaboration, 2011). As part of its 2008-2011 SPT-SZ survey, SPT obtained high-resolution and high-sensitivity maps of the CMB that partially overlap with the full DES footprint (Omori et al., 2017). At somewhat lower sensitivity and resolution, Planck has obtained full-sky maps of the CMB that overlap completely with the DES footprint. Together, these CMB maps enable high signal-to-noise estimation of the CMB lensing signal across the entire DES footprint Omori et al. (2017, 2022), presenting an opportunity for cross-correlation studies.
From the measurements of galaxy positions (used to compute the galaxy overdensity, ), galaxy lensing (, or for the tangential shear), and CMB lensing (), it is possible to form six two-point functions: galaxy clustering (), galaxy-galaxy lensing (), cosmic shear (), galaxy density-CMB lensing cross-correlation (), galaxy shear-CMB lensing cross-correlation (), and the CMB lensing auto-correlation (). All six of the above will be considered here (hereafter, we refer to this combination as ). The five two-point functions excluding the CMB lensing auto-correlation (referred to as ) all probe structure below about at , and are highly correlated. This combination, which we measure using DES, SPT and Planck data is the primary focus of this work. The CMB lensing autocorrelation measurements used in this analysis are derived from all-sky Planck data, and are minimally correlated with the measurements owing to their small (fractional) sky overlap and sensitivity to higher redshifts Omori et al. (2022). We therefore treat the CMB lensing autocorrelation as an external probe, and combine it with at the likelihood level.
As highlighted above, one of the key reasons to consider cross-correlations of galaxy surveys with CMB lensing is to improve robustness to systematic uncertainties. We will therefore also analyze various subsets of the probes for the purposes of testing robustness and exploring sensitivity to possible systematic errors. Of particular interest for these tests is the unexpected discovery in DES Y3 data of discrepancies in the galaxy bias values preferred by the clustering and lensing measurements. The DES Y3 analysis considered two galaxy samples for the purposes of measuring : MagLim and redMaGiC. The MagLim galaxies at were used for the baseline cosmological results presented in DES Collaboration et al. (2022). Surprisingly, the galaxy bias values inferred for redMaGiC galaxies from their clustering were found to be roughly 10% lower than the bias values inferred from lensing Pandey et al. (2021), with this discrepancy increasing for the highest-redshift galaxies. There is no known physical explanation for this discrepancy, but tests in Pandey et al. (2021) suggest that it may be connected to observational systematics imparting additional clustering power. MagLim galaxies at high redshift () also showed a discrepancy between clustering and lensing Porredon et al. (2021). Further investigating these discrepancies is one of the main goals of the present analysis.
The analysis presented here makes several significant improvements relative to previous cross-correlation analyses between DES and SPT/Planck measurements of CMB lensing Baxter et al. (2016); Kirk et al. (2016); Giannantonio et al. (2016); Baxter et al. (2018); Omori et al. (2019a, b); Abbott et al. (2019). First, the DES data have significantly expanded in going from Y1 observations to Y3, covering roughly a factor of three larger area. Second, the CMB lensing maps from SPT/Planck have been remade with several improvements (described in more detail in Omori et al. (2022)). Foremost among these is that we have used the CMB lensing estimator from Madhavacheril and Hill (2018) to reduce contamination in the lensing maps from the thermal Sunyaev-Zel’dovich (tSZ) effect. This contamination was the dominant source of systematic uncertainty for the analysis of (Abbott et al., 2019), and required us to remove a significant fraction of the small-scale measurements from our analysis to ensure that our results were unbiased. As a result, the total signal-to-noise of the CMB lensing cross-correlations was significantly reduced. Using a CMB lensing map that is immune to contamination from the tSZ effect allows us to extract signal from a wider range of angular scales and hence improve our signal-to-noise ratio. Finally, we have also implemented several improvements to the modeling of the correlation functions, which are described in more detail in Omori et al. (2022).
The analysis presented here is the last in a series of three papers: In (Omori et al., 2022, hereafter Paper I) we described the construction of the combined, tSZ-cleaned SPT+Planck CMB lensing map and the methodology of the cosmological analysis. In (Chang et al., 2022, hereafter Paper II), we presented the measurements of the cross-correlation probes , a series of diagnostic tests of the measurements, and cosmological constraints from this cross-correlation combination. In this paper (Paper III), we present the joint cosmological constraints from all the probes, and tests of consistency between various combinations of two-point functions.
The plan of the paper is as follows. In §II we describe the data sets from DES, SPT and Planck that we use in this analysis, and in §III we provide an abridged summary of our model for the correlation function measurements. In §IV, we present cosmological constraints from the joint analysis of cross-correlations between DES and CMB lensing measurements from SPT and Planck, and discuss several tests of the robustness of these constraints enabled by the cross-correlation measurements. We conclude in §V.
II Data from DES, SPT and Planck
DES (Flaugher, 2005) is a photometric survey in five broadband filters (), with a footprint of nearly of the southern sky, imaging hundreds of millions of galaxies. It employs the 570-megapixel Dark Energy Camera (DECam, Flaugher et al., 2015) on the Cerro Tololo Inter-American Observatory (CTIO) 4m Blanco telescope in Chile. We use data from the first three years (Y3) of DES observations. The foundation of the various DES Y3 data products is the Y3 Gold catalog described in Sevilla-Noarbe et al. (2021), which achieves a depth of S/N10 for extended objects up to i23.0 over an unmasked area of . In this work, we consider two types of galaxy samples: lens galaxies that are used as biased tracers of the underlying density field, and source galaxies which are used to measure the shape-distorting effects of gravitational lensing. We use the same galaxy samples as in the DES analysis (DES Collaboration et al., 2022). That is, the lens galaxies are taken from the four-redshift bin MagLim sample described in (Porredon et al., 2021), and the source galaxy shapes are taken from the four-redshift bin MetaCalibration sample described in (Gatti et al., 2021). We will additionally consider lens galaxies from the redMaGiC sample described in (Pandey et al., 2021). In particular, we will investigate the potential systematic biases that led to that sample not being used as the baseline cosmology sample in DES Collaboration et al. (2022). The redshift distributions for the MagLim, redMaGiC, Metacalibration samples are shown in Fig. 1.
As mentioned above, we use two CMB lensing maps in this work: one covering the SPT-SZ footprint that uses data from SPT-SZ and Planck (with an overlapping area of deg2), and a second that covers the northern part of the DES survey that uses only Planck data (with an overlapping area of deg2). Together, these two CMB lensing maps cover the full DES Y3 survey region. Since the noise levels and beam sizes of SPT-SZ and Planck are different, the resulting CMB lensing maps must be treated separately in our analysis. In Paper II we tested the consistency between the cosmological constraints from these two patches, finding good agreement.

III Modeling and Measurements
The theoretical framework we use in this analysis is laid out in Paper I and Krause et al. (2021). The full data vector consists of six two-point functions. Since there is little correlation between and the CMB lensing autocorrelation measurements from Planck, we combine the corresponding constraints at the likelihood level; this approximation is validated in Paper I.
We fit the data to two different cosmological models: a spatially flat, cosmological constant and cold dark matter model, and a cosmological model where the equation of state parameter of dark energy, , is additionally allowed to vary. Following the DES convention, we will refer to these models as CDM and CDM; note, though, that we allow the sum of the neutrino masses to vary in both of these analyses.
The modeling of the correlations begins with the auto and cross-power spectrum of the three fields (, , ). For the correlation functions other than galaxy clustering, we use the Limber approximation Limber (1953):
(1) |
where , labels the redshift bin, is the non-linear matter power spectrum, which we compute using CAMB and Halofit (Lewis et al., 2000; Takahashi et al., 2012), is the comoving distance from the observer, and is the redshift corresponding to . The weighting functions, , describe how the different probes respond to large-scale structure at different distances, and are given by
(2) |
(3) |
(4) |
where and are the Hubble constant and matter density parameters, respectively, is the scale factor corresponding to comoving distance , is galaxy bias as a function of scale () and redshift, are the normalized redshift distributions of the lens/source galaxies in bin . denotes the comoving distance to the CMB last scattering surface. The model for galaxy clustering is computed without using the Limber approximation, as described in Krause et al. (2021). The angular-space correlation functions are then computed from the auto- and cross-spectra as described in Krause et al. (2021); Omori et al. (2022).
In addition to the basic modeling, described above, we also consider several other physical and observational effects. We list these below but refer the readers to Paper I and Krause et al. (2021) for details.
- •
- •
-
•
Lens magnification: Gravitational lensing by foreground mass changes the observed projected number density of lens galaxies as a result of geometric dilution and modulation of galaxy flux and size. We model this effect based on measurements in simulations as described in (Elvin-Poole et al., 2022; Krause et al., 2021).
-
•
Redshift uncertainties: There are uncertainties associated with the estimation of the redshift distributions of the lens and the source sample, which we model as described in (Myles et al., 2021; Gatti et al., 2022; Cawthon et al., 2020). In (Cordero et al., 2022), an alternate approach to marginalizing over uncertainties in the redshift distributions was also considered, which we explore in Appendix C.
- •
-
•
CMB map filtering: In order to suppress very small-scale noise in the CMB lensing cross-correlations, we apply filtering to the CMB lensing maps. This filtering is included in the model as described in Paper I.
-
•
Point mass marginalization: The correlation functions at small scales are impacted by baryonic effects that are challenging to model, such as galaxy formation. This is particularly problematic for : changes in e.g. the masses of the lens galaxies at very small scales can impact the large-scale because tangential shear is a non-local quantity. To reduce sensitivity of our analysis to small-scale effects in , we therefore adopt the point mass marginalization approach of (MacCrann et al., 2020), which involves modifying the covariance matrix of .
We measure the two-point angular correlation functions of the data using the fast tree-based algorithm TreeCorr Jarvis et al. (2004) as described in (Rodríguez-Monroy et al., 2022; Prat et al., 2022; Amon et al., 2022; Secco et al., 2022; Chang et al., 2022). The shear measurements define a spin-2 field on the sky, and there are several ways of decomposing this field for the purposes of measuring two-point functions. For measuring , we use the and decomposition, while for measuring , we consider the correlation only with tangential shear, (Bartelmann and Schneider, 2001). The covariance matrix associated with these measurements is constructed by combining an analytical halo model covariance, analytical lognormal covariance, and empirical noise estimation from simulations (Friedrich et al., 2021; Omori et al., 2022).
For the final parameter inference, we assume a Gaussian likelihood.111See e.g. (Lin et al., 2020) for tests of the validity of this assumption in the context of cosmic shear, which would also apply here. The priors imposed on the model parameters are shown in Table 2 in Appendix A. The modeling and likelihood framework is built within the CosmoSIS package Zuntz et al. (2015). We generate parameter samples using the nested sampler PolyChord Handley et al. (2015).
Due to uncertainties in the modeling of the correlation functions on small scales (e.g., nonlinear galaxy bias and baryonic effects on the matter power spectrum), in our likelihood analysis we remove the small-scale measurements that could potentially bias our cosmological constraints. The procedure of determining these “scale cuts” is described in (Krause et al., 2021) and Paper I. Note that the choice of angular scales used in the analysis varies somewhat depending on whether we assume a linear or nonlinear galaxy bias model. We focus on the results with linear bias, but consider the results from the nonlinear bias analysis in Appendix B.
In each of the cosmological analyses performed in this work, we include a separate likelihood constructed using a set of ratios of galaxy-galaxy lensing measurements on small scales (Sánchez et al., 2022). These lensing ratios are found to primarily constrain parameters describing the intrinsic alignment model and redshift biases, and are effectively independent of the data vector.
We utilize two different statistical metrics to assess the consistency of the DES and CMB-lensing cross-correlation measurements, both internally (i.e. between the different two point functions that we measure) and with other cosmological probes. To assess internal consistency, we primarily rely on the posterior predictive distribution (PPD) methods described in Doux et al. (2021). For these assessments, we will quote -values, with taken as significant evidence of inconsistency. To assess external consistency, we rely on the parameter difference methods developed in Raveri and Doux (2021). For this metric, we will quote differences between parameter constraints in terms of effective values, corresponding to the probability values obtained from the non-Gaussian parameter difference metric. When computing the goodness of fit of our measurements to a particular model, we again rely on the PPD methodology, as discussed in Doux et al. (2021). In this case, the associated -values can be thought of as a generalization of the classical -value computed from the statistic that correctly marginalizes over parameter uncertainty.
IV Cosmological constraints
Probe | GoF -value | Comments | |||
---|---|---|---|---|---|
0.023 | DES Collaboration et al. (2022) | ||||
0.50 | CMB lensing cross-correlations, Paper II | ||||
++ | 0.063 | All cross-correlations, §IV.2.2 | |||
0.062 | §IV.1 | ||||
— | §IV.1 |
IV.1 Baseline cosmological constraints
IV.1.1 CDM
We first present constraints on CDM from the joint analysis of two-point functions involving DES galaxy position and lensing measurements, and measurements of CMB lensing from SPT and Planck. Following DES Collaboration et al. (2022), all of the results in this subsection use the four redshift bin MagLim lens galaxy sample.
Fig. 2 shows how the constraints from the CMB lensing cross-correlations compare to those from . The resulting 68% credible intervals on , , and computed from the marginalized and + posteriors are summarized in Table 1. In the same table, we list the goodness of fit -values for and + , computed using the PPD formalism. As noted in DES Collaboration et al. (2022), the goodness of fit for alone is not particularly high, but is still above our threshold of . The goodness of fit for + is acceptable, as described in Chang et al. (2022). While the cross-correlations prefer somewhat lower and higher , they are statistically consistent with . Using the PPD formalism, we find when comparing the two data subsets, indicating acceptable consistency. We are therefore justified in combining the constraints to form , shown with the teal contours in the figure.
Given the weaker constraining power of relative to , the constraints are not much tighter than the constraints: we find an improvement of roughly 10% in the precision of the marginalized constraints on and (see Table 1). The goodness of fit for the full data vector is , indicating an acceptable fit.
In Fig. 3 we compare the constraints from with those from the CMB lensing autospectrum . Owing to the high redshift of the CMB source plane, the CMB lensing-only contour has a different degeneracy direction than , resulting in a weaker constraint when projecting to the direction, but a comparable constraint in the direction. While the CMB lensing autospectrum prefers somewhat higher than , the constraints are generally consistent. Because the CMB lensing autospectrum measurements are treated as an independent probe, we quantify the tension between these measurements and using the parameter shift metric, finding a difference of , indicating no evidence of significant tension. We therefore combine the two to generate constraints from all six two-point functions, , shown with the orange contour in the figure. Due to degeneracy breaking, the joint analysis leads to notably tighter constraints on both and . The 1D posterior constraints on these parameters from are summarized in Table 1. Fig. 3 also shows constraints from Planck measurements of CMB temperature and polarization fluctuations Aghanim et al. (2020). We will assess consistency between our measurements and the Planck measurements in §IV.3.



IV.1.2 CDM
We now consider constraints on CDM, the cosmological model with a constant equation-of-state parameter of dark energy, . The constraints from , and are shown in Fig. 4. We find that there is little improvement in constraining power on CDM when adding the CMB lensing cross-correlations to . Adding the correlation, however, significantly impacts the constraints, presumably because this correlation function adds additional information about structure at . The analysis yields , , and . Therefore, the constraints on the dark-energy equation of state parameter are largely consistent with the cosmological-constant scenario of , and the constraints on and are consistent with those obtained assuming CDM.
IV.2 Robustness tests
In addition to improving cosmological constraints relative to the DES-only analysis, a significant motivation for cross-correlating DES with CMB lensing is to test the robustness of the DES-only constraints. The cross-correlations probe the same large-scale structure as the DES analysis, but with sensitivity to different potential sources of systematic bias, making them powerful cross-checks on the DES results. In this section, we subject the data vector to several tests of internal consistency.

IV.2.1 vs.
We first assess the internal consistency of the combination of probes by comparing constraints from to the other three two-point functions making up , which we call (i.e. ++). This comparison is shown in Fig. 5. We find that the constraining power from is very similar to that of . Because does not constrain galaxy bias or intrinsic alignment parameters very well, applying the PPD methodology to test consistency between and is not well motivated. However, we note that we have already tested the consistency of with + (i.e. part of ), finding acceptable agreement ().
Fig. 5 makes it clear why prefers a somewhat higher value of than . It is not the case that prefers a higher value of than ; indeed, the opposite is true. Rather, the slightly high value of found for is caused by the fact that and have somewhat different degeneracy directions, and intersect at a high value of for both probes.
IV.2.2 Cross-correlations
Cross-correlations between different observables are generally expected to be more robust to systematic biases than auto-correlations of those observables. Additive systematics that impact a single observable are expected to drop out of cross-correlations with another observable that has uncorrelated systematics. In Fig. 6 we compare the cosmological constraints obtained from only cross-correlations to those from the full . It is clear that removing the information from the auto-correlations — particularly cosmic shear — degrades the constraints somewhat. However, we find that the value of inferred only from cross-correlations is consistent with that inferred from the full analysis. This suggests that additive biases are unlikely to be having a major impact on the DES cosmology results. Using the PPD formalism to evaluate the goodness of fit of the cross-correlations conditioned on the posterior from , we find , indicating an acceptable level of consistency between the constraints and the cross-correlations measurements.

IV.2.3 Lensing only

The relationship between galaxy overdensity and the underlying matter field — galaxy bias — presents a significant challenge for analyses of the galaxy distribution. The baseline results presented in DES Collaboration et al. (2022) and the baseline cross-correlation results presented here assume a linear galaxy bias relation when modeling the galaxy field. This model is known to break down at small scales, as investigated for the DES galaxy samples in (Pandey et al., 2021). More complex bias models, such as the perturbation theory-motivated model developed in Pandey et al. (2020), are also expected to have a limited range of validity. There is therefore value in performing analyses that use only lensing information.
Another motivation to consider lensing-only analyses is that the DES galaxy overdensity measurements made with the redMaGiC and high-redshift MagLim galaxies show evidence of systematic biases (i.e. the samples shown with dashed lines in Fig. 1). Measurements of galaxy-galaxy lensing with the redMaGiC galaxies were shown to be inconsistent with clustering measurements using those galaxies Pandey et al. (2021). This inconsistency suggests a potential problem with the redMaGiC overdensity measurements, although it is not clear whether such issues could be impacting the galaxy-galaxy lensing measurements, clustering measurements, or both. Similarly, galaxy-galaxy lensing and clustering measurements with the high-redshift MagLim galaxies were also found to be mutually inconsistent, contributing to a very poor goodness of fit to any of the cosmological models considered. For these reasons, the high-redshift MagLim galaxies were removed from the cosmological analysis in DES Collaboration et al. (2022). These issues, which we investigate further in §IV.2.6, further motivate a cosmological analysis that does not rely on galaxy overdensity measurements.
In Fig. 7, we present cosmological constraints from gravitational lensing only, namely the two-point functions of galaxy lensing and CMB lensing, and their cross-correlation. The lensing-only analysis obtains cosmological constraints that are of comparable precision to those from the full analysis. We find that the lensing-only analysis yields a constraint on that is in excellent agreement with the baseline analysis.
IV.2.4 No galaxy lensing

We also consider the constraints that result from those probes that do not involve galaxy lensing. The galaxy lensing measurements could in principle be biased by systematic errors in photometric redshifts of the source galaxies, shear calibration, or an incorrect intrinsic alignment model. Such issues could bias constraints involving galaxy lensing, but would not impact the galaxy overdensity or CMB lensing measurements. Fig. 8 shows the constraints that result only from probes that do not include galaxy lensing (i.e. and ). Again, we find that the results are consistent with those of . Fig. 8 also shows the + constraints for comparison (i.e. lensing only, but excluding , which receives contributions from higher redshifts than the other two-point functions). We find that the constraints involving lens galaxy overdensities are consistent with the lensing-only constraints.


IV.2.5 Shear calibration
A potentially significant source of systematic uncertainty impacting cosmological constraints from cosmic shear is biases in shear estimation (Hirata and Seljak, 2003). Typically, estimators of lensing shear are calibrated via application to simulated lensed galaxy images. For the DES Year 3 cosmological analysis, calibration of shear biases is described in (MacCrann et al., 2022). While this approach can be used to place tight constraints on shear biases, it has the disadvantage of relying on simulated data. A mismatch between the simulated galaxies used to calibrate the shear estimators and real galaxies could potentially introduce systematic bias.
As pointed out in Vallinotto (2012); Baxter et al. (2016); Schaan et al. (2017), joint analyses of cross-correlations between galaxy surveys and CMB lensing measurements offer the potential of constraining shear calibration biases using only the data. To explore this idea, we repeat our analysis of the and data vectors using very wide, flat priors on the shear calibration parameters: .
The results of this analysis are shown in Fig. 9. Removing the tight priors on the significantly weakens the cosmological constraints from , especially the constraint on . This is because both and impact the amplitude of the lensing correlation functions, leading to strong degeneracy between the two. The shear calibration parameters are also very poorly constrained without the tight priors. However, when the CMB lensing cross-correlations are analyzed jointly with (i.e. forming ), the analysis becomes significantly more robust to shear calibration. Removing the priors on weakens the cosmological constraints, but not nearly as much as for : Removing the priors degrades the constraints on by a factor of 4.7 for , but only by a factor of 2.3 for (see right panel of Fig. 9). The resulting cosmological constraints are consistent with those in the baseline analysis, providing evidence that the DES Y3 and constraints are robust to shear calibration biases. We also find that the data vector achieves constraints on at roughly the 5–10% level depending on the redshift bin, roughly a factor of two improvement over the Y1 analysis presented in (Abbott et al., 2019).


IV.2.6 Investigating the systematic
As noted previously, analyses of and measured with the DES Y3 MagLim Porredon et al. (2021) and redMaGiC Pandey et al. (2021) galaxy samples uncovered discrepancies between the values of galaxy bias preferred by these two correlation functions. The measurements with MagLim galaxies in the two highest redshift bins (i.e. those shown with the dashed lines in the top panel of Fig. 1) prefer higher bias values than by roughly 40 to 60%. Measurements of with the redMaGiC galaxies, on the other hand, show roughly 10% higher bias values than the measurements for the first four redshift bins, with this discrepancy increasing to roughly 40% for the highest redshift bin. In principle, some difference between the bias values inferred from and could result from stochastic biasing (e.g. Baumann et al., 2013). However, the amplitude of the difference seen for the redMaGiC galaxies and the high-redshift MagLim galaxies (roughly 10 to 40% percent) is significantly larger than expected from stochasticity (a few percent) Desjacques et al. (2018); Pandey et al. (2020). In Pandey et al. (2021), a new parameter, , was introduced to explore this effect:
(5) |
where () is the bias parameter for () in lens galaxy redshift bin . The finding that some of the clustering measurements prefer a higher value of galaxy bias than the galaxy-galaxy lensing measurements amounts to a preference for when we expect .
The galaxy-CMB lensing cross-correlations also constrain galaxy bias, providing another handle on the anomalous values of the parameter seen with the redMaGiC and high-redshift MagLim galaxies. We show constraints on the galaxy bias parameters of the MagLim and redMaGiC galaxies from three combinations of probes in Fig. 10 and Fig. 11, respectively. Each of the plotted constraints uses the combination of and — which are effectively independent of the lens galaxies — to constrain the cosmology. The remaining probe is then chosen to be , , or , and this probe is used to constrain the galaxy bias.222This analysis is similar to that presented in Chang et al. (2022), but differs in that we have allowed cosmological parameters to vary, and have included the and measurements (in effect letting the data constrain the cosmological model).
For the two highest redshift bins of MagLim galaxies, we see from Fig. 10 that the measurements prefer higher values of galaxy bias than the measurements, consistent with the preference for described above. Interestingly, it appears that the measurements prefer galaxy bias values more in line with the measurements. This suggest that the preference for is likely driven by . This is perhaps not surprising, given the large residuals of the model fits to seen in Porredon et al. (2021). However, note that there is no obvious reason for a possible failure of the baseline model to fit . As a cross-correlation, the measurements are expected to be quite robust to many observational systematics. Moreover, any systematic impacting would likely show up even more strongly in , and any systematic impacting would likely show up more strongly in . Another possibility is a failure in modeling some physical effect. One such effect is lens magnification, which is known to have a significant impact on the correlations at high redshifts Porredon et al. (2021).
Fig. 11 shows the analogous bias constraints for redMaGiC galaxies. In this case, we see that and measurements both prefer consistently lower values of galaxy bias than , with this difference particularly pronounced in the last redshift bin. This suggests that a possible cause of the redMaGiC preference for is in the measurements. In the case of , it is possible that some observational systematic is modulating the redMaGiC galaxy overdensity field, resulting in a higher than expected clustering amplitude and thus a preference for higher galaxy bias. Such a systematic in the measurements would be expected to have a less noticeable impact on . At the same time, it should be emphasized that the analysis of Rodríguez-Monroy et al. (2022) extensively tested the redMaGiC sample for possible contamination by various observational systematics. While some correlation of known systematics with galaxy density is detected, this correlation is corrected using galaxy re-weighting. It therefore appears to be difficult to explain the anomalous values with any known observational systematic.

The interpretation of the redMaGiC preference for in terms of a systematic impacting redMaGiC measurements is supported by tests with a modified redMaGiC galaxy sample presented in Pandey et al. (2021). The nominal redMaGiC galaxy sample is selected by requiring that galaxies match a red sequence template, as measured by . In Pandey et al. (2021), an alternative, “broad ” sample of galaxies was selected by relaxing the threshold for selection. One would expect that if an observational systematic is modulating the photometry of galaxies, it should have a smaller impact on the “broad ” sample than on the nominal sample. Indeed, it was found that for this alternate sample, the preference for seen for the first four redshift bins disappears. While it might seem surprising that the preference for is possibly driven by two different factors for MagLim and redMaGiC galaxies, this interpretation seems consistent with the observed redshift trends. It may be that observational systematics in redMaGiC galaxy selection are impacting the bias values inferred from at low redshift, while problems in modeling are impacting the bias values inferred for MagLim from at high redshift. The redMaGiC galaxies may be less affected by this latter systematic, as they do not extend to the high redshifts probed by the last two redshift bins of the MagLim sample. We note, though, that even for redMaGiC, the CMB lensing cross-correlations prefer higher galaxy bias than in the highest redshift bin; this could be suggesting that the same problem impacting the high-redshift MagLim galaxies is impacting the high-redshift redMaGiC galaxies. This interpretation would be consistent with mismodeling of at high redshift.
The impact of the apparent systematic in the redMaGiC sample is also noticeable when the cosmological model is changed from CDM to CDM. While the redMaGiC constraints on CDM are quite robust to allowing the parameter to vary, the constraints on CDM shift significantly when this additional freedom is introduced. This is perhaps not surprising given that the systematic biases with redMaGiC appear to be redshift-dependent, and might therefore be somewhat degenerate with the effects of .

Since our analysis above suggests that the problems with redMaGiC may be isolated to the clustering measurements, in Fig. 12 we present constraints on CDM from the combination of probes without the clustering measurements. Interestingly, we see that there is a significant shift in the constraints on relative to the analysis. The constraints without the clustering measurements are in good agreement with the MagLim constraints. This lends additional support to the idea that the redMaGiC clustering measurements may be systematically biased.
To summarize the above discussion, our analysis with CMB lensing cross-correlations suggests that there may be two different sources for the systematic seen with redMaGiC and MagLim galaxies. For redMaGiC galaxies, our analysis suggests a possible bias in the clustering measurements across all redshift bins. Such a bias could conceivably be caused by some observational systematic impacting the redMaGiC selection, which would be consistent with tests performed in Pandey et al. (2021). At the same time, high-redshift MagLim galaxies (and possibly high-redshift redMaGiC galaxies as well) show evidence of a potentially different systematic error that favors a problem with the fits. Such an issue could conceivably be caused by a problem with the modeling, such as an incorrect prescription for magnification effects, which become more pronounced at high redshifts.
IV.3 Consistency with Planck primary CMB measurements
As seen in Fig. 3, we find that the cosmological constraints on CDM from and are not in significant tension with the constraints from the primary CMB measurements of Planck. In particular, we compare our constraints to those from the combination of Planck , , and low- -mode polarization measurements (Planck TT+TE+EE+lowE) Aghanim et al. (2020). Note that we do not include Planck measurements of the CMB lensing power spectrum in this combination. Using the tension metric of Raveri and Doux (2021), we find that the , , and constraints are in agreement with Planck at the level of , , and , respectively. The fact that and are roughly equally consistent with Planck is not surprising, given that the constraints are quite close to those of . Interestingly, while the constraints are significantly tighter than , the level of consistency with Planck remains roughly the same. This results from the preference by for somewhat higher values of , as seen in Fig. 3. Fig. 13 directly compares the and constraints from these and other two-point function combinations, assuming CDM. We note that for consistency with our analysis, we vary the sum of the neutrino masses and impose the priors shown in Table 2 when generating the Planck primary CMB constraints shown in this figure.
V Summary
We have presented cosmological constraints from an analysis of two-point correlation functions between measurements of galaxy positions and galaxy lensing from DES Y3 data, and CMB lensing measurements from SPT and Planck. Our main cosmological constraints are summarized in Table 1.
The high signal-to-noise of the CMB lensing cross-correlation measurements using DES Y3, SPT-SZ and Planck data enables powerful robustness tests of our cosmological constraints. The results of several of these tests are shown in Fig. 13. We summarize the main findings of these tests below:
-
•
The goodness of fit of CDM to the data vector is acceptable (), and the corresponding parameter constraints are consistent with those from measurements by Planck.
-
•
Using only cross-correlations between DES and CMB lensing, we obtain constraints on that are comparable in precision and consistent with the baseline results. This result suggests that additive systematics are not significantly impacting the cosmological constraints.
-
•
Using only gravitational lensing (i.e. no information from galaxy overdensities) yields constraints in agreement with the baseline results. This result suggests that potential systematics impacting the DES galaxy samples, as well as modeling of galaxy bias, are not significantly biasing the cosmological constraints.
-
•
The cosmological constraints from two-point functions of MagLim galaxy overdensity measurements and CMB lensing are generally consistent with the baseline analysis. This result suggests that shear systematics and modeling of galaxy lensing are not significantly biasing the cosmological constraints. We do, however, observe a low-significance increase in when considering only those two-point functions that do not involve galaxy lensing. This shift is driven by the intersection of the + and constraints, and is not present when considering + alone.
-
•
Without priors on shear calibration, the cosmological constraints on from are in good agreement with the baseline results. The data calibrate the shear bias parameters at the 5–10% level, and yield constraints consistent with our nominal priors. These results suggest that shear calibration biases are not significantly impacting the cosmological constraints.
-
•
The constraints on from the different analysis variations are generally consistent. Although the analysis of ++ prefers a somewhat lower value of , this combination of probes is statistically consistent with .
The cosmological constraints from the , , and analyses therefore appear remarkably robust to possible systematic biases.
Assessing the consistency between our constraint on CDM and those of Planck, we find that the and constraints are statistically consistent with Planck at the level, as assessed using the full, multi-dimensional posteriors from these measurements. As seen in Fig. 13, however, essentially all combinations of two point functions that we consider prefer lower values than Planck. Note, though, that there is significant covariance between some of these measurements.
We have also investigated possible issues with the analysis of alternate lens galaxy samples, namely the high-redshift MagLim galaxies and the redMaGiC galaxies. Evidence for biases when analyzing correlation functions measured with these samples was found previously in Porredon et al. (2021), Pandey et al. (2021), Elvin-Poole et al. (2022), and DES Collaboration et al. (2022). The CMB lensing cross-correlations considered here provide a powerful way to probe the sources of these biases. In the context of CDM, our analysis of CMB lensing cross-correlations suggests a possible problem in the modeling of at high redshift for the MagLim galaxies, and possibly the redMaGiC galaxies as well. At the same time, the measurements with redMaGiC suggest a possible observational systematic that impacts redMaGiC galaxy clustering across all redshifts. This interpretation is supported by tests with an alternate redMaGiC galaxy sample in Pandey et al. (2021). In the context of CDM, the measurements with redMaGiC have previously shown to yield constraints inconsistent with the MagLim analysis, and a preference for surprisingly less negative . We show that analysis of +++ (i.e. two-point functions between DES and CMB lensing, excluding galaxy clustering) measured with redMaGiC yields cosmological constraints that are in better agreement with MagLim, and do not show a strong preference for . Finally, we note that while the analyses presented here suggest possible interpretations of the bias, more work with current and future DES data is needed to clarify the true source of this systematic uncertainty.
As the data volume and quality from cosmological surveys continue to improve, we expect similar cross-correlation analyses between galaxy surveys and CMB lensing measurements to play an important role in constraining late-time large scale structure. Excitingly, we expect constraints from such measurements to improve dramatically in the very near future with Year 6 data from DES and new CMB lensing maps from SPT-3G Sobrin et al. (2022) and AdvACT Henderson et al. (2016). These measurements should help to provide a clearer picture of any possible tension. Looking farther forward, cross-correlations between surveys such as the Vera Rubin Observatory Legacy Survey of Space and Time Ivezić et al. (2019); LSST Science Collaboration et al. (2009), the Nancy Grace Roman Space Telescope Doré et al. (2019) , the ESA Euclid mission Laureijs et al. (2011), Simons Observatory Ade et al. (2019), and CMB-S4 Abazajian et al. (2016) will enable significantly more powerful cross-correlation studies that will deliver some of the most precise and accurate cosmological constraints, and that will allow us to continue stress-testing the concordance CDM model.
Acknowledgements.
The South Pole Telescope program is supported by the National Science Foundation (NSF) through the grant OPP-1852617. Partial support is also provided by the Kavli Institute of Cosmological Physics at the University of Chicago. Argonne National Laboratory’s work was supported by the U.S. Department of Energy, Office of Science, Office of High Energy Physics, under contract DE-AC02- 06CH11357. Work at Fermi National Accelerator Laboratory, a DOE-OS, HEP User Facility managed by the Fermi Research Alliance, LLC, was supported under Contract No. DE-AC02- 07CH11359. The Melbourne authors acknowledge support from the Australian Research Council’s Discovery Projects scheme (DP210102386). The McGill authors acknowledge funding from the Natural Sciences and Engineering Research Council of Canada, Canadian Institute for Advanced research, and the Fonds de recherche du Quúbec Nature et technologies. The CU Boulder group acknowledges support from NSF AST-0956135. The Munich group acknowledges the support by the ORIGINS Cluster (funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy – EXC-2094 – 390783311), the MaxPlanck-Gesellschaft Faculty Fellowship Program, and the Ludwig-Maximilians-Universität München. JV acknowledges support from the Sloan Foundation. 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, NFS’s NOIRLab, 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, Texas A&M University, and the OzDES Membership Consortium. Based in part on observations at Cerro Tololo Inter-American Observatory at NSF’s NOIRLab (NOIRLab Prop. ID 2012B-0001; PI: J. Frieman), which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation. The DES data management system is supported by the National Science Foundation under Grant Numbers AST-1138766 and AST-1536171. The DES participants from Spanish institutions are partially supported by MICINN under grants ESP2017-89838, PGC2018-094773, PGC2018-102021, SEV-2016-0588, SEV-2016-0597, and MDM-2015-0509, some of which include ERDF funds from the European Union. IFAE is partially funded by the CERCA program of the Generalitat de Catalunya. Research leading to these results has received funding from the European Research Council under the European Union’s Seventh Framework Program (FP7/2007-2013) including ERC grant agreements 240672, 291329, and 306478. We acknowledge support from the Brazilian Instituto Nacional de Ciência e Tecnologia (INCT) do e-Universo (CNPq grant 465376/2014-2). This manuscript has been authored by Fermi Research Alliance, LLC under Contract No. DE-AC02-07CH11359 with the U.S. Department of Energy, Office of Science, Office of High Energy Physics. We gratefully acknowledge the computing resources provided on Crossover ( and/or Bebop and/or Swing and/or Blues), a high-performance computing cluster operated by the Laboratory Computing Resource Center at Argonne National Laboratory.Appendix A Parameter priors
In Table 2 we list the priors used in our analysis.
Parameter | Prior |
---|---|
MagLim | |
, , , , , | |
, , , , , | |
, , , , , | |
redMaGiC | |
, , , , | |
, , , , | |
, , , , | |
MetaCalibration | |
, , , | |
, , , |
Appendix B Adding small-scale information with nonlinear galaxy bias
Our baseline analysis adopts a linear galaxy bias model to describe the relationship between the galaxy overdensity and the underlying matter field. At small scales, this description of galaxy biasing is known to break down. The breakdown in linear galaxy bias drives our choice of angular scales used to analyzing the correlation, as described in Paper I. By adopting a higher-order bias model, it is possible to include smaller angular scales in the cosmological analysis and potentially improve parameter constraints. At the same time, a more complex bias model necessitates more free parameters, which degrades the parameter constraints to some extent. We now consider the parameter constraints from using the nonlinear galaxy bias model described in Pandey et al. (2021).
The constraints from this analysis are presented in Fig. 14. We find that adopting a nonlinear description of galaxy bias improves the precision of the constraint on by roughly 10%, and the precision of the constraints on both and by roughly 10%.


Appendix C Alternative redshift calibration and IA model
Our baseline analysis assumes that uncertainties in the source galaxy redshift distributions are characterized by shift and stretch parameters, as described in (Myles et al., 2021; Gatti et al., 2022). An alternative approach to characterizing the uncertainties in the redshift distributions is hyperrank, described in Cordero et al. (2022). Rather than attempt to parameterize biases in the redshift distributions, hyperrank provides a way to sample over realizations of the full posteriors on these distributions. Repeating our analysis of the data using this alternative redshift uncertainty prescription yields the constraints shown in Fig. 15. Although there is a small shift in , it is well within our uncertainties.
The intrinsic alignment (IA) model that we adopt in our baseline analysis is TATT (TATT, Blazek et al., 2019). In Fig. 15, we show the results of instead adopting the nonlinear alignment model (NLA, Bridle and King, 2007). The NLA model is more restrictive than TATT in the sense that the latter becomes equivalent to the former in the limit that . We find that switching to NLA results in minimal changes to the parameter constraints from .
References
- Heymans et al. (2021) C. Heymans, T. Tröster, M. Asgari et al., A&A 646, A140 (2021), arXiv:2007.15632 [astro-ph.CO] .
- Hamana et al. (2020) T. Hamana, M. Shirasaki, S. Miyazaki et al., PASJ 72, 16 (2020), arXiv:1906.06041 [astro-ph.CO] .
- DES Collaboration et al. (2018) DES Collaboration et al., Phys. Rev. D 98, 043526 (2018), arXiv:1708.01530 [astro-ph.CO] .
- DES Collaboration et al. (2022) DES Collaboration et al., Phys. Rev. D 105, 023520 (2022), arXiv:2105.13549 [astro-ph.CO] .
- Battye et al. (2015) R. A. Battye, T. Charnock and A. Moss, Phys. Rev. D 91, 103508 (2015), arXiv:1409.2769 [astro-ph.CO] .
- MacCrann et al. (2015) N. MacCrann, J. Zuntz, S. Bridle et al., MNRAS 451, 2877 (2015), arXiv:1408.4742 [astro-ph.CO] .
- Raveri (2016) M. Raveri, Phys. Rev. D 93, 043522 (2016), arXiv:1510.00688 [astro-ph.CO] .
- Raveri and Hu (2019) M. Raveri and W. Hu, Phys. Rev. D 99, 043506 (2019), arXiv:1806.04649 [astro-ph.CO] .
- Aghanim et al. (2020) N. Aghanim et al. (Planck), Astron. Astrophys. 641, A6 (2020), arXiv:1807.06209 [astro-ph.CO] .
- Krolewski et al. (2021) A. Krolewski, S. Ferraro and M. White, JCAP 2021, 028 (2021), arXiv:2105.03421 [astro-ph.CO] .
- Robertson et al. (2021) N. C. Robertson, D. Alonso, J. Harnois-Déraps et al., A&A 649, A146 (2021), arXiv:2011.11613 [astro-ph.CO] .
- Chang et al. (2022) C. Chang, Y. Omori, E. J. Baxter et al., arXiv e-prints , arXiv:2203.12440 (2022), arXiv:2203.12440 [astro-ph.CO] .
- White et al. (2022) M. White, R. Zhou, J. DeRose et al., JCAP 2022, 007 (2022), arXiv:2111.09898 [astro-ph.CO] .
- Chang et al. (2019) C. Chang, M. Wang, S. Dodelson et al., MNRAS 482, 3696 (2019), arXiv:1808.07335 [astro-ph.CO] .
- Carlstrom et al. (2011) J. E. Carlstrom, P. A. R. Ade, K. A. Aird et al., PASP 123, 568 (2011), arXiv:0907.4445 [astro-ph.IM] .
- Planck Collaboration (2011) Planck Collaboration, A&A 536, A1 (2011), arXiv:1101.2022 [astro-ph.IM] .
- Omori et al. (2017) Y. Omori, R. Chown, G. Simard et al., ApJ 849, 124 (2017), arXiv:1705.00743 .
- Omori et al. (2022) Y. Omori, E. J. Baxter, C. Chang et al., arXiv e-prints , arXiv:2203.12439 (2022), arXiv:2203.12439 [astro-ph.CO] .
- Pandey et al. (2021) S. Pandey, E. Krause, J. DeRose et al., arXiv e-prints , arXiv:2105.13545 (2021), arXiv:2105.13545 [astro-ph.CO] .
- Porredon et al. (2021) A. Porredon, M. Crocce, J. Elvin-Poole et al., arXiv e-prints , arXiv:2105.13546 (2021), arXiv:2105.13546 [astro-ph.CO] .
- Baxter et al. (2016) E. Baxter, J. Clampitt, T. Giannantonio et al., MNRAS 461, 4099 (2016), arXiv:1602.07384 [astro-ph.CO] .
- Kirk et al. (2016) D. Kirk, Y. Omori, A. Benoit-Lévy et al., MNRAS 459, 21 (2016), arXiv:1512.04535 .
- Giannantonio et al. (2016) T. Giannantonio, P. Fosalba, R. Cawthon et al., MNRAS 456, 3213 (2016), arXiv:1507.05551 .
- Baxter et al. (2018) E. J. Baxter, S. Raghunathan, T. M. Crawford et al., MNRAS 476, 2674 (2018), arXiv:1708.01360 [astro-ph.CO] .
- Omori et al. (2019a) Y. Omori, T. Giannantonio, A. Porredon et al., Phys. Rev. D 100, 043501 (2019a), arXiv:1810.02342 [astro-ph.CO] .
- Omori et al. (2019b) Y. Omori, E. J. Baxter, C. Chang et al., Phys. Rev. D 100, 043517 (2019b), arXiv:1810.02441 [astro-ph.CO] .
- Abbott et al. (2019) T. M. C. Abbott, F. B. Abdalla, A. Alarcon et al., Phys. Rev. D 100, 023541 (2019), arXiv:1810.02322 [astro-ph.CO] .
- Madhavacheril and Hill (2018) M. S. Madhavacheril and J. C. Hill, Phys. Rev. D 98, 023534 (2018).
- Flaugher (2005) B. Flaugher, International Journal of Modern Physics A 20, 3121 (2005).
- Flaugher et al. (2015) B. Flaugher, H. T. Diehl, K. Honscheid et al. (DES Collaboration), AJ 150, 150 (2015), arXiv:1504.02900 [astro-ph.IM] .
- Sevilla-Noarbe et al. (2021) I. Sevilla-Noarbe, K. Bechtol, M. Carrasco Kind et al., ApJS 254, 24 (2021), arXiv:2011.03407 [astro-ph.CO] .
- Porredon et al. (2021) A. Porredon et al. (DES), Phys. Rev. D 103, 043503 (2021), arXiv:2011.03411 [astro-ph.CO] .
- Gatti et al. (2021) M. Gatti, E. Sheldon, A. Amon et al., MNRAS 504, 4312 (2021), arXiv:2011.03408 [astro-ph.CO] .
- Krause et al. (2021) E. Krause, X. Fang, S. Pandey et al., arXiv e-prints , arXiv:2105.13548 (2021), arXiv:2105.13548 [astro-ph.CO] .
- Limber (1953) D. N. Limber, ApJ 117, 134 (1953).
- Lewis et al. (2000) A. Lewis, A. Challinor and A. Lasenby, ApJ 538, 473 (2000), arXiv:astro-ph/9911177 [astro-ph] .
- Takahashi et al. (2012) R. Takahashi, M. Sato, T. Nishimichi et al., ApJ 761, 152 (2012), arXiv:1208.2701 .
- Blazek et al. (2019) J. A. Blazek, N. MacCrann, M. A. Troxel et al., Phys. Rev. D 100, 103506 (2019), arXiv:1708.09247 [astro-ph.CO] .
- Elvin-Poole et al. (2022) J. Elvin-Poole, N. MacCrann et al., To be submitted to MNRAS (2022).
- Myles et al. (2021) J. Myles, A. Alarcon, A. Amon et al., MNRAS 505, 4249 (2021), arXiv:2012.08566 [astro-ph.CO] .
- Gatti et al. (2022) M. Gatti, G. Giannini, G. M. Bernstein et al., MNRAS 510, 1223 (2022), arXiv:2012.08569 [astro-ph.CO] .
- Cawthon et al. (2020) R. Cawthon et al. (DES), Submitted to MNRAS (2020), arXiv:2012.12826 [astro-ph.CO] .
- Cordero et al. (2022) J. P. Cordero, I. Harrison, R. P. Rollins et al., MNRAS 511, 2170 (2022), arXiv:2109.09636 [astro-ph.CO] .
- MacCrann et al. (2022) N. MacCrann, M. R. Becker, J. McCullough et al., MNRAS 509, 3371 (2022), arXiv:2012.08567 [astro-ph.CO] .
- MacCrann et al. (2020) N. MacCrann, J. Blazek, B. Jain et al., MNRAS 491, 5498 (2020), arXiv:1903.07101 [astro-ph.CO] .
- Jarvis et al. (2004) M. Jarvis, G. Bernstein and B. Jain, MNRAS 352, 338 (2004), arXiv:astro-ph/0307393 [astro-ph] .
- Rodríguez-Monroy et al. (2022) M. Rodríguez-Monroy, N. Weaverdyck, J. Elvin-Poole et al., MNRAS 511, 2665 (2022), arXiv:2105.13540 [astro-ph.CO] .
- Prat et al. (2022) J. Prat, J. Blazek, C. Sánchez et al., Phys. Rev. D 105, 083528 (2022), arXiv:2105.13541 [astro-ph.CO] .
- Amon et al. (2022) A. Amon, D. Gruen, M. A. Troxel et al., Phys. Rev. D 105, 023514 (2022), arXiv:2105.13543 [astro-ph.CO] .
- Secco et al. (2022) L. F. Secco, S. Samuroff, E. Krause et al., Phys. Rev. D 105, 023515 (2022), arXiv:2105.13544 [astro-ph.CO] .
- Bartelmann and Schneider (2001) M. Bartelmann and P. Schneider, Phys. Rep. 340, 291 (2001), astro-ph/9912508 .
- Friedrich et al. (2021) O. Friedrich, F. Andrade-Oliveira, H. Camacho et al., MNRAS 508, 3125 (2021), arXiv:2012.08568 [astro-ph.CO] .
- Lin et al. (2020) C.-H. Lin, J. Harnois-Déraps, T. Eifler et al., MNRAS 499, 2977 (2020), arXiv:1905.03779 [astro-ph.CO] .
- Zuntz et al. (2015) J. Zuntz, M. Paterno, E. Jennings et al., Astronomy and Computing 12, 45 (2015), arXiv:1409.3409 [astro-ph.CO] .
- Handley et al. (2015) W. J. Handley, M. P. Hobson and A. N. Lasenby, MNRAS 453, 4384 (2015), arXiv:1506.00171 [astro-ph.IM] .
- Sánchez et al. (2022) C. Sánchez, J. Prat, G. Zacharegkas et al., Phys. Rev. D 105, 083529 (2022), arXiv:2105.13542 [astro-ph.CO] .
- Doux et al. (2021) C. Doux, E. Baxter, P. Lemos et al., MNRAS 503, 2688 (2021), arXiv:2011.03410 [astro-ph.CO] .
- Raveri and Doux (2021) M. Raveri and C. Doux, Phys. Rev. D 104, 043504 (2021), arXiv:2105.03324 [astro-ph.CO] .
- Pandey et al. (2020) S. Pandey, E. Krause, B. Jain et al., Phys. Rev. D 102, 123522 (2020), arXiv:2008.05991 [astro-ph.CO] .
- Hirata and Seljak (2003) C. Hirata and U. Seljak, MNRAS 343, 459 (2003), arXiv:astro-ph/0301054 [astro-ph] .
- Vallinotto (2012) A. Vallinotto, ApJ 759, 32 (2012), arXiv:1110.5339 [astro-ph.CO] .
- Schaan et al. (2017) E. Schaan, E. Krause, T. Eifler et al., Phys. Rev. D 95, 123512 (2017), arXiv:1607.01761 .
- Baumann et al. (2013) D. Baumann, S. Ferraro, D. Green et al., JCAP 2013, 001 (2013), arXiv:1209.2173 [astro-ph.CO] .
- Desjacques et al. (2018) V. Desjacques, D. Jeong and F. Schmidt, Phys. Rep. 733, 1 (2018), arXiv:1611.09787 [astro-ph.CO] .
- Sobrin et al. (2022) J. A. Sobrin, A. J. Anderson, A. N. Bender et al., ApJS 258, 42 (2022), arXiv:2106.11202 [astro-ph.IM] .
- Henderson et al. (2016) S. W. Henderson, R. Allison, J. Austermann et al., Journal of Low Temperature Physics 184, 772 (2016), arXiv:1510.02809 [astro-ph.IM] .
- Ivezić et al. (2019) Ž. Ivezić, S. M. Kahn, J. A. Tyson et al., ApJ 873, 111 (2019), arXiv:0805.2366 .
- LSST Science Collaboration et al. (2009) LSST Science Collaboration, P. A. Abell, J. Allison et al., “LSST Science Book, Version 2.0,” (2009), arXiv:0912.0201 [astro-ph.IM] .
- Doré et al. (2019) O. Doré, C. Hirata, Y. Wang et al., eprint arXiv:1904.01174 (2019), 1904.01174 .
- Laureijs et al. (2011) R. Laureijs, J. Amiaux, S. Arduini et al., “Euclid definition study report,” (2011), arXiv:1110.3193 [astro-ph.CO] .
- Ade et al. (2019) P. Ade, J. Aguirre, Z. Ahmed et al., Journal of Cosmology and Astroparticle Physics 2019, 056–056 (2019).
- Abazajian et al. (2016) K. N. Abazajian, P. Adshead, Z. Ahmed et al., “Cmb-s4 science book, first edition,” (2016), arXiv:1610.02743 [astro-ph.CO] .
- Bridle and King (2007) S. Bridle and L. King, New Journal of Physics 9, 444 (2007), arXiv:0705.0166 .