We propose a new mechanism for the formation of dark matter clumps in the radiation era. We assume that a light scalar field oscillates harmonically. We include self-interactions and consider the nonrelativistic regime, where the scalar dynamics are akin to the one of a fluid whose pressure depends on both quantum and self-interaction effects. When the squared speed of sound of the scalar fluid becomes negative, an instability arises and the fluctuations of the scalar energy-density field start growing and eventually reach the nonlinear regime, where clumps form. Subsequently, the clumps aggregate and reach a universal regime. We apply this mechanism to a model with a negative quartic term stabilised by a positive self-interaction of order six, and then to axion monodromy, where a subdominant cosine potential corrects a mass term. In the first case, the speed of sound squared becomes negative when the quartic term dominates, leading to a tachyonic instability. For axion monodromy, the instability starts very slowly. First, the density perturbations perform acoustic oscillations due to the quantum pressure. Eventually, they start growing exponentially due to a parametric resonance. The mass and size of the clumps depend crucially on their formation mechanisms. When a tachyonic phase takes place, the core density of the clumps is uniquely determined by the energy density at the beginning of the instability. For axion monodromy, the core density scales with the soliton mass and radius. This difference comes from the role that the quantum pressure plays in the parametric resonance and in balancing the self-interactions and gravity to form static scalar solitons. In both scenarios, we find that the formation redshift of the scalar clumps can span a large range in the radiation era from $10^{-2}$ to $10^4 \, {\rm GeV}$ and the scalar-field mass can be from $10^{-18}$ GeV to $100$ GeV.

We study the predictions for the matter redshift-space power spectrum and correlation function of a Lagrangian-space Gaussian ansatz introduced in a previous work. This model is a natural extension of the Zeldovich approximation, where the displacement and velocity power spectra are determined by the equations of motion, instead of being set equal to the linear power spectrum. It does not contain any free parameter. As for the real-space statistics, we find that this Lagrangian-space approach is much more efficient for the correlation functions than for the power spectra. The damping of the BAO oscillations is well recovered but there is a large smooth drift from the simulations in the power spectra. The multipoles of the correlation functions are well recovered on BAO scales, with an accuracy of $2\%$ for $\xi^s_0$ down to $10 h^{-1}$ Mpc, and of $3\%$ for $\xi^s_2$ down to $26 h^{-1}$ Mpc, at $z \geq 0.35$.

We present a new approach to model the gravitational dynamics of large-scale structures. Instead of solving the equations of motion up to a finite perturbative order or building phenomenological models, we follow the evolution of the probability distribution of the displacement and velocity fields within an approximation subspace. Keeping the exact equations of motion with their full nonlinearity, this provides a nonperturbative scheme that goes beyond shell crossing. Focusing on the simplest case of a curl-free Gaussian ansatz for the displacement and velocity fields, we find that truncations of the power spectra on nonlinear scales directly arise from the equations of motion. This leads to a truncated Zeldovich approximation for the density power spectrum, but with a truncation that is not set a priori and with different power spectra for the displacement and velocity fields. The positivity of their auto power spectra also follows from the equations of motion. Although the density power spectrum is only recovered up to a smooth drift on BAO scales, the predicted density correlation function agrees with numerical simulations within $2\%$ from BAO scales down to $7 h^{-1} {\rm Mpc}$ at $z \geq 0.35$, without any free parameter.

We investigate scalar-tensor theories where matter couples to the scalar field via a kinetically dependent conformal coupling. These models can be seen as the low-energy description of invariant field theories under a global Abelian symmetry. The scalar field is then identified with the Goldstone mode of the broken symmetry. It turns out that the properties of these models are very similar to the ones of ultralocal theories where the scalar-field value is directly determined by the local matter density. This leads to a complete screening of the fifth force in the Solar System and between compact objects, through the ultralocal screening mechanism. On the other hand, the fifth force can have large effects in extended structures with large-scale density gradients, such as galactic halos. Interestingly, it can either amplify or damp Newtonian gravity, depending on the model parameters. We also study the background cosmology and the linear cosmological perturbations. The background cosmology is hardly different from its $\Lambda$-CDM counterpart whilst cosmological perturbations crucially depend on whether the coupling function is convex or concave. For concave functions, growth is hindered by the repulsiveness of the fifth force whilst it is enhanced in the convex case. In both cases, the departures from the $\Lambda$-CDM cosmology increase on smaller scales and peak for galactic structures. For concave functions, the formation of structure is largely altered below some characteristic mass, as smaller structures are delayed and would form later through fragmentation, as in some warm dark matter scenarios. For convex models, small structures form more easily than in the $\Lambda$-CDM scenario.

We consider scalar dark matter models where the theory has a shift symmetry only broken by the scalar mass term. We restrict ourselves to K-essence kinetic terms where the shift symmetric part of the Lagrangian is a function of the first derivatives of the scalar field only. In the low-amplitude and nonrelativistic regime, which applies on large galactic scales, scalar clouds form solitons with a finite core. Close to the center of galaxies, where a supermassive Black Hole (BH) resides, we analyze the scalar field distribution and the fate of the dark matter soliton when subject to the BH gravitational attraction. We show that the scalar field profile around such a central BH can be described by new oscillatory solutions of a modified Klein-Gordon equation, which generalize the harmonic oscillations of free scalar dark matter in a flat environment and the Jacobi elliptic functions of the $\phi^4$ model. Moreover, we find that, depending on the form of the K-essence kinetic term, regular solutions can be constructed or not, which connect the relativistic ingoing wavelike profile of the scalar field at the BH horizon to the nearly static nonrelativistic soliton at large distance. These profiles have a constant flux and represent the slow infall of scalar matter into the BH. We show that this regular behavior is only possible for K-essence functions that satisfy the usual conditions for the absence of ghosts and gradient instabilities, together with a new restriction on the growth of the kinetic function $K(X)$ for large argument. It turns out that the same conditions of stability guarantee that quantum corrections are tamed, provided that the mass of the scalar field is less than $10^{-3}$ eV and the strong coupling scale of the model $\Lambda$ is much larger than the scalar mass.

Context. Scaling relations between cluster properties embody the formation and evolution of cosmic structure. Intrinsic scatters and correlations between X-ray properties are determined from merger history, baryonic processes, and dynamical state. Aims. We look for an unbiased measurement of the scatter covariance matrix between the three main X-ray observable quantities attainable in large X-ray surveys -- temperature, luminosity, and gas mass. This also gives us the cluster property with the lowest conditional intrinsic scatter at fixed mass. Methods. Intrinsic scatters and correlations can be measured under the assumption that the observable properties of the intra-cluster medium hosted in clusters are log-normally distributed around power-law scaling relations. The proposed method is self-consistent, based on minimal assumptions, and requires neither the external calibration by weak lensing, dynamical, or hydrostatic masses nor the knowledge of the mass completeness. Results. We analyzed the 100 brightest clusters detected in the XXL Survey and their X-ray properties measured within a fixed radius of 300 kpc. The gas mass is the less scattered proxy (~8%). The temperature (~20%) is intrinsically less scattered than the luminosity (~30%) but it is measured with a larger observational uncertainty. We found some evidence that gas mass, temperature and luminosity are positively correlated. Time-evolutions are in agreement with the self-similar scenario, but the luminosity-temperature and the gas mass-temperature relations are steeper. Conclusions. Positive correlations between X-ray properties can be determined by the dynamical state and the merger history of the halos. The slopes of the scaling relations are affected by radiative processes.

In scalar-field dark matter scenarios, a scalar-field soliton could form at the center of galactic halos, around the supermassive black holes that sit at the center of galaxies. Focusing on the large scalar-mass limit, where the soliton is formed by the balance between self-gravity and a repulsive self-interaction, we study the infall of the scalar field onto the central Schwarzschild black hole. We derive the scalar-field profile, from the Schwarzschild radius to the large radii dominated by the scalar cloud. We show that the steady state solution selects the maximum allowed flux, with a critical profile that is similar to the transonic solution obtained for the hydrodynamic case. This finite flux, which scales as the inverse of the self-interaction coupling, is small enough to allow the dark matter soliton to survive for many Hubble times.

We propose a simple model that provides a dynamical cancellation mechanism of the vacuum energy density appearing either in the form of a bare cosmological constant, quantum fluctuations of matter fields or the result of phase transitions. This `conformal compensator model' is based on a conformal coupling $A(\varphi)$ between the Einstein and the Jordan frames. This couples a second scalar field $\lambda$ to the trace of the matter energy-momentum tensor, including the bare cosmological constant, and serves as a dynamical Lagrange multiplier. As a result, the scalar $\lambda$ relaxes to a value which cancels the contributions from the vacuum energy density to the Friedmann equation, and adjusts itself to changes of the vacuum energy density after matter phase transitions. This circumvents Weinberg's theorem through the time dependence of the background scalar field $\varphi$. The radiation era, where the vacuum energy is annulled, is recovered in a natural manner. It is also possible to recover the matter era, via a tracking of the matter energy density by the scalar field, as well as the inflationary and dark energy eras, which correspond to regimes where the cancellation mechanism becomes inefficient. This suggests that inflation, dark energy, and the annulation of the vacuum energy density, could be related to the same mechanism. In this setting, the usual fine-tuning of the vacuum energy is avoided, although the onset of the dark energy era at the appropriate time is not explained.

We consider models of scalar dark matter with a generic interaction potential and non-canonical kinetic terms of the K-essence type that are subleading with respect to the canonical term. We analyze the low-energy regime and derive, in the nonrelativistic limit, the effective equations of motions. In the fluid approximation they reduce to the conservation of matter and to the Euler equation for the velocity field. We focus on the case where the scalar field mass $10^{-21} \ll m \lesssim 10^{-4} \, {\rm eV}$ is much larger than for fuzzy dark matter, so that the quantum pressure is negligible on cosmological and galactic scales, while the self-interaction potential and non-canonical kinetic terms generate a significant repulsive pressure. At the level of cosmological perturbations, this provides a dark-matter density-dependent speed of sound. At the nonlinear level, the hydrostatic equilibrium obtained by balancing the gravitational and scalar interactions imply that virialized structures have a solitonic core of finite size depending on the speed of sound of the dark matter fluid. For the most relevant potential in $\lambda_4 \phi^4/4$ or K-essence with a $(\partial \phi)^4$ interaction, the size of such stable cores cannot exceed 60 kpc. Structures with a density contrast larger than $10^6$ can be accommodated with a speed of sound $c_s\lesssim 10^{-6}$. We also consider the case of a cosine self-interaction, as an example of bounded nonpolynomial self-interaction. This gives similar results in low-mass and low-density halos whereas solitonic cores are shown to be absent in massive halos.

We consider the gravitational Wilsonian effective action at low energy when all the particles of the standard model have decoupled. When the ${\cal R}^2$ terms dominate, the theory is equivalent to a scalar-tensor theory with the universal coupling $\beta=1/\sqrt 6$ to matter for which we present strong lower and upper bounds on the scalaron mass $m$ obtained by using results from the E\"ot-Wash experiment on the modification of the inverse-square law, the observations of the hot gas of galaxy clusters and the Planck satellite data on the neutrino masses. In terms of the range of the scalar interaction mediated over a distance of order $m^{-1}$, this leads to the small interval $4\,\mu m \lesssim m^{-1} \lesssim 68\, \mu m$ within reach of future experimental tests of deviations from Newton's gravitational inverse-square law.

We investigate the impact of modified-gravity models on the Lyman-$\alpha$ power spectrum. Building a simple analytical modeling, based on a truncated Zeldovich approximation, we estimate the intergalactic medium power spectrum and the Lyman-$\alpha$ flux decrement power spectrum along the line of sight. We recover the results of numerical simulations for $f(R)$-gravity models and present new results for K-mouflage scenarios. We find that the shape of the distortion due to the modified gravity depends on the model, through the scale-dependence or not of their growth rate. This is more clearly seen in the three-dimensional power spectrum than in the one-dimensional power spectrum, where the line-of-sight integration smoothes the deviation. Whilst the Lyman-$\alpha$ power spectrum does not provide competitive bounds for $f(R)$ theories, it could provide useful constraints for the K-mouflage models. Thus, the efficiency of the Lyman-$\alpha$ power spectrum as a probe of modified-gravity scenarios depends on the type of screening mechanism and the related scale dependence it induces. The prospect of a full recovery of the three-dimensional Lyman-$\alpha$ power spectrum from data would also lead to stronger constraints and a better understanding of screening mechanisms.

We describe the effects of the cosmological background on the K-mouflage screening properties of an astrophysical structure. We show that the K-mouflage screening of the spatial gradients of the scalar field, i.e. the screening of the fifth force, happens inside a dynamically generated screening radius. This radius is smaller than the location where the quasistatic approximation, i.e. where the spatial gradients exceed the time derivative, holds. Even though this quasistatic radius is much smaller than the size of the matter overdensity, spatial gradients remain well described by the quasistatic approximation up to the horizon. Cosmologically we find that the time derivatives can remain dominant at redshifts $z\gtrsim 2$, when the cosmic web shows a faster growth. Despite the existence of K-mouflage screening, we confirm that the values of the scalar field itself are still dominated by the cosmological background, down to the center of the matter overdensity, and that for instance the time drift of Newton's constant due to the large-scale cosmological evolution highly constrains K-mouflage models.

We describe scalar-bimetric theories where the dynamics of the Universe are governed by two separate metrics, each with an Einstein-Hilbert term. In this setting, the baryonic and dark matter components of the Universe couple to metrics which are constructed as functions of these two gravitational metrics. The scalar field, contrary to dark energy models, does not have a potential whose role is to mimic a late-time cosmological constant. The late-time acceleration of the expansion of the Universe can be easily obtained at the background level in these models by appropriately choosing the coupling functions appearing in the decomposition of the vierbeins for the baryonic and dark matter metrics. We explicitly show how the concordance model can be retrieved with negligible scalar kinetic energy. This requires the scalar coupling functions to show variations of order unity during the accelerated expansion era. This leads in turn to deviations of order unity for the effective Newton constants and a fifth force that is of the same order as Newtonian gravity, with peculiar features. The baryonic and dark matter self-gravities are amplified although the gravitational force between baryons and dark matter is reduced and even becomes repulsive at low redshift. This slows down the growth of baryonic density perturbations on cosmological scales, while dark matter perturbations are enhanced. In our local environment, the upper bound on the time evolution of Newton's constant requires an efficient screening mechanism that both damps the fifth force on small scales and decouples the local value of Newton constant from its cosmological value. This cannot be achieved by a quasi-static chameleon mechanism, and requires going beyond the quasi-static regime and probably using derivative screenings, such as Kmouflage or Vainshtein screening, on small scales.

We investigate cosmological constraints on K-mouflage models of modified gravity. We consider two scenarios: one where the background evolution is free to deviate from $\Lambda$CDM (K-mouflage) and another one which reproduces a $\Lambda$CDM expansion (K-mimic), implementing both of them into the EFTCAMB code. We discuss the main observational signatures of these models and we compare their cosmological predictions to different datasets, including CMB, CMB lensing, SNIa and different galaxy catalogues. We argue about the possibility of relieving the $H_0$ and weak lensing tensions within these models, finding that K-mouflage scenarios effectively ease the tension on the Hubble Constant. Our final 95\% C.L. bounds on the $\epsilon_{2,0}$ parameter that measures the overall departure from $\Lambda$CDM (corresponding to $\epsilon_{2,0}=0$) are $-0.04\leq \epsilon_{2,0} <0$ for K-mouflage and $0< \epsilon_{2,0} <0.002$ for K-mimic. In the former case the main constraining power comes from changes in the background expansion history, while in the latter case the model is strongly constrained by measurements of the amplitude of matter perturbations. The sensitivity of these cosmological constraints closely matches that of solar system probes. We show that these constraints could be significantly tightened with future ideal probes like CORE.

We analyse the role, on large cosmological scales and laboratory experiments, of the leading curvature squared contributions to the low energy effective action of gravity. We argue for a natural relationship $c_0\lambda^2\simeq 1$ at low-energy between the ${\cal R}^2$ coefficients $c_0$ of the Ricci scalar squared term in this expansion and the dark energy scale $\Lambda=(\lambda M_{\rm Pl})^4$ in four dimensional Planck mass units. We show how the compatibility between the acceleration of the expansion rate of the Universe, local tests of gravity and the quantum stability of the model all converge to select such a relationship up to a coefficient which should be determined experimentally. When embedding this low energy theory of gravity into candidates for its ultraviolet completion, we find that the proposed relationship is guaranteed in string-inspired supergravity models with modulus stabilisation and supersymmetry breaking leading to de Sitter compactifications. In this case, the scalar degree of freedom of ${\cal R}^2$ gravity is associated to a volume modulus. Once written in terms of a scalar-tensor theory, the effective theory corresponds to a massive scalar field coupled with the universal strength $\beta=1/\sqrt{6}$ to the matter stress-energy tensor. When the relationship $c_0\lambda^2\simeq 1$ is realised we find that on astrophysical scales and in cosmology the scalar field is ultralocal and therefore no effect arises on such large scales. On the other hand, the scalar field mass is tightly constrained by the non-observation of fifth forces in torsion pendulum experiments such as E\"ot-Wash. It turns out that the observation of the dark energy scale in cosmology implies that the scalar field could be detectable by fifth force experiments in the near future.

Consistency relations of large-scale structures provide exact nonperturbative results for cross-correlations of cosmic fields in the squeezed limit. They only depend on the equivalence principle and the assumption of Gaussian initial conditions, and remain nonzero at equal times for cross-correlations of density fields with velocity or momentum fields, or with the time derivative of density fields. We show how to apply these relations to observational probes that involve the integrated Sachs-Wolfe effect or the kinematic Sunyaev-Zeldovich effect. In the squeezed limit, this allows us to express the three-point cross-correlations, or bispectra, of two galaxy or matter density fields, or weak lensing convergence fields, with the secondary Cosmic Microwave Background (CMB) distortion in terms of products of a linear and a nonlinear power spectrum. In particular, we find that cross-correlations with the integrated Sachs-Wolfe effect show a specific angular dependence. These results could be used to test the equivalence principle and the primordial Gaussianity, or to check the modeling of large-scale structures.

We investigate how the consistency relations of large-scale structures are modified when the initial density field is not Gaussian. We consider both scenarios where the primordial density field can be written as a nonlinear functional of a Gaussian field and more general scenarios where the probability distribution of the primordial density field can be expanded around the Gaussian distribution, up to all orders over $\delta_{L0}$. Working at linear order over the non-Gaussianity parameters $f_{\rm NL}^{(n)}$ or $S_n$, we find that the consistency relations for the matter density fields are modified as they include additional contributions that involve all-order mixed linear-nonlinear correlations $\langle \prod \delta_L \prod \delta \rangle$. We derive the conditions needed to recover the simple Gaussian form of the consistency relations. This corresponds to scenarios that become Gaussian in the squeezed limit. Our results also apply to biased tracers, and velocity or momentum cross-correlations.

The XXL survey currently covers two 25 sq. deg. patches with XMM observations of ~10ks. We summarise the scientific results associated with the first release of the XXL data set, that occurred mid 2016. We review several arguments for increasing the survey depth to 40 ks during the next decade of XMM operations. X-ray (z<2) cluster, (z<4) AGN and cosmic background survey science will then benefit from an extraordinary data reservoir. This, combined with deep multi-$\lambda$ observations, will lead to solid standalone cosmological constraints and provide a wealth of information on the formation and evolution of AGN, clusters and the X-ray background. In particular, it will offer a unique opportunity to pinpoint the z>1 cluster density. It will eventually constitute a reference study and an ideal calibration field for the upcoming eROSITA and Euclid missions.

We consider K-mouflage models which are K-essence theories coupled to matter. We analyse their quantum properties and in particular the quantum corrections to the classical Lagrangian. We setup the renormalisation programme for these models and show that K-mouflage theories involve a recursive construction whereby each set of counter-terms introduces new divergent quantum contributions which in turn must be subtracted by new counter-terms. This tower of counter-terms can be constructed by recursion and allows one to calculate the finite renormalised action of the model. In particular, the classical action is not renormalised and the finite corrections to the renormalised action contain only higher derivative operators. We establish an operational criterion for classicality, where the corrections to the classical action are negligible, and show that this is satisfied in cosmological and astrophysical situations for (healthy) K-mouflage models which pass the solar system tests. We also find that these models are quantum stable around astrophysical and cosmological backgrounds. We then consider the possible embedding of the K-mouflage models in an Ultra-Violet completion. We find that the healthy models which pass the solar system tests all violate the positivity constraint which would follow from the unitarity of the putative UV completion, implying that these healthy K-mouflage theories have no UV completion. We then analyse their behaviour at high energy and we find that the classicality criterion is satisfied in the vicinity of a high energy collision implying that the classical K-mouflage theory can be applied in this context. Moreover, the classical description becomes more accurate as the energy increases, in a way compatible with the classicalisation concept.

We describe K-mouflage models of modified gravity using the effective field theory of dark energy. We show how the Lagrangian density $K$ defining the K-mouflage models appears in the effective field theory framework, at both the exact fully nonlinear level and at the quadratic order of the effective action. We find that K-mouflage scenarios only generate the operator $(\delta g^{00}_{(u)})^n$ at each order $n$. We also reverse engineer K-mouflage models by reconstructing the whole effective field theory, and the full cosmological behaviour, from two functions of the Jordan-frame scale factor in a tomographic manner. This parameterisation is directly related to the implementation of the K-mouflage screening mechanism: screening occurs when $ K'$ is large in a dense environment such as the deep matter and radiation eras. In this way, K-mouflage can be easily implemented as a calculable subclass of models described by the effective field theory of dark energy which could be probed by future surveys.

We present exact kinematic consistency relations for cosmological structures that do not vanish at equal times and can thus be measured in surveys. These rely on cross-correlations between the density and velocity, or momentum, fields. Indeed, the uniform transport of small-scale structures by long wavelength modes, which cannot be detected at equal times by looking at density correlations only, gives rise to a shift in the amplitude of the velocity field that could be measured. These consistency relations only rely on the weak equivalence principle and Gaussian initial conditions. They remain valid in the non-linear regime and for biased galaxy fields. They can be used to constrain non-standard cosmological scenarios or the large-scale galaxy bias.

We investigate the effects of a K-mouflage modification of gravity on the dynamics of clusters of galaxies. We extend the description of K-mouflage to situations where the scalar field responsible for the modification of gravity is coupled to a perfect fluid with pressure. We describe the coupled system at both the background cosmology and cosmological perturbations levels, focusing on cases where the pressure emanates from small-scale nonlinear physics. We derive these properties in both the Einstein and Jordan frames, as these two frames already differ by a few percents at the background level for K-mouflage scenarios, and next compute cluster properties in the Jordan frame that is better suited to these observations. Galaxy clusters are not screened by the K-mouflage mechanism and therefore feel the modification of gravity in a maximal way. This implies that the halo mass function deviates from $\Lambda$-CDM by a factor of order one for masses $M\gtrsim 10^{14} \ h^{-1} M_\odot$. We then consider the hydrostatic equilibrium of gases embedded in galaxy clusters and the consequences of K-mouflage on the X-ray cluster luminosity, the gas temperature, and the Sunyaev-Zel'dovich effect. We find that the cluster temperature function, and more generally number counts, are largely affected by K-mouflage, mainly due to the increased cluster abundance in these models. Other scaling relations such as the mass-temperature and the temperature-luminosity relations are only modified at the percent level due to the constraints on K-mouflage from local Solar System tests.

Super-chameleon models where all types of matter belong to three secluded sectors, i.e. the dark, supersymmetry breaking and matter sectors, are shown to be dynamically equivalent to ultra-local models of modified gravity. In the dark sector, comprising both dark matter and dark energy, the interaction range between the dark energy field and dark matter is constrained to be extremely short, i.e. shorter than the inverse gravitino mass set by supersymmetry breaking. This realises an extreme version of chameleon screening of the dark energy interaction. On the other hand, the baryonic matter sector decouples from the dark energy in a Damour-Polyakov way. These two mechanisms preclude the existence of any modification of gravity locally in the Solar System due to the presence of the super-chameleon field. On larger scales, the super-chameleon can have effects on the growth of structure and the number of dark matter halos. It can also affect the dynamics of galaxies where the fifth force interaction that it induces can have the same order of magnitude as Newton's interaction.

We present a class of modified-gravity theories which we call ultra-local models. We add a scalar field, with negligible kinetic terms, to the Einstein-Hilbert action. We also introduce a conformal coupling to matter. This gives rise to a new screening mechanism which is not entirely due to the non-linearity of the scalar field potential or the coupling function but to the absence of the kinetic term. As a result this removes any fifth force between isolated objects in vacuum. The predictions of these models only depend on a single free function, as the potential and the coupling function are degenerate, with an amplitude given by a parameter $\alpha \lesssim 10^{-6}$, whose magnitude springs from requiring a small modification of Newton's potential astrophysically and cosmologically. This singles out a redshift $z_{\alpha} \sim \alpha^{-1/3} \gtrsim 100$ where the fifth force is the greatest. The cosmological background follows the $\Lambda$-CDM history within a $10^{-6}$ accuracy, while cosmological perturbations are significantly enhanced (or damped) on small scales, $k \gtrsim 2 h {\rm Mpc}^{-1}$ at $z=0$. The spherical collapse and the halo mass function are modified in the same manner. We find that the modifications of gravity are greater for galactic or sub-galactic structures. We also present a thermodynamic analysis of the non-linear and inhomogeneous fifth-force regime where we find that the Universe is not made more inhomogeneous before $z_\alpha$ when the fifth force dominates, and does not lead to the existence of clumped matter on extra small scales inside halos for large masses while this possibility exists for masses $M\lesssim 10^{11} M_\odot$ where the phenomenology of ultra-local models would be most different from $\Lambda$-CDM.

Weak lensing (WL) promises to be a particularly sensitive probe of both the growth of large scale structure (LSS) as well as the fundamental relation between matter density perturbations and metric perturbations, thus providing a powerful tool with which we may constrain modified theories of gravity (MG) on cosmological scales. Future deep, wide-field WL surveys will provide an unprecedented opportunity to constrain deviations from General Relativity (GR). Employing a three-dimensional (3D) analysis based on the spherical Fourier-Bessel (sFB) expansion, we investigate the extent to which MG theories will be constrained by a typical 3D WL survey configuration including noise from the intrinsic ellipticity distribution $\sigma_{\epsilon}$ of source galaxies. Here we focus on two classes of screened theories of gravity: i) $f(R)$ chameleon models and ii) environmentally dependent dilaton models. We use one-loop perturbation theory combined with halo models in order to accurately model the evolution of matter power-spectrum with redshift in these theories. Using a Fisher information matrix based approach, we show that for an all-sky spectroscopic survey, the parameter $f_{R_0}$ can be constrained in the range $f_{R_0}< 5\times 10^{-6}(9\times 10^{-6})$ for $n=1(2)$ with a 3$\sigma$ confidence level. This can be achieved by using relatively low order angular harmonics $\ell<100$. Including higher order harmonics $\ell>100$ can further tighten the constraints, making them comparable to current solar-system constraints. We also employ a Principal Component Analysis (PCA) in order to study the parameter degeneracies in the MG parameters. Our results can trivially be extended to other MG theories, such as the K-mouflage models. The confusion from intrinsic ellipticity correlation and modification of the matter power-spectrum at small scale due to feedback mechanisms is briefly discussed.

Context. The XXL Survey is the largest survey carried out by the XMM-Newton satellite and covers a total area of 50 square degrees distributed over two fields. It primarily aims at investigating the large-scale structures of the Universe using the distribution of galaxy clusters and active galactic nuclei as tracers of the matter distribution. Aims. This article presents the XXL bright cluster sample, a subsample of 100 galaxy clusters selected from the full XXL catalogue by setting a lower limit of $3\times 10^{-14}\,\mathrm{erg \,s^{-1}cm^{-2}}$ on the source flux within a 1$^{\prime}$ aperture. Methods. The selection function was estimated using a mixture of Monte Carlo simulations and analytical recipes that closely reproduce the source selection process. An extensive spectroscopic follow-up provided redshifts for 97 of the 100 clusters. We derived accurate X-ray parameters for all the sources. Scaling relations were self-consistently derived from the same sample in other publications of the series. On this basis, we study the number density, luminosity function, and spatial distribution of the sample. Results. The bright cluster sample consists of systems with masses between $M_{500}=7\times 10^{13}$ and $3\times 10^{14} M_\odot$, mostly located between $z=0.1$ and 0.5. The observed sky density of clusters is slightly below the predictions from the WMAP9 model, and significantly below the predictions from the Planck 2015 cosmology. In general, within the current uncertainties of the cluster mass calibration, models with higher values of $\sigma_8$ and/or $\Omega_m$ appear more difficult to accommodate. We provide tight constraints on the cluster differential luminosity function and find no hint of evolution out to $z\sim1$. We also find strong evidence for the presence of large-scale structures in the XXL bright cluster sample and identify five new superclusters.

We present the XXL Survey, the largest XMM programme totaling some 6.9 Ms to date and involving an international consortium of roughly 100 members. The XXL Survey covers two extragalactic areas of 25 deg2 each at a point-source sensitivity of ~ 5E-15 erg/sec/cm2 in the [0.5-2] keV band (completeness limit). The survey's main goals are to provide constraints on the dark energy equation of state from the space-time distribution of clusters of galaxies and to serve as a pathfinder for future, wide-area X-ray missions. We review science objectives, including cluster studies, AGN evolution, and large-scale structure, that are being conducted with the support of approximately 30 follow-up programmes. We describe the 542 XMM observations along with the associated multi-lambda and numerical simulation programmes. We give a detailed account of the X-ray processing steps and describe innovative tools being developed for the cosmological analysis. The paper provides a thorough evaluation of the X-ray data, including quality controls, photon statistics, exposure and background maps, and sky coverage. Source catalogue construction and multi-lambda associations are briefly described. This material will be the basis for the calculation of the cluster and AGN selection functions, critical elements of the cosmological and science analyses. The XXL multi-lambda data set will have a unique lasting legacy value for cosmological and extragalactic studies and will serve as a calibration resource for future dark energy studies with clusters and other X-ray selected sources. With the present article, we release the XMM XXL photon and smoothed images along with the corresponding exposure maps. The XMM XXL observation list (Table B.1) is available in electronic form at the CDS. The present paper is the first in a series reporting results of the XXL-XMM survey.

We study Modified Gravity (MG) theories by modelling the redshifted matter power spectrum in a spherical Fourier-Bessel (sFB) basis. We use a fully non-linear description of the real-space matter power-spectrum and include the lowest-order redshift-space correction (Kaiser effect), taking into account some additional non-linear contributions. Ignoring relativistic corrections, which are not expected to play an important role for a shallow survey, we analyse two different modified gravity scenarios, namely the generalised Dilaton scalar-tensor theories and the $f({R})$ models in the large curvature regime. We compute the 3D power spectrum ${\cal C}^s_{\ell}(k_1,k_2)$ for various such MG theories with and without redshift space distortions, assuming precise knowledge of background cosmological parameters. Using an all-sky spectroscopic survey with Gaussian selection function $\varphi(r)\propto \exp(-{r^2 / r^2_0})$, $r_0 = 150 \, h^{-1} \, {\textrm{Mpc}}$, and number density of galaxies $\bar {\textrm{N}} =10^{-4}\;{\textrm{Mpc}}^{-3}$, we use a $\chi^2$ analysis, and find that the lower-order $(\ell \leq 25)$ multipoles of ${\cal C}^s_\ell(k,k')$ (with radial modes restricted to $k < 0.2 \, h \,{\textrm{Mpc}}^{-1}$) can constraint the parameter $f_{R_0}$ at a level of $2\times 10^{-5} (3\times 10^{-5})$ with $3 \sigma$ confidence for $n=1(2)$. Combining constraints from higher $\ell > 25$ modes can further reduce the error bars and thus in principle make cosmological gravity constraints competitive with solar system tests. However this will require an accurate modelling of non-linear redshift space distortions. Using a tomographic $\beta(a)$-$m(a)$ parameterization we also derive constraints on specific parameters describing the Dilaton models of modified gravity.

We show that Solar System tests can place very strong constraints on K-mouflage models of gravity, which are coupled scalar field models with nontrivial kinetic terms that screen the fifth force in regions of large gravitational acceleration. In particular, the bounds on the anomalous perihelion of the Moon imposes stringent restrictions on the K-mouflage Lagrangian density, which can be met when the contributions of higher-order operators in the static regime are sufficiently small. The bound on the rate of change of the gravitational strength in the Solar System constrains the coupling strength $\beta$ to be smaller than $0.1$. These two bounds impose tighter constraints than the results from the Cassini satellite and Big Bang Nucleosynthesis. Despite the Solar System restrictions, we show that it is possible to construct viable models with interesting cosmological predictions. In particular, relative to $\Lambda$-CDM, such models predict percent-level deviations for the clustering of matter and the number density of dark matter haloes. This makes these models predictive and testable by forthcoming observational missions.

We use the cosmic shear data from the Canada-France-Hawaii Telescope Lensing Survey to place constraints on $f(R)$ and {\it Generalized Dilaton} models of modified gravity. This is highly complimentary to other probes since the constraints mainly come from the non-linear scales: maximal deviations with respects to the General-Relativity + $\Lambda$CDM scenario occurs at $k\sim1 h \mbox{Mpc}^{-1}$. At these scales, it becomes necessary to account for known degeneracies with baryon feedback and massive neutrinos, hence we place constraints jointly on these three physical effects. To achieve this, we formulate these modified gravity theories within a common tomographic parameterization, we compute their impact on the clustering properties relative to a GR universe, and propagate the observed modifications into the weak lensing $\xi_{\pm}$ quantity. Confronted against the cosmic shear data, we reject the $f(R)$ $\{ |f_{R_0}|=10^{-4}, n=1\}$ model with more than 99.9% confidence interval (CI) when assuming a $\Lambda$CDM dark matter only model. In the presence of baryonic feedback processes and massive neutrinos with total mass up to 0.2eV, the model is disfavoured with at least 94% CI in all different combinations studied. Constraints on the $\{ |f_{R_0}|=10^{-4}, n=2\}$ model are weaker, but nevertheless disfavoured with at least 89% CI. We identify several specific combinations of neutrino mass, baryon feedback and $f(R)$ or Dilaton gravity models that are excluded by the current cosmic shear data. Notably, universes with three massless neutrinos and no baryon feedback are strongly disfavoured in all modified gravity scenarios studied. These results indicate that competitive constraints may be achieved with future cosmic shear data.

The clustering ratio $\eta$, a large-scale structure observable originally designed to constrain the shape of the power spectrum of matter density fluctuations, is shown to provide a sensitive probe of the nature of gravity in the cosmological regime. We apply this analysis to $F(R)$ theories of gravity using the luminous red galaxy (LRG) sample extracted from the spectroscopic Sloan Digital Sky Survey (SDSS) data release 7 and 10 catalogues. We find that General Relativity (GR), complemented with a Friedmann-Robertson-Walker (FRW) cosmological model with parameters fixed by the Planck satellite, describes extremely well the clustering of galaxies up to $z\sim 0.6$. On large cosmic scales, the absolute amplitude of deviations from GR, $|f_{R_0 }|$, is constrained to be smaller than $4.6 \times 10^{-5}$ at the $95\%$ confidence level. This bound makes cosmological probes of gravity almost competitive with the sensitivity of Solar System tests, although still one order of magnitude less effective than astrophysical tests. We also extrapolate our results to future large surveys like Euclid and show that the astrophysical bound will certainly remain out of reach for such a class of modified-gravity models that only differ from $\Lambda$CDM at low redshifts.

We present the redshift-space generalization of the equal-time angular-averaged consistency relations between $(\ell+n)$- and $n$-point polyspectra of the cosmological matter density field. Focusing on the case of $\ell=1$ large-scale mode and $n$ small-scale modes, we use an approximate symmetry of the gravitational dynamics to derive explicit expressions that hold beyond the perturbative regime, including both the large-scale Kaiser effect and the small-scale fingers-of-god effects. We explicitly check these relations, both perturbatively, for the lowest-order version that applies to the bispectrum, and nonperturbatively, for all orders but for the one-dimensional dynamics. Using a large ensemble of $N$-body simulations, we find that our squeezed bispectrum relation is valid to better than $20\%$ up to $1h$Mpc$^{-1}$, for both the monopole and quadrupole at $z=0.35$, in a $\Lambda$CDM cosmology. Additional simulations done for the Einstein-de Sitter background suggest that these discrepancies mainly come from the breakdown of the approximate symmetry of the gravitational dynamics. For practical applications, we introduce a simple ansatz to estimate the new derivative terms in the relation using only observables. Although the relation holds worse after using this ansatz, we can still recover it within $20\%$ up to $1h$Mpc$^{-1}$, at $z=0.35$ for the monopole. On larger scales, $k = 0.2 h\mathrm{Mpc}^{-1}$, it still holds within the statistical accuracy of idealized simulations of volume $\sim8h^{-3}\mathrm{Gpc}^3$ without shot-noise error.

We present a comprehensive derivation of linear perturbation equations for different matter species, including photons, baryons, cold dark matter, scalar fields, massless and massive neutrinos, in the presence of a generic conformal coupling. Starting from the Lagrangians, we show how the conformal transformation affects the dynamics. In particular, we discuss how to incorporate consistently the scalar coupling in the equations of the Boltzmann hierarchy for massive neutrinos and the subsequent fluid approximations. We use the recently proposed K-mouflage model as an example to demonstrate the numerical implementation of our linear perturbation equations. K-mouflage is a new mechanism to suppress the fifth force between matter particles induced by the scalar coupling, but in the linear regime the fifth force is unsuppressed and can change the clustering of different matter species in different ways. We show how the CMB, lensing potential and matter power spectra are affected by the fifth force, and find ranges of K-mouflage parameters whose effects could be seen observationally. We also find that the scalar coupling can have the nontrivial effect of shifting the amplitude of the power spectra of the lensing potential and density fluctuations in opposite directions, although both probe the overall clustering of matter. This paper can serve as a reference for those who work on generic coupled scalar field cosmology, or those who are interested in the cosmological behaviour of the K-mouflage model.

We describe how the kinematic consistency relations satisfied by density correlations of the large-scale structures of the Universe can be derived within the usual Newtonian framework. These relations express a kinematic effect and show how the $(\ell+n)$-density correlation factors in terms of the $n$-point correlation and $\ell$ linear power spectrum factors, in the limit where the $\ell$ soft wave numbers become linear and much smaller than the $n$ other wave numbers. We describe how these relations extend to multifluid cases. In the standard cosmology, these consistency relations derive from the equivalence principle. A detection of their violation would indicate non-Gaussian initial conditions, non-negligible decaying modes, or a modification of gravity that does not converge to General Relativity on large scales.

We investigate the small-scale static configurations of K-mouflage models defined by a general function $K(\chi)$ of the kinetic terms. The fifth force is screened by the nonlinear K-mouflage mechanism if $K'(\chi)$ grows sufficiently fast for large negative $\chi$. In the general non-spherically symmetric case, the fifth force is not aligned with the Newtonian force. For spherically symmetric static matter density profiles, the results depend on the potential function $W_{-}(y) = y K'(-y^2/2)$, which must be monotonically increasing to $+\infty$ for $y \geq 0$ to guarantee the existence of a single solution throughout space for any matter density profile. Small radial perturbations around these static profiles propagate as traveling waves with a velocity greater than the speed of light. Starting from vanishing initial conditions for the scalar field and for a time-dependent matter density corresponding to the formation of an overdensity, we numerically check that the scalar field converges to the static solution. If $W_{-}$ is bounded, for high-density objects there are no static solutions throughout space, but one can still define a static solution restricted to large radii. Our dynamical study shows that the scalar field relaxes to this static solution at large radii, whereas spatial gradients keep growing with time at smaller radii. If $W_{-}$ is not bounded but non-monotonic, there is an infinite number of discontinuous static solutions but they are not physical and these models are not theoretically sound. Such K-mouflage scenarios provide an example of theories that can appear viable at the cosmological level, for the cosmological background and perturbative analysis, while being meaningless at a nonlinear level for small-scale configurations. This shows the importance of small-scale nonlinear analysis of screening models.

We study structure formation in K-mouflage cosmology whose main feature is the absence of screening effect on quasilinear scales. We show that the growth of structure at the linear level is affected by both a new time dependent Newton constant and a friction term which depend on the background evolution. These combine with the modified background evolution to change the growth rate by up to ten percent since $z\sim 2$. At the one loop level, we find that the nonlinearities of the K-mouflage models are mostly due to the matter dynamics and that the scalar perturbations can be treated at tree level. We also study the spherical collapse in K-mouflage models and show that the critical density contrast deviates from its $\Lambda$-CDM value and that, as a result, the halo mass function is modified for large masses by an order one factor. Finally we consider the deviation of the matter spectrum from $\Lambda$-CDM on nonlinear scales where a halo model is utilized. We find that the discrepancy peaks around $1\ h{\rm Mpc}^{-1}$ with a relative difference which can reach fifty percent. Importantly, these features are still true at larger redshifts, contrary to models of the chameleon-$f(R)$ and Galileon types.

We study the cosmology of K-mouflage theories at the background level. We show that the effects of the scalar field are suppressed at high matter density in the early Universe and only play a role in the late time Universe where the deviations of the Hubble rate from its $\Lambda$-CDM counterpart can be of the order five percent for redshifts $1 \lesssim z \lesssim 5$. Similarly, we find that the equation of state can cross the phantom divide in the recent past and even diverge when the effective scalar energy density goes negative and subdominant compared to matter, preserving the positivity of the squared Hubble rate. These features are present in models for which Big Bang Nucleosynthesis is not affected. We analyze the fate of K-mouflage when the nonlinear kinetic terms give rise to ghosts, particle excitations with negative energy. In this case, we find that the K-mouflage theories can only be considered as an effective description of the Universe at low energy below $1$ keV. In the safe ghost-free models, we find that the equation of state always diverges in the past and changes significantly by a few percent since $z\lesssim 1$.

This paper presents 52 X-ray bright galaxy clusters selected within the 11 deg$^2$ XMM-LSS survey. 51 of them have spectroscopic redshifts ($0.05<z<1.06$), one is identified at $z_{\rm phot}=1.9$, and all together make the high-purity "Class 1" (C1) cluster sample of the XMM-LSS, the highest density sample of X-ray selected clusters with a monitored selection function. Their X-ray fluxes, averaged gas temperatures (median $T_X=2$ keV), luminosities (median $L_{X,500}=5\times10^{43}$ ergs/s) and total mass estimates (median $5\times10^{13} h^{-1} M_{\odot}$) are measured, adapting to the specific signal-to-noise regime of XMM-LSS observations. The redshift distribution of clusters shows a deficit of sources when compared to the cosmological expectations, regardless of whether WMAP-9 or Planck-2013 CMB parameters are assumed. This lack of sources is particularly noticeable at $0.4 \lesssim z \lesssim 0.9$. However, after quantifying uncertainties due to small number statistics and sample variance we are not able to put firm (i.e. $>3 \sigma$) constraints on the presence of a large void in the cluster distribution. We work out alternative hypotheses and demonstrate that a negative redshift evolution in the normalization of the $L_{X}-T_X$ relation (with respect to a self-similar evolution) is a plausible explanation for the observed deficit. We confirm this evolutionary trend by directly studying how C1 clusters populate the $L_{X}-T_X-z$ space, properly accounting for selection biases. We point out that a systematically evolving, unresolved, central component in clusters and groups (AGN contamination or cool core) can impact the classification as extended sources and be partly responsible for the observed redshift distribution.[abridged]

We explicitly test the equal-time consistency relation between the angular-averaged bispectrum and the power spectrum of the matter density field, employing a large suite of cosmological $N$-body simulations. This is the lowest-order version of the relations between $(\ell+n)$-point and $n$-point polyspectra, where one averages over the angles of $\ell$ soft modes. This relation depends on two wave numbers, $k'$ in the soft domain and $k$ in the hard domain. We show that it holds up to a good accuracy, when $k'/k\ll 1$ and $k'$ is in the linear regime, while the hard mode $k$ goes from linear ($0.1\,h\mathrm{Mpc}^{-1}$) to nonlinear ($1.0\,h\mathrm{Mpc}^{-1}$) scales. On scales $k\lesssim 0.4\,h\mathrm{Mpc}^{-1}$, we confirm the relation within the statistical error of the simulations (typically a few percent depending on the wave number), even though the bispectrum can already deviate from leading-order perturbation theory by more than $30\%$. We further examine the relation on smaller scales with higher resolution simulations. We find that the relation holds within the statistical error of the simulations at $z=1$, whereas we find deviations as large as $\sim 7\%$ at $k \sim 1.0\,h\mathrm{Mpc}^{-1}$ at $z=0.35$. We show that this can be explained partly by the breakdown of the approximation $\Omega_\mathrm{m}/f^2\simeq1$ with supplemental simulations done in the Einstein-de Sitter background cosmology. We also estimate the impact of this approximation on the power spectrum and bispectrum.

The cosmological dynamics of gravitational clustering satisfies an approximate invariance with respect to the cosmological parameters that is often used to simplify analytical computations. We describe how this approximate symmetry gives rise to angular averaged consistency relations for the matter density correlations. This allows one to write the $(\ell+n)$ density correlation, with $\ell$ large-scale linear wave numbers that are integrated over angles, and $n$ fixed small-scale nonlinear wave numbers, in terms of the small-scale $n$-point density correlation and $\ell$ prefactors that involve the linear power spectra at the large-scale wave numbers. These relations, which do not vanish for equal-time statistics, go beyond the already known kinematic consistency relations. They could be used to detect primordial non-Gaussianities, modifications of gravity, limitations of galaxy biasing schemes, or to help designing analytical models of gravitational clustering.

We investigate the possible accuracy that can be reached by analytical models for the matter density power spectrum and correlation function. Using a realistic description of the power spectrum that combines perturbation theory with a halo model, we study the convergence rate of several perturbative expansion schemes and the impact of nonperturbative effects, as well as the sensitivity to phenomenological halo parameters. We check that the simple reorganization of the standard perturbative expansion, with a Gaussian damping prefactor, provides a well-ordered convergence and a finite correlation function that yields a percent accuracy at the baryon acoustic oscillation peak (as soon as one goes to second order). Lagrangian-space expansions are somewhat more efficient, when truncated at low orders, but may diverge at high orders. We find that whereas the uncertainty on the halo-profile mass-concentration relation is not a strong limitation, the uncertainty on the halo mass function can severely limit the accuracy of theoretical predictions for $P(k)$ (this also applies to the power spectra measured in numerical simulations). The real-space correlation function provides a better separation between perturbative and nonperturbative effects, which are restricted to $x \lesssim 10 h^{-1}$Mpc at all redshifts.

We estimate the amplitude of the source-lens clustering bias and of the intrinsic-alignment bias of weak lensing estimators of the two-point and three-point convergence and cosmic-shear correlation functions. We use a linear galaxy bias model for the galaxy-density correlations, as well as a linear intrinsic-alignment model. For the three-point and four-point density correlations, we use analytical or semi-analytical models, based on a hierarchical ansatz or a combination of one-loop perturbation theory with a halo model. For two-point statistics, we find that the source-lens clustering bias is typically several orders of magnitude below the weak lensing signal, except when we correlate a very low-redshift galaxy ($z_2 \la 0.05$) with a higher redshift galaxy ($z_1 \ga 0.5$), where it can reach $10\%$ of the signal for the shear. For three-point statistics, the source-lens clustering bias is typically of order $10\%$ of the signal, as soon as the three galaxy source redshifts are not identical. The intrinsic-alignment bias is typically about $10\%$ of the signal for both two-point and three-point statistics. Thus, both source-lens clustering bias and intrinsic-alignment bias must be taken into account for three-point estimators aiming at a better than $10\%$ accuracy.

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.

We study the effects of screened modified gravity of the $f(R)$, dilaton and symmetron types on structure formation, from the quasi-linear to the non-linear regime, using semi-analytical methods. For such models, where the range of the new scalar field is typically within the Mpc range and below in the cosmological context, non-linear techniques are required to understand the deviations of the power spectrum of the matter density contrast compared to the $\Lambda$-CDM template. This is nowadays commonly tackled using extensive N-body simulations. Here we present new results combining exact perturbation theory at the one loop level (and a partial resummation of the perturbative series) with a halo model. The former allows one to extend the linear perturbative analysis up to $k\lesssim 0.15{\rm h Mpc}^{-1}$ at the perturbative level while the latter leads to a reasonable, up to a few percent, agreement with numerical simulations for $k\lesssim 3{\rm h Mpc}^{-1}$ for large curvature $f(R)$ models, and $k\lesssim 1{\rm h Mpc}^{-1}$ for dilatons and symmetrons, at $z=0$. We also discuss how the behaviors of the perturbative expansions and of the spherical collapse differ for $f(R)$, dilaton, and symmetron models.

We investigate whether the late-time (at $z\leq 100$) velocity dispersion expected in Warm Dark Matter scenarios could have some effect on the cosmic web (i.e., outside of virialized halos). We consider effective hydrodynamical equations, with a pressurelike term that agrees at the linear level with the analysis of the Vlasov equation. Then, using analytical methods, based on perturbative expansions and the spherical dynamics, we investigate the impact of this term for a 1 keV dark matter particle. We find that the late-time velocity dispersion has a negligible effect on the power spectrum on perturbative scales and on the halo mass function. However, it has a significant impact on the probability distribution function of the density contrast at $z \sim 3$ on scales smaller than $0.1 h^{-1}$Mpc, which correspond to Lyman-$\alpha$ clouds. Finally, we note that numerical simulations should start at $z_i\geq 100$ rather than $z_i \leq 50$ to avoid underestimating gravitational clustering at low redshifts.

Large ongoing and upcoming galaxy cluster surveys in the optical, X-ray and millimetric wavelengths will provide rich samples of galaxy clusters at unprecedented depths. One key observable for constraining cosmological models is the correlation function of these objects, measured through their spectroscopic redshift. We study the redshift-space correlation functions of clusters of galaxies, averaged over finite redshift intervals, and their covariance matrices. Expanding as usual the angular anisotropy of the redshift-space correlation on Legendre polynomials, we consider the redshift-space distortions of the monopole as well as the next two multipoles, $2\ell=2$ and 4. Taking into account the Kaiser effect, we developed an analytical formalism to obtain explicit expressions of all contributions to these mean correlations and covariance matrices. We include shot-noise and sample-variance effects as well as Gaussian and non-Gaussian contributions. We obtain a reasonable agreement with numerical simulations for the mean correlations and covariance matrices on large scales ($r> 10 h^{-1}$Mpc). Redshift-space distortions amplify the monopole correlation by about $10-20%$, depending on the halo mass, but the signal-to-noise ratio remains of the same order as for the real-space correlation. This distortion will be significant for surveys such as DES, Erosita, and Euclid, which should also measure the quadrupole $2\ell=2$. The third multipole, $2\ell=4$, may only be marginally detected by Euclid.

We present a new approach to computing the matter density power spectrum, from large linear scales to small highly nonlinear scales. Instead of explicitly computing a partial series of high-order diagrams, as in perturbative resummation schemes, we embed the standard perturbation theory within a realistic nonlinear Lagrangian-space ansatz. We also point out that an "adhesion-like" regularization of the shell-crossing regime is more realistic than a "Zel'dovich-like" behavior, where particles freely escape to infinity. This provides a "cosmic web" power spectrum with good small-scale properties that provide a good matching with a halo model on mildly nonlinear scales. We obtain a good agreement with numerical simulations on large scales, better than 3% for $k\leq 1 h$Mpc$^{-1}$, and on small scales, better than 10% for $k \leq 10 h$Mpc$^{-1}$, at $z \geq 0.35$, which improves over previous methods.

We study the growth of structures in modified gravity models where the Poisson equation and the relationship between the two Newtonian potentials are modified by explicit functions of space and time. This parameterisation applies to the $f(R)$ models and more generally to screened modified gravity models. We investigate the linear and weakly nonlinear regimes using the "standard" perturbative approach and a resummation technique, while we use the spherical dynamics to go beyond low-order results. This allows us to estimate the matter density power spectrum and bispectrum from linear to highly nonlinear scales, the full probability distribution of the density contrast on weakly nonlinear scales, and the halo mass function. We analyse the impact of modifications of gravity on these quantities for a few realistic models. In particular, we find that the standard one-loop perturbative approach is not sufficiently accurate to probe these effects on the power spectrum and it is necessary to use resummation methods even on weakly nonlinear scales which provide the best observational window for modified gravity as relative deviations from General Relativity do not grow significantly on smaller scales where theoretical predictions become increasingly difficult.

We investigate the performance of an analytic model of the 3D matter distribution, which combines perturbation theory with halo models, for weak-lensing configuration-space statistics. We compared our predictions for the weak-lensing convergence two-point and three-point correlation functions with numerical simulations and fitting formulas proposed in previous works. We also considered the second- and third-order moments of the smoothed convergence and of the aperture-mass. As in our previous study of Fourier-space weak-lensing statistics, we find that our model agrees better with simulations than previously published fitting formulas. Moreover, we recover the dependence on cosmology of these weak-lensing statistics and we can describe multi-scale moments. This approach allows us to obtain the quantitative relationship between these integrated weak-lensing statistics and the various contributions to the underlying 3D density fluctuations, decomposed over perturbative, two-halo, or one-halo terms.

We investigate the performance of an analytic model of the 3D matter distribution, which combines perturbation theory with halo models, for weak-lensing statistics. We compare our predictions for the weak-lensing convergence power spectrum and bispectrum with numerical simulations and fitting formulas proposed in previous works. We find that this model provides better agreement with simulations than published fitting formulas. This shows that building on systematic and physically motivated models is a promising approach. Moreover, this makes explicit the link between the weak-lensing statistics and the underlying properties of the 3D matter distribution, as a function of scale $\ell$. Thus, we obtain the contributions to the lensing power spectrum and bispectrum that arise from perturbative terms (complete up to one-loop) and nonperturbative terms (e.g., "1-halo" term). Finally, we show that this approach recovers the dependence on cosmology (for realistic scenarios).

We study the mean number counts and two-point correlation functions, along with their covariance matrices, of cosmological surveys such as for clusters. In particular, we consider correlation functions averaged over finite redshift intervals, which are well suited to cluster surveys or populations of rare objects, where one needs to integrate over nonzero redshift bins to accumulate enough statistics. We develop an analytical formalism to obtain explicit expressions of all contributions to these means and covariance matrices, taking into account both shot-noise and sample-variance effects. We compute low-order as well as high-order (including non-Gaussian) terms. We derive expressions for the number counts per redshift bins both for the general case and for the small window approximation. We estimate the range of validity of Limber's approximation and the amount of correlation between different redshift bins. We also obtain explicit expressions for the integrated 3D correlation function and the 2D angular correlation. We compare the relative importance of shot-noise and sample-variance contributions, and of low-order and high-order terms. We check the validity of our analytical results through a comparison with the Horizon full-sky numerical simulations, and we obtain forecasts for several future cluster surveys.

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.

We investigate how unified models should be built to be able to predict the matter-density bispectrum (and power spectrum) from very large to small scales and that are at the same time consistent with perturbation theory at low $k$ and with halo models at high $k$. We use a Lagrangian framework to decompose the bispectrum into "3-halo", "2-halo", and "1-halo" contributions, related to "perturbative" and "non-perturbative" terms. We describe a simple implementation of this approach and present a detailed comparison with numerical simulations. We show that the 1-halo and 2-halo contributions contain counterterms that ensure their decay at low $k$, as required by physical constraints, and allow a better match to simulations. Contrary to the power spectrum, the standard 1-loop perturbation theory can be used for the perturbative 3-halo contribution because it does not grow too fast at high $k$. Moreover, it is much simpler and more accurate than two resummation schemes investigated in this paper. We obtain a good agreement with numerical simulations on both large and small scales, but the transition scales are poorly described by the simplest implementation. This cannot be amended by simple modifications to the halo parameters, but we show how it can be corrected for the power spectrum and the bispectrum through a simple interpolation scheme that is restricted to this intermediate regime. Then, we reach an accuracy on the order of 10% on mildly and highly nonlinear scales, while an accuracy on the order of 1% is obtained on larger weakly nonlinear scales. This also holds for the real-space two-point correlation function.

We investigate the building of unified models that can predict the matter-density power spectrum and the two-point correlation function from very large to small scales, being consistent with perturbation theory at low $k$ and with halo models at high $k$. We use a Lagrangian framework to re-interpret the halo model and to decompose the power spectrum into "2-halo" and "1-halo" contributions, related to "perturbative" and "non-perturbative" terms. We describe a simple implementation of this model and present a detailed comparison with numerical simulations, from $k \sim 0.02$ up to $100 h$Mpc$^{-1}$, and from $x \sim 0.02$ up to $150 h^{-1}$Mpc. We show that the 1-halo contribution contains a counterterm that ensures a $k^2$ tail at low $k$ and is important not to spoil the predictions on the scales probed by baryon acoustic oscillations, $k \sim 0.02$ to $0.3 h$Mpc$^{-1}$. On the other hand, we show that standard perturbation theory is inadequate for the 2-halo contribution, because higher order terms grow too fast at high $k$, so that resummation schemes must be used. We describe a simple implementation, based on a 1-loop "direct steepest-descent" resummation for the 2-halo contribution that allows fast numerical computations, and we check that we obtain a good match to simulations at low and high $k$. Our simple implementation already fares better than standard 1-loop perturbation theory on large scales and simple fits to the power spectrum at high $k$, with a typical accuracy of 1% on large scales and 10% on small scales. We obtain similar results for the two-point correlation function. However, there remains room for improvement on the transition scale between the 2-halo and 1-halo contributions, which may be the most difficult regime to describe.

We study the effect of nonperturbative corrections associated with the behavior of particles after shell crossing on the matter power spectrum. We compare their amplitude with the perturbative terms that can be obtained within the fluid description of the system, to estimate the range of scales where such perturbative approaches are relevant. We use the simple Zeldovich dynamics as a benchmark, as it allows the exact computation of the full nonlinear power spectrum and of perturbative terms at all orders. Then, we introduce a "sticky model" that coincides with the Zeldovich dynamics before shell crossing but shows a different behavior afterwards. Thus, their power spectra only differ in their nonperturbative terms. We consider both the real-space and redshift-space power spectra. We find that the potential of perturbative schemes is greater at higher redshift for a $\Lambda$CDM cosmology. For the real-space power spectrum, one can go up to order $66$ of perturbation theory at $z=3$, and to order $9$ at $z=0$, before the nonperturbative correction surpasses the perturbative correction of that order. This allows us to increase the upper bound on $k$ where systematic theoretical predictions may be obtained by perturbative schemes, beyond the linear regime, by a factor $\sim 26$ at $z=3$ and $\sim 6.5$ at $z=0$. This provides a strong motivation to study perturbative resummation schemes, especially at high redshifts $z \geq 1$. We find similar results for the redshift-space power spectrum, with characteristic wavenumbers that are shifted to lower values as redshift-space distortions amplify higher order terms of the perturbative expansions while decreasing the resummed nonlinear power at high $k$.

We revisit the derivation of the mass function and the bias of dark matter halos for non-Gaussian initial conditions. We use a steepest-descent approach to point out that exact results can be obtained for the high-mass tail of the halo mass function and the two-point correlation of massive halos. Focusing on primordial non-Gaussianity of the local type, we check that these results agree with numerical simulations. The high-mass cutoff of the halo mass function takes the same form as the one obtained from the Press-Schechter formalism, but with a linear threshold $\delta_L$ that depends on the definition of the halo. We show that a simple formula, which obeys this high-mass asymptotic and uses the fit obtained for Gaussian initial conditions, matches numerical simulations while keeping the mass function normalized to unity. Next, by deriving the real-space halo two-point correlation in the spirit of Kaiser (1984) and taking a Fourier transform, we obtain good agreement with simulations for the correction to the halo bias due to primordial non-Gaussianity. Therefore, neither the halo mass function nor the bias require an ad-hoc parameter $q$ provided one uses the correct linear threshold $\delta_L$ and pays attention to halo displacements. The nonlinear real-space expression can be useful for checking that the "linearized" bias is a valid approximation. Moreover, it clearly shows how the baryon acoustic oscillation at $\sim 100 h^{-1}$Mpc is amplified by the bias of massive halos and modified by primordial non-Gaussianity. On smaller scales, the correction to the real-space bias roughly scales as $\fNL \, b_M(\fNL=0) \, x^2$. The low-$k$ behavior of the halo bias does not imply a divergent real-space correlation, so that one does not need to introduce counterterms that depend on the survey size.

We build a simple analytical model for the bias of dark matter halos that applies to objects defined by an arbitrary density threshold, $200\leq\deltas\leq 1600$, and that provides accurate predictions from low-mass to high-mass halos. We point out that it is possible to build simple and efficient models, with no free parameter for the halo bias, by using integral constraints that govern the behavior of low-mass and typical halos, whereas the properties of rare massive halos are derived through explicit asymptotic approaches. We also describe how to take into account the impact of halo motions on their bias, using their linear displacement field. We obtain a good agreement with numerical simulations for the halo mass functions and large-scale bias at redshifts $0\leq z \leq 2.5$, for halos defined by a nonlinear density threshold $200\leq\deltas\leq 1600$. We also evaluate the impact on the halo bias of two common approximations, i) neglecting halo motions, and ii) linearizing the halo two-point correlation.

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.

We revisit the study of the mass functions and the bias of dark matter halos. Focusing on the limit of rare massive halos, we point out that exact analytical results can be obtained for the large-mass tail of the halo mass function. This is most easily seen from a steepest-descent approach, that becomes asymptotically exact for rare events. We also revisit the traditional derivation of the bias of massive halos, associated with overdense regions in the primordial density field. We check that the theoretical large-mass cutoff agrees with the mass functions measured in numerical simulations. For halos defined by a nonlinear threshold $\delta=200$ this corresponds to using a linear threshold $\delta_L\simeq 1.59$ instead of the traditional value $\simeq 1.686$. We also provide a fitting formula that matches simulations over all mass scales and obeys the exact large-mass tail. Next, paying attention to the Lagrangian-Eulerian mapping (i.e. corrections associated with the motions of halos), we improve the standard analytical formula for the bias of massive halos. We check that our prediction, which contains no free parameter, agrees reasonably well with numerical simulations. In particular, it recovers the steepening of the dependence on scale of the bias that is observed at higher redshifts, which published fitting formulae did not capture. This behavior mostly arises from nonlinear biasing.

We study the one-dimensional Burgers equation in the inviscid limit for Brownian initial velocity (i.e. the initial velocity is a two-sided Brownian motion that starts from the origin x=0). We obtain the one-point distribution of the velocity field in closed analytical form. In the limit where we are far from the origin, we also obtain the two-point and higher-order distributions. We show how they factorize and recover the statistical invariance through translations for the distributions of velocity increments and Lagrangian increments. We also derive the velocity structure functions and we recover the bifractality of the inverse Lagrangian map. Then, for the case where the initial density is uniform, we obtain the distribution of the density field and its $n$-point correlations. In the same limit, we derive the $n-$point distributions of the Lagrangian displacement field and the properties of shocks. We note that both the stable-clustering ansatz and the Press-Schechter mass function, that are widely used in the cosmological context, happen to be exact for this one-dimensional version of the adhesion model.

We revisit the one-dimensional Burgers equation in the inviscid limit for white-noise initial velocity. We derive the probability distributions of velocity and Lagrangian increments, measured on intervals of any length $x$. This also gives the velocity structure functions. Next, for the case where the initial density is uniform, we obtain the distribution of the density, over any scale $x$, and we derive the density two-point correlation and power spectrum. Finally, we consider the Lagrangian displacement field and we derive the distribution of increments of the Lagrangian map. We check that this gives back the well-known mass function of shocks. For all distributions we describe the limiting scaling functions that are obtained in the large-scale and small-scale limits. We also discuss how these results generalize to other initial conditions, or to higher dimensions, and make the connection with a heuristic multifractal formalism. We note that the formation of point-like masses generically leads to a universal small-scale scaling for the density distribution, which is known as the ``stable-clustering ansatz'' in the cosmological context (where the Burgers dynamics is also known as the ``adhesion model'').

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 develop optimised estimators of two sorts of power spectra for fields defined on the sky, in the presence of partial sky coverage. The first is the cross-power spectrum of two fields on the sky; the second is the skew spectrum of three fields. These can probe the Integrated Sachs Wolfe Effect (ISW) at large angular scales and the Sunyaev Z\'eldovich (SZ) effect from hot gas in clusters at small angular scales. The skew spectrum, recently introduced by Munshi & Heavens (2009), is an optimised statistic which can be tuned to study a particular form of non-Gaussianity, such as may arise in the early Universe, but which retains information on the nature of non-Gaussianity. In this paper we develop the mathematical formalism for the skew spectrum of 3 different fields. When applied to the CMB, this allows us to explore the contamination of the skew spectrum by secondary sources of CMB fluctuations. Considering the three-point function, the study of the bispectrum provides valuable information regarding cross-correlation of secondaries with lensing of CMB with much higher significance compared to just the study involving CMB sky alone. We use our analytical models to study specific cases of cosmological interest which include cross-correlating CMB with various large scale tracers to probe ISW and SZ effects for cross spectral analysis and use the formalism to study the signal-to-noise ratio for detection of the weak lensing of the CMB by cross-correlating it with different tracers as well as point sources for CMB experiments such as Planck (abridged).

We study the decaying Burgers dynamics in $d$ dimensions for random Gaussian initial conditions. We focus on power-law initial energy spectra, such that the system shows a self-similar evolution. This is the case of interest for the "adhesion model" in cosmology and a standard framework for "decaying Burgers turbulence". We briefly describe how the system can be studied through perturbative expansions at early time or large scale (quasi-linear regime). Next, we develop a saddle-point method, based on spherical instantons, that allows to obtain the asymptotic probability distributions $\cP(\eta_r)$ and $\cP(\ctheta_r)$, of the density and velocity increment over spherical cells, reached in the quasi-linear regime. Finally, we show how this approach can be extended to take into account the formation of shocks and we derive the rare-event tails of these probability distributions, at any finite time and scale. This also gives the high-mass tail of the mass function of point-like singularities (shocks in the one dimensional case).

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.

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.

We apply various expansion schemes that may be used to study gravitational clustering to the simple case of the Zeldovich dynamics. Using the well-known exact solution of the Zeldovich dynamics we can compare the predictions of these various perturbative methods with the exact nonlinear result and study their convergence properties. We find that most systematic expansions fail to recover the decay of the response function in the highly nonlinear regime. ``Linear methods'' lead to increasingly fast growth in the nonlinear regime for higher orders, except for Pade approximants that give a bounded response at any order. ``Nonlinear methods'' manage to obtain some damping at one-loop order but they fail at higher orders. Although it recovers the exact Gaussian damping, a resummation in the high-k limit is not justified very well as the generation of nonlinear power does not originate from a finite range of wavenumbers (hence there is no simple separation of scales). No method is able to recover the relaxation of the matter power spectrum on highly nonlinear scales. It is possible to impose a Gaussian cutoff in a somewhat ad-hoc fashion to reproduce the behavior of the exact two-point functions for two different times. However, this cutoff is not directly related to the clustering of matter and disappears in exact equal-time statistics such as the matter power spectrum. On a quantitative level, the usual perturbation theory, and the nonlinear scheme to which one adds an ansatz for the response function with such a Gaussian cutoff, are the two most efficient methods. These results should hold for the gravitational dynamics as well (this has been checked at one-loop order), since the structure of the equations of motion is identical for both dynamics.

We develop a path-integral formalism to study the formation of large-scale structures in the universe. Starting from the equations of motion of hydrodynamics (single-stream approximation) we derive the action which describes the statistical properties of the density and velocity fields for Gaussian initial conditions. Then, we present large-N expansions (associated with a generalization to N fields or with a semi-classical expansion) of the path-integral defined by this action. This provides a systematic expansion for two-point functions such as the response function and the usual two-point correlation. We present the results of two such expansions (and related variants) at one-loop order for a SCDM and a LCDM cosmology. We find that the response function exhibits fast oscillations in the non-linear regime with an amplitude which either follows the linear prediction (for the direct steepest-descent scheme) or decays (for the 2PI effective action scheme). On the other hand, the correlation function agrees with the standard one-loop result in the quasi-linear regime and remains well-behaved in the highly non-linear regime. This suggests that these large-N expansions could provide a good framework to study the dynamics of gravitational clustering in the non-linear regime. Moreover, the use of various expansion schemes allows one to estimate their range of validity without the need of N-body simulations and could provide a better accuracy in the weakly non-linear regime.

Weak gravitational lensing is responsible for the shearing and magnification of the images of high-redshift sources due to the presence of intervening matter. The distortions are due to fluctuations in the gravitational potential, and are directly related to the distribution of matter and to the geometry and dynamics of the Universe. As a consequence, weak gravitational lensing offers unique possibilities for probing the Dark Matter and Dark Energy in the Universe. In this review, we summarise the theoretical and observational state of the subject, focussing on the statistical aspects of weak lensing, and consider the prospects for weak lensing surveys in the future. Weak gravitational lensing surveys are complementary to both galaxy surveys and cosmic microwave background (CMB) observations as they probe the unbiased non-linear matter power spectrum at modest redshifts. Most of the cosmological parameters are accurately estimated from CMB and large-scale galaxy surveys, so the focus of attention is shifting to understanding the nature of Dark Matter and Dark Energy. On the theoretical side, recent advances in the use of 3D information of the sources from photometric redshifts promise greater statistical power, and these are further enhanced by the use of statistics beyond two-point quantities such as the power spectrum. The use of 3D information also alleviates difficulties arising from physical effects such as the intrinsic alignment of galaxies, which can mimic weak lensing to some extent. (Abridged)

We study the relaxation towards thermodynamical equilibrium of a 1-D gravitational system. This OSC model shows a series of critical energies $E_{cn}$ where new equilibria appear and we focus on the homogeneous ($n=0$), one-peak ($n=\pm 1$) and two-peak ($n=2$) states. Using numerical simulations we investigate the relaxation to the stable equilibrium $n=\pm 1$ of this $N-$body system starting from initial conditions defined by equilibria $n=0$ and $n=2$. We find that in a fashion similar to other long-range systems the relaxation involves a fast violent relaxation phase followed by a slow collisional phase as the system goes through a series of quasi-stationary states. Moreover, in cases where this slow second stage leads to a dynamically unstable configuration (two peaks with a high mass ratio) it is followed by a new sequence ``violent relaxation/slow collisional relaxation''. We obtain an analytical estimate of the relaxation time $t_{2\to \pm 1}$ through the mean escape time of a particle from its potential well in a bistable system. We find that the diffusion and dissipation coefficients satisfy Einstein's relation and that the relaxation time scales as $N e^{1/T}$ at low temperature, in agreement with numerical simulations.

We describe a one-dimensional self-gravitating system derived from the problem of large-scale structure formation in cosmology. Considering small times so that the expansion can be neglected we present a thermodynamical analysis of this system. We find a second-order phase transition at $T_{c1}$ from an homogeneous equilibrium at high temperature to a clustered phase (with a density peak at one of the boundaries of the system) at low temperature. There also exists an infinite series of unstable equilibria which appear at lower temperatures $T_{cn}$, reflecting the scale-free nature of the gravitational interaction and the usual Jeans instability. We find that, as for the similar HMF (Hamiltonian mean field) model, all three micro-canonical, canonical and grand-canonical ensembles agree with each other, as well as with the stability properties associated with a hydrodynamical approach. On the other hand, the collisionless dynamics governed by the Vlasov equation yields the same results except that at low $T$ the equilibrium associated with two density peaks (one at each boundary) becomes stable.

We show that as many as 4000 SNeIa may be required to detect the effect of weak lensing on their flux distribution with a high level of significance. However, if the intrinsic SNeIa magnitude dispersion is unknown one needs an even higher number of SNeIa (an order of magnitude more) to reach a similar level of statistical significance. Moreover, the ability to separate the lensing contribution from the intrinsic scatter depends sensitively on the amplitude of the latter. Using a Kolmogorov - Smirnov (K-S) test we check how the required number of SNeIa changes with level of significance. Our model incorporates a completely analytical description of weak lensing which has been tested extensively against numerical simulations. Thus, future missions such as SNAP may be able to detect non-Gaussianity at a lower significance level of 10% (through the K-S test) only if the intrinsic scatter is known from external data (e.g. from low redshift observations) whereas ALPACA with 100,000 SNe will definitely detect non-Gaussianity with a very high confidence even if the intrinsic magnitude dispersion is not known a priori.

We study the possibility of using the entire probability distribution function (PDF) of the aperture mass Map and its related cumulative probability distribution function (CPDF) to obtain meaningful constraints on cosmological parameters. Deriving completely analytic expressions for the associated covariance matrices, we construct the Fisher matrix and use it to estimate the accuracy with which various cosmological parameters can be recovered from future surveys using such statistics. This formalism also includes the effect of various noises such as intrinsic ellipticity distribution of galaxies and finite survey size. The estimation errors are then compared with the ones derived from low order moments of the PDF (variance and skewness) to check how efficiently the high Map tail can be used to constrain cosmological parameters such as Omega_m, sigma_8 and dark energy equation of state w_de. We find that for future surveys such as JDEM the full PDF does not bring significant tightening of constraints on cosmology beyond what is already achievable by the joint use of second and third order moments.

Focusing on the well motivated aperture mass statistics $\Map$, we study the possibility of constraining cosmological parameters using future space based SNAP class weak lensing missions. Using completely analytical results we construct the covariance matrix for estimators based on two-point and three-point statistics. Our approach incorporates an accurate modelling of higher-order statistics to describe cosmic variance as well as various sources of discrete noise at small angular scales. These results are then fed into a Fisher matrix based analysis to study cosmological parameter degeneracies. Joint and independent analysis, with or without redshift binning, for various parameter combinations are presented. An analytical modelling of the covariance matrix opens up the possibility of testing various approximations which are often used in derivations of semi-analytical results. These include how inclusion of full non-Gaussian terms in covariance matrix affects parameter estimation. Inclusion of three-point information and how such information can enhance the accuracy with which certain parameters can be estimated is also studied in detail. It is shown that broad correlation structure among various angular scales in such circumstances implies reduction in number of available angular scales which carry completely independent information. On the other hand, the effect of theoretical inaccuracies, in modelling either the power-spectrum or bi-spectrum evolution, onto the measure of cosmological parameters from weak lensing surveys is also considered. Several cosmological parameters, $\Om$, $\sigma_8$, spectral index $n_s$, running of spectral index $\alpha_s$ and equation of state of the dark energy $\wde$ are included in the analysis.

Weak gravitational lensing surveys measure the distortion of the image of distant sources due to the deflections of light rays by the fluctuations of the gravitational potential along the line of sight. Since they probe the non-linear matter power spectrum itself at medium redshift such surveys are complimentary to both galaxy surveys (which follow stellar light) and cosmic microwave background observations (which probe the linear regime at high redshift). Ongoing CMB experiments such as WMAP and the future Planck satellite mission will measure the standard cosmological parameters with unprecedented accuracy. The focus of attention will then shift to understanding the nature of dark matter and vacuum energy: several recent studies suggest that lensing is the best method for constraining the dark energy equation of state. During the next 5 year period ongoing and future weak lensing surveys such as the Joint Dark Energy Mission (JDEM, e.g. SNAP) or the Large-aperture Synoptic Survey Telescope (LSST) will play a major role in advancing our understanding of the universe in this direction. In this review article we describe various aspects of weak lensing surveys and how they can help us in understanding our universe.

We study various measures of weak lensing distortions in future surveys, taking into account the noise arising from the finite survey size and the intrinsic ellipticity of galaxies. We also consider a realistic redshift distribution of the sources, as expected for the SNAP mission. We focus on the low order moments and the full distribution function (pdf) of the aperture-mass $\Map$ and of the smoothed shear component $\gammais$. We also propose new unbiased estimators for low-order cumulants which have less scatter than the usual estimators of non-Gaussianity based on the moments themselves. Then, using an analytical model which has already been seen to provide a good description of weak gravitational lensing through comparison against numerical simulations, we study the statistical measures which can be extracted from future surveys like the SNAP experiment. We recover the fact that at small angular scales ($1'<\theta_s<10'$) the variance can be extracted with a few percent level accuracy. Non-Gaussianity can also be measured from the skewness of the aperture-mass (at a 10% level) while the shear kurtosis is more noisy and cannot be easily measured beyond 6'. On the other hand, we find that the pdf of the estimator associated with the aperture-mass can be distinguished both from the Gaussian and the Edgeworth expansion and could provide useful constraints, while this appears to be difficult to realize with the shear component. Finally, we investigate various survey strategies and the possibility to perform a redshift binning of the sample.

The one-point probability distribution function (pdf) of the large-scale density field is an important tool to follow the evolution of cosmological structures. In this paper we present a new model for this pdf for all regimes and all densities, that is from linear to highly non-linear scales and from rare voids up to rare high densities. This is probably the simplest model one can build which is consistent with normalization constraints and known rigorous results (the quasi-linear regime and the rare void limit). It is fully parameterized by the non-linear variance and skewness. We obtain a good agreement with N-body data from realistic cosmological simulations of the VIRGO consortium and we find that it works significantly better than previous models such as the lognormal model or the Extended Perturbation Theory (EPT). We explain this success as a result of the tight constraints onto the pdf provided by these consistency conditions. We also point out that while the Lagrangian dynamics of typical fluctuations is quite complex the statistical outcome seems rather simple. This simple model should be useful for studies which require a realistic and convenient description of this pdf.

Analytical expressions for covariances of weak lensing statistics related to the aperture mass $\Map$ are derived for realistic survey geometries such as SNAP for a range of smoothing angles and redshift bins. We incorporate the contributions to the noise due to the intrinsic ellipticity distribution and the effects of finite size of the catalogue. Extending previous results to the most general case where the overlap of source populations is included in a complete analysis of error estimates, we study how various angular scales in various redshifts are correlated and how the estimation scatter changes with survey parameters. Dependence on cosmological parameters and source redshift distributions are studied in detail. Numerical simulations are used to test the validity of various ingredients to our calculations. Correlation coefficients are defined in a way that makes them practically independent of cosmology. They can provide important tools to cross-correlate one or more different surveys, as well as various redshift bins within the same survey or various angular scales from same or different surveys. Dependence of these coefficients on various models of underlying mass correlation hierarchy is also studied. Generalisations of these coefficients at the level of three-point statistics have the potential to probe the complete shape dependence of the underlying bi-spectrum of the matter distribution. A complete error analysis incorporating all sources of errors suggest encouraging results for studies using future space based weak lensing surveys such as SNAP.

We show that the formation of large-scale structures through gravitational instability in the expanding universe can be fully described through a path-integral formalism. We derive the action S[f] which gives the statistical weight associated with any phase-space distribution function f(x,p,t). This action S describes both the average over the Gaussian initial conditions and the Vlasov-Poisson dynamics. Next, applying a standard method borrowed from field theory we generalize our problem to an N-field system and we look for an expansion over powers of 1/N. We describe three such methods and we derive the corresponding equations of motion at the lowest non-trivial order for the case of gravitational clustering. This yields a set of non-linear equations for the mean $\fb$ and the two-point correlation G of the phase-space distribution f, as well as for the response function R. These systematic schemes match the usual perturbative expansion on quasi-linear scales but should also be able to handle the non-linear regime. Our approach can also be extended to non-Gaussian initial conditions and may serve as a basis for other tools borrowed from field theory.

The present generation of weak lensing surveys will be superseded by surveys run from space with much better sky coverage and high level of signal to noise ratio, such as SNAP. However, removal of any systematics or noise will remain a major cause of concern for any weak lensing survey. One of the best ways of spotting any undetected source of systematic noise is to compare surveys which probe the same part of the sky. In this paper we study various measures which are useful in cross correlating weak lensing surveys with diverse survey strategies. Using two different statistics - the shear components and the aperture mass - we construct a class of estimators which encode such cross-correlations. These techniques will also be useful in studies where the entire source population from a specific survey can be divided into various redshift bins to study cross correlations among them. We perform a detailed study of the angular size dependence and redshift dependence of these observables and of their sensitivity to the background cosmology. We find that one-point and two-point statistics provide complementary tools which allow one to constrain cosmological parameters and to obtain a simple estimate of the noise of the survey.

In this paper we describe the predictions for the smoothed weak lensing shear and aperture-mass of two simple analytical models of the density field: the minimal tree-model and the stellar model. Both models give identical results for the statistics of the 3-d density contrast smoothed over spherical cells and only differ by the detailed angular dependence of the many-body density correlations. We have shown in previous work that they also yield almost identical results for the pdf of the smoothed convergence, $\kappa_s$. We find that both models give rather close results for both the shear and the positive tail of the aperture-mass. However, we note that at small angular scales ($\theta_s \la 2'$) the tail of the pdf $\cP(\Map)$ for negative $\Map$ shows a strong variation between the two models and the stellar model actually breaks down for $\theta_s \la 0.4'$ and $\Map<0$. This shows that the statistics of the aperture-mass provides a very precise probe of the detailed structure of the density field, as it is sensitive to both the amplitude and the detailed angular behaviour of the many-body correlations. On the other hand, the minimal tree-model shows good agreement with numerical simulations over all scales and redshifts of interest, while both models provide a good description of the pdf $\cP(\gamma_{is})$ of the smoothed shear components. Therefore, the shear and the aperture-mass provide robust and complimentary tools to measure the cosmological parameters as well as the detailed statistical properties of the density field.

Weak lensing convergence can be used directly to map and probe the dark mass distribution in the universe. Building on earlier studies, we recall how the statistics of the convergence field are related to the statistics of the underlying mass distribution, in particular to the many-body density correlations. We describe two model-independent approximations which provide two simple methods to compute the probability distribution function, pdf, of the convergence. We apply one of these to the case where the density field can be described by a log-normal pdf. Next, we discuss two hierarchical models for the high-order correlations which allow one to perform exact calculations and evaluate the previous approximations in such specific cases. Finally, we apply these methods to a very simple model for the evolution of the density field from linear to highly non-linear scales. Comparisons with the results obtained from numerical simulations, obtained from a number of different realizations, show excellent agreement with our theoretical predictions. We have probed various angular scales in the numerical work and considered sources at 14 different redshifts in each of two different cosmological scenarios, an open cosmology and a flat cosmology with non-zero cosmological constant. Our simulation technique employs computations of the full 3-d shear matrices along the line of sight from the source redshift to the observer and is complementary to more popular ray-tracing algorithms. Our results therefore provide a valuable cross-check for such complementary simulation techniques, as well as for our simple analytical model, from the linear to the highly non-linear regime.

Weak gravitational lensing surveys are rapidly becoming important tools to probe directly the mass density fluctuations in the universe and its background dynamics. Earlier studies have shown that it is possible to model the statistics of the convergence field on small angular scales by suitably modeling the statistics of the underlying density field in the highly non-linear regime. We extend such methods to model the complete probability distribution function, PDF, of the shear as a function of smoothing angle. Our model relies on a simple hierarchical Ansatz for the behaviour of the higher-order correlations in the density field. We compare our predictions with the results of numerical simulations and find excellent agreement for different cosmological scenarios. Our method provides a new way to study the evolution of non-Gaussianity in gravitational clustering and should help to break the degeneracies in parameter estimation based on analysis of the power spectrum alone.

We point out that two problems of observational cosmology, the facts i) that > 60% of the baryonic content of the universe is not observed at z=0 and ii) that the properties of small clusters do not agree with simple expectations, could be closely related. As shown by recent studies, the shock-heating associated with the formation of large-scale structures heats the intergalactic medium (IGM) and leads to a ``warm IGM'' component for the gas. In the same spirit, we suggest the intracluster medium (ICM) to be a mixture of galaxy-recycled, metal enriched gas and intergalactic gas, shock-heated by the collapsing much larger scales. This could be obtained through two processes: 1) the late infalling gas from the external warm IGM is efficiently mixed within the halo and brings some additional entropy, or 2) the shocks generated by larger non-linear scales are also present within clusters and can heat the ICM. We show that if assumption (1) holds, the entropy brought by the warm IGM is sufficient to explain the observed properties of clusters, in particular the entropy floor and the LX-T relation. On the other hand, we briefly note that the scenario (2) would require a stronger shock-heating because of the larger density of the ICM as compared with filaments. Our scenario of clusters being "born warm" can be checked through the predicted redshift evolution of the entropy floor.

We investigate the error implied by the use of the Zel'dovich approximation to set up the initial conditions at a finite redshift zi in numerical simulations. Using a steepest-descent method developed in a previous work we derive the probability distribution P(delta_R) of the density contrast in the quasi-linear regime. This also provides its dependence on the redshift zi at which the simulation is started. Thus, we find that the discrepancy with the exact pdf (defined by the limit zi->infinity) is negligible after the scale factor has grown by a factor a/a_i>5, for scales which were initially within the linear regime with sigma_i>0.1. This shows that the use of the Zel'dovich approximation to implement the initial conditions is sufficient for practical purposes since these are not very severe constraints.

We investigate the behaviour of cosmological baryons at low redshifts $z<5$ after reionization through analytic means. In particular, we study the density-temperature phase-diagram which describes the history of the gas. We show how the location of the matter in this $(\rho,T)$ diagram expresses the various constraints implied by usual hierarchical scenarios. This yields robust model-independent results which agree with numerical simulations. The IGM is seen to be formed via two phases: a ``cool'' photo-ionized component and a ``warm'' component governed by shock-heating. We also briefly describe how the remainder of the matter is distributed over galaxies, groups and clusters. We recover the fraction of matter and the spatial clustering computed by numerical simulations. We also check that the soft X-ray background due to the ``warm'' IGM component is consistent with observations. We find in the present universe a baryon fraction of 7% in hot gas, 24% in the warm IGM, 38% in the cool IGM, 9% within star-like objects and, as a still un-observed component, 22% of dark baryons associated with collapsed structures, with a relative uncertainty no larger than 30% on these numbers.

We investigate the properties of the standard perturbative expansions which describe the early stages of the dynamics of gravitational clustering. We show that for hierarchical scenarios with no small-scale cutoff perturbation theory always breaks down beyond a finite order $q_+$. Besides, the degree of divergence increases with the order of the perturbative terms so that renormalization procedures cannot be applied. Nevertheless, we explain that despite the divergence of these subleading terms the results of perturbation theory are correct at leading order because they can be recovered through a steepest-descent method which does not use such perturbative expansions. Finally, we investigate the simpler cases of the Zel'dovich and Burgers dynamics. In particular, we show that the standard Burgers equation exhibits similar properties. This analogy suggests that the results of the standard perturbative expansions are valid up to the order $q_+$ (i.e. until they are finite). Moreover, the first ``non-regular'' term of a large-scale expansion of the two-point correlation function should be of the form $R^{-2} \sigma^2(R)$. At higher orders the large-scale expansion should no longer be over powers of $\sigma^2$ but over a different combination of powers of 1/R. However, its calculation requires new non-perturbative methods.

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).

Using a non-perturbative method developed in a previous work (paper II), we derive the probability distribution $P(\delta_R)$ of the density contrast within spherical cells in the quasi-linear regime for some non-Gaussian initial conditions. We describe three such models. The first one is a straightforward generalization of the Gaussian scenario. It can be seen as a phenomenological description of a density field where the tails of the linear density contrast distribution would be of the form $P_L(\delta_L) \sim e^{-|\delta_L|^{-\alpha}}$, where $\alpha$ is no longer restricted to 2 (as in the Gaussian case). We derive exact results for $P(\delta_R)$ in the quasi-linear limit. The second model is a physically motivated isocurvature CDM scenario. Our approach needs to be adapted to this specific case and in order to get convenient analytical results we introduce a simple approximation (which is not related to the gravitational dynamics but to the initial conditions). Then, we find a good agreement with the available results from numerical simulations for the pdf of the linear density contrast for $\delta_{L,R} \ga 0$. We can expect a similar accuracy for the non-linear pdf $P(\delta_R)$. Finally, the third model corresponds to the small deviations from Gaussianity which arise in standard slow-roll inflation. We obtain exact results for the pdf of the density field in the quasi-linear limit, to first-order over the primordial deviations from Gaussianity.

We develop a non-perturbative method to derive the probability distribution $P(\delta_R)$ of the density contrast within spherical cells in the quasi-linear regime. Indeed, since this corresponds to a rare-event limit a steepest-descent approximation can yield asymptotically exact results. We check that this is the case for Gaussian initial density fluctuations, where we recover most of the results obtained by perturbative methods from a hydrodynamical description. Moreover, we correct an error which was introduced in previous works for the high-density tail of the pdf. This feature, which appears for power-spectra with a slope $n<0$, points out the limitations of perturbative approaches which cannot describe the pdf $P(\delta_R)$ for $\delta_R \ga 3$ even in the limit $\sigma \to 0$. This break-up does not involve shell-crossing and it is naturally explained within our framework. Thus, our approach provides a rigorous treatment of the quasi-linear regime, which does not rely on the hydrodynamical approximation for the equations of motion. Besides, it is actually simpler and more intuitive than previous methods. Our approach can also be applied to non-Gaussian initial conditions.

We develop a systematic method to obtain the solution of the collisionless Boltzmann equation which describes the growth of large-scale structures as a perturbative series over the initial density perturbations. We give an explicit calculation of the second-order terms which are shown to agree with the results obtained from the hydrodynamical description of the system. Then, we explain that this identity extends to all orders of perturbation theory and that the perturbative series actually diverge for hierarchical scenarios. However, since the collisionless Boltzmann equation provides the exact description of the dynamics (including the non-linear regime) these results may serve as a basis for a study of the non-linear regime. In particular, we derive a non-perturbative quadratic integral equation which explicitly relates the actual non-linear distribution function to the initial conditions (more precisely, to the linear growing mode). This allows us to write an explicit path-integral expression for the probability distribution of the exact non-linear density field.

We present an analytic model, based on previous works which tackled the reionization history of the universe, which allows us to describe the generation of the secondary CMB anisotropies due to the kinetic SZ effect. We take into account the "patchy pattern" of reionization (HII bubbles), the cross-correlations of these ionized regions, the small-scale fluctuations of the matter density field and the contribution from collapsed objects. For an open universe, we find that the power-spectrum $l(l+1)C_l/(2\pi)$ exhibits a plateau of height $10^{-13}$ in the range $10^3 < l < 10^6$. We find that for large wavenumbers $l > 10^4$ the signal is dominated by the contribution from collapsed halos while for $l < 10^4$ it is governed by the large-scale correlations of HII bubbles. This implies that one cannot discriminate reionization by stars from a quasar-driven scenario since the size of ionized regions never dominates the behaviour of the anisotropies. Moreover, the secondary CMB anisotropies arise from a broad range of redshifts (7.5<z<10 for the IGM and 0<z<7 for galactic halos). The signal expected in our model might bias the cosmological parameter estimation from CMB experiments such as Planck and could be detected by future mm-wavelength interferometers (e.g., ALMA).

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 study the distribution of baryonic and luminous matter within the framework of a hierarchical scenario. Using an analytical model for structure formation which has already been checked against observations for galaxies, Lyman-$\alpha$ clouds, clusters and reionization processes, we present its predictions for the bias of these objects. We describe its dependence on the luminosity (for galaxies or quasars) or the column density (for Lyman-$\alpha$ absorbers) of the considered objects. We also study its redshift evolution, which can exhibit an intricate behaviour. These astrophysical objects do not trace the dark matter density field, the Lyman-$\alpha$ forest clouds being undercorrelated and the bright galaxies overcorrelated, while the intermediate class of Lyman-limit systems is seen to sample the matter field quite well. We also present the distribution of baryonic matter over these various objects. We show that light does not trace baryonic mass, since bright galaxies which contain most of the stars only form a small fraction of the mass associated with virialized and cooled halos. We consider two cosmologies: a critical density universe and an open universe. In both cases, our results agree with observations and show that hierarchical scenarios provide a good model for structure formation and can describe a wide range of objects which spans at least the seven orders of magnitude in mass for which data exist. More detailed observations, in particular of the clustering evolution of galaxies, will constrain the astrophysical models involved.

We use a unified analytical formulation for the multiplicity functions of clusters and galaxies which is free from the cloud-in-cloud problem encountered in earlier approaches and well adapted to the description of the non-linear clustering features. It is especially suited to simultaneously describe rich clusters, groups and galaxies, consistently with the hierarchical picture of gravitational clustering and their evolution in time. Using a simple model for the X-ray luminosity (taking into account entropy considerations), we obtain the X-ray luminosity distribution of groups and clusters. Then, using the same formalism we derive the galaxy and quasar multiplicity functions. In particular, we show that the use of the standard Press-Schechter prescription leads to erroneous conclusions at low redshifts while our approach provides a reasonable agreement with observations in a natural fashion because it is able to distinguish galactic halos from groups or clusters. Thus, we obtain a global and consistent picture of the X-ray emissions from all structures. In particular, we show that future observations (e.g., from AXAF) could provide interesting information on galaxy evolution. Indeed, they will constrain the importance of a possible hot diffuse gaseous phase in galactic halos and they could reveal massive galaxies which are just being formed, through the X-ray emission of their cooling gas.

Density fluctuations in the matter distribution lead to distortions of the images of distant galaxies through weak gravitational lensing effects. This provides an efficient probe of the cosmological parameters and of the density field. In this article, we investigate the statistical properties of the convergence due to weak gravitational lensing by non-linear structures (i.e. we consider small angular windows $\theta \la 1'$). Previous studies have shown how to relate the second and third order moments of the convergence to those of the density contrast while models based on the Press-Schechter prescription provide an estimate of the tail of $P(\kappa)$. Here we present a method to obtain an estimate of the full p.d.f. of the convergence $P(\kappa)$. It is based on a realistic description of the density field which applies to overdense as well as underdense regions. We show that our predictions agree very well with the results of N-body simulations for the convergence. This could allow one to derive the cosmological parameters $(\Omega_m,\Omega_{\Lambda})$ as well as the full p.d.f. $P(\delta_R)$ of the density contrast itself in the non-linear regime from observations. Hence this gives a very powerfull tool to constrain scenarios of structure formation.

In this article we present an analytic calculation of the probability distribution of the magnification of distant sources due to weak gravitational lensing from non-linear scales. We use a realistic description of the non-linear density field, which has already been compared with numerical simulations. Then, we can directly express the p.d.f. $P(\mu)$ of the magnification in terms of the p.d.f. of the density contrast realized on non-linear scales (typical of galaxies) where the local slope of the initial linear power-spectrum is $n=-2$. We recover the behaviour seen by numerical simulations: $P(\mu)$ peaks at a value slightly smaller than the mean $<\mu>=1$ and it shows an extended large $\mu$ tail (as described in another article our predictions also show a good quantitative agreement with results from N-body simulations for a finite smoothing angle). Then, we study the effects of weak lensing on the derivation of the cosmological parameters from SNeIa. We show that the inaccuracy introduced by weak lensing is not negligible: $\Delta \Omega_m \ga 0.3$ for two observations at $z_s=0.5$ and $z_s=1$. However, observations can unambiguously discriminate between $\Omega_m=0.3$ and $\Omega_m=1$. Moreover, in the case of a low-density universe one can clearly distinguish an open model from a flat cosmology (besides, the error decreases as the number of observed SNeIa increases). Since distant sources are more likely to be ``demagnified'' the most probable value of the observed density parameter $\Omega_m$ is slightly smaller than its actual value. On the other hand, one may obtain some valuable information on the properties of the underlying non-linear density field from the measure of weak lensing distortions.

We present in this article an analysis of some of the properties of the density field realized in numerical simulations for power-law initial power-spectra in the case of a critical density universe. We compare our numerical results in the non-linear regime with the predictions of a specific scaling model, focusing on its much wider range of applicability, which is one of its main advantages over the standard Press-Schechter approximation. We first check that the two-point correlation functions agree with the stable-clustering ansatz. Next we show that the statistics of the counts-in-cells obey the scaling law predicted by our scaling model. Then, we turn to mass functions of overdense and underdense regions. We first consider the mass function of "just collapsed" objects defined by a density threshold $\Delta~177$. We note that the usual Press-Schechter prescription agrees reasonably well with the simulations (although there are some discrepancies) while the numerical results are also consistent with the scaling model. Then, we consider more general mass functions defined by different density thresholds which can even be negative. This is out of reach of the Press-Schechter approach while our scaling model can handle these mass functions and it shows a reasonably good agreement with numerical results. Finally, we consider objects defined by a constant radius condition. Thus, we find that the scaling model allows one to study many different classes of objects and it clarifies the links between various statistical tools.

Using a unified analytic model for quasars, galaxies, Lyman-alpha absorbers and the IGM, we obtain the redshift evolution of the temperature and the entropy of the gas and the corresponding cluster temperature - X-ray luminosity relation. We show that although quasars can easily reheat the IGM and raise its entropy up to the level required by current cluster observations the energy provided by supernovae is unlikely to be sufficient. Indeed, the efficiency factor needed for the supernova scenario is of order unity while for quasars we get a value ~ 0.008. Thus the IGM is more likely to have been reheated by quasars. Moreover, we find that if both scenarios are normalized to present observations the reheating due to quasars occurs somewhat earlier (z ~ 2) than for supernovae (z ~ 0.4) because of the sharp drop at low z of the quasar luminosity function. We also show that the Compton parameter y induced by the IGM is well below the observed upper limit in all cases. Finally, we note that such a reheating process may partly account for the decline at low redshift of the comoving star formation rate and of the quasar luminosity function. In particular, we show that the contradictory requirements arising from clusters and galaxies lead to a reheating temperature $T \sim 5 10^5$ K. On the other hand, the reionization process of the universe is almost not modified by these entropy sources.

Within the framework of hierarchical clustering scenarios, we investigate the consequences for the properties of virialized halos of the constraints provided by numerical simulations on the first few correlation functions. Thus, we show that the density field cannot be described by a collection of smooth halos with a universal density profile. This implies that substructures within larger objects play an important role (but a mean spherically averaged density profile may exist). In particular, a possible interpretation is that collapsed objects can be divided into an infinite hierarchy of smaller objects with increasingly large densities (these substructures might also be continuously destroyed and created by the long-range action of gravity). Finally, we present multifractal models (restricted to non-linear scales) which can describe in a natural way such non-linear density fields with increasingly large fluctuations at smaller scales. We relate their properties to the correlation functions and present a few constraints they are expected to satisfy, using theoretical considerations as well as constraints from numerical simulations. Thus, the simplest realistic model is the bifractal model described in Balian & Schaeffer (1989a). Moreover, we show that it should provide (at least) a very good approximation of the multifractal properties of the actual non-linear density field, hence of the probability distribution of the density contrast. The implications of this model (e.g. for galaxies) are detailed in other studies.

We incorporate quasars into an analytic model to describe the reheating and reionization of the universe. In combination with a previous study of galaxies and Lyman-$\alpha$ clouds, we are able to provide a unified description of structure formation, verified against a large variety of observations. We also take into account the clumping of the baryonic gas in addition to the presence of collapsed objects. We consider two cosmologies: a critical universe with a CDM power-spectrum and an open universe with $\Omega_0=0.3$, $\Lambda=0$. The derived quasar luminosity function agrees reasonably well with observations at $z<4.5$ and with constraints over larger redshifts from the HDF. The radiation produced by these objects at $z \sim 16$ slowly reheats the universe which gets suddenly reionized at $z_{ri}=6.8$ for the open universe ($z_{ri}=5.6$ for the critical density universe). The UV background radiation simultaneously increases sharply to reach a maximum of $J_{21} = 0.18$ at $z=2.6$, but shows strong ionization edges until $z \leq 1$. The metallicity of the gas increases quickly at high $z$ and is already larger than $0.01 Z_{\odot}$ at $z=10$. The QSO number counts and the helium opacity constrain the reionization redshift to be $z_{ri} \sim 6$. We confirm that a population of faint quasars is needed in order to satisfy the observations. Due to the low reionization redshift, the damping of CMB fluctuations is quite small, but future observations (e.g. with the NGST) of the multiplicity functions of radiation sources and of the HI and HeII opacities will strongly constrain scenarios in which reionization is due to QSOs.

We present a model for the Lyman-alpha absorbers that treats all objects (from the low-density forest clouds to the dense damped systems) in a unified description. This approach is consistent with an earlier model of galaxies (luminosity function, metallicity) but also with the known description of the density field in the small-scale non-linear regime. We consider two cosmological models: a critical universe $\Omega=1$ with a CDM power-spectrum, and an open CDM universe with $\Omega_0=0.3$, $\Lambda=0$. We reproduce the available data on column density distribution as a function of redshift, the value of the main new parameter, the background ionizing UV flux, being consistent with the observed limits. This allows a quantitatively trustable analytical description of the opacity, mass, size, velocity dispersion and metallicity of these absorbers, over a range of column densities spanning 10 orders of magnitude. Moreover, together with an earlier model of galaxy formation this draws a unified picture of the redshift evolution of structures in the universe, from underdense clouds to massive high density galaxies, from weak to very deep potential wells.

We set up a model for the evolution of the galaxy luminosity function, taking advantage of recent work that brought in some better understanding of the mass function for gravitationally condensed objects. The physics of cooling of the baryonic component allows us to distinguish halos that become groups or clusters from those that eventually form galaxies (possibly within the former objects). With our new mass function and our new application of the cooling criteria - which motivated this paper - we get a satisfactory and natural cutoff at the bright end of the luminosity function, the needed flat slope for faint magnitudes and the correct trend in colors (brighter galaxies are redder ) within the framework of the hierarchical clustering picture. This infirms earlier claims that the latter was inadequate to reproduce the former observations. We find the velocity dispersion to be a much better parameter than mass or radius to characterize galaxies. We find that bright galaxies form at z~2 from mergers with a rather quiet evolution afterwards, whereas small galaxies are the result of a continuous merging process active up to the present epoch. The transition is found to occur at the observed transition between bright spirals and small dwarf ellipticals or irregulars. The quasar multiplicity as a function of redshift is also discussed.

Within the framework of hierarchical clustering we show that a simple Press-Schechter-like approximation, based on spherical dynamics, provides a good estimate of the evolution of the density field in the quasi-linear regime up to $\Sigma \sim 1$. Moreover, it allows one to recover the exact series of the cumulants of the probability distribution of the density contrast in the limit $\Sigma \to 0$ which sheds some light on the rigorous result and on ``filtering''. We also obtain similar results for the divergence of the velocity field. Next, we extend this prescription to the highly non-linear regime, using a stable-clustering approximation. Then we recover a specific scaling of the counts-in-cells which is indeed seen in numerical simulations, over a well-defined range. To this order we also introduce an explicit treatment of the behaviour of underdensities, which takes care of the normalization and is linked to the low-density bubbles and the walls one can see in numerical simulations. We compare this to a 1-dimensional adhesion model, and we present the consequences of our prescription for the power-law tail and the cutoff of the density distribution.

We calculate the multiplicity function of matter condensations by directly considering the actual, deeply non-linear density field, which we compare to the popular Press-Schechter approximation. We show the mass function is a function of a unique parameter that contains all the dependence on the mass and radius of the object. It enables us to derive the mass function of general astrophysical objects (in addition to just-virialized halos) which may be defined by any density threshold, that may even vary with the mass of the object. This is beyond the reach of usual formulations based on the initial gaussian field and gives a clear illustration of the advantages of our approach. We explain why numerical tests seem to favor both Press-Schechter and non-linear prescriptions even though the two approximations differ by their scaling as a function of mass as well as a function of redshift. The difference between the two is due to the fact that the Press-Schechter prescription assumes that present-day mass fluctuations can be recognized in the early linear universe, and that their number is conserved, while the non-linear approach takes into account their evolution which leads to an increase of the number of highly non-linear objects (very large or small masses). This difference is seen to be of the same magnitude as the difference obtained by varying the initial spectrum of density fluctuations. Earlier conclusions drawn about the relevance of the latter using analytical approximations to the mass function must thus be reexamined.