Getting leverage on inflation with a large photometric redshift survey
Abstract
We assess the potential of a future largevolume photometric redshift survey to constrain observational inflationary parameters using three largescale structure observables: the angular shear and galaxy power spectra, and the cluster mass function measured through weak lensing. When used in combination with Plancklike CMB measurements, we find that the spectral index can be constrained to a precision of up to . The sensitivity to the running of the spectral index can potentially improve to , roughly a factor of five better than the present constraint from Planck and auxiliary CMB data, allowing us to test the assumptions of the slowroll scenario with unprecedented accuracy. Interestingly, neither CMB+shear nor CMB+galaxy nor CMB+clusters alone can achieve this level of sensitivity; it is the combined power of all three probes that conspires to break the different parameter degeneracies inherent in each type of observations. We make our forecast software publicly available via download or upon request from the authors.
a]Tobias Basse b]Jan Hamann a,c]Steen Hannestad d]Yvonne Y. Y. Wong
[a]Department of Physics and Astronomy, University of Aarhus, DK8000 Aarhus C, Denmark \affiliation[b]Theory Division, Physics Department, CERN, CH1211 Geneva 23, Switzerland \affiliation[c]Aarhus Institute of Advanced Studies, Aarhus University, DK8000 Aarhus C, Denmark \affiliation[d]School of Physics, The University of New South Wales, Sydney NSW 2052, Australia
basse@phys.au.dk, jan.hamann@cern.ch, sth@phys.au.dk, yvonne.y.wong@unsw.edu.au
1 Introduction
Cosmic inflation is currently the most compelling mechanism to explain the generation of the primordial fluctuations. Precision measurements of anisotropies in the cosmic microwave background (CMB) by the Planck [1] and WMAP [2] missions are consistent with adiabatic, Gaussian primordial scalar perturbations that are nearly scaleinvariant [3], as predicted by the simplest slowroll models of inflation. This picture is likewise compatible with the claimed detection of primordial gravitational waves via the mode signal by the BICEP2 experiment [4]. Interestingly, the data do show a weak preference for a negative running of the scalar spectral index, at [5], which if true would imply the presence of nonnegligible third and higherorder derivatives of the inflaton potential [6]. These are not a generic inflationary prediction, and, if ever established at high statistical significance, would rule out a large number of simple candidate models.
To further the constraints on and/or the detection of , it would be extremely helpful to increase the lever arm in wavenumbers by measuring cosmic perturbations on even smaller length scales. Unfortunately, the CMB’s sensitivity here is limited by a combination of the primary signal’s suppression due to Silk damping, and the increasing importance of foreground anisotropies. Density perturbations of the matter field at later times are however not subject to these limitations, and with several largevolume photometric redshift surveys such as the Large Synoptic Survey Telescope (LSST) [7] and the Euclid mission [8] mapping out a significant fraction of the local universe in the coming decade, it is timely to address the question of to what extent the information about smaller scale perturbations gained by these surveys will further our understanding of the physics governing the inflationary era.
To this end, we consider the three most important cosmological observables of a future largevolume photometric redshift surveys—the angular shear power spectrum, the angular galaxy power spectrum, and the cluster mass function—and their sensitivity to inflationary parameters.
We generate mock data for the said cosmological observables given the survey specifications, and then analyse the mock data using Markov Chain Monte Carlo (MCMC) methods. Our main analysis tool is based on the publicly available CosmoMC [9] code,
with modifications first described in [10] and later extended in [11].
Besides reporting our forecast of the expected sensitivity of future photometric redshift surveys to inflationary parameters,
we would also like to take this opportunity to publicly release the code for general use.
The paper is organised as follows. In section 2 we briefly describe the parameterisation of primordial fluctuations used in the analysis. Section 3 reviews the cosmological observables and recapitulates the technical details of the mock data generation from [10, 11], while section 4 describes the forecast model parameter space and the likelihood functions. Our results are presented in section 5, followed by more indepth discussions pertaining to inflationary physics in section 6. Section 7 contains our conclusions.
2 Parametrisation of the primordial power spectra
We parameterise the primordial scalar and tensor power spectra in the usual way. Starting with
(1) 
we Taylorexpand to second order around the pivot scale , so that
(2) 
The choice of pivot scale is arbitrary. In this work, we use Mpc following common convention.
It is also common to express the tensor fluctuation amplitude in terms of the tensortoscalar ratio , defined as the ratio of tensor to scalar power at the pivot scale,
(3) 
Then, using the inflationary slowroll consistency relation , and neglecting any running that might be present in , we are left with four parameters
(4) 
to describe the primordial scalar and tensor fluctuations.
3 Cosmological observables and mock data generation
We consider four cosmological observables, CMB (temperature and polarisation), the angular cosmic shear power spectra, the angular galaxy power spectra, and the cluster mass function. We briefly review below these observables and the generation of synthetic data thereof, but refer the reader to more detailed descriptions in [10] (CMB, shear, and galaxy) and [11] (clusters).
3.1 Shear and galaxy power spectra
Surveys such as the LSST and Euclid are dedicated cosmic shear surveys that will measure the shapes and the positions of millions of galaxies. The cosmic shear power spectrum will be a prime product of these measurements, but the galaxy power spectrum and its crosscorrelation with cosmic shear will also come for free. Because photometric surveys are relatively poor at determining the radial positions of galaxies, it is more advantageous to work with angular power spectra , instead of the full threedimensional power spectra . (See, e.g., [12] for a parameter sensitivity forecast using the threedimensional matter power spectrum.) However, even in this case it is possible to group the observed galaxies into broad redshift bins. Thus, the final observables are , where refer to either shear (s) or galaxies (g), and to the redshift bin number of and respectively.
Once the cosmological model has been specified, the angular power spectra can be computed from theory using
(5) 
Here, is the dimensionless power spectrum of the primordial curvature perturbations , and the source functions for shear and galaxies read, respectively,
(6) 
where the integration variable is the comoving distance, is the conformal time today, and are the spherical Bessel functions.
The source functions incorporate different transfer functions of the metric perturbations and the matter density fluctuations , respectively [10]. The window functions likewise differ, and are given by
(7) 
with
(8) 
where is a surveyspecific source galaxy distribution taken here to be of the form
(9) 
and normalised to the surface density of source galaxies in bin via
(10) 
Distribution (9) is a universal form for magnitude limited surveys [13, 14, 15]. Here we use and , so as to reproduce Euclid’s projected median survey redshift of .
Finally, the galaxy window function also depends on the galaxy bias relating the galaxy number density fluctuations to the underlying matter density perturbations.
Measurement errors
Even if observations were perfect, the construction of the galaxy and shear power spectra still suffer from a sample (cosmic) variance of . This irreducible uncertainty will be accounted for at the level of the likelihood function (section 4.2). Measurementrelated uncertainties, however, must be factored into the mock data generation. We discuss these briefly below.
Shot noise The source density of galaxies determines the smallest angular separation (and thus the highest ) on the sky at which can be reliably measured. This finite surface density appears as a shotnoise term directly in the power spectrum
(11) 
which should be added to the primary power spectrum to obtain t the total observed power. Here, for the galaxy power spectrum measurement, corresponds to the rootmeansquare ellipticity of the source galaxies ( for Euclid [11]), and the Kronecker deltas and signify that the noise term contributes only to the autospectra. The source galaxy surface densities in the redshift bin can be computed using equation (10).
Photometric redshift uncertainties We model the photometric redshift uncertainty as a simple Gaussian error of standard deviation [8] and with no bias. (We note that the Dark Energy Survey has already achieved reliably and with very little bias [16].) This error can be folded into the window functions as described in reference [17]. However, we stress that, using no more than about ten redshift bins, the bin widths are typically much larger than ; increasing the number of redshift bins to beyond ten (and consequently decreasing the bin widths) also does not substantially improve the parameter sensitivies [10]. Therefore, even if finally could only be controlled to , it would have little impact on our forecasted parameter sensitivities.
In principle, redshift space distortions should also be taken into account. However, because these are of order , corresponding to , they will be subdominant to the effect of the photometric redshift error at all but the smallest redshifts. We therefore ignore them in our analysis. Likewise we ignore possible errors in the shear measurements arising from intrinsic alignment and the shape measurement of source galaxies.
Effective range of the shear and galaxy data
As we see from equation (5), the angular power spectrum is formed by integrating over two source functions and the primordial curvature perturbation power spectrum . The role of the source functions, therefore, is to pick out those values of that are relevant for the multipole in the redshift bins concerned.
Figure 1 shows the source functions for shear and galaxy clustering at the highest accessible multipole in three redshift bins. Clearly, in every redshift bin, the shear source function consistently picks out higher values than the galaxy source function, even though the highest accessible in the latter observable generally exceeds the former case. This is because shear is a cumulative lineofsight effect; while the galaxy source function selects predominantly , where is, say, the mean redshift of the bin, for the same the shear source function also finds contributions at larger values from all lower redshifts where .
We therefore expect the shear angular power spectrum to have more constraining power for simply because of the longer lever arm. Later we shall see that this is indeed the case.
3.2 Cluster mass function
As our cluster observable we use the number of clusters in the redshift bin and the mass bin (depending on ), defined as,
(12) 
where is the solid angle covered by the survey, the comoving volume element at redshift , is the window function defining the redshift (index ) and redshiftdependent mass bin (index ), and
(13) 
is the Sheth–Tormen mass function [18]. Here, , , and , while
(14) 
where is the linear threshold density established from the spherical collapse model (see [11] for details), and
(15) 
is the variance of the linear matter density field smoothed with a realspace tophat filter on a length scale corresponding to the cluster mass . The function is the Fourier transform of the realspace tophat filter, and is the linear dimensionless matter power spectrum.
Measurement errors
Photometric cluster observations are subject in principle to errors in both the determination of the cluster redshift and the cluster mass. However, spectroscopic followup on the detected clusters is expected to reduce the redshift uncertainty to a negligible fraction of the redshift bin widths employed in this study. We therefore treat the redshift as infinitely welldetermined.
Weak lensing mass determination A Euclidlike survey will detect clusters through their lensing signal. The minimum cluster mass that can be detected, or the mass detection threshold , is a function of redshift, and can be estimated following the procedure of [19, 20] outlined in [11].
In addition, the weak lensing inferred mass of a cluster is subject to scatter and bias with respect to the true cluster mass . We assume that the bias can be controlled to the required level of accuracy and model the scatter in the observed mass with a lognormal distribution according to [21, 22],
(16) 
where is the probability of inferring mass under the condition that the true cluster mass is , and so that the mean of the distribution corresponds to the true mass . The mass scatter is determined by body simulations [21]. Integrating this distribution over in the observed mass interval of the th mass bin and accounting for the mass detection threshold and redshift binning, the window function for the cluster survey can be constructed as
(17) 
where is the Heaviside step function.
Completeness and efficiency Cluster survey also suffer from incomplete detection of clusters and false cluster detections. The errors associated with these effects are characterised respectively by the completeness , defined as the fraction of clusters actually detected as peaks by the cluster finding algorithm, and the efficiency , defined as the fraction of detected peaks that correspond to real clusters. In general these quantities can be established precisely only with the help of mock cluster catalogues generated from body simulations (see, e.g., [23, 24, 25, 26, 27, 20]). Here, we follow the simplistic approach of [28], and assume both and to be mass and redshiftindependent. Then, both and can be folded into the likelihood function in a simple way (see section 4.2).
Effective range of the cluster data
As in section 3.1.2, we can ask here what is the effective range probed by the cluster mass function (13). At first glance does not appear to depend on . However, upon closer inspection, we see that the variance appearing in the variable is in fact constructed from integrating the linear dimensionless matter power spectrum filtered by a function in equation (15) for each value of . The combination , then, is the counterpart of the source functions we examined previously in section 3.1.2.
Figure 2 shows at the cluster mass detection threshold in three redshift bins. Defined this way, we see that the maximum reach of a cluster survey is comparable to that of galaxy clustering.
3.3 CMB temperature and polarisation power spectra
In order to realistically assess future parameter sensitivities, we must consider CMB anisotropy measurements in addition to the largescale structure observables discussed above. Indeed, some cosmological parameters, notably the tensortoscalar ratio , are accessible almost exclusively to CMB measurements alone.
Information from the CMB can, to some extent, be introduced into a parameter sensitivity forecast in the form of external priors on the relevant cosmological parameters. However, we caution against this shortcut for the following reasons. Firstly, such an approach can only make sense if the external priors have been derived under the same model assumptions as the cosmological model we wish to analyse. Such priors may be easy to find in the literature for simple, standard cosmological models, but are rarely available for more complicated parameter spaces. Secondly, even if the first condition can be satisfied, the marginalisation and interval construction procedures employed in the construction of a conventional onedimensional Gaussian prior remove important information about the shape of the likelihood function and the correlations between parameters. Highly nonGaussian likelihoods and/or significant correlations between parameters can strongly impact on the outcome of a parameter sensitivity forecast.
The fundamental observable we work with are the CMB power spectra , where , representing, respectively, temperature, polarisation, and polarisation. Parity conservation implies . We will therefore only work with the three autocorrelation power spectra and , and the temperature–polarisation crosscorrelation power spectrum. The use of synthetic CMB data for parameter sensitivity forecast has been discussed extensively in [29]. We follow the same approach in this work.
3.4 Synthetic data sets
Having gone through the different types of data we briefly summarise how we construct mock data for each type of observable as well as their possible crosscorrelations.

The cosmic shear autospectrum is , with , independently of redshift. As discussed in [10], shotnoise always dominates well before reaches , so that the value of the actual cutoff is immaterial. The indices label the redshift bin, where the redshift slicing is such that all bins contain similar numbers of source galaxies and so suffer the same amount of shot noise. In [10] it was found that, for most cosmological parameters, slicing the shear data beyond two redshift bins offers no significant improvement on the sensitivities. We shall test this assumption for the inflationary parameters studied here.

The galaxy autospectrum uses multipole moments running from to in redshift bins , where is set following the prescription of [10] so as to keep the nonlinear contamination to a minimum. It was found in [10] that most cosmological parameter sensitivities saturate at around . Again, we shall test whether or not this limit holds true also for inflationary parameters.
Regarding the galaxy bias, a realistic bias model will likely be a continuous and monotonic function of redshift, probably described by two or three parameters, calibrated, within errors, against simulations. In the absence of a concrete model, however, we opt to work in the optimistic extreme and assume the galaxy bias in each redshift bin to be perfectly known. The opposite extreme would be to assume the bias to be unknown and uncorrelated between redshift bins. However, this equally unrealistic assumption would lead to parameter constraints that are completely uncompetitive if the galaxy data were used on their own, and, if analysed in combination with shear data, hardly better than the shearonly constraints [10]. We therefore do not consider the case of an unknown galaxy bias.

The shear–galaxy crossspectrum in the shear redshift bin and galaxy redshift bin runs from to , where is determined by the galaxy redshift binning.

The cluster data set contains clusters in the redshift range and with masses in the range , where denotes the redshiftdependent mass detection threshold, and is defined as the redshift at which exceeds . Following [11], the redshift bins are divided such that for the fiducial model all redshift bins contain equal numbers of clusters. Clusters in each redshift bin are further subdivided into mass bins , where, again, the bin boundaries are chosen so that the same number of clusters falls into each bin [11]. We shall determine in this work how many bins are required in order to extract the maximum possible information from the cluster data.

A mock data set from a Plancklike CMB measurement is generated according to the procedure of [29]. This simplistic framework obviously does not allow us to reproduce the full complexity of the published Planck temperature likelihood [30]. Nonetheless, it can be adapted so as to capture some of the gross features of the Planck likelihood, such as the fact that its effective sky coverage is scaledependent. We outline below our adaptation and justification.
The Planck likelihood consists of a lowmultipole part () using information from 91% of the sky (but includes additional component separation uncertainties), and a highmultipole part, which has a sky coverage of 58% in the 100 GHz and 37% in the 143 GHz and 217 GHz channels. Up to around , the 100 GHz channel is limited only by cosmic variance and therefore dominates the information. On smaller scales, the 100 GHz channel becomes noisedominated; information on these scales therefore comes mainly from the 143 GHz and 217 GHz channels and the effective sky coverage drops to 37%.
We approximate these behaviours by generating our mock data in two components: one component running in the range with and the noise properties of the 100 GHz channel (see table 1), and a second in the range with and the noise properties of the 143 GHz and 217 GHz channels. We generate mock power spectra for ,  and polarisation autocorrelation, as well as crosscorrelation, assuming that temperature and polarisation are measured on identical parts of the sky.
The primary CMB signal becomes swamped by foreground anisotropies on the smallest scales. Rather than modelling explicitly these foregrounds, we opt to mimic their effects on parameter estimation by introducing a small scale cutoff, . Using simulated temperature data plus a Gaussian prior on , we find that we can recover the parameter sensitivities of real Planck data with .
Frequency/GHz 100 9.59 11.0 10.9 143 7.18 6.0 11.4 217 4.87 12.0 26.7 Table 1: Experimental specifications assumed in our Planck data simulations. The beam width and temperature pixel noise correspond to Planck’s actual inflight performance [1], while the polarisation pixel noise is set to the design performance goal given in [31].
4 Forecasting
We now describe our parameter sensitivity forecast for a Euclidlike photometric survey including a measurement of the cluster mass function. The forecast is based on the construction of a likelihood function for the mock data, whereby the survey’s sensitivities to cosmological parameters can be explored using Bayesian inference techniques.
4.1 Model parameter space
The forecast code has previously been used to probe extensions of the CDM model, either in the neutrino sector [10] or in dark energy [11]. Here, we take a similar approach, and extend the standard CDM model in the inflation sector. The noninflationary parameters of the extended model are the physical baryon density , the physical cold dark matter density , the dimensionless Hubble parameter , and the reionisation redshift ; the inflation sector of the model is described by the amplitude of the scalar fluctuation spectrum , the spectral index , the running of the spectral index , and the tensortoscalar ratio , all defined at the pivot scale. Thus, in total the model has eight free parameters:
(18) 
For the noninflationary part of the parameter space, our fiducial model is defined by the parameter values
(19) 
For the inflationary parameter space we will study a variety of fiducial models, but
(20) 
will be our most common choice.
4.2 Likelihood functions
Cmb The effective from CMB can be constructed using the method outlined in [29], and is given by
(21) 
where is the dimensional mock data covariance matrix, the model covariance matrix, and the expression has been normalised such that when .
We set when only and data are used, and 3 when data are also included.
Shear and galaxies The effective for shear and galaxies has exactly the same form as equation (21), but with and identified with and respectively, where . The dimension parameter is determined by the number of redshift bins used in the analysis: for shearonly (s), for galaxiesonly (g), and if we are interested in shear–galaxy crosscorrelation (x) as well.
Clusters The cluster likelihood function in one redshift and mass bin can be modelled at the most basic level as a Poisson distribution in the observed number of clusters , with the theoretical prediction as the mean. However, the imperfect completeness and efficiency of the survey necessitate that we rescale the uncertainty on each data point by an amount (see [11]). We accomplish this by defining an effective number of observed clusters , and likewise an effective theoretical prediction . Then, the effective probability distribution is
(22) 
In a real survey, the effective observed number of clusters in any one bin is necessarily an integer so that equation (22) applies directly. In our forecast, however, corresponds to the theoretical expectation value of the fiducial model which generally does not evaluate to an integer. To circumvent this inconvenience, we generalise the likelihood function (22) by linearly interpolating the logarithm of the discrete distribution in the interval , i.e.,
(23) 
The total cluster loglikelihood function is then obtained straightforwardly by summing over all redshift and mass bins.
5 Results
5.1 Binning
Before proceeding with a more detailed discussion of our results, we first test how our choice of redshift and mass binning affects parameter constraints on and for each of the three largescale structure observables. The study is performed in the same manner as in [10] and [11], i.e., we combine each of the three data sets with CMB data and test how parameter constraints vary with the number of reshift (shear, galaxy, and cluster data) and mass (cluster data) bins.
Figure 3 shows the parameter constraints when varying the number of redshift bins from 1 to 13. In the case of clusters, we also fix, for simplicity, the number of mass bins to be the same as the number of redshift bins . For an example of how the choice of disparate and values affects neutrino and dark energy constraints, see [11].
The results in figure 3 confirm the trends found in [10, 11], namely, that parameter constraints do not improve significantly beyond two redshift bins for shear data, about ten bins for galaxies, and for clusters. In fact, because we are considering parameters related to the initial conditions in this study, rather than parameters that affect the latetime evolution of the observables (e.g., the dark energy equation of state), the redshift binning is even less important here than was found in [10, 11]. For ease of comparison, however, we shall from now on use the same binning choices as [11], i.e., , , and . We emphasise again, however, that none of our results depends strongly on the binning choice.
5.2 Breaking degeneracies
Results from our main runs are shown in table 2, where we have derived parameter sensitivities from a Plancklike CMB experiment plus the three largescale structure observables, individually as well as in combination. Clearly, the sensitivities to and reflect strongly on the effective range available to the observables used: while adding galaxy or cluster data to CMB offer some improvement over the CMBonly sensitivities, both are outperformed by CMB+shear, which, as shown in figure 1, has a much longer lever arm.
Data  

c  0.011  0.0052  0.0074  0.028 
cs  0.0091  0.0030  0.0030  0.027 
cg  0.0046  0.0035  0.0048  0.027 
ccl  0.0068  0.0034  0.0064  0.026 
cscl  0.0066  0.0028  0.0029  0.026 
csgxcl  0.0032  0.0025  0.0017  0.026 
Interestingly, even though the sensitivities of galaxy and cluster data per se are generally inferior to what can be gained from shear measurements, when all used in combination the former two probes in fact contribute to breaking important parameter degeneracies in the shear power spectrum, leading to significant improvements in the the parameter sensitivities over using shear measurements alone. This effect has already been observed in [10] and [11] in relation to neutrino mass constraints. Our results in table 2 show that the same thing is happening also to the inflationary parameters (see also figure 4 for the corresponding twodimensional posteriors). For example, the sensitivity to improves from for CMB+shear, to for CMB+shear+galaxy+clusters+crosscorrelation.
A caveat is in place here. As discussed in section 3.4, in deriving the our parameter sensitivities we have assumed the linear galaxy bias to be perfectly known in every redshift bin. This is likely an overly optimistic assumption. For this reason, in addition to the “csgxcl” data combination, we consider also “cscl”, i.e., without galaxy data. The parameter sensitivities derived from this data combination would correspond roughly to those in the (equally unrealistic) case in which the galaxy bias is completely unknown and uncorrelated between the eleven redshift bins. Under this condition, we see from table 2 that most parameter sensitivities essentially revert back to their CMB+shear values, while mimics its CMB+cluster counterpart. Thus, how well the angular galaxy power spectrum breaks the parameter degeneracies in other largescale structure observables does depend strongly on how well the galaxy bias is understood. Should a viable bias model become available, the true parameter sensitivities would lie somewhere between the “csgxcl” and “cscl” extremes.
5.3 Dependence on the fiducial model
So far we have assumed the fiducial model (20) to be the true model of the universe. Here, we consider how the parameter sensitivities might change with respect to this assumption.
Table 3 shows the sensitivities of the data combination“csgxcl” to the inflationary parameters with respect to variations in the fiducial values of and . Evidently, with the exception of , the parameter sensitivities are quite independent of our choice of fiducial model. In contrast, the sensitivity to is almost a factor of three better when the fiducial model has , compared with the choice of . This difference can be put down to two effects:

If we assume the posterior to be a Gaussian with width parameter centred on the fiducial value, then the fact that the distribution is cut off at alone will cause the standard deviation of to depend on the fiducial value. The extreme example is = 0, in which case is a halfGaussian, and the corresponding standard deviation is . This is in contrast with the case where , for which the standard deviation approaches . By itself, however, this effect cannot entirely account for the observed difference in the sensitivities.

The cosmic variance contribution to the effective uncertainty of the CMB spectra depends on the chosen . This is generally a negligible effect for the and spectra (even at low ), because contributions from scalar perturbations always dominate over tensor contributions. For the polarisation power spectrum, however, which in our case is sourced only by tensor perturbations, the larger cosmic variance following from a larger value can have an important effect on the sensitivity to .
To test this effect, we repeat the analysis for the fiducial models and , but this time using only and data in the CMB component. We find and , respectively, a much less dramatic disparity than that between and we found earlier with the mode included. In fact, the ratio between and is approximately , which is very close to the ratio of expected from the effect discussed above in point 1. This confirms that the larger mode error bars accompanying a larger value are indeed the main culprit for the reduced sensitivity to .
Data  

csgxcl  0  0  0.0031  0.0025  0.0017  0.0090 
csgxcl  0.008  0  0.0031  0.0025  0.0018  0.0088 
csgxcl  0  0.1  0.0032  0.0025  0.0017  0.026 
csgxcl  0.008  0.1  0.0032  0.0025  0.0018  0.026 
5.4 Comparison with other works
Reference [34] (also quoted in [35]) found that Plancklike CMB data combined with data from a Euclidlike spectroscopic survey could attain a sensitivity of and , comparable to our 0.0025 and 0.0021 from the “csgxcl” data combination. There are however some subtle differences between their and our analyses.
Firstly, reference [34] made much more optimistic assumptions about their CMB data than we do in our study. The sensitivities from Plancklike CMB data alone were assumed to be and , compared with our 0.0053 and 0.0075. Secondly, the addition of largescale structure observables in our case improves the constraining power by a factor of 2.12 for and 3.57 for , whereas for the spectroscopic survey studied in [34] the corresponding factors are 1.76 and 1.7.
That photometric data carry more weight than do spectroscope data is most likely due to the former’s larger effective survey volume and hence statistical power. After all, parameters like and do not produce any redshiftdependent effects, and therefore benefit little from highaccuracy redshift measurements. Rather, they induce slow variations with , and are hence better constrained by broadband power spectrum measurements with high statistics.
6 Implications for single field slow roll inflation
Prima facie, the running of the spectral index is simply a phenomenological parameter. However, as it turns out, the large improvements in sensitivity to afforded by Euclidlike data can potentially have a profound impact on our understanding of the dynamics of slowroll inflation. In the slowroll picture, the four parameters describing the primordial perturbations can be expressed in terms of the Hubble parameter and its derivatives during inflation. In addition, as discussed for instance in [6], the running of the spectral index is intimately related to the total amount of inflation, with large negative values implying that inflation ends prematurely. In the following section, we briefly remind the reader of the theoretical background.
6.1 Slowroll inflation
The evolution of the Hubble parameter during inflation is determined by an infinite hierarchy of the Hubble slowroll parameters [36]:
(24) 
where denotes a derivative with respect to the scalar field , and the usual slow roll parameters are and . Inflation continues as long as .
If the hierarchy is truncated at (which is in fact a very good approximation for the majority of generic singlefield models of inflation), the solution for the Hubble parameter as a function of then takes the form [37]
(25) 
where the subscript “” refers to the value of the respective quantity at the moment the mode leaves the horizon. Without loss of generality we can set .
Given a set of values for the phenomenological observables , and , the corresponding slowroll parameters, up to and including and up to second order, can be established via the relations [38, 39],
(26) 
where , and is the Euler–Mascheroni constant. Solving these relations for , , and then allows for a determination of the number of foldings until the end of inflation, , via
(27) 
where is the field value at the end of inflation, i.e., .
For successful inflation, one typically requires , with the exact number depending on the details of the reheating process (see, e.g., [3]). As a consequence, only a certain region of (, , )space is actually consistent with our assumptions about inflation. If observational data should decisively favour parameter combinations outside of this region, then some of the assumptions going into expression (25) must be wrong. One possibility in that case would be the presence of additional, sizeable higherorder terms from the hierarchy in order to describe the evolution of during inflation; the mere presence of these higherorder terms would rule out most simple models of singlefield slowroll inflation. Alternatively, one could invoke multifield dynamics, particularly if data should prefer the region. Either way, such a discovery would completely change the current picture of inflation. Let us now turn to exploring the potential of a future Euclidlike survey in this regard.
6.2 Slowroll analysis and results
We consider first a fiducial model without tensor perturbations. Figure 5 shows that, given our assumption of for , a negative running of is required in order to remain in the desired region of foldings (i.e., ). However, the projected constraints for a fiducial model with still comfortably overlap with this region, so even if the true running of the spectral index was zero (and thus our assumption invalid due to inflation lasting for too long), the data would not be able to distinguish the two cases.
If on the other hand the true running should be large and negative, inflation does not continue for long enough to achieve the required number of foldings. For , the entire 95%credible region lies outside the band of sufficient inflation, which means that the data become able to discriminate between models that are consistent with equation (25) and those that are not.
If the underlying model has , the discriminatory power of the data will further improve. This is demonstrated in figure 6 for a tensortoscalar ratio of : the viable region has now migrated closer to zero running and is significantly compressed with respect to the case. Since implies , the second order terms in the expression for in equation (26) become of the same order as ; for small enough and close to zero, these terms can even dominate, and force to become negative. This can cause the inflaton potential, , to develop a local minimum, thereby trapping the field; inflation continues indefinitely under these circumstances, and the lines of constant stack up near the bottom right corner of the figure.
7 Conclusions
In this work we have studied how much the sensitivity to the observational parameters related to inflation can be improved by combining Plancklike CMB measurements with largescale structure observables from a future largevolume photometric survey such as LSST or Euclid. Largescale structure observables constrain length scales significantly smaller than the reach of CMB measurements, and provide an extended lever arm that is especially useful for measuring the scalar fluctuation spectral index and testing its possible running.
Using a combination of angular shear power spectrum, angular galaxy power spectrum and cluster mass function, together with Plancklike measurements of the CMB temperature and  and polarisation, we find that running of the spectral index can potentially be probed to a precision better than at , about a factor of five improvement over current measurements. The sensitivity to the spectral index likewise sees significant, although less dramatic, improvements: the corresponding precision is , roughly a factor of three tighter than current constraints. Not surprisingly, largescale structure observables offer no direct information on the tensortoscalar ratio . There is likewise no significant correlation between and those parameters constrainable by these observables, as indicated by the fact that the sensitivity to hardly depends on whether or which largescale structure data have been included in the analysis.
Importantly, the greatly enhanced constraining power with respect to the running of the spectral index will open a new window into the physics of inflation and test the paradigm of minimal singlefield slowroll inflation. If , as hinted at by Planck CMB measurements, then the requirement that inflation lasts long enough will necessitate the presence of nonnegligible higherorder derivatives, or multifield dynamics, which are absent in the simplest inflationary models.
The next generation of largescale structure surveys will thus truly represent a new milestone in observational cosmology, with the potential to unravel the mysteries of not only the neutrino and dark energy sectors, but also possibly revolutionising our knowledge of the first moments of the universe.
Acknowledgements
We acknowledge computing resources obtained via a grant from the Danish eInfrastructure Cooperation (DeIC).
Footnotes
 A copy of the forecast software based on the July 2014 version of CosmoMC can be downloaded from http://jhamann.web.cern.ch/jhamann/simdata/simdata.tar.gz.
 Our modelling of the CMB likelihood function assumes the signal to be Gaussiandistributed, and clearly does not capture the large nonGaussian covariance of the lensed mode [32]. This nonGaussianity is very important for future highresolution, lownoise polarisation experiments [33], but can be neglected for experiments of Plancklike sensitivity, such as considered in this work.
 Note that with current data, the bestfit is about 50% larger than this threshold.
References
 Planck Collaboration, P. Ade et al., “Planck 2013 results. I. Overview of products and scientific results,” arXiv:1303.5062 [astroph.CO].
 WMAP Collaboration, C. Bennett et al., “NineYear Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Final Maps and Results,” Astrophys.J.Suppl. 208 (2013) 20, arXiv:1212.5225 [astroph.CO].
 Planck Collaboration, P. Ade et al., “Planck 2013 results. XXII. Constraints on inflation,” arXiv:1303.5082 [astroph.CO].
 BICEP2 Collaboration, P. Ade et al., “Detection of BMode Polarization at Degree Angular Scales by BICEP2,” Phys.Rev.Lett. 112 (2014) 241101, arXiv:1403.3985 [astroph.CO].
 Planck Collaboration, P. Ade et al., “Planck 2013 results. XVI. Cosmological parameters,” arXiv:1303.5076 [astroph.CO].
 R. Easther and H. Peiris, “Implications of a Running Spectral Index for Slow Roll Inflation,” JCAP 0609 (2006) 010, arXiv:astroph/0604214 [astroph].
 LSST Science Collaborations, LSST Project Collaboration, P. A. Abell et al., “LSST Science Book, Version 2.0,” arXiv:0912.0201 [astroph.IM].
 R. Laureijs, J. Amiaux, S. Arduini, J.L. Augueres, J. Brinchmann, et al., “Euclid Definition Study Report,” arXiv:1110.3193 [astroph.CO].
 A. Lewis and S. Bridle, “Cosmological parameters from CMB and other data: A Monte Carlo approach,” Phys.Rev. D66 (2002) 103511, arXiv:astroph/0205436 [astroph].
 J. Hamann, S. Hannestad, and Y. Y. Y. Wong, “Measuring neutrino masses with a future galaxy survey,” JCAP 1211 (2012) 052, arXiv:1209.1043 [astroph.CO].
 T. Basse, O. E. Bjælde, J. Hamann, S. Hannestad, and Y. Y. Y. Wong, “Dark energy properties from large future galaxy surveys,” JCAP 1405 (2014) 021, arXiv:1304.2321 [astroph.CO].
 B. Audren, J. Lesgourgues, S. Bird, M. G. Haehnelt, and M. Viel, “Neutrino masses and cosmological parameters from a Euclidlike survey: Markov Chain Monte Carlo forecasts including theoretical errors,” JCAP 1301 (2013) 026, arXiv:1210.2194 [astroph.CO].
 N. Kaiser, “Weak lensing and cosmology,” Astrophys.J. 498 (1998) 26, arXiv:astroph/9610120 [astroph].
 D. C. Koo, N. P. Vogt, A. C. Phillips, R. Guzman, K. Wu, et al., “Redshift field galaxies observed with the Keck Telescope and the HST,” Astrophys.J. 469 (1996) 535, arXiv:astroph/9604113 [astroph].
 W. Hu, “Power spectrum tomography with weak lensing,” Astrophys.J. 522 (1999) L21–L24, arXiv:astroph/9904153 [astroph].
 C. Sánchez, M. C. Kind, H. Lin, R. Miquel, F. Abdalla, et al., “Photometric redshift analysis in the Dark Energy Survey Science Verification data,” arXiv:1406.4407 [astroph.IM].
 Z.M. Ma, W. Hu, and D. Huterer, “Effect of photometric redshift uncertainties on weak lensing tomography,” Astrophys.J. 636 (2005) 21–29, arXiv:astroph/0506614 [astroph].
 R. K. Sheth and G. Tormen, “Large scale bias and the peak background split,” Mon.Not.Roy.Astron.Soc. 308 (1999) 119, arXiv:astroph/9901122 [astroph].
 S. Wang, J. Khoury, Z. Haiman, and M. May, “Constraining the evolution of dark energy with a combination of galaxy cluster observables,” Phys.Rev. D70 (2004) 123008, arXiv:astroph/0406331 [astroph].
 T. Hamana, M. Takada, and N. Yoshida, “Searching for massive clusters in weak lensing surveys,” Mon.Not.Roy.Astron.Soc. 350 (2004) 893, arXiv:astroph/0310607 [astroph].
 Y. M. Bahe, I. G. McCarthy, and L. J. King, “Mock weak lensing analysis of simulated galaxy clusters: bias and scatter in mass and concentration,” Mon.Not.Roy.Astron.Soc. 421 (2012) 1073–1088, arXiv:1106.2046 [astroph.CO].
 M. R. Becker and A. V. Kravtsov, “On the Accuracy of Weak Lensing Cluster Mass Reconstructions,” Astrophys.J. 740 (2011) 25, arXiv:1011.1681 [astroph.CO].
 M. J. White, L. van Waerbeke, and J. Mackey, “Completeness in weak lensing searches for clusters,” Astrophys.J. 575 (2002) 640–649, arXiv:astroph/0111490 [astroph].
 F. Feroz, P. Marshall, and M. Hobson, “Cluster detection in weak lensing surveys,” arXiv:0810.0781 [astroph].
 F. Pace, M. Maturi, M. Meneghetti, M. Bartelmann, L. Moscardini, et al., “Testing the reliability of weak lensing cluster detections,” arXiv:astroph/0702031 [ASTROPH].
 J. F. Hennawi and D. N. Spergel, “Mass selected cluster cosmology. 1: Tomography and optimal filtering,” Astrophys.J. (2004) , arXiv:astroph/0404349 [astroph].
 M. Lima and W. Hu, “Selfcalibration of cluster dark energy studies: Observablemass distribution,” Phys.Rev. D72 (2005) 043006, arXiv:astroph/0503363 [astroph].
 S. Wang, Z. Haiman, W. Hu, J. Khoury, and M. May, “Weighing neutrinos with galaxy cluster surveys,” Phys.Rev.Lett. 95 (2005) 011302, arXiv:astroph/0505390 [astroph].
 L. Perotto, J. Lesgourgues, S. Hannestad, H. Tu, and Y. Y. Y. Wong, “Probing cosmological parameters with the CMB: Forecasts from full Monte Carlo simulations,” JCAP 0610 (2006) 013, arXiv:astroph/0606227 [astroph].
 Planck Collaboration, P. Ade et al., “Planck 2013 results. XV. CMB power spectra and likelihood,” arXiv:1303.5075 [astroph.CO].
 Planck Collaboration, J. Tauber et al., “The Scientific programme of Planck,” arXiv:astroph/0604069 [astroph].
 K. M. Smith, W. Hu, and M. Kaplinghat, “Weak lensing of the CMB: Sampling errors on Bmodes,” Phys.Rev. D70 (2004) 043002, arXiv:astroph/0402442 [astroph].
 C. Li, T. L. Smith, and A. Cooray, “NonGaussian Covariance of CMB Bmodes of Polarization and Parameter Degradation,” Phys.Rev. D75 (2007) 083501, arXiv:astroph/0607494 [astroph].
 Z. Huang, L. Verde, and F. Vernizzi, “Constraining inflation with future galaxy redshift surveys,” JCAP 1204 (2012) 005, arXiv:1201.5955 [astroph.CO].
 Euclid Theory Working Group Collaboration, L. Amendola et al., “Cosmology and fundamental physics with the Euclid satellite,” Living Rev.Rel. 16 (2013) 6, arXiv:1206.1225 [astroph.CO].
 A. R. Liddle, P. Parsons, and J. D. Barrow, “Formalizing the slow roll approximation in inflation,” Phys.Rev. D50 (1994) 7222–7232, arXiv:astroph/9408015 [astroph].
 A. R. Liddle, “Inflationary flow equations,” Phys.Rev. D68 (2003) 103504, arXiv:astroph/0307286 [astroph].
 W. H. Kinney, “Inflation: Flow, fixed points and observables to arbitrary order in slow roll,” Phys.Rev. D66 (2002) 083508, arXiv:astroph/0206032 [astroph].
 H. Peiris and R. Easther, “Recovering the Inflationary Potential and Primordial Power Spectrum With a Slow Roll Prior: Methodology and Application to WMAP 3 Year Data,” JCAP 0607 (2006) 002, arXiv:astroph/0603587 [astroph].