We present a comparison of major methodologies of fast generating mock halo or galaxy catalogues. The comparison is done for two-point and the three-point clustering statistics. The reference catalogues are drawn from the BigMultiDark N-body simulation. Both friend-of-friends (including distinct halos only) and spherical overdensity (including distinct halos and subhalos) catalogs have been used with the typical number density of a large-volume galaxy surveys. We demonstrate that a proper biasing model is essential for reproducing the power spectrum at quasilinear and even smaller scales. With respect to various clustering statistics a methodology based on perturbation theory and a realistic biasing model leads to very good agreement with N-body simulations. However, for the quadrupole of the correlation function or the power spectrum, only the method based on semi-N-body simulation could reach high accuracy (1% level) at small scales, i.e., r<25 Mpc/h or k>0.15 h/Mpc. Full N-body solutions will remain indispensable to produce reference catalogues. Nevertheless, we have demonstrated that the far more efficient approximate solvers can reach a few percent accuracy in terms of clustering statistics at the scales interesting for the large-scale structure analysis after calibration with a few reference N-body calculations. This makes them useful for massive production aimed at covariance studies, to scan large parameter spaces, and to estimate uncertainties in data analysis techniques, such as baryon acoustic oscillation reconstruction, redshift distortion measurements, etc.

We present the online MultiDark Database -- a Virtual Observatory-oriented, relational database for hosting various cosmological simulations. The data is accessible via an SQL (Structured Query Language) query interface, which also allows users to directly pose scientific questions, as shown in a number of examples in this paper. Further examples for the usage of the database are given in its extensive online documentation (www.multidark.org). The database is based on the same technology as the Millennium Database, a fact that will greatly facilitate the usage of both suites of cosmological simulations. The first release of the MultiDark Database hosts two 8.6 billion particle cosmological N-body simulations: the Bolshoi (250/h Mpc simulation box, 1/h kpc resolution) and MultiDark Run1 simulation (MDR1, or BigBolshoi, 1000/h Mpc simulation box, 7/h kpc resolution). The extraction methods for halos/subhalos from the raw simulation data, and how this data is structured in the database are explained in this paper. With the first data release, users get full access to halo/subhalo catalogs, various profiles of the halos at redshifts z=0-15, and raw dark matter data for one time-step of the Bolshoi and four time-steps of the MultiDark simulation. Later releases will also include galaxy mock catalogs and additional merging trees for both simulations as well as new large volume simulations with high resolution. This project is further proof of the viability to store and present complex data using relational database technology. We encourage other simulators to publish their results in a similar manner.

[abridged] We present a detailed comparison of fundamental dark matter halo properties retrieved by a substantial number of different halo finders. These codes span a wide range of techniques including friends-of-friends (FOF), spherical-overdensity (SO) and phase-space based algorithms. We further introduce a robust (and publicly available) suite of test scenarios that allows halo finder developers to compare the performance of their codes against those presented here. This set includes mock haloes containing various levels and distributions of substructure at a range of resolutions as well as a cosmological simulation of the large-scale structure of the universe. All the halo finding codes tested could successfully recover the spatial location of our mock haloes. They further returned lists of particles (potentially) belonging to the object that led to coinciding values for the maximum of the circular velocity profile and the radius where it is reached. All the finders based in configuration space struggled to recover substructure that was located close to the centre of the host halo and the radial dependence of the mass recovered varies from finder to finder. Those finders based in phase space could resolve central substructure although they found difficulties in accurately recovering its properties. Via a resolution study we found that most of the finders could not reliably recover substructure containing fewer than 30-40 particles. However, also here the phase space finders excelled by resolving substructure down to 10-20 particles. By comparing the halo finders using a high resolution cosmological volume we found that they agree remarkably well on fundamental properties of astrophysical significance (e.g. mass, position, velocity, and peak of the rotation curve).

We present developing of method of the numerical analysis of polarization in the Gauss--Legendre Sky Pixelization (GLESP) scheme for the CMB maps. This incorporation of the polarization transforms in the pixelization scheme GLESP completes the creation of our new method for the numerical analysis of CMB maps. The comparison of GLESP and HEALPix calculations is done.

We report some results from one of the largest hydrodynamical cosmological simulations of large scale structures that has been done up to date. The MareNostrum Universe SPH simulation consists of 2 billion particles (2 times 1024^3) in a cubic box of 500 h^-1 Mpc on a side. This simulation has been done in the MareNostrum parallel supercomputer at the Barcelona SuperComputer Center. Due to the large simulated volume and good mass resolution, our simulated catalog of dark matter halos comprises more than half a million objects with masses larger than a typical Milky Way galaxy halo. From this dataset we have studied several statistical properties such as the evolution of the halo mass function, the void distribution, the shapes of dark and gas halos and the large scale distribution of baryons.

Properties of $\sim$ 6 000 Ly-$\alpha$ absorbers observed in 19 high resolution spectra of QSOs are investigated using the model of formation and evolution of DM structure elements based on the Zel'dovich theory. This model asserts that absorbers are formed in the course of both linear and nonlinear adiabatic or shock compression of dark matter (DM) and gaseous matter. It allows us to link the column density and overdensity of DM and gaseous components with the observed column density of neutral hydrogen, redshifts and Doppler parameters of absorbers and demonstrates that at high redshifts we observe a self similar period of structure evolution with the Gaussian initial perturbations. We show that the colder absorbers are associated with rapidly expanded regions of a galactic scale which represent large amplitude negative density perturbations. We extend and improve the method of measuring the power spectrum of initial perturbations proposed in Demia\'nski & Doroshkevich (2003b). Our method links the observed separations and the DM column density of absorbers with the correlation function of the initial velocity field. We recover the cold dark matter (CDM) like power spectrum at scales 10> D > 0.15Mpc/h with a precision of ~15%. However at scales $\sim 3 - 150 h^{-1}$kpc the measured and CDM--like spectra are different. This result suggests a possible complex inflation with generation of excess power at small scales.

Using high resolution DM simulations we study the shape of dark matter halos. Halos become more spherical with decreasing mass. This trend is even more pronounced for the inner part of the halo. Angular momentum and shape are correlated. The angular momenta of neighboring halos are correlated.

A new scheme of sky pixelization is developed for CMB maps. The scheme is based on the Gauss--Legendre polynomials zeros and allows one to create strict orthogonal expansion of the map. A corresponding code has been implemented and comparison with other methods has been done.

We report the release of the Gauss--Legendre Sky Pixelization (GLESP) software package version 1.0. In this report we present the main features and functions for processing and manipulation of sky signals. Features for CMB polarization is underway and to be incorporated in a future release. Interested readers can visit http://www.glesp.nbi.dk (www.glesp.nbi.dk) and register for receiving the package.

Properties of $\sim$ 5000 observed Ly-$\alpha$ absorbers are investigated using the model of formation and evolution of DM structure elements based on the Zel'dovich theory. This model is generally consistent with simulations of absorbers formation, accurately describes the Large Scale Structure observed in the galaxy distribution at small redshifts and emphasizes the generic similarity of the LSS and absorbers. The simple physical model of absorbers asserts that they are composed of DM and gaseous matter and it allows us to estimate the column density and overdensity of DM and gaseous components and the entropy of the gas trapped within the DM potential wells. The parameters of DM component are found to be consistent with theoretical expectations for the Gaussian initial perturbations with the WDM--like power spectrum. We demonstrate the influence of the main physical factors responsible for the absorbers evolution. The analysis of redshift distribution of absorbers confirms the self consistence of the assumed physical model and allows to estimate the shape of the initial power spectrum at small scales what in turn restricts the mass of dominant fraction of DM particles to $M_{DM}\approx 0.6 - 2$ keV. Our results verify the redshift variations of intensity of the UV background by about $4 - 6$ times and indicate that, perhaps, the available observations underestimate the intensity of this background.

The main observed properties of Ly-$\alpha$ absorbers are investigated on the basis of theoretical model of formation and evolution of DM structure elements. This model is generally consistent with simulations of absorbers formation and with statistical description of structure evolution based on the Zel'dovich theory. The analysis of redshift variations of comoving linear number density of absorbers was performed in our previous paper. We show that the observed characteristics of Doppler parameter can be related to the size of DM structure elements what allows us to explain the observed distribution of Doppler parameter. This distribution is found to be consistent with the Gaussian initial perturbations. The observed characteristics of entropy and column density, $N_{HI}$, confirm that merging of pancakes is the main evolutionary process at redshifts $z\geq 2$. The observed sample of absorbers characterizes mainly the matter distribution within large low density regions and therefore it is difficult to reconstruct the density field from the distribution of absorbers.

The observed redshift distribution of Ly-$\alpha$ lines and metal systems is examined in order to discriminate and to trace the evolution of structure elements observed in the galaxy distribution, at small redshifts, and to test the theoretical description of structure evolution. We show that the expected evolution of filamentary component of structure describes quite well the redshift distribution of metal systems and stronger Ly-$\alpha$ lines with $\log(N_{HI})\geq$14, at $z\leq$ 3. The redshift distribution of weaker Ly-$\alpha$ lines can be attributed to the population of poorer structure elements (Zel'dovich pancakes), which were formed at high redshifts from the invisible DM and non luminous baryonic matter, and at lower redshifts they mainly merged and dispersed.

The large scale matter distribution in three different simulations of CDM models is investigated and compared with corresponding results of the Zel'dovich theory of nonlinear gravitational instability. We show that the basic characteristics of wall-like structure elements are well described by this theory, and that they can be expressed by the cosmological parameters and a few spectral moments of the perturbation spectrum. Therefore the characteristics of such elements provide reasonable estimates of these parameters. We show that the compressed matter is relaxed and gravitationally confined, what manifests itself in the existence of walls as (quasi)stationary structure elements with life time restricted by their disruption into high density clouds. The matter distribution is investigated both in the real and redshift spaces. In both cases almost the same particles form the walls, and we estimate differences in corresponding wall characteristics. The same methods are applied to several mock catalogues of 'galaxies' what allows us to characterize a large scale bias between the spatial distribution of dark matter and of simulated `galaxies'.

Using high resolution N-body simulations we have studied the formation and evolution of clusters and groups in a LCDM cosmological model. Groups of galaxies have been formed already before z=4. The total number of small bound systems and the total number of galaxies in these small systems rapidly decreases after z=1.5. The fraction of isolated galaxies remains approximately constant after z =1, whereas the fraction of galaxies in groups decreases. Merging of groups and accretion leads to cluster formation at z < 2. Some of the groups merge into large isolated halos.

The simulated matter distribution on large scales is studied using core-sampling, cluster analysis, inertia tensor analysis, and minimal spanning tree techniques. Seven simulations in large boxes for five cosmological models with COBE normalized CDM-like power spectra are studied. The wall-like Super Large Scale Structure with parameters similar to the observed one is found for the OCDM and LCDM models with Omega_m h = 0.3 & 0.245. In these simulations, the rich structure elements with a typical value for the largest extension of \~(30 - 50) h^-1 Mpc incorporate ~40% of matter with overdensity of about 10 above the mean. These rich elements are formed due to the anisotropic nonlinear compression of sheets with original size of ~(15 - 25) h^{-1} Mpc. They surround low-density regions with a typical diameter ~(50 - 70) h^{-1} Mpc. The statistical characteristics of these structures are found to be approximately consistent with observations and theoretical expectations. The cosmological models with higher matter density \Omega_m=1 in CDM with Harrison-Zeldovich or tilted power spectra cannot reproduce the characteristics of the observed galaxy distribution due to the very strong disruption of the rich structure elements. Another model with a broken scale invariant initial power spectrum (BCDM) shows not enough matter concentration in the rich structure elements.

The large supercluster structures obvious in recent galaxy redshift surveys are quantified using an one-dimensional cluster analysis (core sampling) and a three-dimensional cluster analysis based on the minimal spanning tree. The comparison with the LCRS reveals promising stable results. At a mean overdensity of about ten, the supercluster systems form huge wall-like structures comprising about 40% of all galaxies. The overdense clusters have a low mean transverse velocity dispersion of about 400 km/s, i.e. they look quite narrow in redshift space. We performed N-body simulations with large box sizes for six cosmological scenarios. The quantitative analysis shows that the observed structures can be understood best in low density models with $\Omega_m \le 0.5$ with or without a cosmological constant.