We present a new method, called x-cut cosmic shear, which optimally removes sensitivity to poorly modeled scales from the two-point cosmic shear signal. We show that the x-cut cosmic shear covariance matrix can be computed from the correlation function covariance matrix in a few minutes, enabling a likelihood analysis at virtually no additional computational cost. Further we show how to generalize x-cut cosmic shear to galaxy-galaxy lensing. Performing an x-cut cosmic shear analysis of the Dark Energy Survey Year 1 (DESY1) shear data, we reduce the error on S8 by 20% relative to a correlation function analysis, while showing our constraints are robust to different baryonic feedback models. Largely driven by information at small angular scales, our result, S8 = 0.734 plus/minus 0.038, yields a 2.6 sigma tension with the Planck Legacy analysis of the cosmic microwave background. As well as alleviating baryonic modelling uncertainties, our method can be used to optimally constrain a large number of theories of modified gravity where computational limitations make it infeasible to model the power spectrum down to extremely small scales. The key parts of our code are made publicly available.

We show that it is possible to build effective matter density power spectra in tomographic cosmic shear observations that exhibit the Baryonic Acoustic Oscillations (BAO) features once a nulling transformation has been applied to the data. The precision with which the amplitude and position of these features can be reconstructed is quantified in terms of sky coverage, intrinsic shape noise, median source redshift and number density of sources. BAO detection in Euclid or LSST like wide surveys will be possible with a modest signal-to-noise ratio. It would improve dramatically for slightly deeper surveys.

We address key points for an efficient implementation of likelihood codes for modern weak lensing large-scale structure surveys. Specifically, we focus on the joint weak lensing convergence power spectrum-bispectrum probe and we tackle the numerical challenges required by a realistic analysis. Under the assumption of (multivariate) Gaussian likelihoods, we have developed a high performance code that allows highly parallelised prediction of the binned tomographic observables and of their joint non-Gaussian covariance matrix accounting for terms up to the 6-point correlation function and super-sample effects. This performance allows us to qualitatively address several interesting scientific questions. We find that the bispectrum provides an improvement in terms of signal-to-noise ratio (S/N) of about 10% on top of the power spectrum, making it a non-negligible source of information for future surveys. Furthermore, we are capable to test the impact of theoretical uncertainties in the halo model used to build our observables; with presently allowed variations we conclude that the impact is negligible on the S/N. Finally, we consider data compression possibilities to optimise future analyses of the weak lensing bispectrum. We find that, ignoring systematics, 5 equipopulated redshift bins are enough to recover the information content of a Euclid-like survey, with negligible improvement when increasing to 10 bins. We also explore principal component analysis and dependence on the triangle shapes as ways to reduce the numerical complexity of the problem.

Weak lensing of galaxies and CMB photons through the large-scale structure of the Universe is one of the most promising cosmological probes with upcoming experiments dedicated to its measurements such as Euclid/LSST and CMB Stage 4 experiments. With increasingly precise measurements, there is a dire need for accurate theoretical predictions. In this work, we focus on higher order statistics of the weak lensing convergence field, namely its cumulants such as skewness and kurtosis and its one-point probability distribution (PDF), and we quantify using perturbation theory the corrections coming from post-Born effects, meaning beyond the straight-line and independent lenses approximations. At first order, two such corrections arise: lens-lens couplings and geodesic deviation. Though the corrections are small for low source redshifts (below a few percents) and therefore for galaxy lensing, they become important at higher redshifts, notably in the context of CMB lensing, where the non-gaussianities computed from tree-order perturbation theory are found to be of the same order as the signal itself. We include these post-Born corrections on the skewness into a prediction for the one-point convergence PDF obtained with large deviation theory and successfully test these results against numerical simulations. The modelled PDF is indeed shown to perform better than the percent for apertures above ~ 10 arcminutes and typically in the three sigmas region around the mean.

The power spectrum response function of the large-scale structure of the Universe describes how the evolved power spectrum is modified by a small change in initial power through non-linear mode coupling of gravitational evolution. It was previously found that the response function for the coupling from small to large scales is strongly suppressed in amplitude, especially at late times, compared to predictions from perturbation theory (PT) based on the single-stream approximation. One obvious explanation for this is that PT fails to describe the dynamics beyond shell-crossing. We test this idea by comparing measurements in $N$-body simulations to prescriptions based on PT but augmented with adaptive smoothing to account for the formation of non-linear structures of various sizes in the multi-stream regime. We first start with one-dimensional (1D) cosmology, where the Zel'dovich approximation provides the exact solution in the single stream regime. Similarly to the 3D case, the response function of the large-scale modes exhibits a strong suppression in amplitude at small scales which cannot be explained by the Zel'dovich solution alone. However, by performing adaptive smoothing of initial conditions to identify haloes of different sizes and solving approximately post-collapse dynamics in the 3-streams regime, agreement between theory and simulations drastically improves. We extend our analyses to the 3D case using PINOCCHIO algorithm, in which similar adaptive smoothing is implemented on the Lagrangian PT fields to identify haloes and is combined with a spherical halo prescription to account for post-collapse dynamics. Again, a suppression is found in the coupling between small- and large-scale modes and the agreement with simulations is improved.

The distribution of the cosmic convergence field is modeled using a large-deviation principle where all non-Gaussian contributions are computed from first principles. The geometry of the past light cone is accounted for by constructing the total weak-lensing signal from contributions of the matter density in thin disk slices. The prediction of this model is successfully tested against numerical simulation with ray tracing, and found to be accurate within at least 5 per cent in the tails at redshift 1 and opening angle of 10 arcmin and even more so with increasing source redshift and opening angle. An accurate analytical approximation to the theory is also provided for practical implementation. The lensing kernel that mixes physical scales along the line-of-sight tends to reduce the domain of validity of this theoretical approach compared to the three dimensional case of cosmic densities in spherical cells. This effect is shown to be avoidable if a nulling procedure is implemented in order to localise the lensing line-of-sight integrations in a tomographic analysis. Accuracy in the tails is thus achieved within a percent for source redshifts between 0.5 and 1.5 and an opening angle of 10 arcmin. Applications to future weak-lensing surveys like Euclid and the specific issue of shape noise are discussed.

This white paper is the result of the Tri-Agency Working Group (TAG) appointed to develop synergies between missions and is intended to clarify what LSST observations are needed in order to maximally enhance the combined science output of LSST and Euclid. To facilitate LSST planning we provide a range of possible LSST surveys with clear metrics based on the improvement in the Dark Energy figure of merit (FOM). To provide a quantifiable metric we present five survey options using only between 0.3 and 3.8% of the LSST 10 year survey. We also provide information so that the LSST DDF cadence can possibly be matched to those of \emph{Euclid} in common deep fields, SXDS, COSMOS, CDFS, and a proposed new LSST deep field (near the Akari Deep Field South). Co-coordination of observations from the Large Synoptic Survey Telescope (LSST) and Euclid will lead to a significant number of synergies. The combination of optical multi-band imaging from LSST with high resolution optical and near-infrared photometry and spectroscopy from \emph{Euclid} will not only improve constraints on Dark Energy, but provide a wealth of science on the Milky Way, local group, local large scale structure, and even on first galaxies during the epoch of reionization. A detailed paper has been published on the Dark Energy science case (Rhodes et al.) by a joint LSST/Euclid working group as well as a white paper describing LSST/Euclid/WFIRST synergies (Jain et al.), and we will briefly describe other science cases here. A companion white paper argues the general science case for an extension of the LSST footprint to the north at airmass < 1.8, and we support the white papers for southern extensions of the LSST survey.

We study the accuracy with which cosmological parameters can be determined from real space power spectrum of matter density contrast at weakly nonlinear scales using analytical approaches. From power spectra measured in $N$-body simulations and using Markov chain Monte-Carlo technique, the best-fitting cosmological input parameters are determined with several analytical methods as a theoretical template, such as the standard perturbation theory, the regularized perturbation theory, and the effective field theory. We show that at redshift 1, all two-loop level calculations can fit the measured power spectrum down to scales $k \sim 0.2 \, h \, \mathrm{Mpc}^{-1}$ and cosmological parameters are successfully estimated in an unbiased way. Introducing the Figure of bias (FoB) and Figure of merit (FoM) parameter, we determine the validity range of those models and then evaluate their relative performances. With one free parameter, namely the damping scale, the regularized perturbation theory is found to be able to provide the largest FoM parameter while keeping the FoB in the acceptance range.

We use a simple cosmological model with two parameters $(A_s, n_s)$ to illustrate the impact of using Halofit on error forecast based on Fisher information matrix for a $h^{-3} Gpc^3$ volume survey. We show that Halofit fails to reproduce well the derivatives of the power spectrum with respect to the cosmological parameters despite the good fit produced for its amplitude. We argue that the poor performance on the derivatives prediction is a general feature of this model and we exhibit the response function for the Halofit to show how it compares with the same quantity measured on simulations. The analytic structure of the Halofit response function points towards the origin of its weak performance at reproducing the derivatives of the non-linear power spectrum, which translates into unreliable Fisher information matrices.

If left unchecked modeling uncertainties at small scales, due to poorly understood baryonic physics and non-linear structure formation, will significantly bias Stage IV cosmic shear two-point statistic parameter constraints. While it is perhaps possible to run N-body or hydrodynamical simulations to determine the impact of these effects this approach is computationally expensive; especially to test a large number of theories of gravity. Instead we propose directly removing sensitivity to small-scale structure from the lensing spectrum, creating a statistic that is robust to these uncertainties. We do this by taking a redshift-dependent l-cut after applying the Bernardeau-Nishimichi-Taruya (BNT) nulling scheme. This reorganizes the information in the lensing spectrum to make the relationship between the angular scale, l, and the structure scale, k, much clearer compared to standard cosmic shear power spectra -- for which no direct relationship exists. We quantify the effectiveness of this method at removing sensitivity to small scales and compute the predicted Fisher error on the dark energy equation of state, w0, for different k-cuts in the matter power spectrum.

We present numerical measurements of the power spectrum response function of the gravitational growth of cosmic structures, defined as the functional derivative of the nonlinear spectrum with respect to the linear counterpart, based on $1,400$ cosmological simulations. We develop a simple analytical model based on a regularization of the standard perturbative calculation. Using the model prediction, we show that this function gives a natural way to interpolate the nonlinear power spectrum over cosmological parameter space from single or multi-step interpolations. We demonstrate that once an accurate numerical spectrum template is available for one (or a small number of) cosmological model(s), it doubles the range in $k$ for which percent level accuracy can be obtained even for large change in the cosmological parameters. The python package RESPRESSO we developed to make those predictions is publicly available.

Large deviation statistics is implemented to predict the statistics of cosmic densities in cylinders applicable to photometric surveys. It yields few percent accurate analytical predictions for the one-point probability distribution function (PDF) of densities in concentric or compensated cylinders; and also captures the density-dependence of their angular clustering (cylinder bias). All predictions are found to be in excellent agreement with the cosmological simulation Horizon Run 4 in the quasi-linear regime where standard perturbation theory normally breaks down. These results are combined with a simple local bias model that relates dark matter and tracer densities in cylinders and validated on simulated halo catalogues. This formalism can be used to probe cosmology with existing and upcoming photometric surveys like DES, Euclid or WFIRST containing billions of galaxies.

Non-Gaussianities of dynamical origin are disentangled from primordial ones using the formalism of large deviation statistics with spherical collapse dynamics. This is achieved by relying on accurate analytical predictions for the one-point probability distribution function (PDF) and the two-point clustering of spherically-averaged cosmic densities (sphere bias). Sphere bias extends the idea of halo bias to intermediate density environments and voids as underdense regions. In the presence of primordial non-Gaussianity, sphere bias displays a strong scale dependence relevant for both high and low density regions, which is predicted analytically. The statistics of densities in spheres are built to model primordial non-Gaussianity via an initial skewness with a scale-dependence that depends on the bispectrum of the underlying model. The analytical formulas with the measured nonlinear dark matter variance as input are successfully tested against numerical simulations. For local non-Gaussianity with a range from $f_{\rm NL}=-100$ to $+100$ they are found to agree within 2\% or better for densities $\rho\in[0.5,3]$ in spheres of radius 15 Mpc$/h$ down to $z=0.35$. The validity of the large deviation statistics formalism is thereby established for all observationally relevant local-type departures from perfectly Gaussian initial conditions. The corresponding estimators for the amplitude of the nonlinear variance $\sigma_8$ and primordial skewness $f_{\rm NL}$ are validated using a fiducial joint maximum likelihood experiment. The influence of observational effects and the prospects for a future detection of primordial non-Gaussianity from joint one- and two-point densities-in-spheres statistics are discussed.

This paper presents a general formalism that allows the derivation of the cumulant generating function and one-point Probability Distribution Function (PDF) of the aperture mass ($\hat{M}_{ap}$), a common observable for cosmic shear observations. Our formalism is based on the Large Deviation Principle (LDP) applied, in such cosmological context, to an arbitrary set of densities in concentric cells. We show here that the LDP can indeed be used for a much larger family of observables than previously envisioned, such as those built from continuous and nonlinear functionals of the density profiles. The general expression of the observable aperture mass depends on reduced shear profile making it a rather involved function of the projected density field. Because of this difficulty, an approximation that is commonly employed consists in replacing the reduced shear by the shear in such a construction neglecting therefore non-linear effects. We were precisely able to quantify how this approximation affects the $\hat{M}_{ap}$ statistical properties. In particular we derive the corrective term for the skewness of the $\hat{M}_{ap}$ and reconstruct its one-point PDF.

Starting from a very accurate model for density-in-cells statistics of dark matter based on large deviation theory, a bias model for the tracer density in spheres is formulated. It adopts a mean bias relation based on a quadratic bias model to relate the log-densities of dark matter to those of mass-weighted dark haloes in real and redshift space. The validity of the parametrised bias model is established using a parametrisation-independent extraction of the bias function. This average bias model is then combined with the dark matter PDF, neglecting any scatter around it: it nevertheless yields an excellent model for densities-in-cells statistics of mass tracers that is parametrised in terms of the underlying dark matter variance and three bias parameters. The procedure is validated on measurements of both the one and two point statistics of subhalo densities in the state-of-the-art Horizon Run 4 simulation showing excellent agreement for measured dark matter variance and bias parameters. Finally, it is demonstrated that this formalism allows for a joint estimation of the nonlinear dark matter variance and the bias parameters using solely the statistics of subhaloes. Having verified that galaxy counts in hydrodynamical simulations sampled on a scale of 10 Mpc/h closely resemble those of subhaloes, this work provides important steps towards making theoretical predictions for density-in-cells statistics applicable to upcoming galaxy surveys like Euclid or WFIRST.

Simple parameter-free analytic bias functions for the two-point correlation of densities in spheres at large separation are presented. These bias functions generalize the so-called Kaiser bias to the mildly non-linear regime for arbitrary density contrasts. The derivation is carried out in the context of large deviation statistics while relying on the spherical collapse model. A logarithmic transformation provides a saddle approximation which is valid for the whole range of densities and shown to be accurate against the 30 Gpc cube state-of-the-art Horizon Run 4 simulation. Special configurations of two concentric spheres that allow to identify peaks are employed to obtain the conditional bias and a proxy to BBKS extrema correlation functions. These analytic bias functions should be used jointly with extended perturbation theory to predict two-point clustering statistics as they capture the non-linear regime of structure formation at the percent level down to scales of about 10 Mpc/h at redshift 0. Conversely, the joint statistics also provide us with optimal dark matter two-point correlation estimates which can be applied either universally to all spheres or to a restricted set of biased (over- or underdense) pairs. Based on a simple fiducial survey, this estimator is shown to perform five times better than usual two-point function estimators. Extracting more information from correlations of different types of objects should prove essential in the context of upcoming surveys like Euclid, DESI, PFS or LSST.

The analytical formalism to obtain the probability distribution functions (PDFs) of spherically-averaged cosmic densities and velocity divergences in the mildly non-linear regime is presented. A large-deviation principle is applied to those cosmic fields assuming their most likely dynamics in spheres is set by the spherical collapse model. We validate our analytical results using state-of-the-art dark matter simulations with a phase-space resolved velocity field finding a 2% percent level agreement for a wide range of velocity divergences and densities in the mildly nonlinear regime (~10Mpc/h at redshift zero), usually inaccessible to perturbation theory. From the joint PDF of densities and velocity divergences measured in two concentric spheres, we extract with the same accuracy velocity profiles and conditional velocity PDF subject to a given over/under-density which are of interest to understand the non-linear evolution of velocity flows. Both PDFs are used to build a simple but accurate maximum likelihood estimators for the redshift evolution of the variance of both the density and velocity divergence fields, which have smaller relative errors than their sample variances when non-linearities appear. Given the dependence of the velocity divergence on the growth rate, there is a significant gain in using the full knowledge of both PDFs to derive constraints on the equation of state of dark energy. Thanks to the insensitivity of the velocity divergence to bias, its PDF can be used to obtain unbiased constraints on the growth of structures ($\sigma_8$,f) or it can be combined with the galaxy density PDF to extract bias parameters.

In order to infer the impact of the small-scale physics to the large-scale properties of the universe, we use a series of cosmological $N$-body simulations of self-gravitating matter inhomogeneities to measure, for the first time, the response function of such a system defined as a functional derivative of the nonlinear power spectrum with respect to its linear counterpart. Its measured shape and amplitude are found to be in good agreement with perturbation theory predictions except for the coupling from small to large-scale perturbations. The latter is found to be significantly damped, following a Lorentzian form. These results shed light on validity regime of perturbation theory calculations giving a useful guideline for regularization of small scale effects in analytical modeling. Most importantly our result indicates that the statistical properties of the large-scale structure of the universe are remarkably insensitive to the details of the small-scale physics, astrophysical or gravitational, paving the way for the derivation of robust estimates of theoretical uncertainties on the determination of cosmological parameters from large-scale survey observations.

We present a first principle approach to obtain analytical predictions for spherically-averaged cosmic densities in the mildly non-linear regime that go well beyond what is usually achieved by standard perturbation theory. A large deviation principle allows us to compute the leading-order cumulants of average densities in concentric cells. In this symmetry, the spherical collapse model leads to cumulant generating functions that are robust for finite variances and free of critical points when logarithmic density transformations are implemented. They yield in turn accurate density probability distribution functions (PDFs) from a straightforward saddle-point approximation valid for all density values. Based on this easy-to-implement modification, explicit analytic formulas for the evaluation of the one- and two-cell PDF are provided. The theoretical predictions obtained for the PDFs are accurate to a few percent compared to the numerical integration, regardless of the density under consideration and in excellent agreement with N-body simulations for a wide range of densities. This formalism should prove valuable for accurately probing the quasi-linear scales of low redshift surveys for arbitrary primordial power spectra.

In order to quantify the error budget in the measured probability distribution functions of cell densities, the two-point statistics of cosmic densities in concentric spheres is investigated. Bias functions are introduced as the ratio of their two-point correlation function to the two-point correlation of the underlying dark matter distribution. They describe how cell densities are spatially correlated. They are computed here via the so-called large deviation principle in the quasi-linear regime. Their large-separation limit is presented and successfully compared to simulations for density and density slopes: this regime is shown to be rapidly reached allowing to get sub-percent precision for a wide range of densities and variances. The corresponding asymptotic limit provides an estimate of the cosmic variance of standard concentric cell statistics applied to finite surveys. More generally, no assumption on the separation is required for some specific moments of the two-point statistics, for instance when predicting the generating function of cumulants containing any powers of concentric densities in one location and one power of density at some arbitrary distance from the rest. This exact "one external leg" cumulant generating function is used in particular to probe the rate of convergence of the large-separation approximation.

The recently published analytic probability density function for the mildly non-linear cosmic density field within spherical cells is used to build a simple but accurate maximum likelihood estimate for the redshift evolution of the variance of the density, which, as expected, is shown to have smaller relative error than the sample variance. This estimator provides a competitive probe for the equation of state of dark energy, reaching a few percent accuracy on wp and wa for a Euclid-like survey. The corresponding likelihood function can take into account the configuration of the cells via their relative separations. A code to compute one-cell density probability density functions for arbitrary initial power spectrum, top-hat smoothing and various spherical collapse dynamics is made available online so as to provide straightforward means of testing the effect of alternative dark energy models and initial power-spectra on the low-redshift matter distribution.

Redshift-space distortions are generally considered in the plane parallel limit, where the angular separation between the two sources can be neglected. Given that galaxy catalogues now cover large fractions of the sky, it becomes necessary to consider them in a formalism which takes into account the wide angle separations. In this article we derive an operational formula for the matter correlators in the Newtonian limit to be used in actual data sets, both in configuration and in Fourier spaces without relying on a plane-parallel approximation. We then recover the plane-parallel limit not only in configuration space where the geometry is simpler, but also in Fourier space, and we exhibit the first corrections that should be included in large surveys as a perturbative expansion over the plane-parallel results. We finally compare our results to existing literature, and show explicitly how they are related.

We present an application of Large Deviation Theory to the problem of structure growth on large-scale structure cosmology. Starting from gaussian distributed overdensities on concentric spherical shells, we show that a Large Deviation Principle holds for the densities on the corresponding shells after gravitational evolution if no shell-crossing happens. As consequences of the Large Deviation Principle we obtain the cumulant generating function for the non-linear densities, and present formulae to compute the cumulant generating function for general window functions.

We analyse the clustering of cosmic large scale structure using a consistent modified gravity perturbation theory, accounting for anisotropic effects along and transverse to the line of sight. The growth factor has a particular scale dependence in f(R) gravity and we fit for the shape parameter f_{R0} simultaneously with the distance and the large scale (general relativity) limit of the growth function. Using more than 690,000 galaxies in the Baryon Oscillation Spectroscopy Survey Data Release 11, we find no evidence for extra scale dependence, with the 95\% confidence upper limit |f_{R0}| <8 \times 10^{-4}. Future clustering data, such as from the Dark Energy Spectroscopic Instrument, can use this consistent methodology to impose tighter constraints.

We propose a procedure to evaluate the impact of nonlinear couplings on the evolution of massive neutrino streams in the context of large-scale structure growth. Such streams can be described by general nonlinear conservation equations, derived from a multiple-flow perspective, which generalize the conservation equations of non-relativistic pressureless fluids. The relevance of the nonlinear couplings is quantified with the help of the eikonal approximation applied to the subhorizon limit of this system. It highlights the role played by the relative displacements of different cosmic streams and it specifies, for each flow, the spatial scales at which the growth of structure is affected by nonlinear couplings. We found that, at redshift zero, such couplings can be significant for wavenumbers as small as $k=0.2\,h$/Mpc for most of the neutrino streams.

In the context of count-in-cells statistics, the joint probability distribution of the density in two concentric spherical shells is predicted from first first principle for sigmas of the order of one. The agreement with simulation is found to be excellent. This statistics allows us to deduce the conditional one dimensional probability distribution function of the slope within under dense (resp. overdense) regions, or of the density for positive or negative slopes. The former conditional distribution is likely to be more robust in constraining the cosmological parameters as the underlying dynamics is less evolved in such regions. A fiducial dark energy experiment is implemented on such counts derived from Lambda-CDM simulations.

Motion equations describing streams of relativistic particles and their properties are explored in detail in the framework of Cosmological Perturbation Theory. Those equations, derived in any metric both in the linear and nonlinear regimes, express the matter and momentum conservation. In this context we extend the setup of adiabatic initial conditions - that was initially performed in the Conformal Newtonian gauge - to the Synchronous gauge. The subhorizon limit of the nonlinear motion equations written in a generic perturbed Friedmann-Lema\^{i}tre metric is then derived and analyzed. We show in particular that the momentum field $P_{i}(x)$ is always potential in the linear regime and remains so at subhorizon scales in the nonlinear regime. Finally the equivalence principle is exploited to highlight invariance properties satisfied by such a system of equations, extending that known for streams of non-relativistic particles, namely the extended Galilean invariance.

An explicit full nulling scheme for cosmic shear observations is presented. It makes possible the construction of shear maps from extended source distributions for which the lens distance distribution is restricted to a definite interval. Such a construction allows to build totally independent shear maps, at all scales, that can be taken advantage of to constrain background cosmological parameters and systematics using the full statistical power of cosmic shear observations. Another advantage of such construction is that, as the lens redshift distribution can be made arbitrarily narrow, scale mixing due to projection effects can be limited allowing controlled predictions on the large scale shear power spectrum from perturbation theory calculations.

Based on the multi-point propagator expansion, we present resummed perturbative calculations for cosmological power spectra and correlation functions in the context of modified gravity. In a wide class of modified gravity models that have a screening mechanism to recover general relativity (GR) on small scales, we apply the eikonal approximation to derive the governing equation for resummed propagator that partly includes the non-perturbative effect in the high-k limit. The resultant propagator in the high-k limit contains the new corrections arising from the screening mechanism as well as the standard exponential damping. We explicitly derive the expression for new high-k contributions in specific modified gravity models, and find that in the case of f(R) gravity for a currently constrained model parameter, the corrections are basically of the sub-leading order and can be neglected. Thus, in f(R) gravity, similarly to the GR case, we can analytically construct the regularized propagator that reproduces both the resummed high-k behavior and the low-k results computed with standard perturbation theory, consistently taking account of the nonlinear modification of gravity valid at large scales. With the regularized multi-point propagators, we give predictions for power spectrum and correlation function at one-loop order, and compare those with N-body simulations in f(R) gravity model. As an important application, we also discuss the redshift-space distortions and compute the anisotropic power spectra and correlation functions.

In the context of upcoming large-scale structure surveys such as Euclid, it is of prime importance to quantify the effect of peculiar velocities on geometric probes. Hence the formalism to compute in redshift space the geometrical and topological one-point statistics of mildly non-Gaussian 2D and 3D cosmic fields is developed. Leveraging the partial isotropy of the target statistics, the Gram-Charlier expansion of the joint probability distribution of the field and its derivatives is reformulated in terms of the corresponding anisotropic variables. In particular, the cosmic non-linear evolution of the Minkowski functionals, together with the statistics of extrema are investigated in turn for 3D catalogues and 2D slabs. The amplitude of the non-Gaussian redshift distortion correction is estimated for these geometric probes. In 3D, gravitational perturbation theory is implemented in redshift space to predict the cosmic evolution of all relevant Gram-Charlier coefficients. Applications to the estimation of the cosmic parameters sigma(z) and beta=f/b1 from upcoming surveys is discussed. Such statistics are of interest for anisotropic fields beyond cosmology.

In dimension 2 and above, the Burgers dynamics, the so-called "adhesion model" in cosmology, can actually give rise to several dynamics in the inviscid limit. We investigate here the statistical properties of the density field when it is defined by a "geometrical model" associated with this Burgers velocity field and where the matter distribution is fully determined, at each time step, by geometrical constructions. Our investigations are based on a set of numerical experiments that make use of an improved algorithm, for which the geometrical constructions are efficient and robust. In this work we focus on Gaussian initial conditions with power-law power spectra of slope $n$ in the range $-3<n<1$, where a self-similar evolution develops, and we compute the behavior of power spectra, density probability distributions and mass functions. As expected for such dynamics, the density power spectra show universal high-$k$ tails that are governed by the formation of pointlike masses. The two other statistical indicators however show the same qualitative properties as those observed for 3D gravitational clustering. In particular, the mass functions obey a Press-Schechter like scaling up to a very good accuracy in 1D, and to a lesser extent in 2D. Our results suggest that the "geometrical adhesion model", whose solution is fully known at all times, provides a precious tool to understand some of the statistical constructions frequently used to study the development of mass halos in gravitational clustering.

These lecture notes introduce analytical tools, methods and results describing the growth of cosmological structure beyond the linear regime. The presentation is focused on the single flow regime of the Vlasov-Poisson equation describing the development of gravitational instabilities in a pressureless fluid. These notes include the introduction of diagrammatic representations of the standard perturbation theory with applications to the calculation of the so-called loop contributions to the power spectra. A large part of these notes is devoted to the exploration of the convergence properties of these terms from the contribution of both the long-wave modes and the short-wave modes. The resulting performances of the two-loop corrections of the power spectra are then presented. Finally other avenues that use different methods are explored. In particular it is shown how joint density and profile probability distribution functions can be constructed out of the multiple-variable cumulant generating functions computed at tree order.

A new analytical approach allowing to account for massive neutrinos in the non-linear description of the growth of the large-scale structure of the universe is proposed. Unlike the standard approach in which neutrinos are described as a unique hot fluid, it is shown that the overall neutrino fluid can be equivalently decomposed into a collection of independent flows. Starting either from elementary conservation equations or from the evolution equation of the phase-space distribution function, we derive the two non-linear motion equations that each of these flows satisfies. Those fluid equations describe the evolution of macroscopic fields. We explain in detail the connection between the collection of flows we defined and the standard massive neutrino fluid. Then, in the particular case of adiabatic initial conditions, we explicitly check that, at linear order, the resolution of this new system of equations reproduces the results obtained in the standard approach based on the collisionless Boltzmann hierarchy. Besides, the approach advocated in this paper allows to show how each neutrino flow settles into the cold dark matter flow depending on initial velocities. It opens the way to a fully non-linear treatment of the dynamical evolution of neutrinos in the framework of large-scale structure growth.

The joint probability distribution function (PDF) of the density within multiple concentric spherical cells is considered. It is shown how its cumulant generating function can be obtained at tree order in perturbation theory as the Legendre transform of a function directly built in terms of the initial moments. In the context of the upcoming generation of large-scale structure surveys, it is conjectured that this result correctly models such a function for finite values of the variance. Detailed consequences of this assumption are explored. In particular the corresponding one-cell density probability distribution at finite variance is computed for realistic power spectra, taking into account its scale variation. It is found to be in agreement with $\Lambda$-CDM simulations at the few percent level for a wide range of density values and parameters. Related explicit analytic expansions at the low and high density tails are given. The conditional (at fixed density) and marginal probability of the slope -- the density difference between adjacent cells -- and its fluctuations is also computed from the two-cells joint PDF; it also compares very well to simulations, in particular in under-dense regions, with a significant reduced cosmic scatter compared to over-dense regions. It is emphasized that this could prove useful when studying the statistical properties of voids as it can serve as a statistical indicator to test gravity models and/or probe key cosmological parameters.

We use the so-called eikonal approximation, recently introduced in the context of cosmological perturbation theory, to compute power spectra for multi-component fluids. We demonstrate that, at any given order in standard perturbation theory, multipoint power spectra do not depend on the large-scale adiabatic modes. Moreover, we employ perturbation theories to decipher how nonadiabatic modes, such as a relative velocity between two different components, damp the small-scale matter power spectrum, a mechanism recently described in the literature. In particular, we do an explicit calculation at 1-loop order of this effect. While the 1-loop result eventually breaks down, we show how the damping effect can be fully captured by the help of the eikonal approximation. A relative velocity not only induces mode damping but also creates large-scale anisotropic modulations of the matter power spectrum amplitude. We illustrate this for the Local Group environment.

We introduce the eikonal approximation to study the effect of the large-scale motion of cosmic fluids on their small-scale evolution. This approach consists in collecting the impact of the long-wavelength displacement field into a single or finite number of random variables, whose statistical properties can be computed from the initial conditions. For a single dark matter fluid, we show that we can recover the nonlinear propagators of renormalized perturbation theory. These are obtained with no need to assume that the displacement field follows the linear theory. Then we extend the eikonal approximation to many fluids. In particular, we study the case of two non-relativistic components and we derive their resummed propagators in the presence of isodensity modes. Unlike the adiabatic case, where only the phase of small-scale modes is affected by the large-scale advection field, the isodensity modes change also the amplitude on small scales. We explicitly solve the case of cold dark matter-baryon mixing and find that the isodensity modes induce only very small corrections to the resummed propagators.

Using a full implementation of resummed perturbation theory (PT) from a multi-point propagator expansion, we put forward new theoretical predictions for the two-point statistics of matter fluctuations in redshift space. The predictions consistently include PT corrections up to the two-loop order and are based on an improved prescription of the redshift-space distortions (RSD) that properly takes into account their non-Gaussian impact from a systematic low-k expansion. In contrast to the previous studies that partly used standard PT calculations, the present treatment is able to provide a consistent prediction for both power spectra and correlation functions. These results are compared with $N$-body simulations with which a very good agreement is found up to the quadrupole moment. The theoretical predictions for the hexadecapole moment of the power spectra are however found to significantly departs from the numerical results at low redshift. We examine this issue and found it to be likely related to an improper modeling of the RSD damping effects on which this moment shows large dependence.

We explore the properties of two-point cosmic propagators when Perturbation Theory (PT) loop corrections are consistently taken into account. We show in particular how the interpolation scheme proposed in arXiv:1112.3895 can be explicitly used up to two-loop order introducing the notion of regular parts for the contributing terms. Extending the one-loop results, we then derive and give semi analytical forms of the two-loop contributions for both the cosmic density and velocity propagators. These results are tested against numerical simulations and shown to significantly improve upon one-loop results for redshifts above 0.5. We found however that at lower redshift two-loop order corrections are too large partly due to a strong sensitivity of those terms to the small scale modes. We show that this dependence is expected to be even larger for higher order loop corrections both from theoretical investigations and numerical tests, the latter obtained with Monte Carlo evaluations of the three-loop contributions. This makes small-scale regularization schemes necessary for the use of such higher order corrections.

We put forward and test a simple description of multi-point propagators (MP), which serve as building-blocks to calculate the nonlinear matter power spectrum. On large scales these propagators reduce to the well-known kernels in standard perturbation theory, while at smaller scales they are suppresed due to nonlinear couplings. Through extensive testing with numerical simulations we find that this decay is characterized by the same damping scale for both two and three-point propagators. In turn this transition can be well modeled with resummation results that exponentiate one-loop computations. For the first time, we measure the four components of the non-linear (two-point) propagator using dedicated simulations started from two independent random Gaussian fields for positions and velocities, verifying in detail the fundamentals of propagator resummation. We use these results to develop an implementation of the MP-expansion for the nonlinear power spectrum that only requires seconds to evaluate at BAO scales. To test it we construct six suites of large numerical simulations with different cosmologies. From these and LasDamas runs we show that the nonlinear power spectrum can be described at the ~ 2% level at BAO scales for redshifts in the range [0-2.5]. We make a public release of the MPTbreeze code with the hope that it can be useful to the community.

We present a specific prescription for the calculation of cosmological power spectra, exploited here at two-loop order in perturbation theory (PT), based on the multi-point propagator expansion. In this approach power spectra are constructed from the regularized expressions of the propagators that reproduce both the resummed behavior in the high-k limit and the standard PT results at low-k. With the help of N-body simulations, we show that such a construction gives robust and accurate predictions for both the density power spectrum and the correlation function at percent-level in the weakly non-linear regime. We then present an algorithm that allows accelerated evaluations of all the required diagrams by reducing the computational tasks to one-dimensional integrals. This is achieved by means of pre-computed kernel sets defined for appropriately chosen fiducial models. The computational time for two-loop results is then reduced from a few minutes, with the direct method, to a few seconds with the fast one. The robustness and applicability of this method are tested against the power spectrum cosmic emulator from which a wide variety of cosmological models can be explored. The fortran program with which direct and fast calculations of power spectra can be done, RegPT, is publicly released as part of this paper.

We present a new scheme for the general computation of cosmic propagators that allow to interpolate between standard perturbative results at low-k and their expected large-k resummed behavior. This scheme is applicable to any multi-point propagator and allows the matching of perturbative low-k calculations to any number of loops to their large-k behavior, and can potentially be applied in case of non-standard cosmological scenarios such as those with non-Gaussian initial conditions. The validity of our proposal is checked against previous prescriptions and measurements in numerical simulations showing a remarkably good agreement. Such a generic prescription for multi-point propagators provides the necessary building blocks for the computation of polyspectra in the context of the so-called Gamma-expansion introduced by Bernardeau et al. (2008). As a concrete application we present a consistent calculation of the matter bispectrum at one-loop order.

Future lensing surveys will be nearly full-sky and reach an unprecedented depth, probing scales closer and closer to the Hubble radius. This motivates the study of the cosmic shear beyond the small-angle approximation and including general relativistic corrections that are usually suppressed on sub-Hubble scales. The complete expression of the reduced cosmic shear at second order including all relativistic effects was derived in [1]. In the present paper we compute the resulting cosmic shear bispectrum when all these effects are properly taken into account and we compare it to primordial non-Gaussianity of the local type. The new general relativistic effects are generically smaller than the standard non-linear couplings. However, their relative importance increases at small multipoles and for small redshifts of the sources. The dominant effect among these non standard corrections is due to the inhomogeneity of the source redshift. In the squeezed limit, its amplitude can become of the order of the standard couplings when the redshift of the sources is below 0.5. Moreover, while the standard non-linear couplings depend on the angle between the short and long mode, the relativistic corrections do not and overlap almost totally with local type non-Gaussianity. We find that they can contaminate the search for a primordial local signal by f_NL>10.

The Geometrical Adhesion Model (GAM) we described in previous papers provides a fully solved model for the nonlinear evolution of fields that mimic the cosmological evolution of pressureless fluids. In this context we explore the expected late time properties of the cosmic propagators once halos have formed, in a regime beyond the domain of application of perturbation theories. Whereas propagators in Eulerian coordinates are closely related to the velocity field we show here that propagators defined in Lagrangian coordinates are intimately related to the halo mass function. Exact results can be obtained in the 1D case. In higher dimensions, the computations are more intricate because of to the dependence of the propagators on the detailed shape of the underlying Lagrangian-space tessellations, that is, on the geometry of the regions that eventually collapse to form halos. We illustrate these results for both the 1D and the 2D dynamics. In particular we give here the expected asymptotic behaviors obtained for power-law initial power spectra. These analytical results are compared with the results obtained with dedicated numerical simulations.

The non-linear evolution of the energy density of the radiation is shown to induce spectral distortions of the cosmic microwave background both at recombination and during the reionization era. This distortion has the same spectral signature as the one produced by the re-scattering of photons by non-relativistic hot electrons, the thermal Sunyaev-Zeldovich effect, whose amplitude is quantified by a Compton y parameter. A diffuse y-sky is then expected to emerge from mode couplings in the non-linear evolution of the cosmological perturbations and to superimpose to the point source contributions of galaxy clusters. The equations describing the evolution of the y field and a hierarchy governing its angular multipoles are derived from the second order Boltzmann equation. These equations are then integrated numerically to obtain the first predicted power spectrum of the diffuse y-sky. It is found to be a remarkable tracer of the reionization history of the Universe.

We consider the effect of modified gravity on the growth of large-scale structures at second order in perturbation theory. We show that modified gravity models changing the linear growth rate of fluctuations are also bound to change, although mildly, the mode coupling amplitude in the density and reduced velocity fields. We present explicit formulae which describe this effect. We then focus on models of modified gravity involving a scalar field coupled to matter, in particular chameleons and dilatons, where it is shown that there exists a transition scale around which the existence of an extra scalar degree of freedom induces significant changes in the coupling properties of the cosmic fields. We obtain the amplitude of this effect for realistic dilaton models at the tree-order level for the bispectrum, finding them to be comparable in amplitude to those obtained in the DGP model.

This article constructs flat-sky approximations in a controlled way in the context of the cosmic microwave background observations for the computation of both spectra and bispectra. For angular spectra, it is explicitly shown that there exists a whole family of flat-sky approximations of similar accuracy for which the expression and amplitude of next to leading order terms can be explicitly computed. It is noted that in this context two limiting cases can be encountered for which the expressions can be further simplified. They correspond to cases where either the sources are localized in a narrow region (thin-shell approximation) or are slowly varying over a large distance (which leads to the so-called Limber approximation). Applying this to the calculation of the spectra it is shown that, as long as the late integrated Sachs-Wolfe contribution is neglected, the flat-sky approximation at leading order is accurate at 1% level for any multipole. Generalization of this construction scheme to the bispectra led to the introduction of an alternative description of the bispectra for which the flat-sky approximation is well controlled. This is not the case for the usual description of the bispectrum in terms of reduced bispectrum for which a flat-sky approximation is proposed but the next-to-leading order terms of which remain obscure.

The evolution of high order correlation functions of a test scalar field in arbitrary inflationary backgrounds is computed. Whenever possible, exact results are derived from quantum field theory calculations. Taking advantage of the fact that such calculations can be mapped, for super-horizon scales, into those of a classical system, we express the expected correlation functions in terms of classical quantities, power spectra, Green functions, that can be easily computed in the long-wavelength limit. Explicit results are presented that extend those already known for a de Sitter background. In particular the expressions of the late time amplitude of bispectrum and trispectrum, as well as the whole high-order correlation structure, are given in terms of the expansion factor behavior. When compared to the case of a de Sitter background, power law inflation and chaotic inflation induced by a massive field are found to induce high order correlation functions the amplitudes of which are amplified by almost one order of magnitude. These results indicate that the dependence of the related non-Gaussian parameters - such as f_NL - on the wave-modes is at percent level.

We explore the noiseless Burgers dynamics in the inviscid limit, the so-called ``adhesion model'' in cosmology, in a regime where (almost) all the fluid particles are embedded within point-like massive halos. Following previous works, we focus our investigations on a ``geometrical'' model, where the matter evolution within the shock manifold is defined from a geometrical construction. This hypothesis is at variance with the assumption that the usual continuity equation holds but, in the inviscid limit, both models agree in the regular regions. Taking advantage of the formulation of the dynamics of this ``geometrical model'' in terms of Legendre transforms and convex hulls, we study the evolution with time of the distribution of matter and the associated partitions of the Lagrangian and Eulerian spaces. We describe how the halo mass distribution derives from a triangulation in Lagrangian space, while the dual Voronoi-like tessellation in Eulerian space gives the boundaries of empty regions with shock nodes at their vertices. We then emphasize that this dynamics actually leads to halo fragmentations for space dimensions greater or equal to 2 (for the inviscid limit studied in this article). This is most easily seen from the properties of the Lagrangian-space triangulation and we illustrate this process in the two-dimensional (2D) case. In particular, we explain how point-like halos only merge through three-body collisions while two-body collisions always give rise to two new massive shock nodes (in 2D). This generalizes to higher dimensions and we briefly illustrate the three-dimensional (3D) case. This leads to a specific picture for the continuous formation of massive halos through successive halo fragmentations and mergings.

This article presents the first computation of the complete bispectrum of the cosmic microwave background temperature anisotropies arising from the evolution of all cosmic fluids up to second order, including neutrinos. Gravitational couplings, electron density fluctuations and the second order Boltzmann equation are fully taken into account. Comparison to limiting cases that appeared previously in the literature are provided. These are regimes for which analytical insights can be given. The final results are expressed in terms of equivalent fNL for different configurations. It is found that for moments up to lmax=2000, the signal generated by non-linear effects is equivalent to fNL~5 for both local-type and equilateral-type primordial non-Gaussianity.

We show here how Renormalized Perturbation Theory (RPT) calculations applied to the quasi-linear growth of the large-scale structure can be carried on in presence of primordial non-Gaussian (PNG) initial conditions. It is explicitly demonstrated that the series reordering scheme proposed in Bernardeau, Crocce and Scoccimarro (2008) is preserved for non-Gaussian initial conditions. This scheme applies to the power spectrum and higher order spectra and is based on a reorganization of the contributing terms into sum of products of multi-point propagators. In case of PNG new contributing terms appear, the importance of which is discussed in the context of current PNG models. The properties of the building blocks of such resummation schemes, the multi-point propagators, are then investigated. It is first remarked that their expressions are left unchanged at one-loop order irrespectively of statistical properties of the initial field. We furthermore show that the high-momemtum limit of each of these propagators can be explicitly computed even for arbitrary initial conditions. They are found to be damped by an exponential cutoff whose expression is directly related to the moment generating function of the one-dimensional displacement field. This extends what had been established for multi-point propagators for Gaussian initial conditions. Numerical forms of the cut-off are shown for the so-called local model of PNG.

Mechanisms for the generation of primordial non-Gaussian metric fluctuations in the context of multiple-field inflation are reviewed. As long as kinetic terms remain canonical, it appears that nonlinear couplings inducing non-gaussianities can be split into two types. The extension of the one-field results to multiple degrees of freedom leads to gravity mediated couplings that are ubiquitous but generally modest. Multiple-field inflation offers however the possibility of generating non-gravity mediated coupling in isocurvature directions that can eventually induce large non-Gaussianities in the metric fluctuations. The robustness of the predictions of such models is eventually examined in view of a case study derived from a high-energy physics construction.

We compute the reduced cosmic shear up to second order in the gravitational potential without relying on the small angle or thin-lens approximation. This is obtained by solving the Sachs equation which describes the deformation of the infinitesimal cross-section of light bundle in the optical limit, and maps galaxy intrinsic shapes into their angular images. The calculation is done in the Poisson gauge without a specific matter content, including vector and tensor perturbations generated at second order and taking account of the inhomogeneities of a fixed redshift source plane. Our final result is expressed in terms of spin-2 operators on the sphere and is valid on the full sky. Beside the well known lens-lens and Born corrections that dominate on small angular scales, we find new non-linear couplings. These are a purely general relativistic intrinsic contribution, a coupling between the gravitational potential at the source with the lens, couplings between the time delay with the lens, couplings between two photon deflections, as well as non-linear couplings due to the second-order vector and tensor components. The inhomogeneity in the redshift of the source induces a coupling between the photon redshift with the lens. All these corrections become important on large angular scales and should thus be included when computing higher-order observables such as the bispectrum, in full or partially full-sky surveys.

Motivated by theoretical studies of gravitational clustering in the Universe, we compute propagators (response functions) in the adhesion model. This model, which is able to reproduce the skeleton of the cosmic web and includes nonlinear effects in both Eulerian and Lagrangian frameworks, also corresponds to the Burgers equation of hydrodynamics. Focusing on the one-dimensional case with power-law initial conditions, we obtain exact results for Eulerian and Lagrangian propagators. We find that Eulerian propagators can be expressed in terms of the one-point velocity probability distribution and show a strong decay at late times and high wavenumbers, interpreted as a "sweeping effect" but not a genuine damping of small-scale structures. By contrast, Lagrangian propagators can be written in terms of the shock mass function -- which would correspond to the halo mass function in cosmology -- and saturate to a constant value at late times. Moreover, they show a power-law dependence on scale or wavenumber which depends on the initial power-spectrum index and is directly related to the low-mass tail of the shock mass function. These results strongly suggest that Lagrangian propagators are much more sensitive probes of nonlinear structures in the underlying density field and of relaxation processes than their Eulerian counterparts.

We introduce the concept of multi-point propagators between linear cosmic fields and their nonlinear counterparts in the context of cosmological perturbation theory. Such functions express how a non-linearly evolved Fourier mode depends on the full ensemble of modes in the initial density field. We identify and resum the dominant diagrams in the large-$k$ limit, showing explicitly that multi-point propagators decay into the nonlinear regime at the same rate as the two-point propagator. These analytic results generalize the large-$k$ limit behavior of the two-point propagator to arbitrary order. We measure the three-point propagator as a function of triangle shape in numerical simulations and confirm the results of our high-$k$ resummation. We show that any $n-$point spectrum can be reconstructed from multi-point propagators, which leads to a physical connection between nonlinear corrections to the power spectrum at small scales and higher-order correlations at large scales. As a first application of these results, we calculate the reduced bispectrum at one-loop in renormalized perturbation theory and show that we can predict the decrease in its dependence on triangle shape at redshift zero, when standard perturbation theory is least successful.

This article investigates the non-linear evolution of cosmological perturbations on sub-Hubble scales in order to evaluate the unavoidable deviations from Gaussianity that arise from the non-linear dynamics. It shows that the dominant contribution to modes coupling in the cosmic microwave background temperature anisotropies on small angular scales is driven by the sub-Hubble non-linear evolution of the dark matter component. The perturbation equations, involving in particular the first moments of the Boltzmann equation for the photons, are integrated up to second order in perturbations. An analytical analysis of the solutions gives a physical understanding of the result as well as an estimation of its order of magnitude. This allows to quantify the expected deviation from Gaussianity of the cosmic microwave background temperature anisotropy and, in particular, to compute its bispectrum on small angular scales. Restricting to equilateral configurations, we show that the non-linear evolution accounts for a contribution that would be equivalent to a constant primordial non-Gaussianity of order fNL~25 on scales ranging approximately from l~1000 to l~3000.

It has been found recently that propagators, e.g. the cross-correlation spectra of the cosmic fields with the initial density field, decay exponentially at large-k in an Eulerian description of the dynamics. We explore here similar quantities defined for a Lagrangian space description. We find that propagators in Lagrangian space do not exhibit the same properties: they are found not to be monotonic functions of time, and to track back the linear growth rate at late time (but with a renormalized amplitude). These results have been obtained with a novel method which we describe alongside. It allows the formal resummation of the same set of diagrams as those that led to the known results in Eulerian space. We provide a tentative explanation for the marked differences seen between the Eulerian and the Lagrangian cases, and we point out the role played by the vorticity degrees of freedom that are specific to the Lagrangian formalism. This provides us with new insights into the late-time behavior of the propagators.

The contribution of cosmological perturbations to the time drift of the cosmological redshift is derived. It is shown that the dominant correction arises from the local acceleration of both the emitter and the observer. The amplitude of this effect is estimated to be of the order of 1% of the drift signal at z=2-4, but can easily be lowered down to 0.1% by using many absorption lines and quasars.

Lagrangian reconstruction of large-scale peculiar velocity fields can be strongly affected by observational biases. We develop a thorough analysis of these systematic effects by relying on specially selected mock catalogues. For the purpose of this paper, we use the MAK reconstruction method, although any other Lagrangian reconstruction method should be sensitive to the same problems. We extensively study the uncertainty in the mass-to-light assignment due to luminosity incompleteness, and the poorly-determined relation between mass and luminosity. The impact of redshift distortion corrections is analyzed in the context of MAK and we check the importance of edge and finite-volume effects on the reconstructed velocities. Using three mock catalogues with different average densities, we also study the effect of cosmic variance. In particular, one of them presents the same global features as found in observational catalogues that extend to 80 Mpc/h scales. We give recipes, checked using the aforementioned mock catalogues, to handle these particular observational effects, after having introduced them into the mock catalogues so as to quantitatively mimic the most densely sampled currently available galaxy catalogue of the nearby universe. Once biases have been taken care of, the typical resulting error in reconstructed velocities is typically about a quarter of the overall velocity dispersion, and without significant bias. We finally model our reconstruction errors to propose an improved Bayesian approach to measure Omega_m in an unbiased way by comparing the reconstructed velocities to the measured ones in distance space, even though they may be plagued by large errors. We show that, in the context of observational data, a nearly unbiased estimator of Omega_m may be built using MAK reconstruction.

We explore extensions of hybrid inflationary models in the context of supersymmetric D-term inflation. We point out that a large variety of inflationary scenarios can be encountered when the field content is extended. It is not only possible to get curvaton type models but also scenarios in which different fields, with nontrivial statistical properties, contribute to the primordial curvature fluctuations. We explore more particularly the parameter space of these multiple field inflationary models. It is shown that there exists a large domain in which significant primordial non-Gaussianities can be produced while preserving a scale free power spectrum for the metric fluctuations. In particular we explicitly compute the expected bi- and trispectrum for such models and compared the results to the current and expected observational constraints. It is shown that it is necessary to use both the bi- and tri-spectra of CMB anisotropies to efficiently reduce their parameter space.

We present a class of models in which the primordial metric fluctuations do not necessarily obey Gaussian statistics. These models are realizations of mechanisms in which non-Gaussianity is first generated by a light scalar field and then transferred into curvature fluctuations during or at the end of inflation. For this class of models we present generic results for the probability distribution functions of the metric perturbation at the end of inflation. It is stressed that finite volume effects can induce non trivial effects that we sketch.

We present a phenomenological investigation of non-Gaussian effects that could be seen on CMB temperature maps. Explicit expressions for the temperature correlation functions are given for different types of primordial mode couplings. We argue that a simplified description of the radial transfer function for the temperature anisotropies allows to get insights into the general properties of the bi and tri-spectra. The accuracy of these results is explored together with the use of the small scale approximation to get explicit expressions of high order spectra. The bi-spectrum is found to have alternate signs for the successive acoustic peaks. Sign patterns for the trispectra are more complicated and depend specifically on the type of metric couplings. Local primordial couplings are found to give patterns that are different from those expected from weak lensing effects.

We present a method for calculating explicit expressions of the shear three-point function for various cosmological models. The method is applied here to the one-halo model in case of power law density profiles for which results are detailed. The three-point functions are found to reproduce to a large extent patterns in the shear correlations obtained in numerical simulations and may serve as a guideline to implement optimized methods for detecting the shear three-point function. In principle, the general method presented here can also be applied for other models of matter correlation.

We compute the radiative corrections to the mass of a test boson field in an inflating space-time. The calculations are carried out in case of a boson part of a supersymmetric chiral multiplet. We show that its mass is preserved up to logarithmic divergences both in ultraviolet and infrared domains. Consequences of these results for inflationary models are discussed.

With the developments of large galaxy surveys or cosmic shear surveys it is now possible to map the dark matter distribution at truly cosmological scales. Detailed examinations of the statistical properties of the dark matter distribution reveal the detail of the large-scale structure growth of the Universe. In particular it is shown here that the behavior of the density field bi-spectrum is sensitive to departure from normal gravity in a way which depends only weakly on the background evolution. The cosmic shear bispectrum appears to be particularly sensitive to changes in the Poisson equation: we show that the current cosmic shear data can already be used to infer constraints on the scale of a possible higher-dimensional gravity, above 2 h^{-1}Mpc.

We investigate a possible asymmetry in the statistical properties of the cosmic microwave background temperature field and to do so we construct an estimator aiming at detecting a dipolar modulation. Such a modulation is found to induce correlations between multipoles with $\Delta\ell=1$. Applying this estimator, to the V and W bands of the WMAP data, we found a significant detection in the V band. We argue however that foregrounds and in particular point sources are the origin of this signal.

Inflation universally produces classical almost scale free Gaussian inhomogeneities of any light scalars. Assuming the coupling constants at the time of inflation depend on some light moduli fields, we encounter the generation of modulated cosmological fluctuations from (p)reheating. This is an alternative mechanism to generate observable (almost) scale free adiabatic metric perturbations. We extend this idea to the class of hybrid inflation, where the bifurcation value of the inflaton is modulated by the spatial inhomogeneities of the couplings. As a result, the symmetry breaking after inflation occurs not simultaneously in space but with the time laps in different Hubble patches inherited from the long-wavelength moduli inhomogeneities. To calculate modulated fluctuations we introduce techniques of general relativistic matching conditions for metric perturbations at the time hypersurface where the equation of state after inflation undergoes a jump, without evoking the detailed microscopic physics, as far as it justifies the jump. We apply this theory to the modulated fluctuations from the hybrid and chaotic inflations. We discuss what distinguish the modulated from the inflation-driven fluctuations, in particular, their spectral index, modification of the consistency relation and the issue of weak non-Gaussianity.

We present the expressions of the three- and four-point correlation functions of a self interacting light scalar field in a de Sitter spacetime at tree order respectively for a cubic and a quartic potential. Exact expressions are derived and their limiting behaviour on super-horizon scales are presented. Their essential features are shown to be similar to those obtained in a classical approach.

Models of multi-field inflation exhibiting primordial non-Gaussianity have recently been introduced. This is the case in particular if the fluctuations of a light field scalar field, transverse to the inflaton direction, with quartic coupling can be transferred to the metric fluctuations. So far in those calculations only the ensemble statistical properties have been considered. We explore here how finite volume effects could affect those properties. We show that the expected non-Gaussian properties survive at a similar level when the finite volume effects are taken into account and also find that they can skew the metric distribution even though the ensemble distribution is symmetric.

We present an overview of the cosmic shear results obtained so far. We focus on the 2-point statistics only. Evidences supporting the cosmological origin of the measured signal are reviewed, and issues related to various systematics are discussed.

We present the current status of cosmic shear based on all surveys done so far. Taken together, they cover more about 70 deg$^2$ and concern more than 3 million galaxies with accurate shape measurement. Theoretical expectations, observational results and their cosmological interpretations are discussed in the framework of standard cosmology and CDM scenarios. The potentials of the next generation cosmic shear surveys are discussed.

The present status of weak lensing analyses of clusters of galaxies and of cosmic shear surveys are presented and discussed. We focus on the impact of very large telescopes on present-day and future surveys and compare their potential with HST or wide field 4 meter telescopes.

We construct explicit models of multi-field inflation in which the primordial metric fluctuations do not necessarily obey Gaussian statistics. These models are realizations of mechanisms in which non-Gaussianity is first generated by a light scalar field and then transferred into curvature fluctuations. The probability distribution functions of the metric perturbation at the end of inflation are computed. This provides a guideline for designing strategies to search for non-Gaussian signals in future CMB and large scale structure surveys.

We have detected non-Gaussian signatures in the VIRMOS-DESCART weak lensing survey from a measurement of the three-point shear correlation function, following the method developed by Bernardeau et al. (2002). We obtain a 2.4-sigma signal over four independent angular bins, or equivalently, a 4.9-sigma confidence level detection with respect to measurements errors on scale of about 2 to 4 arcmin. The amplitude and the shape of the signal are consistent with theoretical expectations obtained from ray-tracing simulations. This result supports the idea that the measure corresponds to a cosmological signal due to the gravitational instability dynamics. Its properties could be used to put constraints on the cosmological parameters, in particular on the density parameter of the Universe, but the error level as well as the cosmic variance are still too large to permit secure conclusions.

We explore the scale and angular dependence of the cosmic shear three-point correlation function. The shear field is found to have a much more complex three-point function than the convergence field, but it also exhibits specific motifs which draw signatures of gravitational lensing. Exact shear patterns are inferred analytically for some geometrical shear triplets configurations, when simple interpretations can be derived. A more complete description of their geometrical properties has then been carried out from ray-tracing numerical simulations and simple models. These patterns reveal typical features that can be utilized for non-gaussian signal detection in cosmic shear surveys. We test the robustness of these properties against complex noise statistics and non-trivial survey topologies. From these conclusive checks, we predict that the VIRMOS-DESCART survey should allow a detection of a non-gaussian signal with a comfortable significance for a low matter density Universe.

This article investigates the generation of non-Gaussianity during inflation. In the context of multi-field inflation, we detail a mechanism that can create significant primordial non-Gaussianities in the adiabatic mode while preserving the scale invariance of the power spectrum. This mechanism is based on the generation of non-Gaussian isocurvature fluctuations which are then transfered to the adiabatic modes through a bend in the classical inflaton trajectory. Natural realizations involve quartic self-interaction terms for which a full computation can be performed. The expected statistical properties of the resulting metric fluctuations are shown to be the superposition of a Gaussian and a non-Gaussian contribution of the same variance. The relative weight of these two contributions is related to the total bending in field space. We explicit the non-Gaussian probability distribution function which appears to be described by a single new parameter. Only two new parameters therefore suffice in describing the non-Gaussianity.

We present the current status of cosmic shear studies and their implications on cosmological models. Theoretical expectations and observational results are discussed in the framework of standard cosmology and CDM scenarios. The potentials of the next generation cosmic shear surveys are discussed.

We review the formalism and applications of non-linear perturbation theory (PT) to understanding the large-scale structure of the Universe. We first discuss the dynamics of gravitational instability, from the linear to the non-linear regime. This includes Eulerian and Lagrangian PT, non-linear approximations, and a brief description of numerical simulation techniques. We then cover the basic statistical tools used in cosmology to describe cosmic fields, such as correlations functions in real and Fourier space, probability distribution functions, cumulants and generating functions. In subsequent sections we review the use of PT to make quantitative predictions about these statistics according to initial conditions, including effects of possible non Gaussianity of the primordial fields. Results are illustrated by detailed comparisons of PT predictions with numerical simulations. The last sections deal with applications to observations. First we review in detail practical estimators of statistics in galaxy catalogs and related errors, including traditional approaches and more recent developments. Then, we consider the effects of the bias between the galaxy distribution and the matter distribution, the treatment of redshift distortions in three-dimensional surveys and of projection effects in angular catalogs, and some applications to weak gravitational lensing. We finally review the current observational situation regarding statistics in galaxy catalogs and what the future generation of galaxy surveys promises to deliver.

The correlation between source galaxies and lensing potentials causes a systematic effect on measurements of cosmic shear statistics, known as the source-lens clustering (SLC) effect. The SLC effect on the skewness of lensing convergence, $S_3$, is examined using a nonlinear semi-analytic approach and is checked against numerical simulations. The semi-analytic calculations have been performed in a wide variety of generic models for the redshift distribution of source galaxies and power-law models for the bias parameter between the galaxy and dark matter distributions. The semi-analytic predictions are tested successfully against numerical simulations. We find the relative amplitude of the SLC effect on $S_3$ to be of the order of five to forty per cent. It depends significantly on the redshift distribution of sources and on the way the bias parameter evolves. We discuss possible measurement strategies to minimize the SLC effects.

We explore the possibility of putting constraints on quintessence models with large-scale structure observations. In particular we compute the linear and second order growth rate of the fluctuations in different flavors of quintessence scenarios. We show that effective models of quintessence (e.g. with a constant equation of state) do not account of the results found in more realistic scenarios. The impact of these results on observational quantities such as the shape of the non-linear power spectrum in weak lensing surveys or the skewness of the convergence field are investigated. It appears that the observational signature of quintessence models are specific and rather large. They clearly cannot be mistaken with a change of (\Omega_{0}).

We present a general method for the reconstruction of the one-point Probability Distribution Function of the local aperture mass in weak lensing maps. Exact results, that neglect the lens-lens coupling and departure form the Born approximation, are derived for both the quasilinear regime at leading order and the strongly nonlinear regime assuming the tree hierarchical model is valid. We describe in details the projection effects on the properties of the PDF and the associated generating functions. In particular, we show how the generic features which are common to both the quasilinear and nonlinear regimes lead to two exponential tails for $P(\Map)$. We briefly investigate the dependence of the PDF with cosmology and with the shape of the angular filter. Our predictions are seen to agree reasonably well with the results of numerical simulations and should be able to serve as foundations for alternative methods to measure the cosmological parameters that take advantage of the full shape of the PDF.

We report a measurement of cosmic shear correlations using an effective area of 6.5 sq. deg. of the VIRMOS deep imaging survey in progress at the Canada-France-Hawaii Telescope. We measured various shear correlation functions, the aperture mass statistic and the top-hat smoothed variance of the shear with a detection significance exceeding 12 sigma for each of them. We present results on angular scales from 3 arc-seconds to half a degree. The consistency of different statistical measures is demonstrated and confirms the lensing origin of the signal through tests that rely on the scalar nature of the gravitational potential. For Cold Dark Matter models we find $\sigma_8 \Omega_0^{0.6}=0.43^{+0.04}_{-0.05}$ at the 95% confidence level. The measurement over almost three decades of scale allows to discuss the effect of the shape of the power spectrum on the cosmological parameter estimation. The degeneracy on sigma_8-Omega_0 can be broken if priors on the shape of the linear power spectrum (that can be parameterized by Gamma) are assumed. For instance, with Gamma=0.21 and at the 95% confidence level, we obtain 0.6<sigma_8<1.1 and 0.2<Omega_0<0.5 for open models, and sigma_8>0.65 and Omega_0<0.4 for flat (Lambda-CDM) models. From the tangential/radial modes decomposition we can set an upper limit on the intrinsic shape alignment, which was recently suggested as a possible contribution to the lensing signal. Within the error bars, there is no detection of intrinsic shape alignment for scales larger than 1'.

Gravitational weak shear produced by large-scale structures of the universe induces a correlated ellipticity distribution of distant galaxies. The amplitude and evolution with angular scale of the signal depend on cosmological models and can be inverted in order to constrain the power spectrum and the cosmological parameters. We present our recent analysis of 50 uncorrelated VLT fields and the very first constrains on ($\Omega_m,\sigma_8$) and the nature of primordial fluctuations based on the join analysis of present-day cosmic shear surveys.

The recent detections of cosmic shear signal announced by several groups have demonstrated the feasibility of this challenging program and convinced astro- nomers of its potential for cosmology. Cosmic shear analysis demands to handle Gigabytes of data in order to probe several square degrees in subarcsecond deep imaging mode. The success of these surveys is sensitive to the designs of the observation strategy, the organization of the data reduction pipelines and the links of the data base with surveys like X-ray or spectroscopic follow up. We describe the cosmic shear surveys we have carried out at the VLT and at CFHT and the way we handle this huge data set in a more general context including the VIRMOS and the XMM-LSS surveys, and the future CMB surveys.

Recent developments in gravitational lensing astronomy have paved the way to genuine mappings of the gravitational potential at cosmological scales. We stress that comparing these data with traditional large scale structure surveys will provide us with a test of gravity at such scales. These constraints could be of great importance in the framework of higher dimensional cosmological models.

We observed with the camera FORS1 on the VLT (UT1, ANTU) 50 randomly selected fields and analyzed the cosmic shear inside circular apertures with diameter ranging from 0.5 to 5.0 arc-minutes. The images were obtained in optimal conditions by using the Service Observing proposed by ESO on the VLT, which enabled us to get a well-defined and homogeneous set of data. The 50 fields cover a 0.65 square-degrees area spread over more than 1000 square-degrees which provides a sample ideal for minimizing the cosmic variance. Using the same techniques as in Van Waerbeke et al. 2000, we measured the cosmic shear signal and investigated the systematics of the VLT sample. We find a significant excess of correlations between galaxy ellipticities on those angular scales. The amplitude and the shape of the correlation as function of angular scale are remarkably similar to those reported so far. Using our combined VLT and CFHT data and adding the results published by other teams we put the first joint constraints on Omega_0 and sigma_8 using cosmic shear surveys. From a deduced average of the redshift of the sources the combined data are consistent with sigma_8\simeq (0.59 \pm 0.03) Omega_0^{-0.47} (for a CDM power spectrum and $\Gamma=0.21$), in excellent agreement with the cosmological constraints obtained from the local cluster abundance.

The first detections of cosmic shear signal reported recently by 4 independent groups cover angular scales between one and 10 arcmin. On those scales, the cosmic shear is a signature of non-linear perturbations, like groups and clusters of galaxies. We present the results obtained by our team on CFHT and on the VLT and discuss its impact for the analysis of cluster abundances and cosmology

As pointed out in previous studies, the measurement of the skewness of the convergence field $\kappa$ will be useful in breaking the degeneracy among the cosmological parameters constrained from weak lensing observations. The combination of shot noise and finite survey volume implies that such a measurement is likely to be done in a range of intermediate scales ($0.5'$ to 20') where neither perturbation theory nor the hierarchical ansatz apply. Here we explore the behavior of the skewness of $\kappa$ at these intermediate scales, based on results for the non-linear evolution of the mass bispectrum. We combined different ray-tracing simulations to test our predictions, and we find that our calculations describe accurately the transition from the weakly non-linear to the strongly non-linear regime. We show that the single lens-plane approximation remains accurate even in the non-linear regime, and we explicitly calculate the corrections to this approximation. We also discuss the prospects of measuring the skewness in upcoming weak lensing surveys.

We explore the stability of the variance and skewness of the cosmic gravitational convergence field, using two different approaches: first we simulate a whole MEGACAM survey (100 sq. degrees). The reconstructed mass map, obtained from a shear map, shows that the state-of-the-art data analysis methods can measure weak-lensing statistics at angular scales ranging from 2.5' to 25'. We looked also at the influence of a varying signal-to-noise ratio over the shear map (due to local variations of source density) on the mass reconstruction, by means of Monte-Carlo simulation. The effect at small scales can easily be corrected-for in most of the relevant cases. These results enhance the confidence in the capability of future large surveys to measure accurately cosmologically interesting quantities.

We report on the detection of cosmic shear on angular scales of 1.3-6.5 arcmin using 45 independent empty fields observed with the Very Large Telescope (VLT). This result confirms previous measurements obtained with the CFH Telescope at the same angular scales. We present the data analysis and the first preliminary results.

We present results of our cosmic shear survey in progress at the Canada- France-Hawaii Telescope (CFHT). We analysed 1.7 sq. degrees of data, out of which we measured a significant correlation of galaxy shape-orientation over several arcmin scale. We present measurements of the variance of the shear $<\gamma^2>$ and of the correlation functions $<e_r(\theta)e_r(0)>$, $<e_t(\theta)e_t(0)>$. The signal is consistent with gravitational lensing by large scale structures predictions. The level of residual systematics after PSF correction is small compared to our measured signal. We outline several evolutions of our work by using additional data covering 9 sq. degrees.

We present investigations of lens phenomenological properties of cosmic strings for deep galaxy surveys. General results that have obtained for lineic energy distribution are presented first. We stress that generically the local convergence always vanishes in presence of strings although there might be some significant distortions. We then propose a simplified model of strings, we call ``Poisson strings'', for which exhaustive investigations can be done either numerically or analytically.

We reconsider the general properties of gravitational lensing effects induced by cosmic string systems, taking into account their equation of state and motion equations. We explicitly show that the deflection patterns induced by a string is equivalent to the one of a lineic and static distribution of matter. We then rigorously show that the convergence part of the deformation field is always zero, except on the intersection of the string worldsheet and the observer past light cone, extending previous results obtained in peculiar cases. Phenomenological consequences of this result on multiple image systems are investigated.

We present a novel approach for investigating lens phenomenology of cosmic strings in order to elaborate detection strategies in galaxy deep field images. To account for the complexity of the projected energy distribution of string networks we assume their lens effects to be similar to those of a straight string carrying a {\em random} lineic energy distribution. In such a model we show that, unlike the case of uniform strings, critical phenomena naturally appear. We explore the properties of the critical lines and caustics. In particular, assuming that the energy coherence length along the string is much smaller than the observation scale, we succeeded in computing the total length of critical lines per unit string length and found it to be $4/\sqrt{3} {\bf E}(3/4)$. The length of the associated caustic lines can also be computed to be $16/(\pi \sqrt{3}) {\bf E}(3/4)$. The picture we obtain here for the phenomenology of cosmic string detection is clearly at variance with common lore.

We report the detection of a significant (5.5 sigma) excess of correlations between galaxy ellipticities at scales ranging from 0.5 to 3.5 arc-minutes. This detection of a gravitational lensing signal by large-scale structure was made using a composite high quality imaging survey of 6300 arcmin^2 obtained at the Canada France Hawaii Telescope (CFHT) with the UH8K and CFH12K panoramic CCD cameras. The amplitude of the excess correlation is 2.2\pm 0.2 % at 1 arcmin scale, in agreement with theoretical predictions of the lensing effect induced by large-scale structure.We provide a quantitative analysis of systematics which could contribute to the signal and show that the net effect is small and can be corrected for. We show that the measured ellipticity correlations behave as expected for a gravitational shear signal. The relatively small size of our survey precludes tight constraints on cosmological models. However the data are in favor of cluster normalized cosmological models, and marginally reject Cold Dark Matter models with (Omega=0.3, sigma_8<0.6) or (Omega=1, sigma_8=1). The detection of cosmic shear demonstrates the technical feasibility of using weak lensing surveys to measure dark matter clustering and the potential for cosmological parameter measurements, in particular with upcoming wide field CCD cameras.

We study heterogeneous quantities that efficiently cross-correlate the lensing information encoded in CMB polarization and large-scale structures recovered from weak lensing galaxy surveys. These quantities allow us to take advantage of the special features of weak lensing effect on CMB B-polarization and of the high (40%) cross-correlation between the two data sets. We show that these objects are robust to filtering effects, have a low intrinsic cosmic variance (around 8% for small 100 square degrees surveys) and can be used as an original constraint on the vacuum energy density.

We explore the possibility of using the B-type polarization of the CMB to map the large-scale structures of the Universe taking advantage of the lens effects on the CMB polarization. The functional relation between the B component with the primordial CMB polarization and the line-of-sight mass distribution is explicited. Noting that a sizeable fraction (at least 40%) of the dark halo population which is responsible of this effect can also be detected in galaxy weak lensing survey, we present statistical quantities that should exhibit a strong sensitivity to this overlapping. We stress that it would be a sound test of the gravitational instability picture, independent on many systematic effects that may hamper lensing detection in CMB or galaxy survey alone. Moreover we estimate the intrinsic cosmic variance of the amplitude of this effect to be less than 8% for a 100, deg^2 survey with a 10' CMB beam. Its measurement would then provide us with an original mean for constraining the cosmological parameters, more particularly, as it turns out, the cosmological constant Lambda.

Extending the Kaiser-Stebbins mechanism we propose here a method for detecting relics of topological defects such as cosmic strings based on lens-induced small-scale B-type polarization in the CMB. Models of inflation, in which large scale structures of the Universe emerge from the inflaton fluctuations,do not exclude the formation of topological defects at the end of the inflationary phase. In such case, we show that the lens effect of a string on the small scale E-type polarization of the CMB induces a significant amount of B-type polarization along the line-of-sight. The amplitude of the effect is estimated for different resolutions of CMB experiments.

The errors on statistics measured in finite galaxy catalogs are exhaustively investigated. The theory of errors on factorial moments by Szapudi & Colombi (1996) is applied to cumulants via a series expansion method. All results are subsequently extended to the weakly non-linear regime. Together with previous investigations this yields an analytic theory of the errors for moments and connected moments of counts in cells from highly nonlinear to weakly nonlinear scales. The final analytic formulae representing the full theory are explicit but somewhat complicated. Therefore as a companion to this paper we supply a FORTRAN program capable of calculating the described quantities numerically (abridged).

The low redshift structures of the Universe act as lenses in a similar way on the Cosmic Microwave Background light and on the distant galaxies (say at redshift about unity). As a consequence, the CMB temperature distortions are expected to be statistically correlated with the galaxy shear, exhibiting a non-uniform distribution of the relative angle between the CMB and the galactic ellipticities. Investigating this effect we find that its amplitude is as high as a 10% excess of alignement between CMB and the galactic ellipticities relative to the uniform distribution. The relatively high signal-to-noise ratio we found should makes possible a detection with the planned CMB data sets, provided that a galaxy survey follow up can be done on a sufficiently large area. It would provide a complementary bias-independent constraint on the cosmological parameters.

We present results of the investigations of the statistical properties of a joint density and velocity divergence probability distribution function (PDF) in the mildly non-linear regime. For that purpose we use both perturbation theory results, extended here for a top-hat filter, and numerical simulations. In particular we derive the quantitative (complete as possible up to third order terms) and qualitative predictions for constrained averages and constrained dispersions -- which describe the nonlinearities and the stochasticity properties beyond the linear regime -- and compare them against numerical simulations. We find overall a good agreement for constrained averages; however, the agreement for constrained dispersions is only qualitative. Scaling relations for the Omega-dependence of these quantities are satisfactory reproduced. Guided by our analytical and numerical results, we finally construct a robust phenomenological description of the joint PDF in a closed analytic form. The good agreement of our formula with results of N-body simulations for a number of cosmological parameters provides a sound validation of the presented approach. Our results provide a basis for a potentially powerful tool with which it is possible to analyze galaxy survey data in order to test the gravitational instability paradigm beyond the linear regime and put useful constraints on cosmological parameters. In particular we show how the nonlinearity in the density--velocity relation can be used to break the so-called Omega-bias degeneracy in cosmic density-velocity comparisons.

The question we address in this paper is the determination of the correlation properties of the dark matter halos appearing in cosmic density fields once they underwent a strongly nonlinear evolution induced by gravitational dynamics. Assuming that the high-order correlation functions of the matter field behave as products of two-body correlation functions, we derive the correlation properties of the halos, that are assumed to represent the correlation properties of galaxies or clusters. The hierarchical pattern originally induced by gravity is shown to be conserved for the halos. The strength of their correlations at any order varies, however, but is found to depend only on their internal properties, namely on the parameter x=m/r^(3-gamma) where m is the mass of the halo, r its size and gammma is the power law index of the two-body correlation function. We were able to derive the explicit form of the generating function of the moments of the halo count probability distribution function. In particular we show explicitely that, generically, S_p(x) -> p^(p-2) in the rare halo limit. Various illustrations of our general results are presented. As a function of the properties of the underlying matter field, we construct the count probabilities for halos and in particular discuss the halo void probability. We also evaluate the dependence of the halo mass function on the environment. We found that within clusters, hierarchical clustering implies the higher masses are favored. We stress that this bias is naturally induced by gravity.

A fundamental hypothesis for the interpretation of the measured large-scale line-of-sight peculiar velocities of galaxies is that the large-scale cosmic flows are irrotational. In order to assess the validity of this assumption, we estimate, within the frame of the gravitational instability scenario, the amount of vorticity generated after the first shell crossings in large-scale caustics. In the Zel'dovich approximation the first emerging singularities form sheet like structures. Here we compute the expectation profile of an initial overdensity under the constraint that it goes through its first shell crossing at the present time. We find that this profile corresponds to rather oblate structures in Lagrangian space. Assuming the Zel'dovich approximation is still adequate not only at the first stages of the evolution but also slightly after the first shell crossing, we calculate the size and shape of those caustics and their vorticity content as a function of time and for different cosmologies. The average vorticity created in these caustics is small: of the order of one (in units of the Hubble constant). To illustrate this point we compute the contribution of such caustics to the probability distribution function of the filtered vorticity at large scales. We find that this contribution that this yields a negligible contribution at the 10 to 15 $h^{-1}$Mpc scales. It becomes significant only at the scales of 3 to 4 $h^{-1}$Mpc, that is, slightly above the galaxy cluster scales.

The purpose of these lecture notes is to describe the gravitational lens effects in different astrophysical contexts. These notes are voluntarily focused on the fundamental mechanisms and the basic concepts that are useful to describe these effects. The observational consequences are presented in more details in accompanying notes by Y. Mellier. The content of these notes is the following. In the first section describe of the basic mechanisms of gravitational lenses, techniques and approximations that are usually employed. The second section is devoted to the case of a very simple deflector, a point-like mass distribution. This corresponds to microlensing events in which the deflectors are compact objects of a fraction of a solar mass that may populate the halo of our Galaxy. The last two sections are devoted to cosmological applications. After a presentation of the geometrical quantities that are specific to cosmology, I will present the various phenomena that can be observed in this context. Finally I describe the weak lensing regime. This is a rapidly developing area that should eventually allow us to map the mass distribution in the Universe.

We apply a mass reconstruction technique to simulated large-scale structure gravitational distortion maps, from 2.5' to 10 degree scales, for different cosmological scenarii. The projected mass is reconstructed using a non-parametric least square method after noise due to the galaxy intrinsic ellipticities has been added on. The shearing of the galaxies is performed using the full lensing equation, without any hypothesis like the weak lensing approximation, or other linearization. It is shown that, in the reconstructed maps the noise acts as a perfect uncorrelated Poisson noise, with no propagation at large scales. The measured power spectrum and first four moments of the convergence can be corrected accurately for this source of noise. Using 60 realizations for each model, we show that a weak lensing survey of 5x5 degrees with a typical sheared background population of 30 gal/arcmin^2 at a redshift z_s=1, is able to probe the amplitude of the power spectrum with a few percents accuracy. Such a survey would lead to a 6 sigma separation between open (Omega=0.3) and flat (Omega=1) model. This separation is shown to be robust against different hypothesis for the normalization or the shape of the power spectrum, and does not require very deep surveys. The observational strategy for an optimal measurement of the power spectrum and the moments of the convergence is discussed.

I investigate the effects of source clustering on the weak lensing statistics, more particularly on the statistical properties of the local convergence, kappa, at large angular scales. The Perturbation Theory approach shows that the variance is not affected by source clustering at leading order but higher order moments such as the third and fourth moments can be. I compute the magnitude of these effects in case of an Einstein-de Sitter Universe for the angular top-hat filtered convergence. In these calculations the so-called Broadhurst and multiple lens coupling effects are neglected. The source clustering effect is found to be particularly important when the redshift distribution is broad enough so that remote background sources can be significantly lensed by closer concentrations of galaxy sources. The source clustering effects are shown to remain negligible, for both the skewness and the kurtosis, when the dispersion of the redshift of the sources is less than about 0.15.

Sujets qui sont plus particulierement abordes dans ce cours: 1. Elements de cosmologie, l'Univers homogene Les principes cosmologiques et leurs verifications observationnelles Histoire thermique de l'Univers Inflation, motivations et principe de base 2. Elements de cosmologie, l'Univers inhomogene La croissance des fluctuations Pourquoi un modele avec de la Matiere Noire Froide ? Evidences de l'existence de matiere noire 3. La dynamique gravitationnelle, les theories lineaires L'approximation Newtonnienne avec un seul flot Description Eulerienne ou description Lagrangienne Vers la dynamique non-lineaire: l'effondrement spherique 4. Le regime quasilineaire Effet du couplage de mode: la skewness La hierarchie des correlations en regimes quasi-lineaire Les termes sous-dominants du developpement perturbatif 5. Application aux proprietes statistiques des champs de distorsions gravitationnelles La relation entre convergence locale et densite projetee Variance et skewness de la convergence locale 6. Vers le regime fortement non-lineaire Les solutions auto-similaires et les modeles hierarchiques Distribution de matiere et distribution de lumiere

We describe two new -- stochastic-geometrical -- methods to obtain reliable velocity field statistics from N-body simulations and from any general density and velocity fluctuation field sampled at a discrete set of locations. These methods, the Voronoi tessellation method and Delaunay tessellation method, are based on the use of the Voronoi and Delaunay tessellations of the point distribution defined by the locations at which the velocity field is sampled. Adjusting themselves automatically to the density of sampling points, they represent the optimal estimator for volume-averaged quantities. They are therefore particularly suited for checking the validity of the predictions of quasi-linear analytical density and velocity field perturbation theory through the results of N-body simulations of structure formation. We illustrate the subsequent succesfull application of the two methods to estimate the bias-independent value of $\Omega$ in the N-body simulations on the basis of the predictions of perturbation theory for the $\Omega$-dependence of the moments and PDF of the velocity divergence in gravitational instability structure formation scenarios with Gaussian initial conditions. We will also shortly discuss practical and complicating issues involved in the obvious extension of the Voronoi and Delaunay method to the analysis of observational samples of galaxy peculiar velocities.

Weak lensing effects are known to introduce non-linear couplings in CMB temperature maps. In inflationary scenario, the primary CMB anisotropies are expected to form a 2D Gaussian map, for which, the probability distribution function of the ellipticity defined from the local temperature curvature matrix has a very specific shape. I show that lenses alter significantly the shape of this PDF, inducing an excess of elongated structures. The precise functional form is computed for both the field points and the temperature extrema. These analytical results are confirmed by numerical experiments on 10x10 square degree maps. These numerical results allow to investigate the effects of smoothing and to estimate the cosmic variance. For the best resolution and sky coverage of the Planck mission the signal to noise ratio for the statistical indicators presented here is about 3 to 6 depending on the cosmological models. A marginal detection should therefore be possible.

The last decade has shown a considerable development of gravitational lensing for cosmology because it probes the amount and the nature of dark matter, and provides information on the density parameter $\Omega$, the cosmological constant $\Lambda$ and the Hubble constant $H_{o}$. Therefore, gravitational lensing can constrain the cosmological scenario which gave birth to the Universe as it appears today. The ongoing programs and future projects which are developing now all over the world show that gravitational lensing is considered as a major cosmological tool for the coming years as well. In this review, we summarize some of the most recent advances in the fields relevant for the dark matter issue. We will focus on the microlensing, the arc(let)s and the weak lensing studies. The possibility to check the existence of a non-zero $\Lambda$ is presented elsewhere.

We study the domain of validity of Perturbation Theory (PT), by comparing its predictions for the reduced skewness, s_3, and kurtosis, s_4, of the projected cosmological density field, with the results of N-body simulations. We investigate models with different linear power spectra and consider as physical applications both angular galaxy catalogues and weak lensing surveys. We first find that the small-angle approximation for the predicted skewness provides a good match to the exact numerical PT results. On the other hand, results from non-linear simulated catalogues agree well with PT results on quasi-linear angular scales, which correspond to scales larger than about 1 degree in the applications we have considered. We also point out that, on smaller scale, the projection effects tend to attenuate the effects of the strong nonlinearities in the angular skewness and kurtosis.

The possibility of detecting weak lensing effects from deep wide field imaging surveys has open new means of probing the large-scale structures of the Universe and measuring cosmological parameters. In this paper we present a systematic study of the expected dependence of the low order moments of the filtered gravitational local convergence with the power spectrum of the density fluctuations and with the cosmological parameters Omega_0 and Lambda. The results show a significant dependence with all these parameters. Though we note that this degeneracy could be partially raised by considering two populations of sources, at different redshifts, computing the third moment is more promising since it is expected, in the quasi-linear regime and for Gaussian initial conditions, to be only $\Omega_0$ dependent (with a slight degeneracy with Lambda) when it is correctly expressed in terms of the second moment. More precisely we show that the variance of the convergence varies approximately as P(k) Omega_0^{1.5} z_s^{1.5}, whereas the skewness varies as Omega_0^{-0.8} z_s^{-1.35}, where P(k) is the projected power spectrum and z_s the redshift of sources. Thus, used jointly they can provide both P(k) and Omega_0. However, the dependence with the redshift of sources is large and could be a major concern for a practical implementation. We have estimated the errors expected for these parameters in realistic scenario and sketched what would be the observational requirements for doing such measurements. A more detailed study of an observational strategy is left for a second paper.

Perturbation Theory (PT) applied to a cosmological density field with Gaussian initial fluctuations suggests a specific hierarchy for the correlation functions when the variance is small. In particular quantitative predictions have been made for the moments and the shape of the one-point probability distribution function (PDF) of the top-hat smoothed density. In this paper we perform a series of systematic checks of these predictions against N-body computations both in 2D and 3D with a wide range of featureless power spectra. In agreement with previous studies, we found that the reconstructed PDF-s work remarkably well down to very low probabilities, even when the variance approaches unity. Our results for 2D reproduce the features for the 3D dynamics. In particular we found that the PT predictions are more accurate for spectra with less power on small scales. The nonlinear regime has been explored with various tools, PDF-s, moments and Void Probability Function (VPF). These studies have been done with unprecedented dynamical range, especially for the 2D case, allowing in particular more robust determinations of the asymptotic behaviour of the VPF. We have also introduced a new method to determine the moments based on the factorial moments. Results using this method and taking into account the finite volume effects are presented.

Analytical studies based on perturbative theory have shown that the moments of the Probability Distribution Function (PDF) of the local smoothed velocity divergence are expected to have a very specific dependence on the density parameter Omega in the quasi-linear regime. This dependence is particularly interesting as it does not involve the possible bias between the galaxy spatial distribution and the underlying mass distribution. This implies a new and promising method for determining a bias-independent value of Omega based on a reliable determination of the velocity divergence PDF. In this paper we study the Omega dependence of the velocity divergence PDF and its first moments in a set of N-body simulations, using the so-called Voronoi and Delaunay methods. We show that this dependence is in agreement with the theoretical prediction, even while the number density of velocity field tracers has been diluted to a value comparable to that available in current galaxy catalogues. In addition, we demonstrate that a sufficiently reliable determination of these statistical quantities is also possible when the measurement of the galaxy peculiar velocities is restricted to the one component along the line-of-sight. Under ideal, noise-free circumstances we can successfully discriminate between low and high Omega.

The weak lensing effects are known to change only weakly the shape of the power spectrum of the Cosmic Microwave Background (CMB) temperature fluctuations. I show here that they nonetheless induce specific non-Gaussian effects that can be detectable with the four-point correlation function of the CMB anisotropies. The magnitude and geometrical dependences of this correlation function are investigated in detail. It is thus found to scale as the square of the derivative of the two-point correlation function and as the angular correlation function of the gravitational displacement field. It also contains specific dependences on the shape of the quadrangle formed by the four directions. When averaged at a given scale, the four-point function, that identifies with the connected part of the fourth moment of the probability distribution function of the local filtered temperature, scales as the square of logarithmic slope of its second moment, and as the variance of the gravitational magnification at the same angular scale. All these effects have been computed for specific cosmological models. It is worth noting that, as the amplitude of the gravitational lens effects has a specific dependence on the cosmological parameters, the detection of the four-point correlation function could provide precious complementary constraints to those brought by the temperature power spectrum.

Perturbation theory makes it possible to calculate the probability distribution function (PDF) of the large scale density field in the small variance limit. For top hat smoothing and scale-free Gaussian initial fluctuations, the result depends only on the linear variance, sigma_linear, and its logarithmic derivative with respect to the filtering scale -(n_linear+3)=dlog sigma_linear^2/dlog L (Bernardeau 1994). In this paper, we measure the PDF and its low-order moments in scale-free simulations evolved well into the nonlinear regime and compare the results with the above predictions, assuming that the spectral index and the variance are adjustable parameters, n_eff and sigma_eff=sigma, where sigma is the true, nonlinear variance. With these additional degrees of freedom, results from perturbation theory provide a good fit of the PDFs, even in the highly nonlinear regime. The value of n_eff is of course equal to n_linear when sigma << 1, and it decreases with increasing sigma. A nearly flat plateau is reached when sigma >> 1. In this regime, the difference between n_eff and n_linear increases when n_linear decreases. For initial power-spectra with n_linear=-2,-1,0,+1, we find n_eff ~ -9,-3,-1,-0.5 when sigma^2 ~ 100.

The gravitational lensing effect is one of the most promising tools for cosmology. Indeed it probes directly the total mass distribution in large-scale structures and can as well provide valuable informations on the values of the density parameter Omega and the cosmological constant lambda. In this review we summarize the most spectacular observational and theoretical advances obtained in this field during the last five years.

The understanding of the large-scale structure formation requires the resolution of coupled nonlinear equations describing the cosmic density and velocity fields. This is a complicated problem that, for the last decade, has been essentially addressed with N-body simulations. There is however a regime, the so-called quasi-linear regime, for which the relative density fluctuations are on average below unity. It is then possible to apply Perturbation Theory techniques where the perturbation expansions are made with respect to the initial fluctuations. I review here the major results that have been obtained in this regime.

It is known that in gravitational instability scenarios the nonlinear dynamics induces non-Gaussian features in cosmological density fields that can be investigated with perturbation theory. Here, I derive the expression of the joint moments of cosmological density fields taken at two different locations. The results are valid when the density fields are filtered with a top-hat filter window function, and when the distance between the two cells is large compared to the smoothing length. In particular I show that it is possible to get the generating function of the coefficients C_{p,q} defined by < delta^p(x_1) delta^q(x_2) >_c= C_{p,q} < delta^2(x) >^{p+q-2} < delta(x_1) delta(x_2) > where delta(x) is the local smoothed density field. It is then possible to reconstruct the joint density probability distribution function (PDF), generalizing for two points what has been obtained previously for the one-point density PDF. I discuss the validity of the large separation approximation in an explicit numerical Monte Carlo integration of the C_{2,1} parameter as a function of |x_1-x_2|. A straightforward application is the calculation of the large-scale `bias' properties of the over-dense (or under-dense) regions. The properties and the shape of the bias function are presented in detail and successfully compared with numerical results obtained in an N-body simulation with CDM initial conditions.

A lot of predictions for the statistical properties of the cosmic velocity field at large-scale have been obtained recently using perturbation theory. In this contribution I report the outcomes of a set of numerical tests that aim to check these results. Using Voronoi and Delaunay tessellations for defining the velocity field by interpolation between the particle velocities in numerical simulations, we have been able to get reliable estimates of the local velocity gradients. Thus, we have been able to show that the properties of the velocity divergence are in very good agreement with the analytical results. In particular we have confirmed the $\Omega$ dependence expected for the shape of its distribution function.

We introduce two new methods to obtain reliable velocity field statistics from N-body simulations, or indeed from any general density and velocity fluctuation field sampled by discrete points. These methods, the {\it Voronoi tessellation method} and {\it Delaunay tessellation method}, are based on the use of the Voronoi and Delaunay tessellations of the point distribution defined by the locations at which the velocity field is sampled. In the Voronoi method the velocity is supposed to be uniform within the Voronoi polyhedra, whereas the Delaunay method constructs a velocity field by linear interpolation between the four velocities at the locations defining each Delaunay tetrahedron. The most important advantage of these methods is that they provide an optimal estimator for determining the statistics of volume-averaged quantities, as opposed to the available numerical methods that mainly concern mass-averaged quantities. As the major share of the related analytical work on velocity field statistics has focussed on volume-averaged quantities, the availability of appropriate numerical estimators is of crucial importance for checking the validity of the analytical perturbation calculations. In addition, it allows us to study the statistics of the velocity field in the highly nonlinear clustering regime. Although in most cases both the Voronoi and the Delaunay method lead to equally good velocity field estimates, the Delaunay method may be slightly preferable. In particular it performs considerably better at small radii. Note that it is more CPU-time intensive while its requirement for memory space is almost a factor 8 lower than the Voronoi method.

For Gaussian initial conditions the perturbation theory predicts a very specific hierarchy for the projected matter $p$-point correlation functions. In the small angle approximation and assuming a power-law spectrum I derive the exact expressions of the coefficients $s_p$ relating the averaged $p$-order angular correlation function, $\omb_p$ to the second one, $\omb_p=s_p\ \omb_2^{p-1}$. These results are valid for any selection function, but for a top-hat angular filter only. These coefficients are found to be significantly higher than their 3D counterparts, $S_p=\xib_p/\xib_2^{p-1}$. For the coefficient $s_3$ I discussed the accuracy of the small angle approximation by computing, for particular examples, its angular dependence with Monte-Carlo numerical integrations. It is found that the accuracy of the small angle approximation for $\theta\approx 1^o$ slightly depends on the selection function. Using the selection function expected for galaxy catalogues the approximation is found to be reasonably good. The measurements of the $s_p$ parameters made in the APM angular survey are found to give systematic lower values than the theoretical predictions. How significant this discrepancy is and what the implications would be for galaxy formation models is discussed in the last section.

We propose a method for measuring the cosmological density parameter $\Omega$ from the statistics of the divergence field, $\theta \equiv H^{-1} \div v$, the divergence of peculiar velocity, expressed in units of the Hubble constant, $H \equiv 100 h km/s/Mpc$. The velocity field is spatially smoothed over $\sim 10 h^{-1} Mpc$ to remove strongly nonlinear effects. Assuming weakly-nonlinear gravitational evolution from Gaussian initial fluctuations, and using second-order perturbative analysis, we show that $ <\theta^3> \propto -\Omega^{-0.6} <\theta^2>^2$. The constant of proportionality depends on the smoothing window. For a top-hat of radius R and volume-weighted smoothing, this constant is $26/7-\gamma$, where $\gamma=-d\log <\theta^2> / d\log R$. If the power spectrum is a power law, $P(k)\propto k^n$, then $\gamma=3+n$. A Gaussian window yields similar results. The resulting method for measuring $\Omega$ is independent of any assumed biasing relation between galaxies and mass. The method has been successfully tested with numerical simulations. A preliminary application to real data, provided by the POTENT recovery procedure from observed velocities favors $\Omega \sim 1$. However, because of an uncertain sampling error, this result should be treated as an assessment of the feasibility of our method rather than a definitive measurement of $\Omega$.

The properties of the probability distribution function of the cosmological continuous density field are studied. We present further developments and compare dynamically motivated methods to derive the PDF. One of them is based on the Zel'dovich approximation (ZA). We extend this method for arbitrary initial conditions, regardless of whether they are Gaussian or not. The other approach is based on perturbation theory with Gaussian initial fluctuations. We include the smoothing effects in the PDFs. We examine the relationships between the shapes of the PDFs and the moments. It is found that formally there are no moments in the ZA, but a way to resolve this issue is proposed, based on the regularization of integrals. A closed form for the generating function of the moments in the ZA is also presented, including the smoothing effects. We suggest the methods to build PDFs out of the whole series of the moments, or out of a limited number of moments -- the Edgeworth expansion. The last approach gives us an alternative method to evaluate the skewness and kurtosis by measuring the PDF around its peak. We note a general connection between the generating function of moments for small r.m.s $\sigma$ and the non-linear evolution of the overdense spherical fluctuation in the dynamical models. All these approaches have been applied in 1D case where the ZA is exact, and simple analytical results are obtained. The 3D case is analyzed in the same manner and we found a mutual agreement in the PDFs derived by different methods in the the quasi-linear regime. Numerical CDM simulation was used to validate the accuracy of considered approximations. We explain the successful log-normal fit of the PDF from that simulation at moderate $\sigma$ as mere fortune, but not as a universal form of density PDF in general.

It has been shown that the large--scale correlation functions of the density field (and velocity divergence field) follow a specific hierarchy in the quasilinear regime and for Gaussian initial conditions (Bernardeau 1992). The exact relationships between the cumulants of the probability distribution functions (the so-called $S_p$ parameters) are however sensitive to the smoothing window function applied to the fields. In this paper, I present a method to derive the whole series of the $S_p$ parameters when the density field is smoothed with a top--hat window function. The results are valid for any power spectrum and any cosmological parameters. Similar calculations are presented for the velocity divergence field. The resulting shapes of the one--point probability distribution functions of the cosmic density and the velocity divergence fields are given as a function of the power spectrum and $\Omega$. Simple analytical fits are proposed when the index of the power spectrum is -1. Comparisons with numerical simulations prove these analytical results to be extremely accurate.

In this paper, I present the calculation of the third and fourth moments of both the distribution function of the large--scale density and the large--scale divergence of the velocity field, $\theta$. These calculations are made by the mean of perturbative calculations assuming Gaussian initial conditions and are expected to be valid in the linear or quasi linear regime. The moments are derived for a top--hat window function and for any cosmological parameters $\Omega$ and $\Lambda$. It turns out that the dependence with $\Lambda$ is always very weak whereas the moments of the distribution function of the divergence are strongly dependent on $\Omega$. A method to measure $\Omega$ using the skewness of this field has already been presented by Bernardeau et al. (1993). I show here that the simultaneous measurement of the skewness and the kurtosis allows to test the validity of the gravitational instability scenario hypothesis. Indeed there is a combination of the first three moments of $\theta$ that is almost independent of the cosmological parameters $\Omega$ and $\Lambda$, $${(<\theta^4>-3<\theta^2>^2) <\theta^2> \over <\theta^3>^2}\approx 1.5,$$ (the value quoted is valid when the index of the power spectrum at the filtering scale is close to -1) so that any cosmic velocity field created by gravitational instabilities should verify such a property.

In this paper I consider the nonlinear evolution of a rare density fluctuation in a random density field with Gaussian fluctuations, and I rigorously show that it follows the spherical collapse dynamics applied to its mean initial profile. This result is valid for any cosmological model and is independent of the shape of the power spectrum. In the early stages of the dynamics the density contrast of the fluctuation is seen to follow with a good accuracy the form $$\delta=(1-\delta_L/1.5)^{-1.5}-1,$$ where $\delta_L$ is the linearly extrapolated overdensity. I then investigate the validity domain of the rare event approximation in terms of the parameter $\nu=\delta_L/\sigma$ giving the initial overdensity scaled by the rms fluctuation at the same mass scale, and find that it depends critically on the shape of the power spectrum. When the power law index $n$ is lower than $-1$ the departure from the spherical collapse is expected to be small, at least in the early stages of the dynamics, and even for moderate values of $\nu$ ($\vert\nu\vert\ge 2$). However, for $n\ge -1$ and whatever the value of $\nu$, the dynamics seems to be dominated by the small--scale fluctuations and the subsequent evolution of the peak may not be necessarily correlated to its initial overdensity. I discuss the implications of these results for the nonlinear dynamics and the formation of astrophysical objects.

Two quasi-linear approximations, the frozen flow approximation (FFA) and the frozen potential approximation (FPA), have been proposed recently for studying the evolution of a collisionless self-gravitating fluid. In the FFA it is assumed that the velocity field remains unchanged from its value obtained from the linear theory whereas in FPA the same approximation is made for the gravitational potential. In this paper we compare these and the older Zel'dovich approximation by calculating the evolution of the density in perturbation theory. In particular we compute the skewness, including the smoothing effects, and the kurtosis for the FFA, FPA and Zel'dovich approximation and compare their relative accuracy.

Velocity dispersion $\sigma$, radius $R$ and luminosity $L$ of elliptical galaxies are known to be related, leaving only two degrees of freedom and defining the so-called ``fundamental plane". In this {\em Letter} we present observational evidence that rich galaxy clusters exhibit a similar behaviour. Assuming a relation $L \propto R^{\alpha}\sigma^{2 \beta}$, the best-fit values of $\alpha$ and $\beta$ are very close to those defined by galaxies. The dispersion of this relation is lower than 10 percent, i.e. significantly smaller than the dispersion observed in the $L-\sigma$ and $L-R$ relations. We briefly suggest some possible implications on the spread of formation times of objects and on peculiar velocities of galaxy clusters.