A numerical study of the effects of primordial nonGaussianities on weak lensing statistics
Abstract
While usually cosmological initial conditions are assumed to be Gaussian, inflationary theories can predict a certain amount of primordial nonGaussianity which can have an impact on the statistical properties of the lensing observables. In order to evaluate this effect, we build a large set of realistic maps of different lensing quantities starting from lightcones extracted from large darkmatter only Nbody simulations with initial conditions corresponding to different levels of primordial local nonGaussianity strength . Considering various statistical quantities (PDF, power spectrum, shear in aperture, skewness and bispectrum) we find that the effect produced by the presence of primordial nonGaussianity is relatively small, being of the order of few per cent for values of compatible with the present CMB constraints and reaching at most 1015 per cent for the most extreme cases with . We also discuss the degeneracy of this effect with the uncertainties due to the power spectrum normalization and matter density parameter , finding that an error in the determination of () of about 3 (10) per cent gives differences comparable with nonGaussian models having . These results suggest that the possible presence of an amount of primordial nonGaussianity corresponding to is not hampering a robust determination of the main cosmological parameters in present and future weak lensing surveys, while a positive detection of deviations from the Gaussian hypothesis is possible only breaking the degeneracy with other cosmological parameters and using data from deep surveys covering a large fraction of the sky.
keywords:
cosmology: theory  gravitational lensing: weak  cosmological parameters  largescale structure of the Universe  Methods: Nbody simulations1 Introduction
In recent years, the interest for an accurate measurement of the amount of nonGaussianity present in the primordial density field has largely increased. The main reason is that this test is now considered not only a general probe of the inflationary paradigm, but also a powerful tool to constrain the plethora of its different variants. Only the most standard slowrolling models based on a single field produce in fact almost uncorrelated fluctuations, which is the motivation of the common assumption (and large simplification) that their distribution is Gaussian. In general, small deviations from Gaussianity are predicted even for the simplest inflationary models, while nonstandard models, like the scenarios based on the curvaton, on the inhomogeneous reheating and on the DiracBornInfeld inflation allow much more significant departures (see Bartolo et al., 2004, and references therein).
It has become common to quantify the level of primordial nonGaussianity adopting the dimensionless nonlinearity parameter (see, e.g., Salopek & Bond, 1990; Gangui, 1994; Verde et al., 2000a; Komatsu & Spergel, 2001), that measures the importance of the quadratic term in a sort of Taylor expansion of the gaugeinvariant Bardeen potential ^{1}^{1}1We recall that on scales smaller than the Hubble radius corresponds to the usual Newtonian peculiar potential (but with changed sign).:
(1) 
here represents a Gaussian random field. Hereafter we will adopt the socalled largescale structure (LSS) convention, where is linearly extrapolated to the present epoch^{2}^{2}2With the cosmological parameters adopted in this paper, this corresponds to values for larger by a factor of with respect to the socalled CMB convention, where is instead extrapolated at .. Moreover we will consider only the socalled local shape for nonGaussianity, in which the bispectrum signal is larger on squeezed triangle configurations. For more details about other possible shapes we refer to Bartolo et al. (2004); Verde & Matarrese (2009); Bartolo & Riotto (2009) and references therein.
At present the most stringent constraints on come from the cosmic microwave background (CMB) data. Their discriminating power derives from the fact that its temperature fluctuations trace the density perturbations before the gravitational nonlinearities modify their original distribution. Whatever is the specific test adopted, all CMB analysis consistently allow only very small deviations from Gaussianity: for instance, analyzing the recent WMAP data, Komatsu et al. (2010) found that varies between 13 and 97, while Smith et al. (2009) found (at 95 per cent confidence level); see also Komatsu et al. (2009).
Alternative and complementary constraints on can in principle be derived analyzing the LSS (for recent reviews, see, e.g. Verde, 2010; Desjacques & Seljak, 2010). As already evident from the first generation of nonGaussian Nbody simulations (Messina et al., 1990; Moscardini et al., 1991; Weinberg & Cole, 1992), the presence of a positive (negative) skewness in the PDF of the primordial density field tends to favor (disfavor) the formation of cosmic structures, inducing a different timing in the whole process of gravitational instability. However, to be fully exploited, this approach needs to be complemented with reliable methods to disentangle from the primordial signal the nonGaussian features introduced by the late nonlinear evolution and by the possible presence of a nonlinear bias factor. While lensing statistics entirely avoid the latter problem, to attack the former one it is necessary to make use of both analytic techniques (like highorder perturbation theory and the TimeRenormalization Group approach; see, e.g., Taruya et al., 2008; Bartolo et al., 2010) and full Nbody simulations (see, e.g., Kang et al., 2007; Grossi et al., 2007; Dalal et al., 2008; Viel et al., 2009; Desjacques et al., 2009; Pillepich et al., 2010; Grossi et al., 2009) to properly calibrate the theoretical predictions in the nonlinear regime. The large amount of theoretical work recently done in this direction allowed to better understand what is the size of the effects on the abundance of nonlinear structures (Matarrese et al., 2000; Verde et al., 2000a; Mathis et al., 2004; Kang et al., 2007; Grossi et al., 2007; Grossi et al., 2009; Maggiore & Riotto, 2010; Roncarelli et al., 2010), on the halo biasing (Dalal et al., 2008; McDonald, 2008; Fedeli et al., 2009; Carbone et al., 2010), on the galaxy bispectrum (Sefusatti & Komatsu, 2007; Jeong & Komatsu, 2009; Nishimichi et al., 2009), on the mass density distribution (Grossi et al., 2008; Lam & Sheth, 2009; Lam et al., 2009), on the topology (Matsubara, 2003; Hikage et al., 2008), on the integrated SachsWolfe effect (Afshordi & Tolley, 2008; Carbone et al., 2008), on the Ly flux from lowdensity intergalactic medium (Viel et al., 2009), on the 21 cm fluctuations (Cooray, 2006; Pillepich et al., 2007) and on the reionization process (Crociani et al., 2009). The first attempts of an application to real observational data gave very encouraging results: Slosar et al. (2008), combining the bias measurements for two samples of luminous red galaxies and quasars, found ; Afshordi & Tolley (2008), studying the integrated SachsWolfe effect (ISW) in the NVSS survey, derived ; all error bars are at 2 level^{3}^{3}3We report the values as revised by Grossi et al. (2009) to include a correction mimicking the ellipsoidal collapse and converted to the LSS convention.. Very recently Xia et al. (2010) found from the analysis of the autocorrelation of the brightest NVSS sources on angular scales of several degrees.
In this paper we will focus on estimating the weak lensing signals in scenarios with primordial nonGaussianity. Being based on the measurement of the shear effect produced by the intervening largescale structure of the Universe on the images of background galaxies, gravitational lensing is a direct probe of the total matter distribution. For this reason it is considered one of the most powerful tools to constrain the main cosmological parameters and many dedicated projects are in progress or under study. A very exciting perspective is certainly opened by the ESA Cosmic Vision project EUCLID (Laureijs, 2009), currently under study: the goal of its wide survey is to obtain the shear measurements for about 40 galaxies per arcmin on the entire extragalactic sky with Galactic latitude (approximately 20000 deg). The possibility of using the weak lensing signals to constrain also the amount of primordial nonGaussianity has been already explored by different authors. Amara & Refregier (2004) used a generalized halo model to study the impact on the estimates of the power spectrum normalization of primordial nonGaussianity, modeled assuming various lognormal distributions for the density field. More recently, Fedeli & Moscardini (2010) computed the power spectrum of the weak cosmic shear for nonGaussian models with different values of . In particular, they improved the halo model including more accurate prescriptions for its ingredients (mass functions, bias and halo profile), calibrated on the last generation of nonGaussian Nbody simulations. The application of this model to a survey having the expected characteristics of the EUCLID project showed the possibility of a significant detection of nonGaussianity at the level of few tens, once the remaining parameters are held fixed.
In this paper we investigate weak lensing statistics in nonGaussian scenarios using numerical rather than analytical tools. Specifically, we will create weak lensing maps performing raytracing simulations through very deep lightcones extracted from highresolution Nbody simulations. The advantage of this approach is twofold. First of all, Nbody simulations permit to fully account for the nonlinear evolution which is usually modeled less accurately by analytical means. Second, numerical experiments allow us to extract a large set of realistic weak lensing maps that can be used for better evaluating the statistical robustness of the results.
The main goal of our numerical work is to figure out what are the observational evidences of the presence of some level of primordial nonGaussianity, as quantified by the parameter. In particular we will compute a large set of weak lensing statistics in models with different and we will quantify the deviations from the corresponding results in the Gaussian scenario. This is important not only to address the possibility of a positive detection with future data, but also to establish at which level an amount of primordial nonGaussianity compatible with the present observational constraints can hamper an accurate measurement of the other cosmological parameters. We must recall that Gaussian initial conditions are virtually always assumed in their practical derivation.
The plan of this paper is as follows. In Section 2 we review the basis of the lensing formalism necessary to the present work. In Section 3 we describe the cosmological simulations and the numerical procedure to build the lensing maps. In Section 4 we present our main results about the statistical properties of the different lensing quantities investigated: the probability distribution function, the thirdorder moment, the power spectrum and the bispectrum. In Section 5 we compare the effects produced by primordial nonGaussianity to the uncertainties related to power spectrum normalization and on the matter density parameter . Finally, in Sect. 6 we draw our conclusions.
2 Lensing theory
In this section we give a short summary of the aspects of the theory of gravitational lensing that will be used throughout this work. For more detail we refer to the review by Bartelmann & Schneider (2001).
We can start by describing how light rays are deflected by the presence of structures in the universe. Since the coherence length of cosmic structures is small compared to the Hubble radius, it is possible to slice the largescale structures into shells and use the socalled thinscreen approximation for them, which allows to consider only the projected mass distribution of the slices. Denoting by the projected mass distribution of the lens at the angular position vector , the convergence can be defined as , where
(2) 
is the critical surface density. In the previous equation , and represent the angulardiameter distances between the observer and the source, between the observer and the lens, and between the lens and the source, respectively.
Thanks to the thinscreen approximation, the object acting as lens can be completely described using its lensing potential , which is related to the convergence through the twodimensional Poisson equation
(3) 
The deflection angle is simply the gradient of the lensing potential, .
Up to second order, lensinginduced image distortions are given by
(4) 
(see, e.g., Goldberg & Bacon, 2005; Bacon et al., 2006), where represent the elements of the Jacobian matrix of the lens equation, are the unlensed coordinates, and the tensor is defined as . In the previous equations . The quantities and can be conveniently expressed as a function of the convergence , of the complex shear term , and of their derivatives:
(5)  
in the previous equations the comma indicates the derivative. The shear derivatives can be combined to construct two new quantities, and , called first and second flexion respectively, defined as
(6)  
(7) 
The previous formalism can be easily generalized to the case in which a continuous distribution of matter is considered. The volume between the observer and the sources can be divided in a sequence of subvolumes having a size along the lineofsight sufficiently small compared to the distances between the observer and the subvolumes, and between those and the sources. The matter of each subvolume can be projected onto a plane and then we are allowed to use again the thinscreen approximation described above. The final quantities, estimated on the source plane, will be the weighted sum of the relevant quantities, where the weight is given by a suitable ratio of the involved distances. In multiple lensplane theory, rotation of light bundles can occur in addition to shear and convergence, but numerical simulations have shown that these are negligibly small (Jain et al., 2000). Then, all the lens properties are contained in the lensing potential. Knowing it on all lens planes allows us to write recursion relations representing the deflection angle, the shear, the effective convergence and the two flexions on the source plane. These relations will be given and discussed in more detail in Section 3.3.
3 The lensing simulations
3.1 The cosmological simulations
To study the effect of nonGaussianity on the weak lensing statistics, we use the outputs of a set of Nbody cosmological simulations following the evolution of dark matter particles only. These simulations have been already used by Grossi et al. (2007), Hikage et al. (2008), Grossi et al. (2008) and Roncarelli et al. (2010). Here we will summarize the information relevant for this paper, referring the interested readers to the original papers for further details.
The simulations were carried out using the publicly available code GADGET2 (Springel, 2005) assuming a CDM model with total matter density , baryon density and cosmological constant density ; the Hubble parameter is set to , while the power spectrum is normalized to . The simulated box has a comoving length of 500 Mpc/h and contains particles with a mass resolution of . The comoving softening scale is . The different outputs are equispaced in comoving space by 250 Mpc/h. The set is composed by seven cosmological simulations considering different values of the nonGaussianity parameter : , , , plus the standard Gaussian case . Notice that in this work we prefer to make use of this set of Nbody simulations in spite of the more recent one presented in Grossi et al. (2009). The reason for this choice is twofold. First, the mass resolution is better and this allows us to have more robust results on small scales, in particular for shear and flexions; second, the assumed values for the parameter span a larger range, covering also cases where the weak lensing signal produced by primordial nonGaussianity is larger.
3.2 Building the mock lightcones
To construct the mock lightcones used to perform the raytracing simulations, we follow the same procedure described in Pace et al. (2007), to which we refer for more details. Since the different outputs of an Nbody simulation represent the redshift evolution of the same initial matter distribution, we need to apply a specific procedure to avoid the introduction of biases related to the fact that the same structures appear approximately at the same positions in different outputs. For this reason we randomly shift and rotate the particle positions exploiting the periodicity of the simulation boxes in the plane perpendicular to the lineofsight. As said before, given our choice for the output redshifts, they overlap for per cent of their comoving sidelength, therefore we can just consider particles in the lower half of the rotated and translated boxes.
Once selected the particles to be used, we project them along the lineofsight on a regular twodimensional grid and compute the projected mass density field using the triangularshapedcloud (TSC) mass assignment (Hockney & Eastwood, 1988). Finally via fast Fourier transform (FFT) techniques it is possible to recover the lensing potential associated with the considered matter distribution. More details on the numerical procedure will be given in the next section.
Our mock lightcones extend along the lineofsight up to redshift , using outputs. As all the distances are in comoving units, the opening angle of the raytracing simulation can be calculated using the comoving distance of the last plane in the stack, resulting in degrees. Using a grid of points, the corresponding angular resolution of the produced maps is 23.7 arcsec.
3.3 Raytracing simulations
Raytracing simulations consist in tracing back a bundle of light rays through the matter distribution of the lightcone, from the observer to the sources, which we place at redshift .
Projected mass maps are converted to projected densitycontrast maps by
(8) 
where is the the area of the grid cell on the th plane, is the mass projected on the grid cell with indices belonging to the th plane and Mpc/h is the depth of the th subvolume used to build mock cone. The density contrast is deliberately defined such as to have the unit of a length (see Hamana & Mellier, 2001). In the previous equation represents the average comoving density of the Universe.
The lensing potential on each plane is related to the projected density contrast through the twodimensional Poisson equation, namely:
(9) 
Exploiting the periodic boundary conditions of the projected maps, Eq. (9) can be solved adopting the FFT techniques. Having the potential on each plane, the lensing quantities can be derived adopting standard finitedifference schemes.
In order to perform the raytracing simulations, we need to apply the multipleplane theory to compute the total effect taking into account the contributions from each single lensing plane. A lightray is deflected on each plane by the amount , thus the total deflection is given by the sum of all contributions. In particular, if the lightcone is sampled into lens planes and the sources are located on the plane, the relation giving the deflection angle on the th plane of a ray with image position reads (in comoving units)
(10) 
here is the unscaled lensing potential, i.e. the Newtonian potential projected along the lineofsight, is a function depending on the cosmology, is the comoving distance and represents the scale factor of the lens plane. Note that in general the lightrays will intercept the lens planes at arbitrary points, while the potential is defined only at the grid points. Thus, it is necessary to use a bilinear interpolation to compute the lensing quantities.
Differentiating Eq. (10) with respect to and defining and the matrix containing the second derivatives of the lensing potential, one obtains
(11) 
In the previous equation, represents the identity matrix. On the source plane, the Jacobian matrix is given by
(12) 
where is now the effective convergence and is the effective shear. The term , called rotation, represents the asymmetry introduced by multiple lenses.
Differentiating Eq. (11) with respect to , a recursive relation for the two flexions can be obtained:
(13) 
where is a tensor containing the third derivatives of the lensing potential.
On the source plane, the tensor reads as
(14)  
where, as for the Jacobian matrix, and represent the asymmetric terms. In Eq. (14) we do not take into account the additional terms, called twist and turn, introduced by Bacon & Schäfer (2009), but several tests performed at different resolutions assure us that their inclusion would not significantly affect our results.
4 Results
In this section we present the results of our analysis of the statistical properties of several lensing quantities extracted from the raytracing simulations described above. In particular we discuss the probability distribution function (PDF) in Sect. 4.1, the thirdorder moment (skewness) in Sect. 4.2, the power spectrum in Sect. 4.3, the shear in aperture in Sect. 4.4 and finally the bispectrum in Sect. 4.5.
4.1 Probability distribution functions
We start considering the PDFs for several important lensing quantities. All our results are grouped in Fig. 1. The panels on the left present the actual PDFs of a given quantity, as obtained averaging over 60 different lightcone realizations. Given the smallness of the differences, we prefer to show only the results for the Gaussian model (red points and lines) and for the most extreme nonGaussian models: (blue lines) and (green lines): models with smaller nonGaussianity are intermediate between the displayed lines. For reference, we also show (for the Gaussian model only) the error bars, referring to the r.m.s. of the corresponding sample of realizations. Finally, the panels on the right present for all models (represented by different colour lines, as described in the label) the ratio between the nonGaussian and Gaussian results.
The top row of Fig. 1 refers to the modulus of the deflection angle . The presence of primordial nonGaussianity affects the PDFs by at most few percent in the most extreme cases with , while the models with primordial nonGaussianity more consistent with current constraints () are almost indistinguishable from the Gaussian case. We also notice that the differences produced by nonGaussianity are more evident in the highvalue tails. These large deflections reflect lensing events associated with rare, large structures. In this sense these results confirm those on the halo abundance: unlikely events (large haloes) are useful probes of primordial nonGaussianity. However, we have to note that in the highvalue tails the statistics can be quite poor and noisy: a good modelization of the strong nonlinear effects acting on the same scales is necessary to allow a positive detection.
In the panels in the second row of Fig. 1 we show the analogous plots for the (modulus of the) shear . Since a negative value of the parameter favors lower values of , we find that the ratio between nonGaussian and Gaussian models is in this case larger than unity; the opposite trend holds for models with positive primordial nonGaussianity. As already pointed out for the deflection angle, the effect of a mild nonGaussianity on the shear distribution is tiny, below 1 per cent, which is comparable to the size of the error bars (shown only in the left panels) but can grow up to 10 per cent for the very extreme tail. Only models with very high values of display a deviation with respect to the Gaussian case that can be as large as the error bars obtained by averaging over the simulated maps.
The PDFs for the first and second flexion are shown in the panels in the third and forth rows of Fig. 1, respectively. The trend is very similar to the case of the shear: compared to the Gaussian model, models with negative have a higher probability of an excess at low values; the opposite trend holds in the highvalue tail. The differences between Gaussian and nonGaussian models are slightly more evident for the second flexion .
Finally the panels in the last row of Fig. 1 refer to the effective convergence, for which we find a slightly different situation. In this case, in the maps there are also negative values, corresponding to underdense regions where the nonGaussian PDFs show large deviations from the Gaussian results: for negative we find less than 1 per cent for , per cent for , and approximately 3 per cent for . Note that this result agrees with the analysis made by Grossi et al. (2008) on the threedimensional density field and its probability distribution: underdense regions tend to maintain the imprint of primordial nonGaussianity, suggesting that statistics based on voids can be a powerful tool to estimate (see also Kamionkowski et al., 2009). The plot shows that also the overdense regions keep the important imprints of the primordial nonGaussianity: for we find 810 per cent for , per cent for , and more than 80 per cent for . However, in this regime, the statistics is quite poor and the ratio between nonGaussian and Gaussian PDFs becomes very noisy. As final comment on PDFs (and on their moments, see the following sections), it is important to stress that the lensing quantities here discussed are observables that directly measure the dark matter distribution. This is not true for the density PDF derived from galaxy surveys, where the complex effect of bias must be corrected for.
4.2 Thirdorder moment
A possible alternative way to detect the signatures produced by the presence of some primordial nonGaussianity is to look at the highorder moments (skewness, kurtosis, etc.) of the distribution of the various lensing quantities. The power of this approach was already evident in the investigation of the corresponding quantities related to the density field (see, e.g., Coles et al., 1993; Grossi et al., 2008; Lam & Sheth, 2009). However, the measurement of these statistics is often affected by large error bars, that increase with the order and hamper their reliable application to the real data. For this reason, in this section we will focus only on the skewness of the twodimensional convergence field , defined as
(15) 
where and represent the mean and the variance, computed on all pixels of the maps. To derive these quantities, we first convolve the convergence map with a Gaussian filter of angular radius , then we subtract the mean and finally we normalize it appropriately.
In the upper panel of Fig. 2 we present the results for the skewness of the effective convergence, as extracted from the Gaussian simulation. For reference we also show the results for the two most extreme nonGaussian models. As expected, (which is computed as average over 60 different realizations) is a decreasing function of the filtering radius : increasing the smoothing reduces the nonGaussian features introduced by the nonlinear evolution. Furthermore the error bars, representing the r.m.s. over the set of different lightcone realizations, are slightly increasing with due to increasing shot noise, since by increasing the smoothing radius we are averaging over a smaller number of circles. This behavior is analogous to what found by Jain et al. (2000) in the determination of the skewness of the effective convergence.
More interesting is the bottom panel of the same figure, where we show the ratio between the thirdorder moments computed in nonGaussian and Gaussian simulations: we find a variation of the order of 12 per cent for , per cent for and only per cent for . It is worth noticing that this effect is almost constant on scales up to 10 arcmin. Comparing these results with the size of the error bars, we can conclude that with relatively small maps, like those analyzed in this paper, we can use the skewness statistics to detect only strong primordial nonGaussianities (). As shown by the results of the analysis made by Fedeli & Moscardini (2010), only with data covering very large area it would be possible to disentangle the effects of nonGaussian models with . Finally we notice that a similar analysis performed on the shear field, which is directly related to the convergence one, would provide comparable results.
4.3 Power Spectrum
An important theoretical quantity, directly related to the matter power spectrum, is the effective convergence power spectrum , which is defined as the squared modulus of the Fourier transform of , averaged on the modes having a given multipole . Starting from , it is possible to derive analytic expressions for the power spectra of the other lensing quantities, such as those for the shear () and the two flexions ( and ), namely:
(16)  
(17) 
In the top panels of Fig. 3 we show, for the Gaussian model only, the comparison between the power spectra extracted from our simulations and the corresponding theoretical predictions including the nonlinear effect. We consider different lensing quantities: effective convergence, shear and reduced shear (upper left panel) and two components of the flexion (top right panel). We notice that the agreement between the numerical results and the theoretical expectations holds, as expected, only in a given range of wavenumbers. On large scales (small ), the departure from theory is due to the low number of modes we can use to average the numerical power spectra. On small scales (large ), the disagreement is produced by the mass resolution (particularly relevant for the flexion power spectrum) and by numerical artifacts produced by the weighting function appearing in the definition of a given power spectrum. A better mass resolution will make the lack of power at high frequencies less severe, but it will not solve it entirely because of the limit of resolution due to the Nyquist frequency intrinsic to every discrete Fourier transformation. Moreover, when numerically performing the raytracing simulations, lowredshift planes are poorly sampled since a few lowredshift pixels contribute to the lensing maps. This poor sampling is amplified by the weighting function appearing in Eqs. 10, 11 and 13. We can use Fig. 3 as a reference to define the range of on which the numerical results can be trusted: typically .
In the bottom panels of Fig. 3 we show the ratio between the power spectra derived from the nonGaussian and the Gaussian simulations. We present here only the results for the effective convergence since the other power spectra can be obtained by introducing the suitable dependence on , and then the corresponding ratios are identical. In this case, we decided to study the effect of primordial nonGaussianity at two different source redshifts, namely and . This is done to verify, with the help of numerical simulations, recent analytical predictions made by Fedeli & Moscardini (2010) using an improved version of the halo model. As expected, the effects of the primordial nonGaussianity become more evident when the absolute value of the parameter increases. At very large scales, where the evolution is still in the linear regime, the Gaussian model and all the nonGaussian ones here considered display approximately the same power spectrum, and the ratio is in practice unity or slightly larger due to the scatter of the power spectrum at low . At smaller scales () the behaviors start to be different: a positive (negative) parameter implies a larger (smaller) power spectrum with respect to the Gaussian case. It is also worth noticing that at very small scales () the ratio tends to an asymptotic value: this is due to the fact that in this range the power spectrum can be very well approximated by a powerlaw relation and all models turn out to have approximately the same slope. We have to recall, however, that some numerical artifacts can affect our results on these small scales. Looking at the plots, we see that, as expected, the differences between Gaussian and nonGaussian models are tiny, for the power spectrum too and increase progressively with the source redshift: for the most extreme models (), the deviations are of the order of 8 per cent at (1012 per cent at ), for the models with they are approximately 45 per cent, while for the models with the deviations are smaller than 1 per cent. From the bottom panels of Fig. 3 we also see that the peak of the ratio between the nonGaussian and the Gaussian models shifts at higher frequencies when the source redshift increases. This is in agreement with what was found recently by Fedeli & Moscardini (2010).
4.4 Shear in aperture
A quantity which is directly related to the shear power spectrum but which can be more easily measured in real data is the socalled shear in aperture . It represents the variance of the shear field within a circular aperture of a given radius :
(18) 
where is the firstorder Bessel function of the first kind.
The top panel of Fig. 4 shows the comparison between the theoretical expectation and the results derived from our Gaussian simulations (averaged over 60 different lightcone realizations): the good agreement is evident up to scales of arcmin. For reference we also plot the results for . The bottom panel of the same figure presents the ratio between the results for nonGaussian and Gaussian models as a function of the aperture angle. Being a weighted integral of the shear power spectrum, has a trend which is quite similar to the one displayed by : the larger the absolute value of the parameter is, the larger are the deviations from the Gaussian expectations. An increment (decrement) in the signal of the order of 78 per cent is expected for the model with (), while the differences reduce to less than 1 per cent for the models with . As expected, we find that the deviations from the Gaussian expectations are larger for small apertures, but rapidly wash out for angles larger than few arcmin, since the average of the shear involves more and more large scale structures. We also notice that the ratio at large smoothing scales does not reach unity since structures are already evolved enough to show high intrinsic nonGaussianity, therefore higher smoothing scales are necessary to compensate the effect.
4.5 Bispectrum
In several works (see e.g. Verde et al., 1998; Verde et al., 2000b; Takada & Jain, 2004; Sefusatti et al., 2006, 2010) it has been shown that a large amount of cosmological information can be inferred from the study of highorder spectra (the socalled polyspectra). Here we just limit our investigation to the simplest highorder spectrum, the bispectrum. We recall that for a Gaussian random field all available information is contained in the power spectrum, while all higherorder spectra vanish or are combinations of the power spectrum. As a consequence the measurement of the bispectrum is quite sensitive to the nonGaussian features of the quantity under investigation.
In general, the bispectrum depends on a triangular configuration and its values are very sensitive to the particular configuration adopted. Since the complete determination of the bispectrum is extremely expensive from a computational point of view, we will restrict our analysis to equilateral triangles only. Notice that, given our assumption of local shape for the primordial nonGaussianity, the nonGaussian signal on squeezed configurations is expected to be larger. In a future work we will use the convergence maps presented in this work to study the full bispectrum configuration and how the signal expected from primordial nonGaussianity is sensitive to the triangular configuration.
In general the bispectrum of a given quantity is defined as:
(19) 
In the previous equation, represents the Fourier transform of the field . Since the wavenumbers at which the bispectrum is computed (, and ) must form a closed triangle in Fourier space, they are related by the constraint . Finally is the angle between and and fixes, together with the triangle condition, .
Similar to Fig. 3, the top panel of Fig. 5 shows, for the Gaussian model only, the comparison between our numerical estimates and the theoretical bispectrum computed for the effective convergence. We show a bispectrum computed for an equilateral configuration, i.e., in the Fourier space, . The results from the simulations represent the average over 60 different lightcone realizations with the corresponding r.m.s. plotted as error bars. The theoretical quantity are derived from the matter density bispectrum , that can be computed according to the perturbation theory. For instance, in the case of the effective convergence , the bispectrum can be obtained by using the Limber’s approximation in an analogous way as for the power spectrum (see e.g. Cooray & Hu, 2001), namely:
(20)  
where is a function depending on the cosmological parameters and contains the distance weight (Bartelmann & Schneider, 2001). For the other lensing quantities, similar relations hold. In order to quantify the effects of primordial nonGaussianity, the panel on the bottom shows the bispectrum computed for the different nonGaussian models, normalized to the corresponding result for the Gaussian one.
As for the case of the power spectrum, we see that the agreement between the theoretical expectations and the numerical predictions holds only for a limited range of angular scales, but in this case the range for which the estimate of bispectrum is robust is much more limited. Focusing on scales where the numerical results can be considered reliable, we find that the contribution of nonGaussianity is not very relevant, reaching up to 12 per cent in the very nonlinear regime, but only for the most extreme nonGaussian models (). For models with in the range allowed by the most recent CMB constraints, we do not expect a difference with respect to the Gaussian case larger than few per cent.
5 The degeneracy between primordial nonGaussianity and the power spectrum normalization and matter density parameter
As shown in the previous sections, the effects due to primordial nonGaussianity on the various weak lensing quantities are small, in particular when we consider models with values of in the range allowed by present data. To better investigate the possibility of their future detection, here we will discuss their degeneracy with the effects produced by the variation of other relevant cosmological parameters. In particular we will focus on the degeneracy between the amount of primordial nonGaussianity and the power spectrum normalization and matter density parameter , which are certainly the most important source of uncertainty.
Since the different timing in the structure formation induced by modifies the power spectrum and its amplitude, nonGaussianity can be confused with a Gaussian model with a different or parameter. In particular a positive (negative) has effects similar to a higher (lower) power spectrum normalization or matter density parameter. In this section we show how much the uncertainty on the value of the normalization affects the possibility of detecting nonGaussianity, keeping all the other parameters fixed. We also show how important is the error in the determination of the matter density parameter. To do so, we vary at the same time both and to have a flat geometry. For the sake of simplicity, we restrict our analysis to the convergence power spectrum, the shear in aperture and the effective convergence bispectrum.
In the left column of Fig. 6 we show a comparison between the contribution of nonGaussianity and the changes produced by a variation of the spectrum normalization while on the right panels we show the effect of varying the matter content, keeping . In particular we consider () in the range between 0.87 and 0.93 (0.27 and 0.33): the size of this variation approximately corresponds to the typical 2 error derived from the most recent CMB analysis. It is clear from the figure that the error on the determination of the parameters has effects more important than the introduction of a weak primordial nonGaussianity. The two most extreme models produce deviations from the Gaussian predictions that are of the same level of uncertainty and . On the contrary, the models with in the range allowed by current constraints have effects well within the range allowed by the present uncertainties, since their effect is smaller than a variation of . This implies that significant constraints on the amount of primordial nonGaussianity can be obtained only reducing the uncertainties on the power spectrum normalization and on the matter density parameter.
6 Conclusions
In this work, we used the outputs of Nbody simulations to create a large set of realistic mock maps for several lensing quantities (deflection angle, effective convergence, shear and the two components of the flexion) in the framework of cosmological models with different amount of primordial nonGaussianity, quantified using the dimensionless parameter . In particular we considered several statistical properties (PDF, power spectrum, shear in aperture, skewness and bispectrum) and compared the results with the corresponding ones obtained in the Gaussian case. Our main results can be summarized as follows.

For all quantities here considered the effect produced by the presence of primordial nonGaussianity is relatively small, amounting to differences of 1 per cent for , 5 per cent for , and 815 per cent for . These results are in good agreement with the analytic predictions presented in Fedeli & Moscardini (2010).

The largest effects are visible on small scales (i.e. for large multipoles ), where, however, also nonlinearity can produce strong effects which have to be accurately modeled.

The most promising statistical tests to search for imprints of primordial nonGaussianity are the (convergence and shear) power spectra and the (convergence) bispectrum, thanks to the smaller size of their error bars at the relevant scales.

The differences of the various PDFs in both rareevent tails are also important, but their discriminating power is reduced by the poor statistics and by the highlevel of noise.

We compared the effects produced by the primordial nonGaussianity with the uncertainties due to the power spectrum normalization and matter density parameter: an error in the determination of of about 3 per cent or of of about 10 per cent gives differences comparable with the nonGaussian models with , while for more realistic nonGaussian models with the effects is smaller than the one induced by .
As said, a significant covariance exists between primordial nonGaussianity and fundamental cosmological parameters, especially with and . The smallness of nonGaussian effects found in our analysis means that one can obtain a precise estimate of these parameters despite their covariance with . As an example, ignoring the possible presence of a primordial nonGaussianity with , consistent with current observational constraints, would induce a mere 0.2 per cent uncertainty in the estimate of the cosmological parameters, i.e. well below the current 1 error. The same argument, in reverse, is telling us that the search for nonGaussian features in the weaklensing statistics is a very challenging task. According to our results, it can be successfully completed only if errorbars can be significantly reduced. For the purpose of discriminating among competing models the reduction needs not to be dramatic. Indeed, as pointed out by Fedeli & Moscardini (2010), the fact that deviations from Gaussianity are small but systematics allows one to estimate by adding the statistical information from different angular bins, this is also shown in our analysis of the power spectrum for sources at different redshifts.
A more significant reduction in the errorbars is required to break the degeneracy between nonGaussianity with cosmological parameters, say or . Indeed, as shown in Fig. 6, removing such degeneracy requires to compare information from a limited number of bins at very different scales. As shown by Fedeli & Moscardini (2010), both tasks could be achieved with nextgeneration, allsky surveys. A good example is certainly represented by the ESA mission EUCLID (Laureijs, 2009). According to the analysis made by Fedeli & Moscardini (2010) and based on analytic predictions for the power spectrum quantitatively confirmed in this numerical work, we expect that the quality and the quantity of the EUCLID data will allow to constrain at the level of few tens, opening the possibility of discriminating between the various inflationary models.
Acknowledgments
Computations have been performed on the IBMSP5 at CINECA (Consorzio Interuniversitario del NordEst per il Calcolo Automatico), Bologna, with CPU time assigned under an INAFCINECA grant, and on the IBMSP4 machine at the “Rechenzentrum der MaxPlanckGesellschaft” at the MaxPlanck Institut fuer Plasmaphysik with CPU time assigned to the “MaxPlanckInstitut für Astrophysik” and at the “LeibnizRechenzentrum” with CPU time assigned to the Project “h0073”. This work was supported by the Deutsche Forschungsgemeinschaft (DFG) under the grants BA 1369/51 and 1369/52 and through the Transregio Sonderforschungsbereich TR 33, as well as by the DAAD and CRUI through their Vigoni program. We also acknowledge partial financial contributions from the ASI contracts I/016/07/0 COFIS, EuclidDune I/064/08/0, ASIINAF I/023/05/0, ASIINAF I/088/06/0, ASI/INAF Agreement I/072/09/0 for the Planck LFI Activity of Phase E2 and from the grant ANR06JCJC0141. K. D. acknowledges the support by the DFG Priority Programme 1177 and additional support by the DFG Cluster of Excellence ”Origin and Structure of the Universe”. We warmly thank Massimo Meneghetti, Cosimo Fedeli and Melita Carbone for useful discussions and comments. We also thank the anonymous referee whose comments helped us to improve the presentation of our results.
References
 Afshordi & Tolley (2008) Afshordi N., Tolley A. J., 2008, Phys. Rev. D, 78, 123507
 Amara & Refregier (2004) Amara A., Refregier A., 2004, MNRAS, 351, 375
 Bacon et al. (2006) Bacon D. J., Goldberg D. M., Rowe B. T. P., Taylor A. N., 2006, MNRAS, 365, 414
 Bacon & Schäfer (2009) Bacon D. J., Schäfer B. M., 2009, MNRAS, 396, 2167
 Bartelmann & Schneider (2001) Bartelmann M., Schneider P., 2001, Physics Reports, 340, 291
 Bartolo et al. (2010) Bartolo N., Beltran Almeida J. P., Matarrese S., Pietroni M., Riotto A., 2010, Journal of Cosmology and AstroParticle Physics, 3, 11
 Bartolo et al. (2004) Bartolo N., Komatsu E., Matarrese S., Riotto A., 2004, Physics Reports, 402, 103
 Bartolo & Riotto (2009) Bartolo N., Riotto A., 2009, Journal of Cosmology and AstroParticle Physics, 3, 17
 Carbone et al. (2010) Carbone C., Mena O., Verde L., 2010, Journal of Cosmology and AstroParticle Physics, 7, 20
 Carbone et al. (2008) Carbone C., Verde L., Matarrese S., 2008, ApJL, 684, L1
 Coles et al. (1993) Coles P., Moscardini L., Lucchin F., Matarrese S., Messina A., 1993, MNRAS, 264, 749
 Cooray (2006) Cooray A., 2006, Physical Review Letters, 97, 261301
 Cooray & Hu (2001) Cooray A., Hu W., 2001, ApJ, 548, 7
 Crociani et al. (2009) Crociani D., Moscardini L., Viel M., Matarrese S., 2009, MNRAS, 394, 133
 Dalal et al. (2008) Dalal N., Doré O., Huterer D., Shirokov A., 2008, Phys. Rev. D, 77, 123514
 Desjacques & Seljak (2010) Desjacques V., Seljak U., 2010, Classical and Quantum Gravity, 27, 12
 Desjacques et al. (2009) Desjacques V., Seljak U., Iliev I. T., 2009, MNRAS, 396, 85
 Fedeli & Moscardini (2010) Fedeli C., Moscardini L., 2010, MNRAS, 405, 681
 Fedeli et al. (2009) Fedeli C., Moscardini L., Matarrese S., 2009, MNRAS, 397, 1125
 Gangui (1994) Gangui A., 1994, Phys. Rev. D, 50, 3684
 Goldberg & Bacon (2005) Goldberg D. M., Bacon D. J., 2005, ApJ, 619, 741
 Grossi et al. (2008) Grossi M., Branchini E., Dolag K., Matarrese S., Moscardini L., 2008, MNRAS, 390, 438
 Grossi et al. (2007) Grossi M., Dolag K., Branchini E., Matarrese S., Moscardini L., 2007, MNRAS, 382, 1261
 Grossi et al. (2009) Grossi M., Verde L., Carbone C., Dolag K., Branchini E., Iannuzzi F., Matarrese S., Moscardini L., 2009, MNRAS, 398, 321
 Hamana & Mellier (2001) Hamana T., Mellier Y., 2001, MNRAS, 327, 169
 Hikage et al. (2008) Hikage C., Coles P., Grossi M., Moscardini L., Dolag K., Branchini E., Matarrese S., 2008, MNRAS, 385, 1613
 Hockney & Eastwood (1988) Hockney R. W., Eastwood J. W., 1988, Computer simulation using particles
 Jain et al. (2000) Jain B., Seljak U., White S., 2000, ApJ, 530, 547
 Jeong & Komatsu (2009) Jeong D., Komatsu E., 2009, ApJ, 703, 1230
 Kamionkowski et al. (2009) Kamionkowski M., Verde L., Jimenez R., 2009, Journal of Cosmology and AstroParticle Physics, 1, 10
 Kang et al. (2007) Kang X., Norberg P., Silk J., 2007, MNRAS, 376, 343
 Komatsu et al. (2009) Komatsu E., Dunkley J., Nolta M. R., et al. 2009, ApJS, 180, 330
 Komatsu et al. (2010) Komatsu E., Smith K. M., Dunkley J., Bennett C. L., Gold et. al 2010, ArXiv eprints, 1001.4538
 Komatsu & Spergel (2001) Komatsu E., Spergel D. N., 2001, Phys. Rev. D, 63, 063002
 Lam et al. (2009) Lam T. Y., Desjacques V., Sheth R. K., 2010, MNRAS, 402, 2397
 Lam & Sheth (2009) Lam T. Y., Sheth R. K., 2009, MNRAS, 395, 1743
 Laureijs (2009) Laureijs R., 2009, ArXiv eprints, 0912.0914
 Maggiore & Riotto (2010) Maggiore M., Riotto A., 2010, ApJ, 717, 526
 Matarrese et al. (2000) Matarrese S., Verde L., Jimenez R., 2000, ApJ, 541, 10
 Mathis et al. (2004) Mathis H., Diego J. M., Silk J., 2004, MNRAS, 353, 681
 Matsubara (2003) Matsubara T., 2003, ApJ, 584, 1
 McDonald (2008) McDonald P., 2008, Phys. Rev. D, 78, 123519
 Messina et al. (1990) Messina A., Moscardini L., Lucchin F., Matarrese S., 1990, MNRAS, 245, 244
 Moscardini et al. (1991) Moscardini L., Matarrese S., Lucchin F., Messina A., 1991, MNRAS, 248, 424
 Nishimichi et al. (2009) Nishimichi T., Taruya A., Koyama K., Sabiu C., 2010, Journal of Cosmology and AstroParticle Physics, 7, 2
 Pace et al. (2007) Pace F., Maturi M., Meneghetti M., Bartelmann M., Moscardini L., Dolag K., 2007, A&A, 471, 731
 Pillepich et al. (2010) Pillepich A., Porciani C., Hahn O., 2010, MNRAS, 402, 191
 Pillepich et al. (2007) Pillepich A., Porciani C., Matarrese S., 2007, ApJ, 662, 1
 Roncarelli et al. (2010) Roncarelli M., Moscardini L., Branchini E., Dolag K., Grossi M., Iannuzzi F., Matarrese S., 2010, MNRAS, 402, 923
 Salopek & Bond (1990) Salopek D. S., Bond J. R., 1990, Phys. Rev. D, 42, 3936
 Sefusatti et al. (2010) Sefusatti E., Crocce M., Desjacques V., 2010, MNRAS, 406, 1014
 Sefusatti et al. (2006) Sefusatti E., Crocce M., Pueblas S., Scoccimarro R., 2006, Phys. Rev. D, 74, 023522
 Sefusatti & Komatsu (2007) Sefusatti E., Komatsu E., 2007, Phys. Rev. D, 76, 083004
 Slosar et al. (2008) Slosar A., Hirata C., Seljak U., Ho S., Padmanabhan N., 2008, Journal of Cosmology and AstroParticle Physics, 8, 31
 Smith et al. (2009) Smith K. M., Senatore L., Zaldarriaga M., 2009, Journal of Cosmology and AstroParticle Physics, 9, 6
 Springel (2005) Springel V., 2005, MNRAS, 364, 1105
 Takada & Jain (2004) Takada M., Jain B., 2004, MNRAS, 348, 897
 Taruya et al. (2008) Taruya A., Koyama K., Matsubara T., 2008, Phys. Rev. D, 78, 123534
 Verde (2010) Verde L., 2010, Advances in Astronomy, 2010, 64
 Verde et al. (2000b) Verde L., Heavens A. F., Matarrese S., 2000b, MNRAS, 318, 584
 Verde et al. (1998) Verde L., Heavens A. F., Matarrese S., Moscardini L., 1998, MNRAS, 300, 747
 Verde & Matarrese (2009) Verde L., Matarrese S., 2009, ApJL, 706, L91
 Verde et al. (2000a) Verde L., Wang L., Heavens A. F., Kamionkowski M., 2000a, MNRAS, 313, 141
 Viel et al. (2009) Viel M., Branchini E., Dolag K., Grossi M., Matarrese S., Moscardini L., 2009, MNRAS, 393, 774
 Weinberg & Cole (1992) Weinberg D. H., Cole S., 1992, MNRAS, 259, 652
 Xia et al. (2010) Xia J., Viel M., Baccigalupi C., De Zotti G., Matarrese S., Verde L., 2010, ApJL, 717, 17