We obtain exact expressions for the effect of primordial non-Gaussianity on the matter density perturbation up to second order in a LambdaCDM cosmology, fully accounting for the general relativistic corrections arising on scales comparable with the Hubble radius. We present our results both in the Poisson gauge and in the comoving and synchronous gauge, which are relevant for comparison to different cosmological observables.

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 study the properties of X-ray galaxy clusters in four cold dark matter models with different baryon fraction $\Omega_{BM}$ ranging from 5 to 20 per cent. By using an original three-dimensional hydrodynamic code based on the piecewise parabolic method, we run simulations on a box with size $64 h^{-1}$ Mpc and we identify the clusters by selecting the peaks in the X-ray luminosity field. We analyse these mock catalogues by computing the mass function, the luminosity function, the temperature distribution and the luminosity-temperature relation. By comparing the predictions of the different models to a series of recent observational results, we find that only the models with low baryonic content agree with the data, while models with larger baryon fraction are well outside the 1-$\sigma$ errorbars. In particular, the analysis of the luminosity functions, both bolometric and in the energy band [0.5--2] keV, requires $\Omega_{BM}\mincir 0.05$ when we fix the values $h=0.5$ and $n=0.8$ for the Hubble parameter and the primordial spectral index, respectively. Moreover we find that, independently of the cosmological scenario, all the considered quantities have a very little redshift evolution, particularly between $z=0.5$ and $z=0$.

We present a hydrodynamical code for cosmological simulations which uses the Piecewise Parabolic Method (PPM) to follow the dynamics of gas component and an N-body Particle-Mesh algorithm for the evolution of collisionless component. The gravitational interaction between the two components is regulated by the Poisson equation which is solved by a standard FFT procedure. In order to simulate cosmological flows we have introduced several modifications to the original PPM scheme which we describe in detail. Various tests of the code are presented including adiabatic expansion, single and multiple pancake formation and three-dimensional cosmological simulations with initial conditions based on the cold dark matter scenario.

We have carried out a study of the hydrodynamics of disconnected quark regions during the final stages of the cosmological quark-hadron transition. A set of relativistic Lagrangian equations is presented for following the evaporation of a single quark drop and results from the numerical solution of this are discussed. A self-similar solution is shown to exist and the formation of baryon number density inhomogeneities at the end of the drop contraction is discussed.

We present the results of a series of one-dimensional N-body and hydrodynamical simulations which have been used for testing the different clustering properties of baryonic and dark matter in an expanding background. Initial Gaussian random density perturbations with a power-law spectrum $P(k) \propto k^n$ are assumed. We analyse the distribution of density fluctuations and thermodynamical quantities for different spectral indices $n$ and discuss the statistical properties of clustering in the corresponding simulations. At large scales the final distribution of the two components is very similar while at small scales the dark matter presents a lumpiness which is not found in the baryonic matter. The amplitude of density fluctuations in each component depends on the spectral index $n$ and only for $n=-1$ the amplitude of baryonic density fluctuations is larger than that in the dark component. This result is also confirmed by the behaviour of the bias factor, defined as the ratio between the r.m.s of baryonic and dark matter fluctuations at different scales: while for $n=1,\ 3$ it is always less than unity except at very large scales where it tends to one, for $n=-1$ it is above 1.4 at all scales. All simulations show also that there is not an exact correspondence between the positions of largest peaks in dark and baryonic components, as confirmed by a cross-correlation analysis. The final temperatures depend on the initial spectral index: the highest values are obtained for $n=-1$ and are in proximity of high density regions.

We present results on the non-linear dynamics of inhomogeneous cosmological models with irrotational dust and a positive cosmological constant, considering, in particular, a wide class with vanishing magnetic Weyl tensor. For those patches of the universe that do not recollapse to a singularity we find a unique attractor, representing a de Sitter vacuum phase. For the (re-)collapsing regions we find a family of (Kasner) attractors, so that generically these regions fall in spindle-like singularities. These results give substantial support to the idea that the universe can be very inhomogeneous on ultra-large, super-horizon scales, whit observers living in those (almost) isotropic regions that emerge from an inflationary phase.

We investigate the local non--linear dynamics of irrotational dust with vanishing magnetic part of the Weyl tensor, $H_{ab}$. Once coded in the initial conditions, this dynamical restriction is respected by the relativistic evolution equations. Thus, the outcome of the latter are {\it exact solutions} for special initial conditions with $H_{ab}=0$, but with no symmetries: they describe inhomogeneous triaxial dynamics generalizing that of a fluid element in a Tolman--Bondi, Kantowski--Sachs or Szekeres geometry. A subset of these solutions may be seen as (special) perturbations of Friedmann models, in the sense that there are trajectories in phase--space that pass arbitrarily close to the isotropic ones. We find that the final fate of ever--expanding configurations is a spherical void, locally corresponding to a Milne universe. For collapsing configurations we find a whole family of triaxial attractors, with vanishing local density parameter $\Omega$. These attractors locally correspond to Kasner vacuum solutions: there is a single physical configuration collapsing to a degenerate {\it pancake}, while the generic configuration collapses to a triaxial {\it spindle} singularity. These {\it silent universe} models may provide a fair representation of the universe on super horizon scales. Moreover, one might conjecture that the non--local information carried by $H_{ab}$ becomes negligible during the late highly non--linear stages of collapse, so that the attractors we find may give all of the relevant expansion or collapse configurations of irrotational dust.

A Lagrangian relativistic approach to the non--linear dynamics of cosmological perturbations of an irrotational collisionless fluid is considered. Solutions are given at second order in perturbation theory for the relevant fluid and geometric quantities and compared with the corresponding ones in the Newtonian approximation. Specifically, we compute the density, the volume expansion scalar, the shear, the ``electric" part, or tide, and the ``magnetic" part of the Weyl tensor. The evolution of the shear and the tide beyond the linear regime strongly depends on the ratio of the characteristic size of the perturbation to the cosmological horizon distance. For perturbations on sub--horizon scales the usual Newtonian approximation applies, at least at the considered perturbative order; on super--horizon scales, instead, a new picture emerges, which we call ``silent universe", as each fluid element evolves independently of the environment, being unable to exchange signals with the surrounding matter through either sound waves or gravitational radiation. For perturbations inside the Hubble radius particular attention is paid in singling out non--local effects during the non--linear evolution of fluid elements. These non--local effects are shown to be carried by a traceless and divergenceless tensor, contained in the magnetic part of the Weyl tensor, which is dynamically generated as soon as the system evolves away from the linear regime.

The non--linear dynamics of cosmological perturbations of an irrotational collisionless fluid is analyzed within General Relativity. Relativistic and Newtonian solutions are compared, stressing the different role of boundary conditions in the two theories. Cosmological implications of relativistic effects, already present at second order in perturbation theory, are studied and the dynamical role of the magnetic part of the Weyl tensor is elucidated.

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.