The clustering of galaxies and galaxy clusters: constraints on primordial nonGaussianity from future widefield surveys
Abstract
We investigate the constraints on primordial nonGaussianity with varied bispectrum shapes that can be derived from the power spectrum of galaxies and clusters of galaxies detected in future wide field optical/nearinfrared surveys. Having in mind the proposed ESA space mission Euclid as a specific example, we combine the spatial distribution of spectroscopically selected galaxies with that of weak lensing selected clusters. We use the physically motivated halo model in order to represent the correlation function of arbitrary tracers of the Large Scale Structure in the Universe. As naively expected, we find that galaxies are much more effective in jointly constrain the level of primordial nonGaussianity and the amplitude of the matter power spectrum than clusters of galaxies, due to the much lower abundance of the latter that is not adequately compensated by the larger effect on the power spectrum. Nevertheless, combination of the galaxy power spectrum with the clustergalaxy cross spectrum can decrease the error on the determination of by up to a factor of . This decrement is particularly evident for the less studied nonGaussian bispectrum shapes, the socalled enfolded and the orthogonal ones. Setting constraints on these models can shed new light on various aspects of the physics of the early Universe, and it is hence of extreme importance. By combining the power spectra of clusters and galaxies with the clustergalaxy cross spectrum we find constraints on primordial nonGaussianity of the order a few, competitive and possibly superior to future CMB experiments.
1 Introduction
It is now commonly accepted that the formation of structures in the Universe originated from seed density fluctuations in the dark matter fluid that were laid down during inflation (Hawking & Moss, 1983). Gravitational instability has the effect of amplifying density perturbations, that at a given point enter the nonlinear regime and collapse to form clumps, voids, and more in general the complex Large Scale Structure (LSS henceforth) that we observe today. For a given distribution of the initial conditions, the statistical properties of the LSS are determined uniquely by the subsequent expansion history of the Universe that is, ultimately, on its matter and energy content. In the past decade much effort has been directed toward understanding the effect of dark energy on the formation of structures (see Cunha, Huterer, & Frieman 2009; Grossi & Springel 2009; Samushia & Ratra 2009; De Boni et al. 2010; Mortonson, Hu, & Huterer 2010 for recent examples) in order to gain insights on its nature, specifically whether it truly is a cosmological constant or it does have some kind of evolution with cosmic time.
After several pioneering works (Messina et al., 1990; Moscardini et al., 1991; Weinberg & Cole, 1992), only recently has the question of the initial conditions gained renewed attention, specifically about the shape of the primordial density fluctuations distribution. While the simplest models of inflation (single slowrolling scalar field) predicts this distribution to be virtually indistinguishable from a Gaussian, a plethora of more complex models have been proposed that, in addition to solving the standard cosmological problems, allow for significant and possibly scale dependent deviations from Gaussianity. The study of nonGaussian cosmologies and the effect they have on the formation and evolution of cosmic structures is thus extremely important in order to rule out inflationary models, and hence to have a better handle on the physics of the early Universe. Moreover, studies on possible detectability of primordial nonGaussianity are very timely, given the recent claims for a positive skewness of the primordial density fluctuations distribution coming from the Cosmic Microwave Background (CMB) and from the angular correlation function of radioselected quasars (Xia et al., 2010).
The problem of constraining deviations from primordial Gaussianity by means different from the CMB intrinsic anisotropies has recently attracted much efforts in the literature, with studies directed towards the abundance of nonlinear structures (Matarrese, Verde, & Jimenez 2000; Verde et al. 2000; Mathis, Diego, & Silk 2004; Grossi et al. 2007, 2009; Maggiore & Riotto 2010c), halo biasing (Dalal et al. 2008; McDonald 2008; Fedeli, Moscardini, & Matarrese 2009), galaxy bispectrum (Sefusatti & Komatsu, 2007; Jeong & Komatsu, 2009), mass density distribution (Grossi et al., 2008) and topology (Matsubara, 2003; Hikage et al., 2008), cosmic shear (Fedeli & Moscardini 2010, Pace et al. 2010), integrated SachsWolfe effect (Afshordi & Tolley, 2008; Carbone, Verde, & Matarrese, 2008), Ly flux from lowdensity intergalactic medium (Viel et al., 2009), cm fluctuations (Cooray, 2006; Pillepich, Porciani, & Matarrese, 2007) and reionization (Crociani et al., 2009).
In this work we focused attention on the spatial distribution of galaxies and galaxy clusters as tracers of the LSS. We were particularly interested in comparing the performances of the power spectra of the two individual tracers in constraining primordial nonGaussianity, and evaluate the improvements in forecasted constraints given by the addition of the clustergalaxy cross power spectrum. Throughout the paper we assumed a fiducial future allsky optical/nearinfrared survey on the model of Euclid (Laureijs, 2009; Beaulieu et al., 2010). In order to fully exploit the potentials of both the imaging and the spectroscopy part of Euclid we considered galaxies as spectroscopically selected according to their H flux, and galaxy clusters as selected to be the high signaltonoise ratio (S/N) peaks in fullsky cosmic shear maps. This approach has the advantage of allowing treatment of the galaxy and cluster samples as independent. The results obtained here are relevant for other planned missions with a similar concept to Euclid, such as the NASA Wide Field Infrared Survey (WFIRST).
The rest of the paper is organized as follows. In Sections 2 and 3 we summarize the nonGaussian models that we have explored in this work, as well as their effect on the mass function and large scale bias of dark matter halos. In Section 4 we describe the halo model, the physically motivated framework that we adopted for modeling the power spectrum of clusters and galaxies as well as the clustergalaxy cross spectrum. In Section 5 we summarize our results and in Section 6 we discuss them, with particular emphasis on alternative survey configurations. Finally, in Section 7 we present our conclusions. Throughout this work we refer to the fiducial cosmological model as the one defined by the parameter set derived by the latest analysis of the WMAP data (Komatsu et al., 2010). This means that density parameters for matter, cosmological constant, and baryons are equal to , , and , respectively, the Hubble constant equals , and the normalization of the matter power spectrum is set by .
2 NonGaussian cosmologies
Extensions of the most standard model of inflation (Starobinskiǐ, 1979; Guth, 1981; Linde, 1982) can produce substantial deviations from a Gaussian distribution of primordial density and potential fluctuations (see Bartolo et al. 2004; Chen 2010; Desjacques & Seljak 2010 for recent reviews). The amount and shape of this deviation depend critically on the kind of nonstandard inflationary model that one has in mind, as will be detailed later on.
A particularly convenient (although not unique) way to describe generic deviations from a Gaussian distribution consists in writing the gaugeinvariant Bardeeen’s potential as the sum of a Gaussian random field and a quadratic correction (Salopek & Bond, 1990; Gangui et al., 1994; Verde et al., 2000; Komatsu & Spergel, 2001), according to
(1) 
The parameter in Eq. (1) determines the amplitude of nonGaussianity, and it is in general dependent on the scale. The symbol denotes convolution between functions, and reduces to standard multiplication upon constancy of . In the following we adopted the largescale structure convention (as opposed to the CMB convention, see Afshordi & Tolley 2008; Carbone et al. 2008; Pillepich, Porciani, & Hahn 2009 and Grossi et al. 2009) for defining the fundamental parameter . According to this, the primordial value of has to be linearly extrapolated at , and as a consequence the constraints given on by the CMB have to be raised by per cent to comply with this paper’s convention (see also Fedeli, Moscardini, & Matarrese 2009 for a concise explanation).
In the case in which the potential is a random field with a nonGaussian probability distribution. Therefore, the field itself cannot be described by the power spectrum alone, rather higherorder moments are needed. The dominant higherorder contribution is generically given by the bispectrum . Only in very peculiar situations the bispectrum vanishes, and one has to resort to the trispectrum or higherorder correlations. The bispectrum is the Fourier transform of the threepoint correlation function and it can hence be implicitly defined as
(2) 
As mentioned above understanding the shape of nonGaussianity is of fundamental importance in order to pinpoint the physics of the early universe and the evolution of the inflaton field in particular. For this reason, in this work we considered four different shapes of the potential bispectrum, arising from different modifications of the standard inflationary scenario. They are all described in the following.
Local shape
The standard singlefield inflationary scenario generates negligibly small deviations from Gaussianity. These deviations are said to be of the local shape, and the related bispectrum of the Bardeen’s potential is maximized for squeezed configurations, where one of the three wavevectors has much smaller magnitude than the other two. In this case the parameter must be a constant (see however Byrnes et al. 2010a, b), and it is expected to be of the same order of the slowroll parameters (Falk, Rangarajan, & Srednicki, 1993), that are very close to zero.
However nonGaussianities of the local shape can also be generated in the case in which an additional light scalar field, different from the inflaton, contributes to the observed curvature perturbations (Babich, Creminelli, & Zaldarriaga, 2004). This happens, for instance, in curvaton models (Sasaki, Väliviita, & Wands, 2006; Assadullahi, Väliviita, & Wands, 2007) or in multifields models (Bartolo, Matarrese, & Riotto, 2002; Bernardeau & Uzan, 2002; Huang, 2009). In this case the parameter is allowed to be substantially different from zero, and the bispectrum of the primordial potential assumes the simple form
(3) 
Due to its simplicity, the local model is the most studied one, especially in terms of cosmological simulations (see however Wagner, Verde, & Boubekeur 2010).
Equilateral shape
In some inflationary models the kinetic term of the inflaton Lagrangian is not standard, containing higherorder derivatives of the field itself. One significant example of this is the DBI model (Alishahiha, Silverstein, & Tong 2004; Silverstein & Tong 2004, see also ArkaniHamed et al. 2004; Seery & Lidsey 2005; Li, Wang, & Wang 2008). In this case the primordial bispectrum is maximized for configurations where the three wavevectors have approximately the same amplitude, and it takes the form (Creminelli et al., 2007)
(4)  
The function in the first line of the previous equation represents a running of the parameter that we have the liberty to insert since this parameter is not forced to be constant in the present case. This running has been considered as an actual part of the equilateral bispectrum in all the calculations that follow. It reads (Chen, 2005; LoVerde et al., 2008; Crociani et al., 2009)
(5) 
We adopted the exponent , that increase the level of nonGaussianity at scales smaller than that corresponding to Mpc. This coincides with the larger multipole used in the CMB analysis by the WMAP team (Komatsu et al., 2009, 2010), . When referring to the equilateral shape in the rest of this paper we always mean the bispectrum given by Eq. (4) with , unless otherwise noted.
Enfolded shape
For deviations from Gaussianity evaluated in the regular BunchDavies vacuum state, the primordial potential bispectrum is of local or equilateral shape, depending on whether or not higherorder derivatives play a significant role in the evolution of the inflaton field. If the BunchDavies vacuum hypothesis is dropped, the resulting bispectrum is maximal for squashed configurations (Chen et al., 2007; Holman & Tolley, 2008).
*ME09.1 found a template that describes very well the properties of this enfoldedshape bispectrum (see also Creminelli et al. 2010), that reads
(6)  
No running of the parameter was introduced for both the enfolded and orthogonal shapes. The reasons are discussed further below.
Orthogonal shape
A shape of the bispectrum can be constructed that is nearly orthogonal to both the local and equilateral forms (Senatore, Smith, & Zaldarriaga, 2010). In this case the potential bispectrum is well approximated by the following template
(7)  
Constraints on the level of nonGaussianity compatible with the CMB in the local, equilateral and orthogonal scenarios were recently given by the WMAP team (Komatsu et al., 2010), while constraints on enfolded nonGaussianity from LSS were given by Verde & Matarrese (2009)
Although there is no theoretical prescription against a running of the parameter with the scale in the enfolded and orthogonal shapes, we decided not to include one. The reason for this is that there is no first principle that can guide one in the choice of a particular kind of running, and until now no work has addressed the problem of a running for these shapes (Fergusson & Shellard, 2009; Fergusson, Liguori, &
Shellard, 2010). As can be noted, in all nonlocal cases except the enfolded one we defined the level of nonGaussianity by equating the bispectrum normalization in the equilateral limit to the same quantity for the local shape (see the discussion in Fergusson & Shellard 2009).
3 Halo mass function and linear bias
Primordial nonGaussianity produces modifications in the statistics of density peaks, resulting in differences in the mass function of cosmic objects and the bias of dark matter halos with respect to the underlying smooth density field. In the following we summarize how these modifications have been taken into account in the present work.
3.1 Mass function
For the nonGaussian modification to the mass function of cosmic objects we adopted the prescription of LoVerde et al. (2008). The main assumption behind it is that the effect of primordial nonGaussianity on the mass function is independent of the prescription adopted to describe the mass function itself. This means that, if and are the nonGaussian and Gaussian mass functions, respectively, computed according to the Press & Schechter (1974) recipe, we can define a correction factor . Then, the nonGaussian mass function computed according to an arbitrary prescription, can be related to its Gaussian counterpart through
(8) 
In order to compute , and hence , LoVerde et al. (2008) performed an Edgeworth expansion (Blinnikov & Moessner, 1998) of the probability distribution for the smoothed density fluctuations field, truncating it at the linear term in . The resulting Press & Schechter (1974)like mass function reads
(9)  
In the previous equation is the comoving matter density in the Universe, is the rms of density fluctuations smoothed on a scale corresponding to the mass , and , where is the growth factor and is the critical linear overdensity for collapse. The function is the reduced skewness of the nonGaussian distribution, and the skewness can be computed as
(10)  
The last thing that remains to be defined is the function , that relates the density fluctuations smoothed on some scale to the respective peculiar potential,
(11) 
where is the matter transfer function and is the Fourier transform of the tophat window function.
In this work we adopted the Bardeen et al. (1986) matter transfer function, with the shape factor correction of Sugiyama (1995). This reproduces fairly well the more sophisticated recipe of Eisenstein & Hu (1998) except for the presence of the baryon acoustic oscillation, that anyway is not of interest here. We additionally adopted as reference mass function the prescription of Sheth & Tormen (2002) (see Jenkins et al. 2001; Warren et al. 2006; Tinker et al. 2008 for alternative prescriptions). Other approaches also exist for computing the nonGaussian correction to the mass function, that give results is broad agreement with those adopted here (Matarrese, Verde, & Jimenez 2000; D’Amico et al. 2010). These semianalytic prescriptions are known to be in good agreement with cosmological simulations having local nonGaussian initial conditions (Grossi et al., 2007, 2009; Desjacques & Seljak, 2010), while more recently the same has been proven to be true also for more generic primordial bispectrum shapes (Wagner, Verde, & Boubekeur, 2010).
In the left panel of Figure 1 we show the mass dependence of the reduced skewness for the four nonGaussian shapes that have been detailed above. It is interesting to note that in the orthogonal case the skewness for positive would be negative, thus bringing to a reduction in the abundance of massive objects. Also, the equilateral shape is the only case in which is not monotonic. In computing the nonGaussian corrections to the mass function we have taken into account the modification to the critical overdensity for collapse suggested by Grossi et al. (2009) (see also Maggiore & Riotto 2010a, b, c), according to which , with .
3.2 Halo bias
Recently much attention has been devoted to the effect of primordial nonGaussianity on halo bias, and the use thereof for constraining (Dalal et al. 2008; Verde & Matarrese 2009; Carbone, Verde, & Matarrese 2008). In particular, Matarrese & Verde (2008) have shown that primordial nonGaussianity introduce a scale dependence on the large scale halo bias. This peculiarity allows to place already stringent constraints from existing data (Slosar et al., 2008; Afshordi & Tolley, 2008).
The nonGaussian halo bias can be written in a relatively straightforward way in terms of its Gaussian counterpart as (Carbone, Mena, & Verde, 2010)
(12) 
where the function encapsulates all the scale dependence of the nonGaussian correction to the bias, and reads
(13)  
where . In the simple case of local bispectrum shape it can be shown that the function should scale as at large scales, so that a substantial boost (if ) in the halo bias is expected at those scales.
In the right panel of Figure 1 we show the scale dependence of the function for a fixed halo mass. In the remainder of this paper we adopted the Sheth, Mo, & Tormen (2001) formula for the Gaussian bias. Thus, since we express the correction to the halo bias as a function of the Gaussian bias itself, we did not need to take into account the ellipsoidal collapse correction suggested by Grossi et al. (2009).
4 Halo model of the largescale structure
The halo model (Ma & Fry, 2000; Seljak, 2000) is a physical framework that allows the description of the correlation function of dark matter as well of different tracers of the LSS such as galaxies ad galaxy clusters. It is based on the assumption that the power spectrum of particles (either dark matter particles or tracers) is given by the sum of two contributions: particle pairs residing in the same structure, and particle pairs residing in two different structures. The implicit hypothesis underlying this assumption is that no particles are found outside bound structures. Accordingly, the cross spectrum of two different kinds of particles x and y (if x = y then is the power spectrum of particle type x) can be written as the sum of two terms,
(14) 
The first one is named halo term, while the second is the halo term. From the discussion above, it is immediately obvious that the first term dominates on very small scales, while the second is dominant at large scales.
The exact form of the two terms depend on the exact kind of particle that we are considering. However, a common feature is that the halo term should depend upon the bias of dark matter halos, since it should represents particle pairs residing in separated structures (Cooray & Sheth, 2002).
4.1 Galaxy spectrum
In this case we set x y g, and the terms contributing to the galaxy power spectrum can be written as
(15) 
and
(16) 
where
(17) 
Note that we have included a scale dependence in the halo bias, in order to account for the effect of primordial nonGaussianity. The functions and are the first and second moment respectively of the halo occupation distribution , which represents the probability that galaxies reside inside a dark matter halo of mass . The function is the linear matter power spectrum.
We observe that the previous integrals do not actually extend down to , since there is a minimum halo mass below which no galaxy formation is possible, i.e., for . Additionally, not all galaxies at all redshifts are accessible to observations, hence should be the minimum galaxy halo mass that, at a given redshift, enters in our fiducial catalogue. We come back on the issue of the minimum mass further below. The function is the mean number density of galaxies, and reads
(18) 
The quantity in Eqs. (15) and (17) represents instead the Fourier transform of the galaxy number density inside dark matter halos of mass and redshift . We set it equal to the number density of dark matter particles, that is . For the dark matter density profile we adopt a Navarro, Frenk, & White (1996) shape (NFW henceforth, see also Navarro, Frenk, & White 1995, 1997) with a concentrationmass relation that reproduces simulated matter power spectra (see Huffenberger & Seljak 2003; Fedeli & Moscardini 2010 for discussions). We discuss more on the adopted density profile in Section 6.
We make the assumption that if a dark matter halo hosts at least one galaxy, then one galaxy should sit at the center of the halo itself. As a consequence, the halo term of the galaxy power spectrum consists of two contributions: galaxy pairs that involve the central objects and all the other pairs. As discussed by Cooray & Sheth (2002) we can selfconsistently assign to the first class () and to the second (). As for the moments of the halo occupation distribution we adopted simple fitting forms that reproduce the results of semianalytic galaxy formation models, that is
(19) 
and
(20) 
(Cooray & Sheth, 2002; Cooray, 2004; Hütsi & Lahav, 2008), where if and otherwise. This corresponds at assuming a Poisson distribution at high halo masses with deviations thereof occurring at low masses.
The three free parameters , and depend on the type of galaxy that is under consideration (Cooray & Sheth, 2002). Typically, for blue galaxies , , for and otherwise. For red galaxies we have instead , , and . Obviously, the average number of galaxies inside a halo of some mass is the sum of the average number of red galaxies and of blue galaxies. Now, the distinction between blue and red galaxies can be rather qualitative, and dependent on the framework we are interested in (e.g. observations vs. semianalytic models). The slitless spectroscopic instrument planned for Euclid will be mainly sensitive to the flux and equivalent width of emission lines, mainly H, and thus will mostly observe starforming galaxies. Although this category includes dustobscured starburst objects, that are not properly blue, we believe the parameter set for blue galaxies to be more suited for the target of Euclid. Hence we shall adopt it in the remainder of this work. In Section 6 we discuss how our results are changed if we consider a multislit spectroscopic instrument, that select a different galaxy population.
4.2 Cluster spectrum
For galaxy clusters the situation is less elaborated, since it is commonly assumed that only one cluster is contained inside each dark matter halo that is massive enough (this is obviously not true for galaxies and dark matter particles). Consequently, the halo term vanishes, , while the halo term reads
(21) 
where
(22) 
In this case as well there is a minimum halo mass below which no cluster is formed (or, below which clusters are not observable), and the integrals are to be evaluated above that mass. The function then reads
(23) 
Further below we detail how the issue of the minimum cluster mass is addressed in this paper.
In Figure 2 we show the power spectrum of galaxy clusters alongside that for galaxies, and compare them to the dark matter power spectrum obtained with the halo model (Fedeli & Moscardini, 2010). At linear scales clusters are more biased than galaxies as is to be expected. Moreover, since the power spectrum of clusters is made only by the halo term, their correlation drops off substantially at small scales as compared to other tracers.
4.3 Clustergalaxy cross spectrum
Extending the previous results, Hütsi & Lahav (2008) derived the two halo model contributions to the cross correlation of clusters and galaxies. They can be written, respectively, as
(24) 
and
(25) 
In Figure 2 we also show the cross spectrum of galaxy clusters and galaxies. As can be seen, at large scales the cross spectrum stays in between the spectra of galaxies and clusters. However, it shows more power on small scales than both clusters and galaxies. The fact that the scale dependence of the cross spectrum is different from that of the galaxy spectrum and cluster spectrum means that it carries a different kind of information on the spatial distribution of LSS tracers, thus the crosscorrelation can make up an important improvement when constraining cosmological parameters.
4.4 Minimum masses
As mentioned above, in order to compute the cluster and galaxy power spectra it is necessary to select a minimum mass for both classes of objects. In general this minimum mass depends on redshift. Let us start with clusters. Bergé, Amara, & Réfrégier (2010) have estimated the weak lensing selection function for galaxy clusters given a Euclidlike survey and different S/N detection thresholds (see also Laureijs 2009). In order to be conservative, we adopted S/N , and fitted the corresponding curve in Figure 1 of Bergé, Amara, & Réfrégier (2010). The result is displayed by the red dashed curve in Figure 3, which represents the we adopted in this work.
For galaxies, we adopted the minimum halo mass which, given the halo occupation moments reported in Eqs. (19) and (20), reproduces in the reference Gaussian cosmology the effective bias predicted for H selected galaxies by Euclid (Orsi et al., 2010). According to Eq. (16), the effective galaxy bias can be estimated within the framework of the halo model as the large scale limit of , i.e. . The effective bias estimates reported in Orsi et al. (2010) have been performed adopting semianalitic models of galaxy formation (Baugh et al., 2005; Bower et al., 2006), and assuming a spectroscopic selection based on H emission. The latter is appropriate for our scope, since we are considering starforming galaxies only. The points reported in Figure 4 are the predictions performed by Orsi et al. (2010) using the Baugh et al. (2005) model for two different H flux thresholds ( and in CGS units). Obviously there is not much difference between the two choices, that are well reproduced by the solid line, representing the outcome of our halo model. In Figure 3 we also report the minimum mass of galaxy halos, compared to the average such mass found again by Orsi et al. (2010).
5 Results
In this Section we summarize our main results. First, we explored what happens to the power and cross spectra described above when the distribution of primordial density fluctuations deviates from a Gaussian. Then, we performed a Fisher matrix analysis in order to forecast the constraints on the main cosmological parameters that can be expected by a Euclidlike allsky survey.
5.1 Effect of nonGaussianity on the power and cross spectra
As a first step of our analysis we computed the Euclidexpected galaxy power spectrum, cluster power spectrum, and clustergalaxy cross spectrum for all the cosmological models considered in this work, and investigated the modifications induced by primordial nonGaussianity on each one. In Figure 5 we show the effect of the different shapes of primordial nonGaussianity examined in the present paper on the power and cross spectra for different tracers of the LSS. We show results for , since this roughly corresponds to the best high confidence limits that are currently placed on this parameter, and for as a comparison. We recall that the power spectrum of dark matter halos in cosmological simulations with local nonGaussian initial conditions has been studied in Grossi et al. (2009) and Desjacques, Seljak, & Iliev (2009). However, here we are interested in a wider variety of bispectrum shapes, and in specific classes of objects that are not necessarily in a onetoone relation with dark matter halos.
The corrections to the power and cross spectra that are due to different kinds of primordial nonGaussianity resemble the corrections to the large scale halo bias that have been shown in Figure 1. This is to be expected, since in the halo model the mass function is normalized over (see Section 4), hence the largest effect is expected to be due to the halo bias. In particular, for the local shape we observe a sharp increase of the power at large scales for all tracers. A similar increase in power at large scales, although not as strong, is also observed for the enfolded and orthogonal shapes. Note also that the trend for the latter is reversed, due to the fact that in this case a positive implies a negative effect on the bias and the mass function. Finally, for the equilateral shape we find a slight increase in the power at intermediate scales, again in agreement with the behavior of the halo bias.
As one might expect, the larger effect is found in correspondence of the tracers that are more biased with respect to the underlying matter density field, that are galaxy clusters. This can be understood by looking at Eq. (12), where it is shown that the difference between the nonGaussian halo bias and its Gaussian counterpart is in fact proportional to the square of the Gaussian bias itself. On the other hand, the effect on the clustergalaxy cross power spectrum is at an intermediate level between the effect on the galaxy power spectrum and that on the cluster power spectrum. Due to their low average bias, starforming galaxies are very poorly affected by primordial nonGaussianity, and this is particularly true for the enfolded and orthogonal bispectrum shapes.
A closer look to the first and third column panels of Figure 5 reveals an interesting fact. While the change in the clustergalaxy cross spectrum due to primordial nonGaussianity is larger than the change in the galaxy power spectrum, it is more so for the enfolded and orthogonal bispectrum shapes, as compared to the local and equilateral ones. This means that the inclusion of the crosscorrelation between clusters and galaxies should bring a better improvement over the constraints on nonGaussianity obtained by using the galaxy spectrum alone in the first two cases than in the second ones. Although it is difficult to gauge what the effect of the additional inclusion of the cluster power spectrum would be, we show in the next subsection that this expectation is to some extent satisfied.
5.2 Fisher matrix analysis
For the Fisher information matrix of the cross correlation between clusters and galaxies we followed the calculations of Hütsi & Lahav (2008), assuming perfect redshift knowledge for both galaxies and clusters. This assumption is justified since galaxies are spectroscopically selected, while clusters, being chosen as the highest S/N cosmic shear peaks, should be very massive, thus making followup confirmations relatively straightforward. In addition, we assumed redshift space distortions to be sufficiently well understood to be modeled away. Accordingly, the Fisher matrix for the redshift bin centered at can be well approximated as
(26) 
where , and
GALAXIES  GALAXY CLUSTERS  

Sky Coverage  sq. deg.  sq. deg. 
Redshift Efficiency  
H Flux Threshold  erg s cm  
Cosmic Shear S/N Threshold  
Redshift Range  
Minimum Scale  Mpc  Mpc 
(27) 
In the above set of equations, the quantities and are the average number densities of clusters and galaxies in the survey at hand (Eqs. 23 and 18, respectively), is the comoving volume contained in the unit redshift around , and is the th parameter of our cosmological model. The two wavenumbers and represent the boundaries of the wavenumber range used in the analysis. The chosen redshift bins should be relatively narrow, in order to treat the spatial number densities of tracers within each bin as constants.The total Fisher matrix is then given by the sum over the adopted redshift bins,
(28) 
In the case in which we are interested only in galaxies or only in clusters, the formalism remains the same, with the only difference being given by the replacements of with (with x c or x g) and of with
(29) 
in Eq. (26). In our analysis we adopted a redshift range typical for H galaxies selected by Euclid, (cf. the Euclid Yellow Book), and a range of scales corresponding to , where matches the largest scale available for a given redshift bin. For completeness, in Table 1 we summarize all the survey parameters that have been adopted in our Fisher matrix analysis including the redshift efficiency, that is the fraction of spectroscopically selected galaxies for which a reliable redshift measurement is expected.
LOCAL  EQUILATERAL  ENFOLDED  ORTHOGONAL  

In Figure 6 we show the joint confidence level contours in the plane that result from the Fisher matrix analysis performed for all the four bispectrum shapes considered in the present paper. All other cosmological parameters are have been kept fixed to the fiducial WMAP7 values summarized in Section 1. We report the results of considering the galaxy power spectrum alone, the cluster power spectrum, the clustergalaxy cross spectrum, and different combinations of the three. For further clarity, in Table 2 we report the errors that are forecasted both on the level of nonGaussianity and on the amplitude of the matter power spectrum . The fiducial model that we adopted has and , as specified in Section 1.
As can be seen, the tracers of the LSS that individually give the best constraints on both parameters considered are galaxies, due to their very large number density which beats down the shot noise in the Fisher matrix. Galaxies alone give constraints on that are more than an order of magnitude better than those of clusters alone, while the performance on the level of primordial nonGaussianity , as expected, depends on the specific choice of the primordial bispectrum shape. Galaxies perform a factor of better than clusters in the equilateral case, but only a factor of better in the orthogonal configuration.
The constraints on the plane that we gather from clusters alone are somewhat different than those reported in Sartoris et al. (2010), which performed a similar Fisher matrix analysis on the cluster sample expected to be detected by the planned Wide Field Xray Telescope (WFXT) (see also Cunha, Huterer, & Doré 2010 for a related work). Particularly, while the forecasted error on is comparable, the error on is much smaller in our case than in their. We ascribe this difference to two factors. First, clusters are selected in different ways in the two cases. Second, Sartoris et al. (2010) introduced nuisance parameters for estimating the uncertainty in the mass calibration, and at the same time adopted priors on several cosmological parameters inspired by Planck forecasts. We infer the second factor to be the dominant one, as rerunning the analysis of Sartoris et al. (2010) without nuisance parameters and priors brings to errors on comparable with ours own (B. Sartoris, private communication). We believe it to be safe to assume that Euclid weak lensing cluster catalogs will have sufficient control over systematics to ignore the contribution of imperfect knowledge of the massobservable relation. This point of view is shared by Bergé, Amara, & Réfrégier (2010).
The behavior of the confidence contours given by the clustergalaxy cross spectrum alone is interesting. Firstly, while the constraining power on is greatly enhanced as compared to the constraining power of clusters alone, the same is not true for . Actually the errors on expected from the clustergalaxy cross correlation are systematically larger than those expected from clusters alone. As a consequence, while galaxies still play a dominant role about constraining , the cross correlation gives comparable constraints on . As a matter of fact, for the enfolded and the orthogonal shapes, the cross spectrum performs even better with respect to , so that combination of all the three probes brings to a significant improvement over galaxies alone. Specifically, while the combined probes decrease the error on only by for the local bispectrum shape, this decrement reaches a factor of in the orthogonal case. The fact that the improvement in constraining power due to the inclusion of the cross spectrum is more pronounced for the enfolded and orthogonal shapes can be deduced from the observation that we made at the end of the last Subsection, justified by the complicated interplay of galaxy and cluster biases modified by primordial nonGaussianity.
A perhaps unexpected result is that the constraints on the parameter in the equilateral case are almost comparable with those for the local case, and much more stringent than for the enfolded and orthogonal shapes. Since the effect of nonGaussianity on the largescale bias is very slight in the equilateral case one would have naively expected the opposite behavior. However, it should be recalled that the correction to the mass function (which also enters the calculations of the power spectra and the Fisher matrices) for the equilateral shape is comparable to that for the local shape. Moreover, in the former case the nonGaussian correction to the bias is concentrated at intermediate scales, and the resulting effect on the power spectra is thus less degenerate with with respect to the other cases. Finally, the integral over wavenumbers in the definition of the Fisher matrix weights more smaller scales. Hence a modification of the power spectrum at a given level has a stronger impact if applied to small scales rather than large ones. In Fedeli, Moscardini, & Matarrese (2009) it is shown how the effect of equilateral primordial nonGaussianity on the mass function and halo bias is reduced when removing the running of . In the next Section we show that this has the effect of loosening the constraints on the level of primordial nonGaussianity.
LOCAL  EQUILATERAL  ENFOLDED  ORTHOGONAL  

0.05  0.02  0.05  0.09  
0.37  0.33  0.27  0.18  
0.13  0.13  0.08  0.04  
0.13  0.12  0.08  0.05  
0.04  0.02  0.03  0.04 
As a final step of our analysis, we computed the correlation coefficients between and for different nonGaussian shapes and for different combinations of LSS probes adopted. The correlation coefficient of the th and the th cosmological parameters in our set is defined as
(30) 
where is the Fisher matrix relative to the probe we are considering. For instance, the Fisher matrix for the combination of the cluster power spectrum with the clustergalaxy cross spectrum is . According to its definition, the correlation coefficient would have value for perfectly correlated parameters, for perfectly anticorrelated ones, and for uncorrelated parameters. The results are summarized in Table 3. The first thing to note is that, while the coefficients are positive for the local, equilateral, and enfolded shapes, they are negative for the orthogonal shape. This reflects the point that, in the latter case, the confidence contours plotted in Figure 6 are tilted in the opposite direction with respect to the former. This is a consequence of the fact that in the orthogonal case, a positive value of corresponds to a negative skewness of the density field and a negative correction to the linear bias, while for the other shapes considered in this work the opposite is true. A closer look to the absolute value of the correlation coefficients reveals that the constraints given by galaxy clusters alone are the most degenerate between and , while the constraints given by galaxies alone are almost undegenerate, in the sense that they pin down much better than .
The positive correlation between and in terms of power spectrum of cosmic tracers can be understood by considering the effect that both parameters have on the halo bias. An increase in has the effect of making structures less rare, and hence less biased with respect to the background matter density field. On the other hand an increase in in the positive direction increases the bias, except for the peculiar orthogonal case, that indeed shows an opposite degeneracy. The resulting degeneracy is however not very large because an increase in also increases the matter power spectrum, thus partially counteracting the decrement in the halo bias.
6 Discussion
In this Section we debate some of the assumptions that have been adopted in the present paper, and discuss overlappings with other works. We begin by noting that Wang et al. (2010) forecasted the constraints on the redshift evolution of the dark energy equation of state parameter expected by the galaxy power spectrum measured by Euclid. The error on that they found is of the order of , thus being much larger than our result. The reasons for this mismatch are the following: i) Wang et al. (2010) let all the cosmological parameters free to vary, while we fixed them all to their WMAP best fits, except for and ; ii) we used the full power spectrum information, while Wang et al. (2010) limited themselves only to the shape and position of the Baryon Acoustic Oscillation.
Recently, Smith, Desjacques, & Marian (2010) have shown by using numerical simulations of structure formation in the presence of local nonGaussianity, that dark matter halos tend to be on average more centrally concentrated than their counterparts in the reference Gaussian cosmology for positive , while the opposite is true for negative . This finding is in agreement with previous work (AvilaReese et al., 2003), and with the naive expectation that a higher efficiency in forming highmass objects implies. Under our assumption that the galaxy number density traces the dark matter density, this would underpin a different distribution of galaxies within dark matter halos as well. However, such an assumption is probably rough, and it is not clear if the galaxy distribution would be modified similarly to the matter distribution. Nevertheless, effects due to the inner structures of clusters show up only on very nonlinear scales, that are important for cosmic shear studies, but have been neglected in the present investigation (see Section 5). Thus, we stick to our choice of a standard NFW profile for the distribution of galaxies within halos for all cosmologies considered here.
The running of the parameter introduced for the equilateral shape in Eq. (4) also deserves further discussion. Namely, different authors tend to use different shapes for this running, or, in some circumstances, no running at all (Verde & Matarrese, 2009). Therefore, we redid our Fisher matrix analysis by setting the running . The results are shown in Figure 7. As can be seen, while the constraints on are almost unchanged those on are visibly loosened. For instance, the forecasted error obtained by galaxy clusters alone goes from to when the running is removed, while that obtained by galaxies alone grows from to . This kind of effect was expected, since the running we introduced had the effect of increasing the level of nonGaussianity on subCMB (cluster) scales. By setting such an increase is removed, and the errors increase. This result stresses the importance of considering the same running when results from different works are compared (see also Shandera, Dalal, & Huterer 2010).
In order to make our analysis more complete, we evaluated constraints on nonGaussianity considering survey configurations different from the fiducial one described in the preceding Sections. As a first step, we changed the S/N threshold for the detection of dark matter halos in cosmic shear maps, going from S/N to S/N . We again fitted the corresponding curve in Figure 1 of Bergé, Amara, & Réfrégier (2010) in order to infer the minimum cluster mass to be adopted. In this case, at is a factor of larger than for the S/N case. Only two bispectrum shapes have been considered here: the local one, since it is the one expected to give the best constraints, and the orthogonal one, since it is the less studied one. In Figure 8 we show the joint confidence levels on and given by this new detection criterion for galaxy clusters, and in Table 4 we report the numerical values of the respective errors. While the constraints given by galaxies alone are obviously unchanged, those given by clusters and the clustergalaxy crosscorrelation are significantly loosened. As a consequence, very little information is added if combining the galaxy power spectrum with power spectra involving clusters, and this remains true for the orthogonal model as well. This result implies that, while selecting cosmic shear peaks with higher S/N we are effectively selecting more massive structures, on whose statistics the effect of nonGaussianity is larger, the respective constraining power is much too suppressed by the exponential decrement in the abundance of such objects.
As a second alternative survey configuration we considered a multislit spectroscopic instrument on the model of the Digital Micromirror Devices (DMDs) conceived for the SPACE mission (Cimatti et al., 2009). This is particularly relevant in view of future space missions alternative to Euclid, such as WFIRST and JEDI (Crotts et al., 2005). In this case selection of galaxies would occur according to their Hband flux. As a consequence Luminous Red Galaxies (LRGs) would be detected as well, thus we had to use the halo occupation distribution relative to a mixture of red and blue galaxies, rather than to blue galaxies only as we did for our fiducial survey. We found that a new constant minimum halo mass of reproduces fairly well the effective bias of galaxies with Hband magnitude brighter than H according to the simulations of Orsi et al. (2010). All the other parameters of the survey are left unchanged, while the redshift range of the analysis is extended to (for both clusters and galaxies, although the former are unaffected by spectroscopy) in order to comply with the different observational specifications. In Table 4 we report the errors forecasted for this alternative configuration. As can be seen, in this case the constraints on both and are significantly improved with respect to our fiducial Euclidlike case. The error on the level of nonGaussianity can for instance reach for the local bispectrum shape. Since the cluster selection is unchanged, we interpret this improvement as due to a combination of the different galaxy selection and of the extended redshift range that is allowed by the multislit configuration. We expand a little bit more on this further below. This improvement is additionally emphasized in Figure 9, where we show the parameter 68% confidence levels on the plane given by the combination of the power spectra of clusters and galaxies and the clustergalaxy cross spectrum for both the multislit configuration we are considering and the fiducial slitless Euclid case. It is interesting to note that, due to the different galaxy population selected in the multislit case, the inclination of the confidence ellipses changes with respect to the fiducial slitless configuration.
LOCAL (S/N )  ORTHOGONAL (S/N )  LOCAL (multislit)  ORTHOGONAL (multislit)  

Finally, we estimated how the constraints on the level of nonGaussianity changed upon modification of the redshift and scale ranges considered in the Fisher matrix analysis, for both the slitless Euclidlike configuration and the multislit case. In Figure 10 we show the results of changing either the minimum scale included in the analysis, , the minimum, and the maximum redshift of both galaxies and clusters. Fiducial vales are Mpc, for Euclid, and for the multislit case. As one could naively expect, the error on decreases by increasing both the minimum scale and the maximum redshift (although in the latter case the trend is quite mild), while it increases by increasing the minimum redshift (because the maximum redshift is held fixed). As mentioned above, the differences between the slitless and the multislit case are due to a combination of the different selection of galaxies with the different redshift range adopted. An exception to this is given by the two middle panels of Figure 10, where the redshift ranges adopted are the same for both configurations, and thus the differences can be ascribed only to the different selection functions.
7 Summary and conclusions
We adopted the physically motivated halo model in order to compute the effect of different kinds of primordial nonGaussianity on the power spectrum of galaxies and galaxy clusters, as well as on the clustergalaxy cross spectrum. Specifically, we considered galaxies selected spectroscopically according to their H flux, and galaxy clusters selected as the largest S/N peaks in cosmic shear maps, having in mind future wide field optical/nearinfrared surveys such as Euclid and WFIRST. We additionally performed a Fisher matrix analysis in order to forecast the errors on the joint determination of the level of nonGaussianity and the amplitude of the matter power spectrum . The main findings of this work can be summarized as follows.

The nonGaussian corrections to the power spectrum of tracers of the LSS resembles the modifications to the largescale bias of dark matter halos, as one might naively expect. The largest effect is seen for the local bispectrum shape, while the smallest appears for the enfolded shape. The power spectrum of massive galaxy clusters is more heavily modified than the spectrum of galaxies, because the former are originally more biased than the latter. The modification to the clustergalaxy cross spectrum lies somewhere in between.

Galaxies have a much higher constraining power on both and as compared to galaxy clusters, due to the much lower abundance of the latter that is not adequately compensated by the larger effect on the relative power spectrum. Assuming a Euclidlike survey, while spectroscopically selected galaxies can constrain to the level of a few and to the level of , errors on parameters derived by clusters alone are at the level of and (with some dependence on the primordial bispectrum shape).

When considering the clustergalaxy cross spectrum alone, the forecasted constraints on are comparable to the clusteronly constraints, while the constraints on are improved by more than one order of magnitude, reaching a predicted error comparable with that coming from galaxies alone, a few. This result highlights the high complementarity of the cluster power spectrum and the clustergalaxy cross spectrum as cosmological probes.

While the constraints on coming from the galaxy power spectrum alone are almost unchanged when combined with the clustergalaxy cross correlation, the constraints on are improved. This is true only to a slight level for the local and equilateral bispectrum shapes, however the error on can be reduced by a factor of for the enfolded and orthogonal cases. The addition of the cluster power spectrum carries only slight change to this conclusion.

As expected the parameters and are degenerate with respect to the power spectra of LSS tracers. This degeneracy is more marked for the cluster power spectrum, with correlation coefficients reaching up to . The degeneracy on the other hand is almost insignificant for the galaxy power spectrum, in the sense that the constraining power weights much more on than on . The degeneracy with respect to the clustergalaxy cross spectrum lies in between the former two, with the exception of the orthogonal bispectrum shape.

We considered several survey configurations alternative to the fiducial Euclidlike one, finding that a multislit spectroscopic instrument would allow more stringent constraints both on and . This improvement is due to a combination of the fact that with this other arrangement the selected galaxy population is different, and that the data analysis can be pushed down to a lower minimum redshift. This result is interesting with respect to future wide field survey concepts alternative to Euclid, such as WFIRST.
The improvement in the constraints on both and when combining the cluster and/or galaxy power spectra with the clustergalaxy cross spectrum is reminescent of the fact that the latter is sensitive to different scales in a different way with respect to the former. The fact that this improvement is more important for the enfolded and orthogonal bispectrum shapes is extremely interesting, since these models are still relatively poorly studied compared to the equilateral and, almost ubiquitous, local shapes. Additionally, the amazing constraining power of the galaxy power spectrum, even when considered alone, stresses the importance of the spectroscopy part for future Euclidlike missions when it comes to restrict cosmology on the basis of the spatial distribution of objects.
The present work reinforces the conclusion that the spatial distribution of tracers of the LSS, especially galaxies, is an incredibly powerful probe for primordial nonGaussianity, thanks to the very strong impact that the latter has on the large scale bias of dark matter halos (Carbone, Verde, & Matarrese, 2008, 2010). The combination of cluster and/or galaxy power spectra with the cross spectrum of clusters and galaxies significantly improves the forecasted constraints, especially for the least studied nonGaussian shapes. Merging all the information from future wide field surveys such as Euclid and WFIRST promise to bring constraints on to the unity level, and constraints on to the level of , in both cases superior to future CMB experiments (Sefusatti et al., 2009; Verde & Matarrese, 2009).
Acknowledgments
We acknowledge financial contributions from contracts ASIINAF I/023/05/0, ASIINAF I/088/06/0, ASI I/016/07/0 ’COFIS’, ASI ’EuclidDUNE’ I/064/08/0, ASIUni BolognaAstronomy Dept. ’EuclidNIS’ I/039/10/0, and PRIN MIUR ’Dark energy and cosmology with large galaxy surveys’. We are grateful to N. Bartolo, A. Gonzalez and B. Sartoris for stimulating discussions about the manuscript, and to A. Orsi for sharing the results of his simulations. We wish to thank the anonymous referee for useful comments that allowed us to improve the presentation of our work.