We study correlations amongst tidal fields originated by the large-scale distribution of matter in the Universe. The two-point tidal correlation is described as a rank-4 tensor, whose elements can be written in terms of four fundamental scalar functions ranging, with respect the spatial separation, from purely transversal to purely longitudinal correlations. Tidal fields, both on galaxy and cluster scales, reveal to be correlated over distances larger than the mass-density correlation lenght, though traceless tidal fields show anti-correlation between diagonal terms along orthogonal directions. The cross-correlation between mass and tidal field is also analyzed. These results are relevant for galaxy formation and the interpretation of large-scale weak lensing phenomena.

We show with analytic models that the assumption of uncorrelated intrinsic ellipticities of target sources that is usually made in searches for weak gravitational lensing due to large-scale mass inhomogeneities (``field lensing'') is unwarranted. If the orientation of the galaxy image is determined either by the angular momentum or the shape of the halo in which it forms, then the image should be aligned preferentially with the component of the tidal gravitational field perpendicular to the line of sight. Long-range correlations in the tidal field will thus lead to long-range ellipticity-ellipticity correlations that mimic the shear correlations due to weak gravitational lensing. We calculate the ellipticity-ellipticity correlation expected if halo shapes determine the observed galaxy shape, and we discuss uncertainties (which are still considerable) in the predicted amplitude of this correlation. The ellipticity-ellipticity correlation induced by angular momenta should be smaller. We consider several methods for discriminating between the weak-lensing (extrinsic) and intrinsic correlations, including the use of redshift information. An ellipticity--tidal-field correlation also implies the existence of an alignment of images of galaxies near clusters. Although the intrinsic alignment may complicate the interpretation of field-lensing results, it is inherently interesting as it may shed light on galaxy formation as well as on structure formation.

We calculate the galaxy bispectrum in both real and redshift space adopting the most common prescriptions for local Eulerian biasing and Lagrangian evolving-bias model. We show that the two biasing schemes make measurably different predictions for these clustering statistics. The Eulerian prescription implies that the galaxy distribution depends only on the present-day local mass distribution, while its Lagrangian counterpart relates the current galaxy distribution to the mass distribution at an earlier epoch when galaxies first formed. Detailed measurement of the galaxy bispectrum (of its reduced amplitude) can help establish whether galaxy positions are determined by the current mass distribution or an earlier mass distribution.

We investigate the clustering of dark matter halos in Lagrangian space in terms of their two-point correlation function, spanning more than 4 orders of magnitudes of halo masses. Analyzing a set of collisionless scale-free 128^3-particle N-body simulations with spectral indices n = -2,-1, we measure the first two Lagrangian bias parameters b_1 and b_2 relating halo and mass correlations. We find that the Mo & White leading-order formula for b_1 describes the clustering of halos with mass M ~ M_* quite accurately, where M_* indicates the characteristic non-linear mass. Smaller halos turn out to be less clustered in Lagrangian space than predicted by Mo & White. Our findings are consistent with the recent results of Jing for the clustering of halo populations in Eulerian space, demonstrating that the discrepancies between the N-body and analytical Mo & White prediction for the bias exist already in Lagrangian space. This shows that a more refined theoretical algorithm for selecting halos in the initial conditions needs to be developed. Finally, we present a very accurate fitting formula for the linear halo bias factor b_1 in Lagrangian space.

We study the spatial distribution of dark matter halos in the Universe in terms of their number density contrast, related to the underlying dark matter fluctuation via a non-local and non-linear bias random field. The description of the matter dynamics is simplified by adopting the `truncated' Zel'dovich approximation to obtain both analytical results and simulated maps. The halo number density field in our maps and its probability distribution reproduce with excellent accuracy those of halos in a high-resolution N-body simulation with the same initial conditions. Our non-linear and non-local bias prescription matches the N-body halo distribution better than any Eulerian linear and local bias.

This paper presents a stochastic approach to the clustering evolution of dark matter haloes in the Universe. Haloes, identified by a Press-Schechter-type algorithm in Lagrangian space, are described in terms of `counting fields', acting as non-linear operators on the underlying Gaussian density fluctuations. By ensemble averaging these counting fields, the standard Press-Schechter mass function as well as analytic expressions for the halo correlation function and corresponding bias factors of linear theory are obtained, thereby extending the recent results by Mo and White. The non-linear evolution of our halo population is then followed by solving the continuity equation, under the sole hypothesis that haloes move by the action of gravity. This leads to an exact and general formula for the bias field of dark matter haloes, defined as the local ratio between their number density contrast and the mass density fluctuation. Besides being a function of position and `observation' redshift, this random field depends upon the mass and formation epoch of the objects and is both non-linear and non-local. The latter features are expected to leave a detectable imprint on the spatial clustering of galaxies, as described, for instance, by statistics like bispectrum and skewness. Our algorithm may have several interesting applications, among which the possibility of generating mock halo catalogues from low-resolution N-body simulations.

We present a stochastic approach to the spatial clustering of dark matter haloes in Lagrangian space. Our formalism is based on a local formulation of the `excursion set' approach by Bond et al., which automatically accounts for the `cloud-in-cloud' problem in the identification of bound systems. Our method allows to calculate correlation functions of haloes in Lagrangian space using either a multi-dimensional Fokker-Planck equation with suitable boundary conditions or an array of Langevin equations with spatially correlated random forces. We compare the results of our method with theoretical predictions for the halo auto-correlation function considered in the literature and find good agreement with the results recently obtained within a treatment of halo clustering in terms of `counting fields' by Catelan et al.. The possible effect of spatial correlations on numerical simulations of halo merger trees is finally discussed.

The formalism that describes the non-linear growth of the angular momentum L of protostructures from tidal torques in a Friedmann Universe, as developed in a previous paper, is extended to include non-Gaussian initial conditions. We restrict our analysis here to a particular class of non-Gaussian primordial distributions, namely multiplicative models. In such models, strongly correlated phases are produced by obtaining the gravitational potential via a nonlinear local transformation of an underlying Gaussian random field. The dynamical evolution of the system is followed by describing the trajectories of fluid particles using second-order Lagrangian perturbation theory. In the Einstein-de Sitter universe, the lowest-order perturbative correction to the variance of the linear angular momentum of collapsing structures grows as t^8/3 for generic non-Gaussian statistics, which contrasts with the t^10/3 growth rate characteristic of Gaussian statistics. This is a consequence of the fact that the lowest-order perturbative spin contribution in the non-Gaussian case arises from the third moment of the gravitational potential, which is identically zero for a Gaussian field. Evaluating these corrections at the maximum expansion time of the collapsing structure, we find that these non-Gaussian and non-linear terms can be as high as the linear estimate, without the degree of non-Gaussianity as quantified by skewness and kurtosis of the density field being unacceptably large. The results suggest that higher-order terms in the perturbative expansion may contribute significantly to galactic spin which contrasts with the straightforward Gaussian case.

We discuss the non-linear evolution of the angular momentum L acquired by protostructures, like protogalaxies and protoclusters, due to tidal interactions with the surrounding matter inhomogeneities. The primordial density distribution is assumed to be Gaussian and the non-linear dynamics of the collisionless mass fluid is followed using Lagrangian perturbation theory. For a Cold Dark Matter spectrum, the inclusion of the leading-order Lagrangian correction terms results in a value of the rms ensemble average <L^2>^{1/2} which is only a factor of 1.3 higher than the corresponding linear estimate, irrespective of the scale. Consequently, the predictions of linear theory are rather accurate in quantifying the evolution of the angular momentum of protostructures before collapse sets in. In the Einstein-de Sitter universe, the initial torque is a good estimate for the tidal torque over the whole period during which the object is spun up.

The growth of the angular momentum L of protogalaxies induced by tidal torques is reconsidered within the Zel'dovich approximation. We obtain a general expression for the ensemble expectation value of the square of L in terms of the first and second invariant of the inertia tensor of the Lagrangian volume enclosing the protoobject's collapsing mass. We then specialize the formalism to the particular case in which this volume is centered on a peak of the smoothed Gaussian density field and approximated by an isodensity ellipsoid. The result is the appropriate analytical estimate for the rms angular momentum of peaks to be compared against simulations that make use of the Hoffman-Ribak algorithm to set up a constrained density field that contains a peak with given shape. Extending the work of Heavens & Peacock, we calculate the joint probability distribution function for several spin parameters and peak mass M using the distribution of peak shapes, for different initial power spectra. The values of observed specific angular momentum versus mass are well fitted by our theoretical isoprobability contours. In contrast, the observed lower values for the specific angular momentum for ellipticals of the same mass cannot be accounted for within our linear regime investigation, highlighting the importance of strongly non-linear phenomena to explain the spin of such objects.

The problem of solving perturbatively the equations describing the evolution of self-gravitating collisionless matter in an expanding universe considerably simplifies when directly formulated in terms of the gravitational and velocity potentials: the problem can be solved {\it exactly}, rather than approximately, even for cosmological models with arbitrary density parameter $\Omega$. The Eulerian approach we present here allows to calculate the higher-order moments of the initially Gaussian density and velocity fields: in particular, we compute the gravitationally induced skewness of the density and velocity-divergence fields for any value of $\Omega$, confirming the extremely weak $\Omega$-dependence of the skewness previously obtained via Lagrangian perturbation theory. Our results show that the separability assumption of higher-order Eulerian perturbative solutions is restricted to the Einstein-de Sitter case only.

We present a new formalism which allows to derive the general Lagrangian dynamical equations for the motion of gravitating particles in a non--flat Friedmann universe with arbitrary density parameter $\Omega$ and no cosmological constant. We treat the set of particles as a Newtonian collisionless fluid. The non--linear dynamical evolution of the fluid element trajectories is then described up to the third--order expansion in Lagrangian coordinates. This work generalizes recent investigations carried out by Bouchet \etal (1992) and Buchert (1994).

We examine the multi--point velocity field for non--Gaussian models as a probe of non--Gaussian behavior. The two--point velocity correlation is not a useful indicator of a non--Gaussian density field, since it depends only on the power spectrum, even for non--Gaussian models. However, we show that the distribution of velocity differences $\bfv_1 - \bfv_2$, where $\bfv_1$ and $\bfv_2$ are measured at the points $\bfr_1$ and $\bfr_2$, respectively, is a good probe of non--Gaussian behavior, in that $p(\bfv_1 - \bfv_2)$ tends to be non--Gaussian whenever the density field is non--Gaussian. As an example, we examine the behavior of $p(\bfv_1 - \bfv_2)$ for non--Gaussian seed models, in which the density field is the convolution of a distribution of points with a set of density profiles. We apply these results to the global texture model.

We discuss the non--linear growth of the kurtosis of the smoothed peculiar velocity field (along an arbitrary direction), in an Einstein--de Sitter universe, induced by Gaussian primordial density fluctuations. Applying the perturbative theory, we show that, for different cosmological models, a departure from the original Gaussian distribution is gravitationally induced only on small scales ($\simlt 20-30$ Mpc). Models with scale--free power spectrum $P(k)\propto k^{\,n}$, with $-1< n \leq 1$, are also considered. When the fluid particles move according to the Zel'dovich approximation the probability distribution of the peculiar velocity field remains unaltered during the evolution.

We consider a weighted biasing scheme for galaxy clustering. This differs from previous treatments in the fact that the biased density field coincides with the background mass--density whenever the latter exceeds a given threshold value. There is some physical motivation for this scheme and it is in better accord with intuitive ideas than models based on the Kaiser (1984) analysis of the clustering of rich clusters. We explain how different classes of object could be biased in different ways with respect to the underlying density distribution but still have $b=1$. We also show that if one applies our scheme consistently a weak dependence of $b$ upon density can be implied. This could also be the reason why the correlation function of galaxies in groups does not differ substantially from the correlation function of all galaxies.

Some recently proposed approximations to follow the non--linear evolution of collisionless matter perturbations in the universe are reviewed. The first one, called frozen--flow approximation, is an Eulerian method within Newtonian theory, and is based on neglecting the role of particle inertia compared to the damping implied by the Hubble drag. The second approach is General Relativistic and Lagrangian; it is based on following the evolution of fluid and geometric observables in the rest frame of each fluid element, under the only approximation of neglecting the back--reaction on the system of the gravitational radiation emitted during non--linear collapse.

We discuss the non--linear growth of the excess kurtosis parameter of the smoothed density fluctuation field $\delta$, $S_4\equiv[\lan\delta^{\,4}\ran-3\lan\delta^{\,2}\ran^2]/ \lan\delta^{\,2}\ran^3$ in an Einstein--de Sitter universe. We assume Gaussian primordial density fluctuations with scale--free power spectrum $P(k)\propto k^{\,n}$ and analyze the dependence of $S_4$ on primordial spectral index $n$, after smoothing with a Gaussian filter. As already known for the skewness ratio $S_3$, the kurtosis parameter is a {\it decreasing function} of $n$, both in exact perturbative theory and in the Zel'dovich approximation. The parameter $S_4$ provides a powerful statistics to test different cosmological scenarios.