Dynamics of gravitational clustering IV. The probability distribution of rare events

Using a non-perturbative method developed in a previous article (paper II) we investigate the tails of the probability distribution $P(\rho_R)$ of the overdensity within spherical cells. We show that our results for the low-density tail of the pdf agree with perturbative results when the latter are finite (up to the first subleading term), that is for power-spectra with $-3<n<-1$. Over the range $-1<n<1$ some shell-crossing occurs (which leads to the break-up of perturbative approaches) but this does not invalidate our approach. In particular, we explain that we can still obtain an approximation for the low-density tail of the pdf. This feature also clearly shows that perturbative results should be viewed with caution (even when they are finite). We point out that our results can be recovered by a simple spherical model but they cannot be derived from the stable-clustering ansatz in the regime $\sigma \gg 1$ since they involve underdense regions which are still expanding. Second, turning to high-density regions we explain that a naive study of the radial spherical dynamics fails. Indeed, a violent radial-orbit instability leads to a fast relaxation of collapsed halos (over one dynamical time) towards a roughly isotropic equilibrium velocity distribution. Then, the transverse velocity dispersion stabilizes the density profile so that almost spherical halos obey the stable-clustering ansatz for $-3<n<1$. We again find that our results for the high-density tail of the pdf agree with a simple spherical model (which takes into account virialization). Moreover, they are consistent with the stable-clustering ansatz in the non-linear regime. Besides, our approach justifies the large-mass cutoff of the Press-Schechter mass function (although the various normalization parameters should be modified).Comment: 27 pages, final version published in A&

Construction of the one-point PDF of the local aperture mass in weak lensing maps

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.Comment: 17 pages, final version published in A&

Expansion schemes for gravitational clustering: computing two-point and three-point functions

We describe various expansion schemes that can be used to study gravitational clustering. Obtained from the equations of motion or their path-integral formulation, they provide several perturbative expansions that are organized in different fashion or involve different partial resummations. We focus on the two-point and three-point correlation functions, but these methods also apply to all higher-order correlation and response functions. We present the general formalism, which holds for the gravitational dynamics as well as for similar models, such as the Zeldovich dynamics, that obey similar hydrodynamical equations of motion with a quadratic nonlinearity. We give our explicit analytical results up to one-loop order for the simpler Zeldovich dynamics. For the gravitational dynamics, we compare our one-loop numerical results with numerical simulations. We check that the standard perturbation theory is recovered from the path integral by expanding over Feynman's diagrams. However, the latter expansion is organized in a different fashion and it contains some UV divergences that cancel out as we sum all diagrams of a given order. Resummation schemes modify the scaling of tree and one-loop diagrams, which exhibit the same scaling over the linear power spectrum (contrary to the standard expansion). However, they do not significantly improve over standard perturbation theory for the bispectrum, unless one uses accurate two-point functions (e.g. a fit to the nonlinear power spectrum from simulations). Extending the range of validity to smaller scales, to reach the range described by phenomenological models, seems to require at least two-loop diagrams.Comment: 24 pages, published in A&