We investigate how long-range scalar interactions affect the properties of dark matter haloes. For doing so we employ the ReBEL model which implements an additional interaction between dark matter particles. On the phenomenological level this is equivalent to a modification of gravity. We analyse the differences between five ReBEL models and $\Lambda$CDM using a series of high resolution cosmological simulations. Emphasis is placed on investigating how halo properties change in the presence of a fifth force. We report that the density profile of ReBEL haloes is well described by the NFW profile but with mean concentrations from $5\%$ to a few times higher than the standard $\Lambda$CDM value. We also find a slight increase of the halo spin for haloes more massive than $5\times10^{11}\M_{\odot}$, reflecting a higher rotational support of those haloes due to scalar forces. In addition, the dark matter haloes in our models are more spherical than their counterparts in $\Lambda$CDM. The ReBEL haloes are also more virialised, with a large difference from $\Lambda$CDM for strong fifth forces and a much smaller change for weak scalar interactions.

In this contribution we present the preliminary results regarding the non-linear BAO signal in higher-order statistics of the cosmic density field. We use ensembles of N-body simulations to show that the non-linear evolution changes the amplitudes of the BAO signal, but has a negligible effect on the scale of the BAO feature. The latter observation accompanied by the fact that the BAO feature amplitude roughly doubles as one moves to higher orders, suggests that the higher-order correlation amplitudes can be used as probe of the BAO signal.

In this study we show that the skewness S_3 of the cosmic density field contains a significant and potentially detectable and clean imprint of Baryonic Acoustic Oscillations. Although the BAO signal in the skewness has a lower amplitude than second order measures like the two-point correlation function and power spectrum, it has the advantage of a considerably lower sensitivity to systematic influences. Because it lacks a direct dependence on bias if this concerns simple linear bias, skewness will be considerably less beset by uncertainties due to galaxy bias. Also, it has a weaker sensitivity to redshift distortion effects. We use perturbation theory to evaluate the magnitude of the effect on the volume-average skewness, for various cosmological models. One important finding of our analysis is that the skewness BAO signal occurs at smaller scales than that in second order statistics. For an LCDM spectrum with WMAP7 normalization, the BAO feature has a maximum wiggle amplitude of ~3% and appears at a scale of ~82Mpc/h. We conclude that the detection of BAO wiggles in future extensive galaxy surveys via the skewness of the observed galaxy distribution may provide us with a useful, and potentially advantageous, measure of the nature of Dark Energy.

In this work we investigate higher order statistics for the $\lcdm$ and ReBEL scalar-interacting dark matter models by analyzing $180\hmpc$ dark matter N-body simulation ensembles. The N-point correlation functions and the related hierarchical amplitudes, such as skewness and kurtosis, are computed using the Count-In-Cells method. Our studies demonstrate that the hierarchical amplitudes $S_n$ of the scalar-interacting dark matter model significantly deviate from the values in the $\lcdm$ cosmology on scales comparable and smaller then the screening length $r_s$ of a given scalar-interacting model. The corresponding additional forces that enhance the total attractive force exerted on dark matter particles at galaxy scales lowers the values of the hierarchical amplitudes $S_n$. We conclude that hypothetical additional exotic interactions in the dark matter sector should leave detectable markers in the higher-order correlation statistics of the density field. We focussed in detail on the redshift evolution of the dark matter field's skewness and kurtosis. From this investigation we find that the deviations from the canonical $\lcdm$ model introduced by the presence of the ``fifth'' force attain a maximum value at redshifts $0.5<z<2$. We therefore conclude that moderate redshift data are better suited for setting observational constraints on the investigated ReBEL models.

The amplitude of cosmological density fluctuations, sigma_8, has been studied and estimated by analysing many cosmological observations. The values of the estimates vary considerably between the various probes. However, different estimators probe the value of sigma_8 in different cosmological scales and do not take into account the nonlinear evolution of the parameter at late times. We show that estimates of the amplitude of cosmological density fluctuations derived from cosmic flows are systematically higher than those inferred at early epochs from the CMB because of nonlinear evolution at later times. We discuss the past and future evolution of linear and nonlinear perturbations, derive corrections to the value of sigma_8 and compare amplitudes after accounting for these differences.

We explore the possibility of improving the $\Lambda$CDM model at megaparsec scales by introducing a scalar interaction which increases the mutual gravitational attraction of dark matter particles. Using N-body simulations, we study the spatial distribution of dark matter particles and halos. We measure the effect of modifications in the Newton's gravity on properties of the two-point correlation function, the dark matter power spectrum, the cumulative halo mass function and density probability distribution functions. The results look promising: the scalar interactions produce desirable features at megaparsec scales without spoiling the $\Lambda$CDM successes at larger scales.

We propose to use spatial correlations of the kinetic Sunyaev-Zeldovich (KSZ) flux as an estimator of the peculiar velocity power spectrum. In contrast with conventional techniques, our new method does not require measurements of the thermal SZ signal or the X-ray temperature. Moreover, this method has the special advantage that the expected systematic errors are always sub-dominant to statistical errors on all scales and redshifts of interest. We show that future large sky coverage KSZ surveys may allow a peculiar velocity power spectrum estimates of an accuracy reaching ~10%.

We show that by directly correlating the cluster kinetic Sunyaev Zeldovich (KSZ) flux, the cluster peculiar velocity power spectrum can be measured to $\sim 10%$ accuracy by future large sky coverage KSZ surveys. This method is almost free of systemics entangled in the usual velocity inversion method. The direct correlation brings extra information of density and velocity clustering. We utilize these information to construct two indicators of the Hubble constant and comoving angular distance and propose a novel method to constrain cosmology.

Using mean relative peculiar velocity measurements for pairs of galaxies, we estimate the cosmological density parameter $\Omega_m$ and the amplitude of density fluctuations $\sigma_8$. Our results suggest that our statistic is a robust and reproducible measure of the mean pairwise velocity and thereby the $\Omega_m$ parameter. We get $\Omega_m = 0.30^{+0.17}_{-0.07}$ and $\sigma_8 = 1.13^{+0.22}_{-0.23}$. These estimates do not depend on prior assumptions on the adiabaticity of the initial density fluctuations, the ionization history, or the values of other cosmological parameters.

We present a new constraint on the biased galaxy formation picture. Gravitational instability theory predicts that the two-point mass density correlation function, \xi(r), has an inflection point at the separation r=r_0, corresponding to the boundary between the linear and nonlinear regime of clustering, \xi = 1. We show how this feature can be used to constrain the square of the biasing parameter, b^2 = \xi_g / \xi on scales r = r_0, where \xi_g is the galaxy-galaxy correlation function, allowed to differ from \xi. We apply our method to real data: the \xi_g(r), estimated from the APM galaxy survey. Our results suggest that the APM galaxies trace the mass at separations r > 5 Mpc/h, where h is the Hubble constant in units of 100 km/s Mpc. The present results agree with earlier studies, based on comparing higher order correlations in the APM with weakly non-linear perturbation theory. Both approaches constrain the "b" factor to be within 20% of unity. If the existence of the feature we identified in the APM \xi_g(r) -- the inflection point near \xi_g = 1 -- is confirmed by more accurate surveys, we may have discovered gravity's smoking gun: the long awaited ``shoulder'' in \xi, predicted by Gott and Rees 25 years ago.

We present two new dynamical tests of the biasing hypothesis. The first is based on the amplitude and the shape of the galaxy-galaxy correlation function, $\xi_g(r)$, where $r$ is the separation of the galaxy pair. The second test uses the mean relative peculiar velocity for galaxy pairs, $\vs(r)$. This quantity is a measure of the rate of growth of clustering and it is related to the two-point correlation function for the matter density fluctuations, $\xi(r)$. Under the assumption that galaxies trace the mass ($\xi_g = \xi$), the expected relative velocity can be calculated directly from the observed galaxy clustering. The above assumption can be tested by confronting the expected $\vs$ with direct measurements from velocity-distance surveys. Both our methods are checked against N-body experiments and then compared with the $\xi_g(r)$ and $\vs$ estimated from the {\sc APM} galaxy survey and the Mark III catalogue, respectively. Our results suggest that cosmological density parameter is low, $\Omega_m \approx 0.3$, and that the {\sc APM} galaxies trace the mass at separations $r \ga 5 \Mlu$, where $h$ is the Hubble constant in units of 100 km s$^{-1}$Mpc. The present results agree with earlier studies, based on comparing higher order correlations in the {\sc APM} with weakly non-linear perturbation theory. Both approaches constrain the linear bias factor to be within 20% of unity. If the existence of the feature we identified in the {\sc APM} $\xi_g(r)$ -- the inflection point near $\xi_g = 1$ -- is confirmed by more accurate surveys, we may have discovered gravity's smoking gun: the long awaited ``shoulder'' in $\xi$, generated by gravitational dynamics and predicted by Gott and Rees 25 years ago.

We compute the skewness of the matter distribution arising from non-linear evolution and from non-Gaussian initial perturbations. We apply our result to a very generic class of models with non-Gaussian initial conditions and we estimate analytically the ratio between the skewness due to non-linear clustering and the part due to the intrinsic non-Gaussianity of the models. We finally extend our estimates to higher moments.

The motions of galaxies can be used to constrain the cosmological density parameter Omega and the clustering amplitude of matter on large scales. The mean relative velocity of galaxy pairs, estimated from the Mark III survey, indicates that Omega = 0.35 +0.35/-0.25. If the clustering of galaxies is unbiased on large scales, Omega = 0.35 +/- 0.15, so that an unbiased Einstein-de Sitter model (Omega = 1) is inconsistent with the data.

It is well known that estimating the pairwise velocity of galaxies, v_{12}, from the redshift space galaxy correlation function is difficult because this method is highly sensitive to the assumed model of the pairwise velocity dispersion. Here we propose an alternative method to estimate v_{12} directly from peculiar velocity samples, which contain redshift-independent distances as well as galaxy redshifts. In contrast to other dynamical measures which determine beta = sigma_8 x Omega^{0.6}, our method can provide an estimate of (sigma_8)^2 x Omega^{0.6} for a range of sigma_8 (here Omega is the cosmological mass density parameter while sigma_8 is the standard normalization parameter for the spectrum of matter density fluctuations). We demonstrate how to measure this quantity from realistic catalogues.

We propose a simple and accurate method for computing analytically the mass correlation function for cold dark matter and scale-free models that fits N-body simulations over a range that extends from the linear to the strongly non-linear regime. The method, based on the dynamical evolution of the pair conservation equation, relies on a universal relation between the pair-wise velocity and the smoothed correlation function valid for high and low density models, as derived empirically from N-body simulations. An intriguing alternative relation, based on the stable-clustering hypothesis, predicts a power-law behavior of the mass correlation function that disagrees with N-body simulations but conforms well to the observed galaxy correlation function if negligible bias is assumed. The method is a useful tool for rapidly exploring a wide span of models and, at the same time, raises new questions about large scale structure formation.

The mean pairwise velocity of galaxies has traditionally been estimated from the redshift space galaxy correlation function. This method is notorious for being highly sensitive to the assumed model of the pairwise velocity dispersion. Here we propose an alternative method to estimate the streaming velocity directly from peculiar velocity samples, which contain redshift-independent distances as well as galaxy redshifts. This method can provide an estimate of $\Omega^{0.6}\sigma_8^2$ for a range of $\sigma_8$ where $\Omega$ is the cosmological density parameter, while $\sigma_8$ is the standard normalization for the power spectrum of density fluctuations. We demonstrate how to measure this quantity from realistic catalogues and identify the main sources of bias and errors

We derive a simple closed-form expression, relating $\vs(r)$ -- the mean relative velocity of pairs of galaxies at fixed separation $r$ -- to the two-point correlation function of mass density fluctuations, $\xi(r)$. We compare our analytic model for $\vs(r)$ with N-body simulations, and find excellent agreement in the entire dynamical range probed by the simulations ($0.1 \lsim \xi \lsim 1000$). Our results can be used to estimate the cosmological density parameter, $\Om$, directly from redshift-distance surveys, like Mark III.

Using an Eulerian perturbative calculation, we show that the distribution of relative pairwise velocities which arises from gravitational instability of Gaussian density fluctuations has asymmetric (skewed) exponential tails. The negative skewness is induced by the negative mean streaming velocity of pairs (the infall prevails over expansion), while the exponential tails arise because the relative pairwise velocity is a number, not volume weighted statistic. The derived probability distribution is compared with N-body simulations and shown to provide a reasonable fit.

If the first stars formed soon after decoupling of baryons from the thermal cosmic background radiation (CBR), the radiation may have been last scattered in a cloudy plasma. We discuss the resulting small-scale anisotropy of the CBR in the limit where the plasma clouds are small compared to the mean distance between clouds along a line of sight. This complements the perturbative analysis valid for mildly nonlinear departures from homogeneity at last scattering. We conclude that reasonable choices for the cloud parameters imply CBR anisotropy consistent with the present experimental limits, in agreement with the perturbative approach. This means the remarkable isotropy of the CBR need not contradict the early small-scale structure formation predicted in some cosmogonies.

Since the appearance of the classical paper of Lifshitz almost half a century ago, linear stability analysis of cosmological models is textbook knowledge. Until recently, however, little was known about the behavior of higher than linear order terms in the perturbative expansion. These terms become important in the weakly nonlinear regime of gravitational clustering, when the rms mass density contrast is only slightly smaller than unity. In the past, theorists showed little interest in studying this regime, and for a good reason: only a decade ago, it would have been an academic excercise - at scales large enough to probe the weakly nonlinear regime, all measures of clustering were dominated by noise. This is no longer the case with present data. The purpose of this talk is to provide a brief summary of recent advances in weakly nonlinear perturbation theory. We present analytical perturbative results together with results of N-body experiments, conducted to test their accuracy. We compare perturbative predictions with measurements from galaxy surveys. Such comparisons can be used to test the gravitational instability theory and to constrain possible deviations from Gaussian statistics in the initial mass distribution; they can be also used to study the nature of physical processes that govern galaxy formation (``biasing''). We also show how future studies of velocity field statistics can provide a new way to determine the density parameter, $\Omega$.

We present a series of N-body experiments which confirm the reality of the previrialization effect. We also use weakly nonlinear perturbative approach to study the phenomenon. These two approaches agree when the rms density contrast, $\sigma$, is small; more surprisingly, they remain in agreement when $\sigma \approx 1$. When the slope of the initial power spectrum is $n>-1$, nonlinear tidal interactions slow down the growth of density fluctuations and the magnitude of the suppression increases when $n$ (i.e. the relative amount of small scale power) is increased. For $n<-1$ we see an opposite effect: the fluctuations grow more rapidly than in linear theory. The transition occurs at $n=-1$ when the weakly nonlinear correction to $\sigma$ is close to zero and the growth rate is close to linear. Our results resolve recent controversy between two N-body studies of previrialization. Peebles (1990) assumed $n=0$ and found strong evidence in support of previrialization, while Evrard \& Crone (1992), who assumed $n=-1$, reached opposite conclusions. As we show here, the initial conditions with $n=-1$ are rather special because the nonlinear effects nearly cancel out for that particular spectrum. In addition to our calculations for scale-free initial spectra, we show results for a more realistic spectrum of Peacock \& Dodds (1994). Its slope near the scale usually adopted for normalization is close to $-1$, so $\sigma$ is close to linear. Our results retroactively justify linear normalization at 8$h^{-1}$ Mpc, while also demonstrating the danger and limitations of this practice.

An attractive and simple hypothesis for the formation of large-scale structure is that it developed by gravitational instability from primordial fluctuations with an initially Gaussian probability distribution. Non-linear gravitational evolution drives the distribution away from the Gaussian form, generating measurable skewness and kurtosis even when the variance of the fluctuations is much smaller than unity. We use perturbation theory to compute the kurtosis of the mass density field and the velocity divergence field that arises during the weakly non-linear evolution of initially Gaussian fluctuations. We adopt an Einstein--de~Sitter universe for the perturbative calculations, and we discuss the generalization to a universe of arbitrary $\Omega$. We obtain semi-analytic results for the case of scale-free, power-law spectra of the initial fluctuations and final smoothing of cosmic fields with a Gaussian filter. We also give an exact analytical formula for the dependence of the skewness of these fields on the power spectrum index. We show that the kurtosis decreases with the power spectrum index, and we compare our more accurate results for the kurtosis to previous estimates from Monte Carlo integrations. We also compare our results to values obtained from cosmological N-body simulations with power-law initial spectra. Measurements of the skewness and kurtosis parameters can be used to test the hypothesis that structure in the universe formed by gravitational instability from Gaussian initial conditions.

We study the effects of peculiar velocities on statistical measures of galaxy clustering. These effects occur when distances to the galaxies are estimated from their redshifts. It is assumed that the clustering pattern results from the gravitational instability of initially Gaussian, small-amplitude perturbations of a Friedman-Lemaitre cosmological model. Explicit expressions are given for an arbitrary density parameter \Omega of the model, both when the cosmological constant, \Lambda, is zero, and when the model is spatially flat, \Omega + \Lambda/3H^2 =1. We show how to compute the skewness, or third moment of the density field, in redshift space, and its Fourier space counterpart, the bispectrum. We rely on a perturbative expansion of particle trajectories in Lagrangian coordinates. This formalism extends to higher orders the Zel'dovich first order solution (1970). We show that a physically consistent and quantitatively accurate analysis of the growth skewness in redshift space can be obtained from second-order Lagrangian theory. We also study the effects of spatial smoothing of the evolved density field. We give analytic expressions for the gravitationally induced skewness as a function of the power spectrum and of \Omega, for a spherical top-hat and a Gaussian smoothing filter. We compare our analytical predictions with measurements performed in numerical simulations, and find good agreement. These results should then prove useful in analyzing large scale redshift surveys. In particular, our results, in conjunction with the recent suggestion of Fry (1994), may help to break the coupling between the bias and the cosmological

This paper deals with the time evolution in the matter era of perturbations in Friedman-Lemaitre models with arbitrary density parameter $\Omega$, with either a zero cosmological constant, $\Lambda = 0$, or with a non-zero cosmological constant in a spatially flat Universe. Unlike the classical Eulerian approach where the density contrast is expanded in a perturbative series, this analysis relies instead on a perturbative expansion of particles trajectories in Lagrangian coordinates. This brings a number of advantages over the classical analysis. In particular, it enables the description of stronger density contrasts. Indeed the linear term is the famous Zel'dovich approximate solution (1970). We present here a systematic and detailed account of this approach. We give analytical results (or fits to numerical results) up to the third order. We then proceed to explore the link between the lagrangian description and statistical measures. We show in particular that Lagrangian perturbation theory provides a natural framework to compute the effect of redshift distortions, using the skewness of the density distribution function as an example. Finally, we show how well the second order theory does as compared to other approximat- ions in the case of spherically symmetric perturbations. We also compare this second order approximation and Zel'dovich solution to N-body simulations in the description of large-scale structure formation starting from a power law (n=-2) power spectrum of Gaussian perturbation. We find that second order theory is both simple and powerful.

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

This paper covers the material of our two talks. We describe a series of projects based upon perturbative expansions to follow the gravitational evolution of the one point probability distribution functions (PDFs) for the density contrast field and for the divergence of the corresponding velocity field. The Edgeworth expansion greatly simplifies the problem. Indeed, the PDF can be described through several low order moments, with the first non-trivial contribution coming from the skewness, or third moment. Gaussian initial conditions imply that the skewness of the smoothed density field is proportional to the square of the variance. The constant of proportionality is then a function of the logarithmic slope of the variance versus scale, and is otherwise independent of scale. This constant is also practically independent of the density parameter $\Omega$, and has very nearly the same value when distances are estimated by redshifts. To show the latter property, we also briefly discuss our Lagrangian approach to perturbation theory. Finally, the observed skewness does depend, of course, on the relation between the distribution of light and mass; we describe the case of local biasing. For the divergence of the velocity field, the skewness of its PDF is also proportional to the square of its variance, but the ratio now strongly depends on the value of $\Omega$. This offers a new method to determine the density parameter which does not involve comparisons

We calculate the cosmological evolution of the 1-point probability distribution function (PDF), using an analytic approximation that combines gravitational perturbation theory with the Edgeworth expansion of the PDF. Our method applies directly to a smoothed mass density field or to the divergence of a smoothed peculiar velocity field, provided that rms fluctuations are small compared to unity on the smoothing scale, and that the primordial fluctuations that seed the growth of structure are Gaussian. We use this `Edgeworth approximation' to compute the evolution of $<\delta|\delta|>$ and $<|delta|>$; these measures are similar to the skewness and kurtosis of the density field, but they are less sensitive to tails of the probability distribution, so they may be more accurately estimated from surveys of limited volume. We compare our analytic calculations to cosmological N-body simulations in order to assess their range of validity. When $\sigma \ll 1$, the numerical simulations and perturbation theory agree precisely, demonstrating that the N-body method can yield accurate results in the regime of weakly non-linear clustering. We show analytically that `biased' galaxy formation preserves the relation $<\delta^3> \propto <\delta^2>^2$ predicted by second-order perturbation theory, provided that the galaxy density is a local function of the underlying mass density.

Large-scale structures, observed today, are generally believed to have grown from random, small-amplitude inhomogeneities, present in the early Universe. We investigate how gravitational instability drives the distribution of these fluctuations away from the initial state, assumed to be Gaussian. Using second order perturbation theory, we calculate the skewness factor, $S_3 \equiv \langle \delta^3 \rangle \, /\, \langle \delta^2 \rangle^2$. Here the brackets, $\langle \ldots \rangle$, denote an ensemble average, and $\delta$ is the density contrast field, smoothed with a low pass spatial filter. We show that $S_3$ decreases with the slope of the fluctuation power spectrum; it depends only weakly on $\Omega$, the cosmological density parameter. We compare perturbative calculations with N-body experiments and find excellent agreement over a wide dynamic range. If galaxies trace the mass, measurements of $S_3$ can be used to distinguish models with Gaussian initial conditions from their non-Gaussian alternatives.