(abridged) We present a model to predict the clustering properties of X-ray clusters in flux-limited surveys. Our technique correctly accounts for past light-cone effects on the observed clustering and follows the non-linear evolution in redshift of the underlying DM correlation function and cluster bias factor. The conversion of the limiting flux of a survey into the corresponding minimum mass of the hosting DM haloes is obtained by using theoretical and empirical relations between mass, temperature and X-ray luminosity of clusters. Finally, our model is calibrated to reproduce the observed cluster counts adopting a temperature-luminosity relation moderately evolving with redshift. We apply our technique to three existing catalogues: BCS, XBACs and REFLEX samples. Moreover, we consider an example of possible future space missions with fainter limiting flux. In general, we find that the amplitude of the spatial correlation function is a decreasing function of the limiting flux and that the EdS models always give smaller correlation amplitudes than open or flat models with low matter density parameter. In the case of XBACs, the comparison with previous estimates of the observational spatial correlation shows that only the predictions of models with Omega_0m=0.3 are in good agreement with the data, while the EdS models have too low a correlation strength. Finally, we use our technique to discuss the best strategy for future surveys. Our results show that the choice of a wide area catalogue, even with a brighter limiting flux, is preferable to a deeper, but with smaller area, survey.

We reconsider the problem of the cosmological reionization due to stellar sources. Using a method similar to that developed by Haiman & Loeb (1997), we investigate the effect of changing the stellar models and the stellar spectra adopted for deriving the ionizing photon production rate. In particular, we study the consequences of adopting zero metallicity stars, which is the natural choice for the first stellar populations. We construct young isochrones representative of Population III stars from existing sets of evolutionary models (Forieri 1982; Cassisi & Castellani 1993) and calculate a suitable library of zero metallicity model atmospheres. The number of ionizing photons emitted by such a zero metal population is about 40% higher than that produced by standard metal poor isochrones. We find that adopting suitable zero metallicity models modifies the reionization epoch. However the latter is still largely affected by current uncertainties in other important physical processes such as the efficiency of the star formation and the fraction of escaping UV photons.

We analyse the spatial clustering properties of the RASS1 Bright Sample, an X-ray flux-limited catalogue of galaxy clusters selected from the southern part of the $ROSAT$ All-Sky Survey. The two-point correlation function $\xi(r)$ of the whole sample is well fitted (in an Einstein-de Sitter model) by the power-law $\xi=(r/r_0)^{-\gamma}$, with $r_0= 21.5^{+3.4}_{-4.4} h^{-1}$ Mpc and $\gamma=2.11^{+0.53}_{-0.56}$ (95.4 per cent confidence level with one fitting parameter). We use the RASS1 Bright Sample as a first application of a theoretical model which aims at predicting the clustering properties of X-ray clusters in flux-limited surveys for different cosmological scenarios. The model uses the theoretical and empirical relations between mass, temperature and X-ray cluster luminosity, and fully accounts for the redshift evolution of the underlying dark matter clustering and cluster bias factor. The comparison between observational results and theoretical predictions shows that the Einstein-de Sitter models display too low a correlation length, while models with a matter density parameter $\Omega_{\rm 0m}=0.3$ (with or without a cosmological constant) are successful in reproducing the observed clustering. The dependence of the correlation length $r_0$ on the X-ray limiting flux and luminosity of the sample is generally consistent with the predictions of all our models. Quantitative agreement is however only reached for $\Omega_{\rm 0m} = 0.3$ models. The model presented here can be reliably applied to future deeper X-ray cluster surveys: the study of their clustering properties will provide a useful complementary tool to the traditional cluster abundance analyses to constrain the cosmological parameters.

We present a theoretical model which aims at predicting the clustering properties of X-ray clusters in flux-limited surveys for different cosmological scenarios. The model uses the theoretical and empirical relations between mass, temperature and X-ray luminosity, and fully accounts for the redshift evolution of the underlying dark matter clustering and cluster bias factor. We apply the model to the RASS1 Bright Sample and to the XBACs catalogue. The results show that the Einstein-de Sitter models display too low a correlation length, while models with a matter density parameter Omega_0m=0.3 (with or without a cosmological constant) are successful in reproducing the observed clustering.

(abridged) The evolution of galaxy clustering from z=0 to z=4.5 is analyzed using the angular correlation function and the photometric redshift distribution of galaxies brighter than I_{AB}\le 28.5 in the HDF North. The reliability of the photometric redshift estimates is discussed on the basis of the available spectroscopic redshifts, comparing different codes and investigating the effects of photometric errors. The redshift bins in which the clustering properties are measured are then optimized to take into account the uncertainties of the photometric redshifts. The results show that the comoving correlation length has a small decrease in the range 0<z<1 followed by an increase at higher z. We compare these results with the theoretical predictions of a variety of cosmological models belonging to the general class of CDM. The comparison with the expected mass clustering evolution indicates that the observed high-redshift galaxies are biased tracers of the dark matter with an effective bias b strongly increasing with redshift. Assuming an Einstein-de Sitter universe, we obtain b\simeq 2 at z=2 and b\simeq 5 at z=4. A comparison of the clustering amplitudes that we measured at z=3 with those reported for LBG suggests that the clustering depends on the abundance of the objects: more abundant objects are less clustered, as expected in the paradigm of hierarchical galaxy formation. The strong clustering and high bias measured at z=3 are consistent with the expected density of massive haloes predicted for the various cosmologies here considered. At z=4, the strong clustering observed in the HDF requires a significant fraction of massive haloes to be already formed by that epoch. This feature could be a discriminant test for the cosmological parameters if confirmed by future observations.

We calculate the statistical clustering of Lyman-break galaxies predicted in a selection of currently-fashionable structure formation scenarios. These models are all based on the cold dark matter model, but vary in the amount of dark matter, the initial perturbation spectrum, the background cosmology and in the presence or absence of a cosmological constant term. If Lyman-break galaxies form as a result of hierarchical merging, the amplitude of clustering depends quite sensitively on the minimum halo mass that can host such a galaxy. Interpretation of the recent observations by Giavalisco et al. (1998) would therefore be considerably clarified by a direct determination of the relevant halo properties. For a typical halo mass around $10^{11} h^{-1} M_\odot$ the observations do not discriminate strongly between cosmological models, but if the appropriate mass is larger, say $10^{12} ~h^{-1} M_\odot$ (which seems likely on theoretical grounds), then the data strongly favour models with a low matter-density.

We discuss the theoretical interpretation of observational data concerning the clustering of galaxies at high redshifts. Building on the theoretical machinery developed by Matarrese et al. (1997), we make detailed quantitative predictions of galaxy clustering statistics for a variety of cosmological models, taking into account differences in spatial geometry and initial fluctuation spectra and exploring the role of bias as a complicating factor in these calculations. We demonstrate that the usual description of evolution (in terms of the parameters $\epsilon$ and $r_0$) is not useful for realistic galaxy clustering models. We compare the detailed predictions of the variation of correlation functions with redshift against current observational data to constrain available models of structure formation. Theories that fit the present-day abundance of rich clusters are generally compatible with the observed redshift evolution of galaxy clustering if galaxies are no more than slightly biased at $z\sim 1$. We also discuss the interpretation of a concentration of Lyman-break galaxies found by Steidel et al. (1998), coming to the conclusion that such concentrations are not unexpected in `standard' models of structure formation.

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.

We discuss how the redshift dependence of the observed two-point correlation function of various classes of objects can be related to theoretical predictions. This relation involves first a calculation of the redshift evolution of the underlying matter correlations. The next step is to relate fluctuations in mass to those of any particular class of cosmic objects; in general terms, this means a model for the bias and how it evolves with cosmic epoch. Only after these two effects have been quantified can one perform an appropriate convolution of the non-linearly evolved two-point correlation function of the objects with their redshift distribution to obtain the `observed' correlation function for a given sample. This convolution in itself tends to mask the effect of evolution by mixing amplitudes at different redshifts. We develop a formalism which incorporates these requirements and, in particular, a set of plausible models for the evolution of the bias factor. We apply this formalism to the spatial, angular and projected correlation functions from different samples of high-redshift objects, assuming a simple phenomenological model for the initial power-spectrum and an Einstein-de Sitter cosmological model. We find that our model is roughly consistent with data on the evolution of QSO and galaxy clustering, but only if the effective degree of biasing is small. We discuss the differences between our analysis and other theoretical studies of clustering evolution and argue that the dominant barrier to making definitive predictions is uncertainty about the appropriate form of the bias and its evolution with cosmic epoch.

We analyse the effect of two relevant physical constraints on the mass multiplicity function of dark matter halos in a Press--Schechter type algorithm. Considering the random--walk of linear Gaussian density fluctuations as a function of the smoothing scale, we simultaneously i) account for mass semi--positivity and ii) avoid the cloud--in--cloud problem. It is shown that the former constraint implies a severe cutoff of low--mass objects, balanced by an increase on larger mass scales. The analysis is performed both for scale--free power--spectra and for the standard cold dark matter model. Our approach shows that the well--known ``infrared" divergence of the standard Press--Schechter mass function is caused by unphysical, negative mass events which inevitably occur in a Gaussian distribution of density fluctuations.

Recent data on the high--redshift abundance of damped Ly$\alpha$ systems are compared with theoretical predictions for `blue' (i.e. $n>1$) Mixed Dark Matter models. The results show that decreasing the hot component fraction $\Omega_\nu$ and/or increasing the primordial spectral index $n$ implies an earlier epoch of cosmic structure formation. However, we also show that varying $\Omega_\nu$ and $n$ in these directions makes the models barely consistent with the observed abundance of galaxy clusters. Therefore, requiring at the same time observational constraints on damped Ly$\alpha$ systems and cluster abundance to be satisfied represents a challenge for the Mixed Dark Matter class of models.

Recently, there has been growing interest on primordial ``blue" ($n>1$) perturbation spectra, motivated both by a composite set of observational data on large scales and from the point of view of theoretical model building. After reviewing the theoretical (inflationary) motivations for these blue spectra, we consider various observational constraints, within both Cold Dark Matter and Mixed Dark Matter scenarios. In particular, using linear theory, we discuss large--scale bulk flows and the X--ray cluster abundance. We also perform various N--body simulations of these models to study the clustering properties of the matter distribution and the peculiar velocity field.

We compute the contribution to the three--point temperature correlation function of the Cosmic Microwave Background coming from the non--linear evolution of Gaussian initial perturbations, as described by the Rees--Sciama (or integrated Sachs--Wolfe) effect. By expressing the collapsed three--point function in terms of multipole amplitudes, we are able to calculate its expectation value for any power spectrum and for any experimental setting on large angular scales. We also give an analytical expression for the {\it rms} collapsed three--point function arising from the cosmic variance of a Gaussian fluctuation field. In the case of {\it COBE} DMR, we find that the predicted signal is about three orders of magnitude below that expected from the cosmic variance.

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 analyse the large--scale velocity field obtained by N--body simulations of cold dark matter (CDM) models with non--Gaussian primordial density fluctuations, considering models with both positive and negative primordial skewness in the density fluctuation distribution. We study the velocity probability distribution and calculate the dependence of the bulk flow, one--point velocity dispersion and Cosmic Mach Number on the filtering size. We find that the sign of the primordial skewness of the density field provides poor discriminatory power on the evolved velocity field. All non--Gaussian models here considered tend to have lower velocity dispersion and bulk flow than the standard Gaussian CDM model, while the Cosmic Mach Number turns out to be a poor statistic in characterizing the models. Next, we compare the large--scale velocity field of a composite sample of optically selected galaxies as described by the Local Group properties, bulk flow, velocity correlation function and Cosmic Mach Number, with the velocity field of mock catalogues extracted from the N--body simulations. The comparison does not clearly permit to single out a best model: the standard Gaussian model is however marginally preferred by the maximum likelihood analysis.

We analyze the cosmic peculiar velocity field as traced by a sample of 1184 spiral, elliptical and S0 galaxies, grouped in 704 objects. We perform a statistical analysis, by calculating the bulk flow, Cosmic Mach Number and velocity correlation function for this sample and for mock catalogs extracted from a set of N--body simulations. We run four cold dark matter (CDM) simulations: two tilted models (with spectral index $n=0.6$ and $n=0.8$), the standard model ($n=1$) and a ``blue" one ($n=1.2$), with different values of the linear bias parameter $b$. By means of a Maximum Likelihood analysis we estimate the ability of our models to fit the observations, as measured by the above statistics, and to reproduce the Local Group properties. On the basis of this analysis we conclude that the best model is the unbiased standard model $(n,~b)=(1,~1)$, even though the overall flatness of the joint likelihood function implies that one cannot strongly discriminate models in the range $0.8 \le n \le 1$, and $1 \le b \le 1.5$. Models with $b \geq 2.5$ are rejected at the $95\%$ confidence level. For $n=0.8$ the values of $b$ preferred by the present analysis, together with the {\em COBE} data, require a negligible contribution to $\Delta T/T$ by gravitational waves. Finally, the blue model, normalized to {\em COBE}, does not provide a good fit to the velocity data.

{}From three quasar samples with a total of 1038 objects in the redshift range $1.0 \div 2.2$ we measure the variance $\sigma^2$ of counts in cells of volume $V_u$. By a maximum likelihood analysis applied separately on these samples we obtain estimates of $\sigma^2(\ell)$, with $\ell \equiv V_u^{1/3}$. The analysis from a single catalog for $\ell = ~40~h^{-1}$ Mpc and from a suitable average over the three catalogs for $\ell = ~60,~80$ and $100~h^{-1}$ Mpc, gives $\sigma^2(\ell) = 0.46^{+0.27}_{-0.27}$, $0.18^{+0.14}_{-0.15}$, $0.05^{+0.14}_{-0.05}$ and $0.12^{+0.13}_{-0.12}$, respectively, where the $70\%$ confidence ranges account for both sampling errors and statistical fluctuations in the counts. This allows a comparison of QSO clustering on large scales with analogous data recently obtained both for optical and IRAS galaxies: QSOs seem to be more clustered than these galaxies by a biasing factor $b_{QSO}/b_{gal} \sim 1.4 - 2.3$.

We analyze the temperature three--point correlation function and the skewness of the Cosmic Microwave Background (CMB), providing general relations in terms of multipole coefficients. We then focus on applications to large angular scale anisotropies, such as those measured by the {\em COBE} DMR, calculating the contribution to these quantities from primordial, inflation generated, scalar perturbations, via the Sachs--Wolfe effect. Using the techniques of stochastic inflation we are able to provide a {\it universal} expression for the ensemble averaged three--point function and for the corresponding skewness, which accounts for all primordial second--order effects. These general expressions would moreover apply to any situation where the bispectrum of the primordial gravitational potential has a {\em hierarchical} form. Our results are then specialized to a number of relevant models: power--law inflation driven by an exponential potential, chaotic inflation with a quartic and quadratic potential and a particular case of hybrid inflation. In all these cases non--Gaussian effects are small: as an example, the {\em mean} skewness is much smaller than the cosmic {\em rms} skewness implied by a Gaussian temperature fluctuation field.

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 investigate inflationary models leading to density perturbations with a spectral index $n>1$ (``blue spectra"). These perturbation spectra may be useful to simultaneously account for both the amount of ultra large-scale power required to fit cosmic microwave background anisotropies, such as those measured by COBE, and that required to give bulk motions and structures on the $\sim 50~h^{-1}$ Mpc scale.

We study the cosmic peculiar velocity field as traced by a sample of 1184 spiral, elliptical and S0 galaxies, grouped in 704 objects. We carry out a statistical analysis, by calculating bulk flows and velocity correlation functions for this sample and for mock catalogs which we extract from N--body simulations. For the simulations we consider tilted (i.e. with spectral index $n\leq 1$) CDM models with different values of the linear bias parameter $b$. By mean of a maximum likelihood analysis we estimate the ability of the models in fitting the observations as measured by the above statistics and in reproducing the Local Group properties.

We investigate the accuracy of the frozen--flow approximation (FFA), recently proposed by Matarrese \etal (1992), for following the nonlinear evolution of cosmological density fluctuations under gravitational instability. We compare a number of statistics between results of the FFA and nbody simulations, including those used by Melott, Pellman \& Shandarin (1993) to test the Zel'dovich approximation. The FFA performs reasonably well in a statistical sense, e.g. in reproducing the counts--in--cell distribution, at small scales, but it does poorly in the crosscorrelation with nbody which means it is generally not moving mass to the right place, especially in models with high small--scale power.

We calculate reduced moments $\overline \xi_q$ of the matter density fluctuations, up to order $q=5$, from counts in cells produced by Particle--Mesh numerical simulations with scale--free Gaussian initial conditions. We use power--law spectra $P(k) \propto k^n$ with indices $n=-3,~-2,~-1,~0,~1$. Due to the supposed absence of characteristic times or scales in our models, all quantities are expected to depend on a single scaling variable. For each model, the moments at all times can be expressed in terms of the variance $\overline \xi_2$, alone. We look for agreement with the hierarchical scaling ansatz, according to which $\overline \xi_q \propto \overline \xi_2^{~q-1}$. For $n\leq -2$ models we find strong deviations from the hierarchy, which are mostly due to the presence of boundary problems in the simulations. A small, residual signal of deviation from the hierarchical scaling is however also found in $n \geq -1$ models. For the first time, due to our large dynamic range and careful checks of scaling and shot--noise effects, we are able to detect evolution away from the perturbation theory result.

We investigate the evolution of the skewness of the distribution of density fluctuations in CDM models with both Gaussian and non--Gaussian initial fluctuations. We show that the method proposed by Coles \& Frenk (1991), which uses the skewness of galaxy counts to test the hypothesis of Gaussian primordial density fluctuations, is a potentially powerful probe of initial conditions. As expected, the mass distribution in models with initially non--Gaussian fluctuations shows systematic departures from the Gaussian behaviour on intermediate to large scales. We investigate the effect of peculiar velocity distortions and normalisation upon the relationship between skewness and variance. These effects are generally small for the models we consider. Comparing our results to the QDOT measurements of the skewness, we find that our initially positive--skew models are clearly excluded by this analysis, but the available data do not rule out the negative--skew models.

We generate artificial Lick maps using $N$--body simulations and compare the angular correlation function, $w(\vartheta)$, measured from the simulations with the APM correlation. For the Gaussian CDM model, neither the standard biassed model nor a more evolved model (as suggested by the COBE data), reproduce the correlations on large angular scales. We come to a similar conclusion about CDM models with positively skewed initial fluctuation distributions. In contrast, models with initially negatively skewed fluctuations produce a $w(\vartheta)$ that declines much more gently on large scales. Such models are therefore in principle capable of reconciling the lack of large--scale power of the CDM spectrum with the observed clustering of APM galaxies.