Cosmological Parameters Degeneracies and Non-Gaussian Halo Bias

We study the impact of the cosmological parameters uncertainties on the measurements of primordial non-Gaussianity through the large-scale non-Gaussian halo bias effect. While this is not expected to be an issue for the standard LCDM model, it may not be the case for more general models that modify the large-scale shape of the power spectrum. We consider the so-called local non-Gaussianity model and forecasts from planned surveys, alone and combined with a Planck CMB prior. In particular, we consider EUCLID- and LSST-like surveys and forecast the correlations among $f_{\rm NL}$ and the running of the spectral index $\alpha_s$, the dark energy equation of state $w$, the effective sound speed of dark energy perturbations $c^2_s$, the total mass of massive neutrinos $M_\nu=\sum m_\nu$, and the number of extra relativistic degrees of freedom $N_\nu^{rel}$. Neglecting CMB information on $f_{\rm NL}$ and scales $k>0.03 h$/Mpc, we find that, if $N_\nu^{\rm rel}$ is assumed to be known, the uncertainty on cosmological parameters increases the error on $f_{\rm NL}$ by 10 to 30% depending on the survey. Thus the $f_{\rm NL}$ constraint is remarkable robust to cosmological model uncertainties. On the other hand, if $N_\nu^{\rm rel}$ is simultaneously constrained from the data, the $f_{\rm NL}$ error increases by $\sim 80%$. Finally, future surveys which provide a large sample of galaxies or galaxy clusters over a volume comparable to the Hubble volume can measure primordial non-Gaussianity of the local form with a marginalized 1--$\sigma$ error of the order $\Delta f_{\rm NL} \sim 2-5$, after combination with CMB priors for the remaining cosmological parameters. These results are competitive with CMB bispectrum constraints achievable with an ideal CMB experiment.


Introduction
Tests of deviations from Gaussian initial conditions offer an important window into the very early Universe and a powerful test for the mechanism which generated primordial perturbations. While standard single-field slow-roll models of inflation lead to small departures from Gaussianity, non-standard scenarios allow for a larger non-Gaussianity (NG) level (e.g. [1,2,3], and refs. therein).
In particular, large NG can be produced if any of the conditions below is violated: a) single field, b) canonical kinetic energy c) slow roll and d) adiabatic (Bunch-Davies) initial vacuum state. The type of NG arising in standard inflation reads [4,5,6,7] where Φ denotes Bardeen's gauge-invariant potential, which, on sub-Hubble scales reduces to the usual Newtonian peculiar gravitational potential up to a minus sign, and φ denotes a Gaussian random field. The NG parameter f NL is often considered to be a constant, yielding NG of the local type with a bispectrum which is maximized for squeezed configurations [8]. NG of the local type is generated in standard inflationary scenarios (where f NL is expected to be of the same order of the slow-roll parameters) as well as in multi-field inflationary scenarios 1 .
The standard observables to constrain NG are the Cosmic Microwave Background (CMB) and the Large-Scale Structure (LSS) of the Universe. Traditionally, the most popular method to detect primordial NG has been to measure the bispectrum or the threepoint correlation function of the CMB [6,9,10], while the LSS bispectrum has been shown to be sensitive to primordial NG only at high redshift [6,11,12,13,14].
Refs. [25] and [26] (hereafter MV08) showed that primordial NG affects the clustering of dark matter halos inducing a scale-dependent large-scale bias. This effect, which goes under the name of non-Gaussian halo bias, is particularly promising, yielding already stringent constraints from existing data [27,17]. Forthcoming constraints on NG exploiting the non-Gaussian halo bias are expected to be similar to those achievable from an ideal CMB survey [16]. These predictions have been confirmed by N-body simulations [25,28,29,30].
Forecasts for f NL constraints from the halo-bias have been carried out so far assuming perfect knowledge of all other cosmological parameters. While for a ΛCDM model this is expected to be a reasonable assumption, for more general models one may expect f NL to be degenerate with other parameters and thus to have a larger marginal error. Here we study the degeneracies among the large-scale non-Gaussian halo bias (for NG of the local type) and the cosmological parameters which affect the large-scale halo power spectrum, focusing on dark energy perturbations, massive neutrinos, number of relativistic species, and running spectral index, which can produce large deviations of the underlying cosmology from the minimal ΛCDM scenario. The paper is organized as follows. In §2 we briefly review the analytic expressions of the non-Gaussian halo power spectrum generalized to redshift dependent transfer functions. The redshift dependence is due to the presence of both dark energy perturbations and massive neutrinos. In §3 we summarize the Fisher matrix formalism applied to the observed halo power spectrum. In §4 we describe the assumed fiducial cosmology. Finally, in §5 and §6 we discuss the results and draw our conclusions.

Non-Gaussian halo bias
Here we summarize the derivation of [26], extending it to the case of a redshift dependent transfer function. In Fourier space, the filtered linear over-density δ R is related to the primordial potential Φ(k) by the Poisson equation: where T (k, z) denotes the matter transfer function (which is redshift dependent in the presence of massive neutrinos and/or dark energy perturbations), W R (k) is the Fourier transform of the top-hat function of width R, H 0 and Ω m,0 are the current values of the Hubble constant and the total matter energy density respectively, and D(z) = (1 + z) −1 g(z)/g(0) is the linear growth-factor of density fluctuations, normalized to D(0) = 1 with g(z) being the growth suppression factor for non Einstein-de Sitter universes. In this context, the non-Gaussian halo power-spectrum takes the form is the Gaussian lagrangian bias 2 for dark matter halos of mass M [34], σ 2 M is the mass variance linearly extrapolated to z = 0, δ c (z) ≡ ∆ c (z)/D(z), being ∆ c (z) the linear overdensity for spherical collapse, which can be considered as a constant ∆ c (z) = 1.686, even in the presence of dark energy [35], q is a factor extracted from N-body simulations ( [28] and references therein) and P δδ (k, z) is the power spectrum of δ R (as it will be shown later in Eq. (4.8)). Finally, β(k, z) is defined as where α = k 2 1 + k 2 + 2k 1 kµ. Here we adopt the so-called CMB f NL normalization where Eq. (1.1) is intended to be deep in the matter-dominated era 3 . Consequently, the non-Gaussian halo Lagrangian bias reads Making the standard assumption that halos move coherently with dark matter, the Eulerian bias is b E = 1 + b L . The halo power spectrum given by Eq. (2.2), is connected directly to the underlying dark matter power spectrum and can be recostructed from the galaxy power spectrum using different techniques (e.g. [39]). It provides important information on the growth of structure, which helps in constraining dark energy and neutrino masses together with primordial non-Gaussianities. We shall exploit information from both the shape and the amplitude of the NG halo power spectrum up to scales k < 0.03h/Mpc, where β(k, z) has a negligible dependence on the halo mass (recall that for local non-Gaussianity β(k, z) ∝ 1/k 2 ). In addition, Eqs. (2.2)-(2.4) are valid only on scales much larger than the Lagrangian radius of the halo [26,16] and have been tested and calibrated on N-body simulations only on scales k < 0.03h/Mpc [28]. We believe that the approach followed here is conservative, since in our parameter forecasts we do not consider scales k > 0.03h/Mpc in the halo power spectrum, which could provide much tighter constraints on cosmology, via e.g. Baryonic Acoustic Oscillations (BAO). In what follows, we shall use Eq. (2.2) to propagate errors of the non-Gaussian halo power spectrum into errors of the cosmological parameters via a Fisher matrix approach, as described below.
2 In general, the Gaussian halo bias may have a non trivial dependence on both the halo formation redshift z f and the observation redshift zo [31,32,33]. However, for objects that did not undergo recent mergers, z f ≫ zo, or in the case of rapid mergers z f ≈ zo for δ 2 c ≫ σ 2 M , i.e. large M and/or high z f , the bias is well approximated by Eq. (2.2).

Methodology
In this paper we adopt the Fisher matrix formalism to make predictions of the cosmological parameter errors including the NG parameter f NL . The Fisher matrix is defined as the second derivative of the likelihood surface about the maximum. As long as the posterior distribution for the parameters is well approximated by a multivariate Gaussian function, its elements are given by [40,41,42] where C = S + N is the total covariance which consists of signal S and noise N terms. The commas in Eq. (3.1) denote derivatives with respect to the cosmological parameters within the assumed fiducial cosmology 5 .
To derive realistic parameter forecasts, we consider future redshift surveys, as Euclid 6and LSST 7 -like galaxy surveys. While for BAO surveys an accurate redshift measurement is crucial [43,44], for our purposes, a precise redshift extraction is not needed, as it will be explained after Eq. (3.3). For the EUCLID-like survey we will assume a redshift coverage of 0.5 < z < 2 and a sky area of f sky = 20000 deg 2 . For the LSST-like survey, we will consider 0.3 < z < 3.6 and f sky = 30000 deg 2 for redshift and area coverages, respectively. In order to explore the cosmological parameters constraints from a given redshift survey, it is mandatory to specify the measurement uncertainties of the halo power spectrum. In general, the statistical error on the measurement of P h (k) at a given wavenumber bin is given by [45] ∆P wheren h is the mean number density of dark matter halos, V is the comoving survey volume of the galaxy survey, and µ is the cosine of the angle between k and the line-ofsight direction. To our purposes it is adequate to perform an angular average over µ. Thus, our Fisher matrix for the large-scale structure data is given by where p α represents the chosen set of cosmological parameters. We divide the surveys in redshift bins of width ∆z = 0.1 (larger than standard photometric and spectroscopic redshift errors), and set k max to be 0.03h/Mpc and k min to be greater than 2π/∆V 1/3 , where ∆V is the volume of the redshift shell. Conservatively, we do not consider here that scales larger than k min can be used by cross-correlating different shells.
The effect of NG alters the broad-band behavior of P h (k) on very large scales, where P h (k) is unaffected by the precision with which the radial positions of the galaxies are measured. Thus, we can treat photometric and spectroscopic surveys on the same footing. Moreover, the requirement of surveying a large volume of the universe and sampling highly biased galaxies to beat shot-noise, which is a key point for BAO surveys, is also a bonus for constraining primordial NG [16].
Note, moreover, that the NG correction of the halo bias is boosted by the Lagrangian Gaussian halo bias factor itself.
In particular, for the value of k min used here, we find that P h (k min ) ≃ P h (k = 0.2h/Mpc), thus, the shot-noise requirement for BAO surveys ofnP (k = 0.2h/Mpc) > 1, implies that for all scales of interest here,nP ≫ 1. We have checked that our results do not change if we imposenP ∼ 3 at all scales.
We compute as well the CMB Fisher matrix to obtain forecasts for the Planck satellite 8 . We follow here the method of [46], considering the likelihood function for a realistic experiment with partial sky coverage, and noisy data and compute its second derivatives to obtain the corresponding Fisher matrix being C XY ℓ the temperature and polarization power spectra (X, Y ≡ {T, E, B}) and N ℓ the noise bias. Finally, f XY sky is the fraction of observed sky which can be different for the T -, E-, and B-modes. In Eq. (3.5) L ≡ ln L, p α and p β denote the cosmological parameters of the assumed model and form the vector p whose fiducial value is given byp.
Combining the Planck and redshift survey Fisher matrices (F αβ = F LSS αβ +F CMB αβ ) we get the joint constraints for p. The 1-σ error on p α marginalized over the other parameters is σ(p α ) = (F −1 ) αα , being F −1 the inverse of the Fisher matrix. We then consider joint constraints in a two-parameter subspace (marginalized over all other cosmological parameters) to study the covariance between f NL and the other cosmological parameters considered in this work.
Furthermore, in order to quantify the level of degeneracy between different parameters and f NL , we estimate the so-called correlation coefficients, given by where p α denotes one of the model parameters. When the coefficient |r| = 1, the two parameters are totally degenerate, while r = 0 means they are uncorrelated.
Here Ω m,0 and Ω b,0 are the total matter and baryon present-day energy densities, respectively, in units of the critical energy density of the Universe, h is given by H 0 = 100h km s −1 Mpc −1 , where H 0 is the Hubble constant, ∆ 2 R (k 0 ) represents the dimensionless amplitude of the primordial curvature perturbations evaluated at a pivot scale k 0 , n s is the scalar spectral index of the primordial matter power spectrum, assumed to be a power-law, and α s is the running of the scalar spectral index. We will consider two different fiducial models matching this ΛCDM cosmology as follows. We adopt the same values for the 7 "base" parameters. We do not consider primordial gravitational waves and assume a flatness prior, Ω K = 0, as predicted by longlasting inflation models, so that Ω de,0 = 1 − Ω m,0 , where Ω K and Ω de,0 are, respectively,  A Gaussian prior of 5% on the present-day Hubble's constant H 0 = 100h km s −1 Mpc −1 is assumed, following the results of [52]. While this uncertainty is comparable to the one achieved by recent WMAP-7yr data 9 in the determination of H 0 for the ΛCDM model [51], this information will improve, as expected, the parameter constraints on models different from the minimal ΛCDM model, such as models with dark energy perturbations, massive neutrinos and non-vanishing running of the spectral index.
We also consider dark energy to be a cosmic fluid with clustering properties on the Gpc scale, described by an equation of state that we assume to be constant where where p de and ρ de represent respectively the pressure and energy density of the dark energy fluid, and we assume a fiducial value w| fid = −0.9 which lies well within the current 95% C.L. limits on a constant dark energy equation of state parameter w. The dark energy fluid will be also described here by an effective sound speed c s which parametrizes the transition between the smooth and clustered dark-energy regimes [47] whit fiducial value c 2 s | fid = 0.9. Dark energy perturbations will arise only if the dark energy equation of state parameter w is different from −1.
We assume the power spectrum of primordial curvature perturbations, P R (k), to be .

(4.3)
where k 0 = 0.002/Mpc and ∆ 2 R (k 0 )| fid = 2.64 × 10 −9 [48]. The matter energy density Ω m,0 includes the neutrino contribution when neutrinos are non-relativistic Ω m,0 = Ω c,0 + Ω b,0 + Ω ν,0 , where Ω ν,0 is related to the sum of neutrino masses M ν ≡ m ν as and the neutrino mass eigenstates are assumed to have a degenerate spectrum, i.e. the three neutrinos have the same mass. We will constrain the following set of "baseline" parameters From here we specify two cosmological models: 9 http://lambda.gsfc.nasa.gov/ • N rel ν -cosmology, where neutrinos are effectively massless but the the number of relativistic species N rel ν can deviate from the standard value N rel ν = 3.04. In this case the fiducial value N rel ν | fid = 3.04 is chosen, fixing M ν = const = 0. N rel ν is given by the energy density associated to radiation where Ω γ,0 = 2.469 × 10 −5 h −2 is the present-day photon energy density parameter for T cmb = 2.725 K [49].
• M ν -cosmology, where N rel ν is fixed at the fiducial value and M ν is allowed to vary. In this case, being still consistent with current data [50,51] we choose a fiducial value M ν | fid = 0.3 eV. This choice is motivated by the fact that for taking two-sided numerical derivatives, the fiducial M ν must be non-zero. It is well known that the error from cosmological observations on the neutrino mass depends somewhat on the fiducial mass chosen; from [53] we estimate that around a fiducial M ν = 0 the error on M ν = 0 would increase by less than 20%. As it will be clear from Sec. 5 the effect of this correction will be negligible on the f NL error estimate. In addition most of the signal to constrain M ν from LSS surveys will come from smaller scales (not considered here).
At the CMB level, if neutrinos are still relativistic at the decoupling epoch, z ≃ 1090, i.e. if the mass of the heaviest neutrino specie is m ν < 0.58 eV, massive neutrinos do not affect the CMB power spectra, except through the gravitational lensing effect [49], and, as a consequence, the dark energy equation of state w is not degenerate with the neutrino mass. However, the limit on the the sum of the neutrino masses degrades significantly when the dark energy equation of state is a function of redshift, if the amplitude of the galaxy spectrum is used for getting constraints on w and M ν , since dark energy and massive neutrinos both affect the growth rate of structures.
In this work, both massive neutrinos and clustering properties of the dark energy perturbations are considered. In this scenario, the growth function of the dark matter perturbations is scale-dependent, even at the linear level. The overall effect induces a redshift-dependent transfer function [54,55,56], and the power spectrum of the linear density field, smoothed on a sphere of radius R, takes the form where D(z) is the scale independent linear growth-factor defined in §2. The redshiftdependent transfer function T (k, z) is directly extracted from CAMB 10 [57] at each redshift z, in order to compute the Fisher matrix, given by Eq. (3.3), within each redshift bin. Our analysis exploit exclusively the linear matter power spectrum, since we restrict ourselves to scales k 0.03h/Mpc, where the details of the halo occupation distribution of galaxies are irrelevant. We focus here on dark matter halos with mass ∼ 10 12 -10 13 (where the lower mass limit is relative to the highest redshift of the survey).
Notice as well that in our Fisher matrix analysis we do not add constraints on the f NL parameter from CMB experiments, so that our forecasts on f NL result exclusively from future redshift survey measurements of the dark matter halo power spectrum on scales k 0.03h/Mpc, i.e. without including information from BAOs which will further reduce the forecasted errors and residual degeneracies. For all the reasons explained above, the results presented in the next section should be considered conservative.

Results
In this Section we present the predicted 1-σ marginalized errors of the f NL parameter, and the f NL covariance with the remaining cosmological parameters considered in our Fisher matrix analysis. We show forecasts both from LSST and EUCLID data only, as well as the expected errors after combining the results from these two experiments with Planck forecasted errors.
First of all, when considering forecasts from the redshift surveys alone, we expect that the f NL parameter will be correlated with all cosmological parameters which affect the amplitude and shape of P h (k) at scales k 0.03h/Mpc. Tables 1-2 show the correlation r (see Eq. (3.6)) among f NL and the other cosmological parameters p α . We expect w and f NL to be correlated. In fact, at scales k 0.03h/Mpc, an increase of w produces a trend on the matter power spectrum which is opposite to the one produced in P h (k) by increasing f NL . So the effect of positively increasing f NL can be mimicked by a larger w (see Fig. 1).
Likewise, f NL is correlated with Ω b,0 h 2 , Ω c,0 h 2 , and M ν , since, on the scales considered, the larger these parameters are, the smaller the matter power spectrum is, and this competes with the rise of P h (k) due to an increasing value of f NL .
On the other hand, f NL is negatively correlated with both the running of the scalar spectral index α s and the effective sound speed of dark energy perturbations c s . In fact, if either α s or c s increase, the matter power spectrum P δδ is modified in a very similar way to the halo power spectrum P h (k) for a larger f NL value. Consequently, the effect of a positively increasing f NL can be mimicked by decreasing either α s or c s . For what concerns the effective number of relativistic species, the f NL -N rel ν correlation is more complicated, as it depends on several factors. Naively, one would expect a positively correlation between these two parameters, since an increase in the number of relativistic particles should suppress the matter power spectrum, which can be compensated by increasing f NL . However, since the redshift evolution of the halo bias is also modified in this situation, the correlation coefficient may change its sign.
When the Planck Fisher matrix information is added to the survey Fisher matrix, all degeneracies are either resolved or drastically reduced. In some cases, the correlation coefficient r can even change sign, see Tabs. 1-2. This change in the behavior of r arises either due to the presence of dominant parameter degeneracies affecting the CMB spectrum, or because of marginalization of a high-dimension parameter space down to two variables. In particular, it is worth noting here that, while f NL and N rel ν are negatively correlated if only galaxy survey data are considered, they are positively correlated after adding Planck priors. In fact, since both N rel ν and Ω c,0 h 2 are strongly, positively correlated at the CMB level via the equality redshift z eq , then, a positive correlation between Ω c,0 h 2 and f NL automatically turns into a correlation of the same sign between N rel ν and f NL . Despite the residual non-zero correlation coefficients, one should bear in mind that the marginalized f NL errors decrease by a factor > 2 when the CMB prior is added.
From our analysis, we conclude that the effective number of relativistic species is the main parameter affecting the constraints on f NL . The 1-σ errors of f NL , for the two fiducial cosmologies considered here, are shown in Tabs. 3 and 4. Let us notice that the marginalized errors are significantly larger than the non-marginalized ones when only the LSS surveys are used for the forecasts. Nonetheless, the marginalized errors become comparable in magnitude to the non-marginalized errors when Planck priors are added, since the CMB mitigates the intrinsic degeneracies between f NL and the other cosmological parameters at the LSS level. It is also worth noting that the non-marginalized 1-σ errors of f NL presented in this work are larger than the corresponding errors presented in [16]. In fact we now consider less highly biased halos, with a fiducial bias parameter more in-line with the expected one for (blue) EUCLID galaxies [58,59]. It is important to note that the f NL effect on the halo bias is modulated by b L,h (z, M ), which depends crucially on the selected halo and its merging history [60,61].
Tabs. 3-4 show as well the effect of each cosmological parameter on the f NL forecasts. Obviously the parameters which have a larger impact on the f NL errors are the ones more degenerated with it, and can be directly inferred from Tabs. 1-2. Moreover, since with the inclusion of the parameters α s , c s , M ν and N rel ν , we have considered cosmologies which deviate substantially from the minimal ΛCDM model, we fix pairs/triplets composed by these parameters to show how much deviations from a ΛCDM cosmology can affect the f NL constraints. For the N rel ν model cosmology there is an important impact on the f NL marginalized errors from Ω c,0 h 2 , h, w and, in particular, from N rel ν . In summary, if N rel ν is assumed to be fixed, the uncertainties on the other cosmological parameters increase the error on f NL only by 10 to 30%, depending on the survey. If N rel ν is considered as an extra parameter to be simultaneously constrained from the data then the uncertainty in the underlying cosmology increases the f NL error by ∼ 80%.
In Figs. 2-3 we show the 2-parameter projected 68% C.L., 95.4% C.L. and 99.73% C.L. contours in the f NL -p α sub-space with p α = w, c 2 s , α s , M ν , N rel ν , obtained after combining LSST and EUCLID data with Planck priors for the two fiducial models considered in this work. The black line shows the 1-parameter confidence level at 1-σ. The orientation of the ellipses reflects the correlations among the parameters shown in Tabs. 1-2.

Conclusions
Deviations from non-Gaussianity, usually parameterized by the parameter f NL , offer a powerful tool to identify the mechanism which generates the seeds for the structures we observe currently in our Universe.
Here we study the impact of the uncertainties of the cosmological parameters on the f NL errors expected for the case of local non-Gaussianity for the large-scale non-Gaussian Figure 3: The same as in Fig. 2, for the fiducial model with massive neutrinos of total mass M ν | fid = 0.3 eV, as described in the text. halo bias effect. We forecast the correlations among f NL and the remaining cosmological parameters (including the running of the spectral index α s , and the dark energy parameters w and c 2 s ) within two possible cosmological models. The first model contains massive neutrinos (hypothesis robustly confirmed by neutrino oscillation data), where the total neutrino mass is a parameter to be constrained by the cosmological data. The second model assumes massless neutrinos (or neutrinos with a mass too small to be relevant for the cosmological observations considered here) and allows for extra relativistic degrees of freedom N rel ν , which could be induced by the presence of sterile neutrinos, non minimally coupled quintessence fields, or even by the violation of the spin statistics theorem in the neutrino sector.
We follow here a conservative approach, assuming that f NL is constrained exclusively from the very large scale halo power spectrum (i.e. we neglect CMB information on f NL ), and restrict ourselves to scales k 0.03h/Mpc, without exploiting information e.g., from BAOs, which will further reduce degeneracies and forecasted errors. We present first the Fisher matrix forecasts for f NL assuming EUCLID-and LSST-like surveys for the two model cosmologies considered here. Then, we add the Planck Fisher forecasts for the remaining cosmological parameters to study the impact on the f NL correlations.
The combined errors on f NL do not change significantly in the presence of a dark energy equation of state, massive neutrinos, running of the spectral index, or clustering of dark energy perturbations, which are the parameters we have particularly focused on, since they are expected to affect the matter power spectrum on large scales, and represent the main deviations from a minimal ΛCDM model. However, the errors on f NL are highly affected in the presence of extra relativistic degrees of freedom N rel ν . We find that if N rel ν is assumed to be fixed, the effect of the uncertainties on the other cosmological parameters increases the error on f NL only by 10 to 30% depending on the survey. If N rel ν is considered as a parameter to be simultaneously constrained from the data, then the uncertainty in the underlying cosmology increases the f NL error by ∼ 80%. We thus conclude that, except for the effect of N rel ν , the halo-bias f NL constraints are remarkably robust to uncertainties in the underlying cosmology.
One important point to discuss is the effect of the (Gaussian) halo bias, as its value boosts the effect of f NL on the halo power spectrum shape, and in our analysis it has been assumed to be known. The (Gaussian) halo bias depends strongly on the type of halos selected by the survey -whether they correspond to extremely high and rare peaks in the initial fluctuation field-, and on their accretion history. Errors on f NL may be improved -at least in principle-by up to a factor of two by optimizing the choice of tracers.
The bias factor itself will need to be estimated from the survey, at the same time as the other cosmological parameters; the signal comes from scales much smaller than those used here, where the NG effect on halo bias is completely negligible. We estimate that the error on the (Gaussian) halo bias will be of the same order (in %) as the error on the linear growth factor f as a function of redshift, which is forecasted to be <∼ 10% [62,63]. Such a residual uncertainly will therefore increase the f NL errors reported here by at most 10%.
Let us recall that the purpose of this work is to show the main correlations between f NL and the other cosmological parameters, and to understand if these degeneracies can degrade dramatically the f NL errors. We have shown that, after the combination with Planck constraints on parameters different from f NL , the degeneracies get mostly broken, independently on the particular cosmological parameter, even without adding information from smaller scales corresponding to k > 0.03h/Mpc. Therefore we conclude that the f NL constraints are very robust against underlying cosmology assumptions.
Finally, future surveys which provide a large sample of galaxies or galaxy clusters over a volume comparable to the Hubble volume (LSST, EUCLID) will measure primordial non-Gaussianity of the local form with a marginalized 1-σ error of the order ∆f NL ∼ 2 − 5, after combination with CMB priors for the remaining cosmological parameters. These results are competitive with CMB bispectrum constraints achievable with an ideal CMB experiment ∆f NL ∼few [64,65].