Fundamental Physics from Future WeakLensing Calibrated SunyaevZel'dovich Galaxy Cluster Counts
Abstract
Future highresolution measurements of the cosmic microwave background (CMB) will produce catalogs of tens of thousands of galaxy clusters through the thermal SunyaevZel’dovich (tSZ) effect. We forecast how well different configurations of a CMB Stage4 experiment can constrain cosmological parameters, in particular the amplitude of structure as a function of redshift , the sum of neutrino masses , and the dark energy equation of state . A key element of this effort is calibrating the tSZ scaling relation by measuring the lensing signal around clusters. We examine how the mass calibration from future optical surveys like the Large Synoptic Survey (LSST) compares with a purely internal calibration using lensing of the CMB itself. We find that, due to its highredshift leverage, internal calibration gives constraints on cosmological parameters comparable to the optical calibration, and can be used as a crosscheck of systematics in the optical measurement. We also show that in contrast to the constraints using the CMB lensing power spectrum, lensingcalibrated tSZ cluster counts can detect a minimal at the 35 level even when the dark energy equation of state is freed up.
I Introduction
The abundance of galaxy clusters is a sensitive probe of the amplitude of density fluctuations that scales strongly with the normalization of the matter power spectrum, , and the matter density, (e.g., Voit, 2005; Allen et al., 2011). Measuring cluster abundance as a function of redshift allows one to probe physics that affects the growth of structure, for example the effect of massive neutrinos and the dark energy equation of state. Recent constraints from measurements of cluster abundances have however been limited by systematic effects (e.g., Vikhlinin et al., 2009; Vanderlinde et al., 2010; Sehgal et al., 2011; Benson et al., 2013; Hasselfield et al., 2013; Planck Collaboration & Ade, 2014; Mantz et al., 2014; Planck Collaboration & Ade, 2016; Mantz et al., 2015; de Haan et al., 2016), the dominant systematic uncertainty being the calibrations of observabletomass relations. Therefore, accurate and precise calibrations of the observabletomass relation is essential for any future cluster cosmological constraint.
Galaxy clusters are observationally identified across the electromagnetic spectrum, from microwaves to Xray energies. The measurements of secondary temperature anisotropies in the CMB that arise from the tSZ effect (Sunyaev & Zeldovich, 1970) are emerging as a powerful tool to find and count clusters. To compare with observational abundances the theoretical abundance predictions are typically forward modeled from cosmological parameters through a mass function (e.g., Press & Schechter, 1974; Tinker et al., 2008) and an observable–mass relation. For example, the recent Planck (Planck Collaboration & Ade, 2016) and South Pole Telescope (SPT, de Haan et al., 2016) cosmological constraints from tSZ cluster counts included weaklensing and Xray masses to calibrate their observable–mass relation, and Atacama Cosmology Telescope (ACT) used velocity dispersions (Hasselfield et al., 2013). In all tSZ cluster count analyses an independent method for measuring and calibrating cluster masses is required and lensing calibrations are expected to be the most unbiased (e.g., Becker & Kravtsov, 2011).
Planned CMB experiments like Advanced ACT, SPT3G, Simons Array, Simons Observatory (SO) and CMB Stage4 will produce catalogs of tens of thousands of galaxy clusters. Overlap with an optical survey like LSST will allow for precise measurements of the shapes of background galaxies behind these galaxy clusters, which will provide accurate mass calibrations for most of the galaxy clusters at low and intermediate redshifts. Another promising and independent way to calibrate the observable–mass relation is to use the socalled CMB halo lensing (Seljak & Zaldarriaga, 2000; Dodelson, 2004; Holder & Kosowsky, 2004; Hu et al., 2007), lensing of the CMB itself (Louis & Alonso, 2017) by dark matter halos. This technique is viable for clusters at any redshift and has completely different systematics from optical weak lensing. The lensing signal from dark matter halos has only recently been detected (Madhavacheril et al., 2015; Baxter et al., 2015; Planck Collaboration & Ade, 2016; Geach & Peacock, 2017; Baxter et al., 2017) and has already been used as a mass calibrator in an tSZ cosmological analysis (Planck Collaboration & Ade, 2016).
In this paper we forecast the constraints that CMB Stage4 can achieve on cosmological parameters using tSZ cluster counts. In the forecasts we include the external calibrations of the tSZ observable–mass relation for clusters from optical weaklensing observations using experiments like LSST and internal CMB Halo lensing calibration from CMB Stage4. The paper is structured as follows: in Section II we describe the assumptions about the experimental setup of CMB Stage4 and its variations that we compare. Section III describes our methodology for modeling cluster detection and cluster abundance. In Section IV, we present how we forecast the ability of an LSSTlike experiment to calibrate the tSZ scaling relation and in Section V, we describe how CMB lensing does the same. Section VI presents our Fisher forecasting assumptions and the cosmological models considered. We conclude with a discussion in Section VII.
Ii Experimental setup
We consider an experimental configuration consisting of a single large telescope with seven bandpasses shown in Table 1. Our baseline configuration has a white noise level of in the 150 GHz and 90 GHz channels. We vary the beam fullwidth halfmaximum (FWHM) in the 150 GHz channel from 1 arcminute to 3 arcminutes, scaling the beam FWHM in the other channels assuming . The noise sensitivities are assumed to correspond to a fraction of sky observed .
In addition to instrumental white noise, we include the effect of atmospheric noise for a groundbased experiment, parameterized through a knee multipole and tilt ,
(1) 
Atmospheric noise can have a considerable impact on the number of clusters detected and consequently on cosmological constraints. In our fiducial analysis, we will assume an of 3500 in temperature and 300 in polarization, and an of 4.5 in temperature and 3.5 in polarization in every frequency bandpass (Louis et al., 2017b). In reality, these parameters may depend on the aperture size of the telescope (amongst other experimental variables) and will vary between bandpasses. Our fiducial values are motivated by the performances of past and ongoing groundbased CMB experiments. A detailed analysis of the dependence of and on aperture size and frequency of observation is beyond the scope of this work. We do however undertake a study of the effect of and on the number of clusters detected in Section III.
As described in Section III, each assumed experimental configuration predicts a certain number of tSZ cluster detections as a function of mass, redshift and signaltonoise. For the sample of tSZ clusters selected this way, we obtain lensing mass calibration either internally using CMB lensing (see Section V), or externally from an LSSTlike optical weak lensing survey configuration (see Section IV). The internal calibration is done either on both temperature and polarization data (T+P) from the 150 GHz channel, or on polarization only (Ponly). The optical lensing calibration is done either for clusters with redshifts or for . We assume that the optical lensing survey provides brightest central galaxies (BCG) which are used as centroids for stacking the optical and CMB lensing signals. We therefore assume that miscentering of the stack and the true mass centroid can be assumed to be negligible (compared to the beam size) for clusters with . For CMB lensing mass calibration of clusters, we do not assume the availability of BCG centers and marginalize over miscentering effects as described in Section V.
The optical survey is also assumed to provide photometric redshifts for at least some member galaxies of each tSZ detected cluster. These redshifts are not required to be very precise since they are needed only for coarse binning of the clusters in redshift. Any tSZ detected cluster that cannot be associated with any possible member galaxies in the optical survey can fairly confidently be assigned to the redshift bin that is calibrated using CMB lensing.

Beam (arcminutes)  

Frequency (GHz)  Noise []  CMBS43.0  CMBS42.5  CMBS42.0  CMBS41.5  CMBS41.0 
21  7.9  21.4  17.9  14.3  10.7  7.1 
29  5.6  15.5  12.9  10.3  7.8  5.2 
40  5.4  11.2  9.4  7.5  5.6  3.8 
95  1.5  4.7  4.0  3.2  2.4  1.6 
150  1.5  3.0  2.5  2.0  1.5  1.0 
220  5.2  2.0  1.7  1.4  1.0  0.7 
270  9.0  1.7  1.4  1.1  0.8  0.6 
We emphasize that this is a strawperson experimental design for CMB Stage4, for example the frequency bands for CMB Stage4 have not yet been determined. These particular bands were chosen to cover the main atmospheric windows around the peak CMB and tSZ sensitivity with extra high and low frequency bands for potential foreground cleaning. The distribution of detectors weights among bands reflects a rough optimization for CMB and tSZ signal and assuming some level of foreground subtraction. In this work we focused on the aperture size of CMB Stage4 and a full optimization of the frequency bands, noise levels, and aperture sizes is beyond the scope of this work. In particular frequency bands optimization requires simulations that include correlated sources of noise like the analyses in Melin et al. (2017).
Iii Methodology
The thermal SZ (tSZ) signal is the observable cluster property that we model onto the theoretical predictions for the abundance of clusters. We use an analytic model for these tSZ selected clusters that accounts for measurement uncertainties in mass calibration and integrated Compton signal. The spectral distortion caused by the tSZ in the observed CMB temperature is a function of frequency
(2) 
here , where , is the Planck constant, and is the Boltzmann constant. Note that we neglected relativistic corrections to the tSZ spectral function (e.g., Nozawa et al., 2006; Chluba et al., 2012). As shown in Equation 2, the amplitude of the tSZ spectral distortion is directly proportional to , which is defined as,
(3) 
Here the physical constants , , and correspond to the speed of light, electron mass, and Thompson crosssection, respectively. The physical properties of the free electron that scatter the CMB photons are: the electron number density and is the electron temperature. Equation 3 is integral along the lineofsight, . For a given spherical pressure profile, the signal for a cluster projected on the sky is,
(4) 
Here , is the angular diameter distance to redshift , and is the 2D angular coordinate on the sky.
For the shape of , we choose the pressure profile from (Arnaud et al., 2010), which was used in the Planck cluster analysis (e.g., Planck Collaboration & Ade, 2016). The parametric form for the profile is a generalized NavarroFrenkWhite profile (Zhao, 1996),
(5) 
where the , and the parameters of the profile have the values, , , , , and . Additionally, we choose the filter scale for each cluster, , such that .
iii.1 Cluster detection
We find tSZ clusters using a matched filter technique that exploits the unique spectral distortion of the tSZ effect Herranz et al. (2002); Melin et al. (2006). We model the maps of the millimeter sky, , as:
(6) 
here, is the amplitude of the tSZ signal for a given halo, is the normalized projected profile, , and is the noise when searching for a tSZ signal. Here the noise is a function of and includes instrumental noise, atmosphere (described in Section II), primary CMB, and other secondary sources.
The estimator we use to measure is a matched filter that is designed to minimize the variance across a given set of frequency bands for an assumed profile
(7) 
Here we sum over and is an unbiased, realspace matched filter that minimizes the variance. In Fourier space this matched filter has the form,
(8) 
Here is the Fourier transform of the normalized projected profile, is the variance, and is the covariance matrix of the noise power spectrum. Note that the Fourier transform of is . The variance is defined as
(9) 
and the noise covariance matrix is defined as,
(10)  
The components of the noise covariance matrix are the CMB crosspower spectra , the secondary cross power spectra , and the debeamed noise that only contributes to the diagonal, where is the Fourier transform of the beam, which we assume to be Gaussian. The FWHM of the beams for corresponding frequency bands are shown in Table 1. The CMB secondary anisotropies that we include are, radio point sources (Poisson term), the cosmic infrared background (CIB, both Poisson and clustered terms), kinetic SunyaevZel’dovich, unresolved tSZ, and the tSZCIB crosscorrelation term (Addison et al., 2012). For the unresolved tSZ contribution, we estimate that half of the total autospectrum power is coming from clusters with masses (e.g., Komatsu & Seljak, 2002; Trac et al., 2011; Battaglia et al., 2012), that will be detected. Therefore, we removed the contribution from these clusters to the autospectrum power for the purposes of additional secondary anisotropy noise. We use the functional forms and parameters for these secondary anisotropies presented in Dunkley et al. (2013).
iii.2 tSZ cluster abundances
For simplicity the abundance of tSZ clusters is modeled as the number of clusters () observed in bins of lensing mass calibration (), tSZ signaltonoise () from the matched filter, and redshift ():
(11)  
Here is the probability distribution function of and given the integrated Comptony, , and halo mass, (see Equation 12), is the differential number of clusters with respect to and (see Equation 13), and is a lognormal distribution of given the mean integrated Comptony (, see Equation 14) and the intrinsic scatter ( see Equation 15).
We model the probability distribution function of and given and as two independent normal distributions (),
(12)  
The measurement errors, , is determined from the matched filter (see Equation 9) and measurement errors, , comes from either the optical weaklensing or CMB halolensing mass calibration (see Sections IV and V for details). The parameter is set to for the main analysis but is allowed to vary with a 1% Gaussian prior when we wish to explore the imposition of a 1% systematic floor on the mass calibration. We apply to and not since we want to impose a systematic floor that is independent of and is irreducible.
The differential numbers counts can be further deconstructed into
(13) 
Here the volume element is and for the mass function, , we use Tinker et al. (2008) while accounting for the neutrino suppression of power in after recombination, where and not (e.g., Ichiki & Takada, 2012; Costanzi et al., 2013; Castorina et al., 2014, 2015).
Following the theoretically motivated selfsimilar evolution of halos (Kaiser, 1986) we model the – scaling relation as a powerlaw that is a function of halo mass and redshift,
(14)  
Where is a constant, the pivot mass is /h, is the mass bias correction, is the first order powerlaw mass dependence, is the second order powerlaw mass dependence, and is an additional redshift dependence beyond the expected selfsimilar scaling. The functions and are the Hubble function and the angular diameter distance, respectively. The fiducial values for the scaling relation parameters in this model are . We model the scatter in this scaling relation as,
(15) 
where is the fiducial scatter and the smooth powerlaw mass and redshift dependence of the scatter are and , respectively. We choose the fiducial values for these scatter parameters to be . All fiducial parameter values and their associated step sizes for our Fisher analyses are shown in Table II.
In Figure 1 we show the number of expected clusters that would be detected as a function of redshift for various choices of aperture size for CMB Stage4. The bottom panel illustrates the ratio of clusters compared to the most pessimistic design of the CMB Stage4 experiment (3 arcminute aperture at 150 GHz) for the purposes of tSZ and secondary anisotropy science. As a function of increasing aperture the gain in detected clusters increases strongly with redshift. Between , we find increases on the order of hundreds when comparing a 1 arcminute to a 3 arcminute aperture. This redshift range of is a new frontier for clusters and protoclusters science, with only a few heterogeneously detected clusters and protoclusters within this redshift range (Overzier, 2016, and references therein). If these clusters contain hot gas, as we have assumed, CMB Stage4 will find them and produce a legacy catalog of uniformly selected, high tSZ clusters that will be ideal to study galaxy formation in high redshift, dense environments.
Atmospheric noise primarily affects large scales in the CMB, but with a sufficiently high knee multipole its effects can degrade scales relevant for cluster finding. We explore the effect of atmospheric noise in Figure 2. There is a strong dependence on , e.g., for the 1arcminute configuration corresponds to detecting 30% fewer clusters than if there were no atmospheric noise ().
Iv Optical Weak Lensing
We model the optical lensing signal as follows. For a given cluster mass , we compute concentration using the the concentrationmass relation derived in Duffy et al. (2008). We then convert to assuming the NavarroFrenkWhite (NFW, Navarro et al., 1997) profile. We assume the optical weak lensing signal at redshift is measured in terms of excess surface density (ESD);
(16)  
where is the galaxycentric transverse comoving distance, is the projected matter density profile along the lineofsight, and is the projected matter density profile averaged over distance . We employ the following halo model for ESD;
(17)  
where is a smoothlytruncated version of the NFW profile proposed in Baltz et al. (2009) with the dimensionless smoothing radius , which is converted from defined against virial radius in Oguri & Hamana (2011). We do not fit for the 2halo term , since we restrict the our analyses to the regime where dominates.
The shape noise of a given radial bin with sources at redshift is estimated as
(18) 
where is the RMS of intrinsic ellipticity
(19) 
where is the angular diameter distance of the lenssource system and comes from our use of comoving coordinates (Mandelbaum et al., 2006). The total shape noise is estimated by assuming Eq. (18) forms the inverse variance weight for given redshift of the lens;
(20) 
We compute the shape noise by asserting that there will be HSClike survey over the entire of the survey. This is not unreasonable considering LSST will be available over a large area if not all of the CMB Stage4 survey area and the parameters of HSC are conservative compared to LSST. For the HSC survey we assumed a source background 20 galaxies per square arcminute with the from Oguri & Hamana (2011);
(21) 
where that corresponds to the mean redshift .
We fit given the derived shape noise errors with an NFW profile over the radial range where the 1halo term is determined using Markov Chain Monte Carlo (MCMC) (ForemanMackey et al., 2013) assuming a fixed concentration mass relation (Duffy et al., 2008). The radial range we use for the fit is in comoving Mpc. We use the width of the inferred weak lensing mass distribution as our weak lensing mass error and take the ratio of this over the median inferred mass as the percent weak lensing mass error . Figure 3 illustrates the S/N per cluster as a function of redshift.
We caution that systematic errors are not taken into account when we forecast the errors bars on . We expect systematic uncertainties to be increasingly important as the redshift of the clusters increases due to photometric redshift and shape measurement biases (e.g., Jarvis et al., 2016; Tanaka et al., 2017). For this reason, we compare constraints from clusters with redshifts and . It is possible that at higher redshifts such systematic uncertainties will be larger than 1%, in this regime CMB halo lensing will become important.
V CMB Halo Lensing
The CMB is lensed by all structure since recombination and hence can in principle be used as a source for any galaxy cluster. In contrast to optical lensing where the sources are distributed in a wide range of uncertain redshifts behind the cluster, the CMB source plane is fixed at a relatively thin and wellmeasured slice at . In addition, the CMB is a diffuse field whose unlensed statistics are well captured by a Gaussian random field specified through a power spectrum. The effect of lensing is to to couple previously independent harmonic modes of the CMB temperature (T) and polarization fields (curlfree E and curllike B). This insight allows one to write a quadratic estimator that sums over pairs of CMB modes optimally to reconstruct the projected lensing potential at any given mode Hu et al. (2007).
The quadratic estimator requires a pair of maps, the first ‘leg’ of the pair effectively serving as a measure of the background gradient at the location of the cluster, and the second leg capturing information about the smallscale fluctuations induced by lensing by the cluster. For example, in the temperatureonly estimator combination (TT), the quadratic estimator can be written as Hu et al. (2007)
where is the lowpass filtered gradient of the temperature map as a function of angle on the sky , is the highpass filtered temperature map, is a normalization in Fourierspace that ensures this estimator is unbiased as a function of angular wavenumber , and and represent 2D Fourier and inverseFourier transforms respectively.
We use the estimator combinations TT, TE, EE, EB and TB, where the first leg in the pair is used in gradient estimation and the second leg is used for smallscale fluctuations. The noise per mode in each estimator is where the estimator normalization is given by Equation 19 in Hu et al. (2007). We consider minimum variance combinations of either all the above estimators (T+P) or ‘Polarization Only’ (P only), i.e., EE and EB. We do not account for the covariance between these estimators since in most cases either TT or EB dominates, and the covariance between TT and EB is zero. We then calculate the total signal and variance of the lensing convergence measured within a radius of when using a matched filter designed to optimally measure the lensing signal.
In this forecast analysis, we only use information for CMB lensing calibration from the CMB Stage4 150 GHz channel, but use multifrequency information from the Planck satellite experiment. Planck has effectively imaged CMB temperature modes at scales to nearly the cosmic variance limit (even if considering a foreground cleaned map like SMICA (Planck Collaboration et al., 2016) or LGMCA (Bobin et al., 2014, 2016)). The CMB halo lensing signal can be contaminated by both noise and systematic biases from cluster foregrounds such as tSZ and CIB. This contamination is to some extent mitigated by the imposition of a lowpass filter on the gradient leg. However, the systematic bias from cluster contaminants can be nearly eliminated by requiring that the gradient leg be foreground free. Since this leg only requires modes , we lose almost no signaltonoise if we use Planck foreground cleaned maps in the temperature gradient leg of the TT, TE and TB estimators while eliminating the largest source of potential bias. The rest of this analysis assumes Planck beam (5 arcminutes) and white noise of 42 arcmin for the first leg of TT, TE, TB (corresponding roughly to the level in a foreground cleaned map) and CMB Stage4 150 GHz (Table 1) in all the other legs of the estimator combinations.
Our model for the CMB lensing convergence signal is an NFW profile with mass , virial radius and concentration :
(22) 
where is angle on the sky, is the speed of light, is the Hubble constant, is the redshift of the cluster, is the scale radius , is the critical density of the Universe today, is the comoving distance to the CMB, is the comoving distance to the lens cluster, is the comoving distance between the CMB and cluster,
(23) 
and
(24) 
The relative error on the lensing mass of a given cluster is
(25) 
where is the integrated convergence within a disk of radius . The matched filter variance in the same region is given by
(26) 
where
(27) 
Here is the power spectrum of the convergence field (a lineofsight integral over the cosmological matter power spectrum including nonlinear corrections from Halofit) which captures fluctuations in the lensing field that are not related to the NFW cluster, and is the lensing reconstruction noise per mode described earlier. We note that at low instrument noise levels, the quadratic estimator is less optimal than maximum likelihood techniques Raghunathan et al. (2017) causing an underestimate of sensitivity. However, in the small lens limit, the approximation in Eq 26 that the noise modes are uncorrelated also breaks down resulting in an overestimate of sensitivity Horowitz et al. (2017). Since most clusters in the cosmological sample are not in this regime of high S/N and since the information on halo masses does not purely come from the small lens limit, we ignore these complications and leave a more complete treatment for later work.
We can now compare the performance of CMB lensing and optical lensing. Using the formalism described above, in Figure 3, we compare the signaltonoiseratio per cluster for clusters of various masses for a 1 arcminute FWHM beam experiment that utilizes both temperature and polarization for CMB lensing (and includes noise from temperature foregrounds). Since shape noise increases with redshift as fewer source galaxies become available, optical weak lensing starts becoming less statistically informative for clusters at redshifts greater than around depending on the mass of the cluster.
When stacking on clusters where we have no optical followup data (for the highestredshift clusters), assuming that the tSZ centroid is the center of the cluster can result in a smearing of the signal. To model this, we convolve the convergence profile with a Rayleigh distribution:
(28) 
where is taken to be half the FWHM of the beam.
While the resolution of the CMB instrument affects the number of clusters detected through the tSZ effect, it also affects the mass sensitivity since a higher resolution experiment images smaller scales in the CMB temperature and polarization field that contribute to the lensing signal. In Figure 4, we show the dependence of the mass sensitivity on beam FWHM for five scenarios. The most optimistic assumes that both temperature and polarization data are used and that there are no sources of noise from foregrounds (discussed in Section IIA). Foregrounds in temperature degrade the mass sensitivity by up to 80%. If one assumes there is no foreground contamination in polarization, utilizing only the EE and EB polarization based estimators results in further degradation of mass sensitivity by around 20%. While galactic foregrounds in polarization uncorrelated with the positions of galaxy clusters are highly uncertain at smallscales, polarized emission from clusters is expected to be well below the 1 level Louis et al. (2017a). We also show in Figure 4 how miscentering enhances the degradation of sensitivity as a function of beam FWHM. Our baseline forecasts assume both temperature and polarization data with foregrounds in temperature, and miscentering only for clusters with .
Parameter  Fiducial  Step  Prior 

0.1194  0.0030  
0.022  0.0008  
67.0  0.5  
2.2  0.1  
0.96  0.01  
0.06  0.02  0.01  
(meV)  60  20  
1  0.05  
0  0.1  
0.8  0.02  
1.79  0.04  
0.127  0.02  
0  0.02  
0  0.02  
0  0.02  
0  0.02  
1  0.1  0.01  
0.75  0.2  beam/2  

Vi Cosmological Constraints
By finding galaxy clusters and calibrating their masses, we are constraining the halo abundance , a function that is sensitive to the amplitude of matter fluctuations and the total matter density . The halo abundance is related to the matter power spectrum (and consequently the growth factor). Physics that affects the matter power spectrum or its growth can lead to differing predictions on the abundance of halos as a function of mass and redshift.
We obtain predicted constraints on cosmological parameters by calculating Fisher matrices for each experimental configuration (e.g., Fisher, 1935; Knox, 1995; Jungman et al., 1996). We proceed as in Shimon et al. (2011); Louis & Alonso (2017), modeling the observed number counts as Poisson distributed about the predicted mean number of clusters (See Equation 11) in each bin,
(29) 
and obtain 68% C.L. constraints by expanding assuming the likelihood is Gaussian in the parameters. This requires us to calculate derivatives of the observed number counts as a function of the parameters varied. The Fisher matrix is
(30) 
where and indicate the cosmological, scaling relation and other nuisance parameters. Marginalized 68 %C.L. constraints on parameter for instance can then be calculated as
(31) 
We vary with respect to the following set of parameters. When including primary CMB information, our parameter set is . When not including primary CMB, we exclude from the Fisher matrix. In addition, we always marginalize over a set of scaling relation parameters described earlier in Section III. The other parameters and external data sets we consider are described below.
Before these derivatives are calculated, the original grids are rebinned such that all clusters with are in a single bin and such that the bins are coarse enough (37 bins between ) given the mass calibration errors.
vi.1 External data sets
The cosmological constraints considered here will benefit from primary CMB information that pins down the amplitude of primordial power. For this purpose, we include a CMB Stage4 Fisher matrix for when varying neutrino mass and when not. For CMB Stage4, Fisher information is only included for the temperature multipole range and the polarization multipole range . In addition, we include a Planck Fisher matrix for temperature and polarization. We avoid including low polarization from Planck and avoid double counting as follows. We include Planck temperature information for , Planck temperature and polarization in the overlapping sky of and Planck temperature and polarization in the nonoverlapping sky of . We use unlensed spectra to effectively exclude information from primary CMB. In lieu of including low polarization from Planck, we impose a flat prior of 0.01 on (unless otherwise specified) whenever including primary CMB data.
Although cluster counts measure , some improvement in cosmological constraints can be obtained by the addition of baryon acoustic oscillations (BAO) measurement through the information it provides on . We consider a BAO experiment like DESI and calculate its Fisher matrix. While not critical for the science targets, we show how much the addition of DESI can improve our constraints. The Fisher formalism for CMB and BAO and the experimental configurations for Planck and DESI follow those used in (Allison et al., 2015).
vi.2 Additional nuisance parameters
We explore the effect of imposing a 1% floor on the systematic uncertainty in the mass calibration since in the case of optical weak lensing a combination of shear multiplicative bias, photoz uncertainties and modeling uncertainties is expected to lead to an overall floor at that level. We impose this floor by rescaling the lensing mass as
(32) 
marginalizing over but with a prior of 0.01. We find that this has very little effect on parameter constraints since the 37 mass bins used in the analysis allow for some selfcalibration (Lima & Hu, 2005) through information in the shape of the mass function. We do not vary in the main results in this work.
For CMB lensing, we marginalize over a miscentering offset discussed in SectionV for clusters with and impose a prior on it that is equal to half the beam FWHM.
vi.3 Amplitude of the matter power spectrum
In order to project the sensitivity to the amplitude of matter fluctuations as a function of redshift, we calculate the following derivatives,
(33)  
where is the step size for the derivative calculation which we take to be (with other values tested to confirm stability). When calculating these derivatives, we fix to their fiducial values varying only the power in each redshift bin that is used in the calculation of halo abundances. These derivatives are then stitched into a Fisher matrix for parameters and a CMB Fisher matrix for is added to it. When reporting constraints on the overall amplitude of power , we involve a single derivative that varies power across all redshifts.
In Figure 5, we show constraints on in redshift bins for the various telescope resolutions in Table 1, where we assume internal CMB lensing calibration with T+P. While the amplitude of the tSZ effect is roughly constant as a function of redshift, a higher resolution experiment is able to find tSZ clusters at higher redshifts that subtend smaller angles on the sky. These clusters are also calibrated better with CMB lensing as the resolution improves as indicated in Figures 3 and 4. The increased resolution thus primarily improves constraints at higher redshifts. This improvement begins to saturate below 1.5 arcminutes due to degeneracies with other parameters (see VII).
vi.4 The sum of neutrino masses
Measurements of neutrino oscillations indicate that neutrinos are massive, but the absolute mass scale (sum of the masses of three neutrino species) is not known. Solar and atmospheric neutrino measurements allow for either a normal or inverted hierarchy of the three species with the minimal possible mass scales in each case being 58 meV and 100 meV respectively. Any cosmological observable sensitive to in principle offers information on the neutrino mass scale since massive neutrinos become nonrelativistic around and affect the growth of structure. In particular, the power spectrum is suppressed on scales smaller than the neutrino freestreaming scale.
The CMB lensing power spectrum (which depends on the integrated lineofsight matter power) is likely the cleanest cosmological probe of since for CMB Stage 4 lensing constraints will be driven by polarization data with potentially fewer astrophysical systematics. However, because of a degeneracy with , CMB lensing will need to be combined with external data from baryon acoustic oscillation (BAO) surveys like DESI in order to approach sensitivities capable of providing evidence for a minimal neutrino mass. the tSZ cluster abundances however are also highly sensitive to the neutrino mass scale since the neutrino free streaming scale below which matter fluctuations are suppressed is larger than the typical virial radii of even the most massive galaxy clusters. A cosmology with a large neutrino mass scale will therefore have fewer clusters. Importantly, cluster abundances also measure , and therefore provide better constraints than CMB lensing alone in the absence of external data.
We allow for a nonzero sum of neutrino masses through its effect on the growth of structure (captured in the matter power spectrum ) and consequently on the number density of detected clusters.
vi.5 The dark energy equation of state
We consider the dark energy equation of state parametrized as
(34) 
and forecast constraints either on (with ) or on the combinations , and .
Measurements of the CMB lensing power spectrum are sensitive to the integrated matter power spectrum and lack redshift resolution. Since the equation of state of dark energy affects the growth of structure at low redshifts () and massive neutrinos suppress power below the freestreaming much earlier on, CMB lensing power spectra measurements suffer a degeneracy between and (Allison et al., 2015). Counting clusters in redshift bins significantly alleviates this degeneracy. In particular, as can be seen in Figure 6 higher resolution telescopes that find more clusters at higher redshifts suffer less degeneracy between and .
In Figure 7, we look at constraints in either a CDM+{,} cosmology (left panel) or a CDM+{,,} cosmology (right panel). In the first case, allowing only clusters when optical weak lensing is used as a mass calibrator (anticipating photometric redshift systematics for source galaxies ) we find that CMB lensing performs comparably or better regardless of whether we restrict ourselves to polarization data only. When is freed up, optical weak lensing performs better relative to CMB lensing. Figure 8 looks at the marginalized constraint on the sum of neutrino masses alone, in cosmologies with fixed, varied with fixed, and both and varied. We find that in contrast to the CMB lensing power spectrum, constraints on neutrino mass are not significantly degraded when the equation of state of dark energy is freed up.
Vii Discussion and Conclusions
Future measurements of CMB secondary anisotropies from experiments like CMB Stage4 will find a very large number of galaxy clusters, , with the precise number depending on the resolution of the experiment. With such large statistics, many of the scaling relation parameters and the mass calibration selfcalibrate. One would expect that the cosmological information would increase roughly with the square root of the number of clusters detected, especially since most of the clusters detected with higher resolution are at high redshifts. However, certain parameter degeneracies prevent there from being significant improvement in constraints beyond a resolution of around 1.5 arcminutes, particularly for the sum of neutrino masses.
Our fiducial forecasts assume knowledge of the optical depth to a precision of 0.01 (realized by the most recent Planck analysis of lowell polarization Planck Collaboration (2016)). Since a measurement of lowredshift amplitude of structure requires the amplitude of scalar fluctuations as a reference for untangling the effects of neutrino mass or dark energy, a degeneracy with the optical depth is introduced when using primary CMB data. This degeneracy is particularly limiting for , and prevents increased cluster counts from higher resolution CMB measurements from commensurately improving constraints. In Figure 8, we show the improvement obtained if the constraint on optical depth were 0.006 (the Planck Blue Book value) and the cosmicvariancelimited value of 0.002. Such improvements would allow the minimal neutrino mass in the normal hierarchy to be detected at the 5sigma level or greater using CMB Stage4 clusters + primary CMB data along, even when freeing up dark energy equation of state parameters, allowing for critical crosschecks of measurements made through other probes such as the CMB lensing power spectrum. The optical depth could be measured to the cosmic variance limit (CVL) by a future spacebased CMB experiment or by a groundbased CMB experiment if significant advances are made in controlling low polarization systematics. It could also be potentially improved beyond CVL with a 21 cm experiment (Liu et al., 2016).
Our knowledge of other parameters such as the Hubble constant is also a limiting factor in these constraints. We explore the contributions of all parameter degeneracies in Figure 9 by starting with our fiducial constraint and cumulatively improving the prior knowledge on every other parameter along the xaxis, switching to the next parameter once the previous parameter saturates the constraint. Perfect knowledge of and significantly improves constraints for and as expected.
Previous forecasts and analyses of cosmological constraints from CMB Stage4 SZ clusters have come to same conclusions as this work with regards to the potential cosmological utility of SZ selected clusters. The constraints on shown in this work, Louis & Alonso (2017), and Melin et al. (2017) are competitive with other largescale structure probes, for example CMB lensing (Abazajian & CMBS4 collaboration, 2016, and references therein). These constraints have the advantage that they do not require an additional constraint on . This work goes beyond previous work (e.g., Louis & Alonso, 2017; Melin et al., 2017) in a few ways. First, we include the atmosphere in the noise modeling, which affects the number of clusters one expects to find (see Figure 2). Second, we include additional degrees of freedom in the Y–M scaling relation modeling, such as mass and redshift dependent intrinsic scatter, and we explore a systematic uncertainty floor of 1% on the lensing mass calibration. Thirdly, we include the unresolved tSZ source background, which is an irreducible source of noise for the match filter, since it has the same frequency dependence as the SZ clusters we are looking for. We consider the optical weaklensing mass calibration from a conservative LSSTlike survey over the entire CMB Stage4 survey area, which will already be completed before CMB Stage4 finishes its surveying. Finally, our cosmology constraints also extend to a modelindependent forecast on the amplitude of matter fluctuations which demonstrates the effectiveness of larger telescopes in distinguishing between differing predictions of the growth of structure at high redshifts.
In this work we have been very conservative in how we implement the observable–mass relation, scatter, and the mass calibrations, however, if tSZ cluster counts are going to drive the design of CMB Stage4 then moving beyond Fisher forecasts is imperative. Important steps forward include the incorporation of correlations between the weaklensing and tSZ mass proxies. As expected these correlations are present (Stanek et al., 2010; White et al., 2010; Shirasaki et al., 2016), and to properly quantify these correlation requires hydrodynamic cosmological simulations. Additionally, correlated foregrounds, in particular the tSZCIB correlation will need to be properly addressed with future analyses of cosmological simulations, similar to the analyses in Melin et al. (2017). A realistic treatment of internal CMB halo lensing requires the full utilization of mutlifrequency information, explicitly projecting out the foregrounds that cause the most bias, a full accounting of the impact of kinetic SZ that cannot be removed through component separation and a comparison with maximum likelihood techniques that takes into account the presence of foregrounds. Finally, we do not consider any uncertainties on the halo mass function that could be as high as 10% when considering the baryonic effects (e.g., Cui et al., 2012, 2014; Bocquet et al., 2016). Fortunately, these baryonic effects can be separated from the impact of on the halo mass function (Mummery et al., 2017). Regardless of these proposed forecasts advancements that will require more sophisticated treatments like simulations, the trends with aperture size and the overall constraints, particularly on , will not be affected substantially, since they are limited by external parameters like .
Acknowledgements.
We thank Steve Allen, David Alonso, Jim Bartlett, Brad Benson, Tom Crawford, Jo Dunkley, Colin Hill, Arthur Kosowsky, Thibaut Louis, Neelima Sehgal and David Spergel for their productive discussions and comments on this work. NB acknowledges the support from the Lyman Spitzer Jr. Fellowship. The Flatiron Institute is supported by the Simons Foundation. HM is supported by the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration.Footnotes
 for a theoretical experiment represents the central frequency of a given frequency band
 In this equation, ellipticity is defined in terms of shear, i.e., , where and is the major and minor axis, respectively
References
 Abazajian, K. N., & CMBS4 collaboration. 2016, ArXiv eprints 1610.02743
 Addison, G. E., Dunkley, J., & Spergel, D. N. 2012, Mon. Not. Roy. Astron. Soc. , 427, 1741
 Allen, S. W., Evrard, A. E., & Mantz, A. B. 2011, ARAA, 49, 409
 Allison, R., Caucal, P., Calabrese, E., Dunkley, J., & Louis, T. 2015, Phys. Rev. D , 92, 123535
 Arnaud, M., Pratt, G. W., Piffaretti, R., Böhringer, H., Croston, J. H., & Pointecouteau, E. 2010, Astron. Astrophys. , 517, A92
 Baltz, E. A., Marshall, P., & Oguri, M. 2009, JCAP, 1, 015
 Battaglia, N., Bond, J. R., Pfrommer, C., & Sievers, J. L. 2012, Astrophys. J. , 758, 75
 Baxter, E. J. et al. 2015, Astrophys. J. , 806, 247
 Baxter, E. J., Raghunathan, S., SPT Collaboration, & DES Collaboration. 2017, ArXiv eprints 1708.01360
 Becker, M. R., & Kravtsov, A. V. 2011, Astrophys. J. , 740, 25
 Benson, B. A. et al. 2013, Astrophys. J. , 763, 147
 Bobin, J., Sureau, F., & Starck, J.L. 2016, Astron. Astrophys. , 591, A50
 Bobin, J., Sureau, F., Starck, J.L., Rassat, A., & Paykari, P. 2014, Astron. Astrophys. , 563, A105
 Bocquet, S., Saro, A., Dolag, K., & Mohr, J. J. 2016, Mon. Not. Roy. Astron. Soc. , 456, 2361
 Castorina, E., Carbone, C., Bel, J., Sefusatti, E., & Dolag, K. 2015, JCAP, 7, 043
 Castorina, E., Sefusatti, E., Sheth, R. K., VillaescusaNavarro, F., & Viel, M. 2014, JCAP, 2, 049
 Chluba, J., Nagai, D., Sazonov, S., & Nelson, K. 2012, Mon. Not. Roy. Astron. Soc. , 426, 510
 Costanzi, M., VillaescusaNavarro, F., Viel, M., Xia, J.Q., Borgani, S., Castorina, E., & Sefusatti, E. 2013, JCAP, 12, 012
 Cui, W., Borgani, S., Dolag, K., Murante, G., & Tornatore, L. 2012, Mon. Not. Roy. Astron. Soc. , 423, 2279
 Cui, W., Borgani, S., & Murante, G. 2014, Mon. Not. Roy. Astron. Soc. , 441, 1769
 de Haan, T. et al. 2016, Astrophys. J. , 832, 95
 Dodelson, S. 2004, Phys. Rev. D , 70, 023009
 Duffy, A. R., Schaye, J., Kay, S. T., & Dalla Vecchia, C. 2008, Mon. Not. Roy. Astron. Soc. , 390, L64
 Dunkley, J. et al. 2013, JCAP, 7, 025
 Fisher, R. A. 1935, Journal of the Royal Statistical Society, 98, pp. 39
 ForemanMackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306
 Geach, J. E., & Peacock, J. A. 2017, Nature Astronomy, 1, 795
 Hasselfield, M. et al. 2013, JCAP, 7, 8
 Herranz, D., Sanz, J. L., Hobson, M. P., Barreiro, R. B., Diego, J. M., MartínezGonzález, E., & Lasenby, A. N. 2002, Mon. Not. Roy. Astron. Soc. , 336, 1057
 Holder, G., & Kosowsky, A. 2004, Astrophys. J. , 616, 8
 Horowitz, B., Ferraro, S., & Sherwin, B. D. 2017, ArXiv eprints 1710.10236
 Hu, W., DeDeo, S., & Vale, C. 2007, New Journal of Physics, 9, 441
 Ichiki, K., & Takada, M. 2012, Phys. Rev. D , 85, 063521
 Jarvis, M. et al. 2016, Mon. Not. Roy. Astron. Soc. , 460, 2245
 Jungman, G., Kamionkowski, M., Kosowsky, A., & Spergel, D. N. 1996, Phys. Rev. D , 54, 1332
 Kaiser, N. 1986, Mon. Not. Roy. Astron. Soc. , 222, 323
 Knox, L. 1995, Phys. Rev. D , 52, 4307
 Komatsu, E., & Seljak, U. 2002, Mon. Not. Roy. Astron. Soc. , 336, 1256
 Lima, M., & Hu, W. 2005, Phys. Rev. D , 72, 043006
 Liu, A., Pritchard, J. R., Allison, R., Parsons, A. R., Seljak, U., & Sherwin, B. D. 2016, Phys. Rev. D , 93, 043013
 Louis, T., & Alonso, D. 2017, Phys. Rev. D , 95, 043517
 Louis, T., Bunn, E. F., Wandelt, B., & Silk, J. 2017a, ArXiv 1707.04102
 Louis, T. et al. 2017b, JCAP, 6, 031
 Madhavacheril, M. et al. 2015, Physical Review Letters, 114, 151302
 Mandelbaum, R., Seljak, U., Cool, R. J., Blanton, M., Hirata, C. M., & Brinkmann, J. 2006, Mon. Not. Roy. Astron. Soc. , 372, 758
 Mantz, A. B., Allen, S. W., Morris, R. G., Rapetti, D. A., Applegate, D. E., Kelly, P. L., von der Linden, A., & Schmidt, R. W. 2014, Mon. Not. Roy. Astron. Soc. , 440, 2077
 Mantz, A. B. et al. 2015, Mon. Not. Roy. Astron. Soc. , 446, 2205
 Melin, J.B., Bartlett, J. G., & Delabrouille, J. 2006, Astron. Astrophys. , 459, 341
 Melin, J.B. et al. 2017, ArXiv 1703.10456
 Mummery, B. O., McCarthy, I. G., Bird, S., & Schaye, J. 2017, Mon. Not. Roy. Astron. Soc. , 471, 227
 Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, Astrophys. J. , 490, 493
 Nozawa, S., Itoh, N., Suda, Y., & Ohhata, Y. 2006, Nuovo Cimento B Series, 121, 487
 Oguri, M., & Hamana, T. 2011, Mon. Not. Roy. Astron. Soc. , 414, 1851
 Overzier, R. A. 2016, Astronomy and Astrophysics Reviews, 24, 14
 Planck Collaboration. 2005, ESA Planck blue book, ESA
 —. 2016, Astron. Astrophys. , 596, A107
 Planck Collaboration et al. 2016, Astron. Astrophys. , 594, A9
 Planck Collaboration, & Ade, P. A. R. 2014, Astron. Astrophys. , 571, A20
 —. 2016, Astron. Astrophys. , 594, A24
 Press, W. H., & Schechter, P. 1974, Astrophys. J. , 187, 425
 Raghunathan, S. et al. 2017, ArXiv 1705.00411
 Sehgal, N. et al. 2011, Astrophys. J. , 732, 44
 Seljak, U., & Zaldarriaga, M. 2000, Astrophys. J. , 538, 57
 Shimon, M., Sadeh, S., & Rephaeli, Y. 2011, Mon. Not. Roy. Astron. Soc. , 412, 1895
 Shirasaki, M., Nagai, D., & Lau, E. T. 2016, Mon. Not. Roy. Astron. Soc. , 460, 3913
 Stanek, R., Rasia, E., Evrard, A. E., Pearce, F., & Gazzola, L. 2010, Astrophys. J. , 715, 1508
 Sunyaev, R. A., & Zeldovich, Y. B. 1970, Comments on Astrophysics and Space Physics, 7, 3
 Tanaka, M. et al. 2017, ArXiv:1704.05988
 Tinker, J., Kravtsov, A. V., Klypin, A., Abazajian, K., Warren, M., Yepes, G., Gottlöber, S., & Holz, D. E. 2008, Astrophys. J. , 688, 709
 Trac, H., Bode, P., & Ostriker, J. P. 2011, Astrophys. J. , 727, 94
 Vanderlinde, K. et al. 2010, Astrophys. J. , 722, 1180
 Vikhlinin, A. et al. 2009, Astrophys. J. , 692, 1060
 Voit, G. M. 2005, Reviews of Modern Physics, 77, 207
 White, M., Cohn, J. D., & Smit, R. 2010, Mon. Not. Roy. Astron. Soc. , 408, 1818
 Zhao, H. 1996, Mon. Not. Roy. Astron. Soc. , 278, 488