In the current paper we consider the prospects of observation of gravitationally lensed sources in far-infrared and submillimeter wavelength ranges by the future space telescopes equipped with actively cooled main mirrors. We consider the possibility of solving some of the important cosmological and astrophysical scientific tasks by observation of gravitationally lensed sources. The number counts of lensed sources for wavelengths from 70$\mu m$ up to 2000$\mu m$ were calculated. We discuss the distribution of lensed sources at different redshifts and magnification coefficients and also mass distribution of lenses. Model sky maps that illustrate the contribution of lensed sources were created.

Recent advances in observational astronomy allow to study various groups of Dark Matter (DM) dominated objects from the dwarf spheroidal (dSph) galaxies to clusters of galaxies that span the mass range from $10^{6} M_{\odot}$ to $10^{15} M_{\odot}$. To analyze data of this divers collection of objects we used a simple toy model of spherical DM halo formation that was initially proposed by Peebles. This model introduced the concept of the epoch or redshift of halo formation. Using this concept we analyzed selected sample of DM dominated objects and we have found empirical correlations between the virial mass, $M_{vir}$, of halos and basic parameters of their cores, namely, the mean DM density, pressure and entropy. These correlations are a natural result of similar evolution of all such objects. It is driven mainly by gravitational interactions what implies a high degree of self similarity of both the process of halos formation and their internal structure. We confirmed the CDM--like shape of both the small and large scale power spectrum. However, our reconstruction of the evolutionary history of observed objects differs from expectations of the standard $\Lambda$CDM cosmology and requires either a multicomponent composition of DM or a more complex primordial power spectrum of density perturbations with significant excess of power at scales of clusters of galaxies and larger. We demonstrated that a model with suitable combination of the heavy DM particles (CDM) and DM particles with large damping scale (HDM) could provide a successful description of the observational data in a wide range of masses.

We analyzed available observational data of a sample of dark matter (DM) dominated galaxies and clusters of galaxies and we have found correlations between the virial mass, $M_{vir}$, of halos and basic parameters of their cores, namely, the mean DM density, pressure and entropy. These correlations are a natural consequence of similar evolution of all such objects. It is driven mainly by gravitational interactions what implies a high degree of self similarity of both the process of halos formation and their internal structure. We confirmed the CDM--like shape of both the small and large scale power spectrum. However, our reconstruction of the evolutionary history of observed objects requires either a multicomponent composition of DM or a more complex primordial power spectrum of density perturbations with significant excess of power at scales of clusters of galaxies and larger. We demonstrated that a model with suitable combination of the heavy DM particles (CDM) and DM particles with large damping scale (HDM) could provide a successful description of the observational data in a wide range of masses.

This paper describes outstanding issues in astrophysics and cosmology that can be solved by astronomical observations in a broad spectral range from far infrared to millimeter wavelengths. The discussed problems related to the formation of stars and planets, galaxies and the interstellar medium, studies of black holes and the development of the cosmological model can be addressed by the planned space observatory Millimetron (the "Spectr-M" project) equipped with a cooled 10-m mirror. Millimetron can operate both as a single-dish telescope and as a part of a space-ground interferometer with very long baseline.

The mass and composition of dark matter (DM) particles and the shape and damping scales of the power spectrum of density perturbations can be estimated from recent observations of the DM dominated relaxed objects -- dwarf galaxies and clusters of galaxies. We confirm that the observed velocity dispersion of dSph galaxies agrees with the possible existence of DM particles with mass $m_w\sim 15 - 20keV$. More complex analysis utilizes the well known semi analytical model of formation of DM halos in order to describe the basic properties of corresponding objects and to estimate their redshifts of formation. For the DM halos this redshift is determined by their masses and the initial power spectrum of density perturbations. This correlation allows us to partly reconstruct the small scale spectrum of perturbations. We consider the available sample of suitable observed objects that includes $\sim 40$ DM dominated galaxies and $\sim 40$ clusters of galaxies and we show that the observed characteristics of these objects are inconsistent with expectations of the standard $\Lambda$CDM cosmological model. However, they are consistent with a more complex DM model with a significant contribution of the hot DM--like power spectrum with relatively large damping scale ($\sim 10 - 30Mpc$). The HDM component of DM decelerates but does not prevent formation of low mass objects. These preliminary inferences require confirmation by a more representative observational data that should include -- if possible -- DM dominated objects with intermediate masses $M\sim 10^{10} - 10^{12} M_\odot$. Comparison of observed properties of such objects with numerical simulations will provide more detailed picture of the process of formation of DM halos.

We apply the well known semi analytical model of formation of DM halos to discuss the process of formation and evolution of the relaxed objects dominated by the DM component (such as the first and dSph galaxies and clusters of galaxies). On the one hand this approach allows us to obtain a simple description of evolution of the first galaxies. On the other hand it reveals links between the observed characteristics of the relaxed DM halos and the initial power spectrum of density perturbations. Results of our analysis of the observed properties of 41 DM dominated dSph galaxies and data of $\sim 300$ clusters of galaxies are consistent with the $\Lambda$CDM like power spectrum of initial perturbations down to the scale of $\sim 10kpc$.

This work investigates a set of cosmological collisionless N-body simulations with featured power spectra of initial perturbations in the context of the core-cusp and satellites problems. On the studied power spectra some scales of fluctuations were suppressed. Such spectral features can be caused by multicomponent dark matter. The density profiles innermost resolved slopes $\alpha\equiv d \log(\rho) /d \log(r) $ of the five largest haloes were measured and its dependence on the parameters of the suppression was traced. In a certain range of the parameters the slopes flatten from initial value of about -1.2 to -0.6 or even to -0.2 in one of the cases. This qualitatively demonstrates that (i) profiles shape depends on initial power spectrum and (ii) this effect may be responsible for the solution of the core-cusp problem. The suppression of some part of the initial power spectrum also leads to the decrease of the number of massive subhaloes.

With help of a set of toy N-body models of dark halo formation we study the impact of small scale initial perturbations on the inner density profiles of haloes. We find a significant flattening of the inner slope $\alpha={d \log \rho \over d \log r}$ to $\alpha=-0.5$ in some range of scales and amplitudes of the perturbations (while in the case of absence of these perturbations the NFW profile with $\alpha=-1$ is reproduced). This effect may be responsible for the formation of cuspless galactic haloes.

We discuss various aspects of the inner structure formation in virialized dark matter (DM) halos that form as primordial density inhomogeneities evolve in the cosmological standard model. The main focus is on the study of central cusps/cores and of the profiles of DM halo rotation curves, problems that reveal disagreements among the theory, numerical simulations, and observations. A method that was developed by the authors to describe equilibrium DM systems is presented, which allows investigating these complex nonlinear structures analytically and relating density distribution profiles within a halo both to the parameters of the initial small-scale inhomogeneity field and to the nonlinear relaxation characteristics of gravitationally compressed matter. It is shown that cosmological random motions of matter `heat up' the DM particles in collapsing halos, suppressing cusp-like density profiles within developing halos, facilitating the formation of DM cores in galaxies, and providing an explanation for the difference between observed and simulated galactic rotation curves. The analytic conclusions obtained within this approach can be confirmed by the N-body model simulation once improved spatial resolution is achieved for central halo regions.

We analyze evolution of the basic properties of simulated large scale structure elements formed by dark matter (DM LSS) and confront it with the observed evolution of the Lyman-$\alpha$ forest. In three high resolution simulations we selected samples of compact DM clouds of moderate overdensity. Clouds are selected at redshifts $0\leq z\leq 3$ with the Minimal Spanning Tree (MST) technique. The main properties of so selected clouds are analyzed in 3D space and with the core sampling approach, what allows us to compare estimates of the DM LSS evolution obtained with two different techniques and to clarify some important aspects of the LSS evolution. In both cases we find that regular redshift variations of the mean characteristics of the DM LSS are accompanied only by small variations of their PDFs, what indicates the self similar character of the DM LSS evolution. The high degree of relaxation of DM particles compressed within the LSS is found along the shortest principal axis of clouds. We see that the internal structure of selected clouds depends upon the mass resolution and scale of perturbations achieved in simulations. It is found that the low mass tail of the PDFs of the LSS characteristics depends upon the procedure of clouds selection.

We show that the popular ILC approach is unstable in respect to the division of the sample of map pixels to the set of ``homogeneous'' subsamples. For suitable choice of such subsamples we can obtain the restored CMB signal with amplitudes ranged from zero to the amplitudes of the observed signal. We propose approach which allows us to obtain reasonable estimates of $C_\ell$ at $\ell\leq 30$ and similar to WMAP $C_\ell$ for larger $\ell$. With this approach we reduce some anomalies of the WMAP results. In particular, our estimate of the quadrupole is well consistent to theoretical one, the effect of the ``axis of evil'' is suppressed and the symmetry of the north and south galactic hemispheres increases. This results can change estimates of quadrupole polarization and the redshift of reionization of the Universe. We propose also new simple approach which can improve WMAP estimates of high $\ell$ power spectrum.

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 analyze basic properties of about 200 metal systems observed in high resolution spectra of 9 quasars by Boksenberg, Sargent & Rauch (2003). The measured Doppler parameters for the hydrogen and carbon lines are found to be different by a factor $\sim\sqrt{m_C/ m_H}$ what indicates the domination of the thermal broadening and high degree of relaxation for majority of the observed metal systems. The analysis of the mean separation of the metal systems confirms that they can be located in regions with the typical proper size $R_{cl}\sim 0.1 - 0.25h^{-1}$ Mpc what corresponds to the typical sizes of protogalaxies with the baryonic masses $M_b\approx 10^9-10^{10}M_\odot$. The metal abundances observed at redshifts $z=3 - 6$ within both IGM and galaxies are found to be quite similar to each other what indicates the strong interaction of young galaxies with the IGM and the important role of dwarf satellites of host galaxies in such interaction. On the other hand, this metal abundance can be considered as the integral measure of nuclear reactions within any stars what in turn restricts the contribution of the stars to the creation of ionizing UV background. Available now estimates of these abundances demonstrate that attempts to explain the reionization of the Universe as a byproduct of the process of synthesis of metals only are problematic and a significant contribution of unobserved and/or non thermal UV sources seems to be required.

We propose a solution of the cusp problem in framework of the standard $\Lambda$CDM cosmology. To do this we describe the linear and nonlinear periods of halo formation by the entropy function of dark matter particles. This approach allows us to take into account together the impact of both the processes of nonlinear relaxation of compressed matter and the small scale initial velocity perturbations in collapsed halos. We show that such random velocities lead to the random variations of the density profile of relaxed halos. As a rule, they suppress the formation of cusp--like halos and favor the creation of core--like ones. This approach allows us to reproduce observed rotation curves, to explain their random scatter and deviations from simulated ones.

We show that the difference between the theoretically expected and measured by WMAP amplitude of the quadrupole fluctuations of CMB can be related to the impact of the anisotropic curvature of the homogeneous universe dominated by the dark energy. In such universe the matter expansion becomes practically isotropic just after the period of inflation and only at small redshifts the anisotropic expansion is generated again by the small curvature $\Omega_K=1- \Omega_m- \Omega_\Lambda \leq 10^{-4}$. For such models the possible deviations from the parameters derived for the standard cosmological model are evidently negligible but the correlations of large scale perturbations and distortions of their Gaussianity are possible. Such models are also compatible with existence of a homogeneous magnetic field and matter rotation which contribute to the low $\ell$ anisotropy and can be considered as ``hidden parameters'' of the model. Their influence can be observed as, for example, special correlations of small scale fluctuations and the Faraday rotation of the CMB and radiation of the farthest quasars. However, both the magnetic field and matter rotation require also modifications of the simple models of isotropic inflation and they change the evolutionary history of the early Universe.

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.

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.

We present circular and linear cross-correlation tests and the "friend-of-friend" analysis for phases of the Internal Linear Combination Map (ILC) and the WMAP foregrounds for all K--W frequency bands at the range of multipoles $\ell\le100$. We compare also Tegmark, de Oliveira-Costa and Hamilton (2003) and Naselsky et al. (2003) cleaned maps with corresponding foregrounds. We have found significant deviations from the expected Poissonian statistics for all the cleaned maps and foregrounds. Our analysis shows that, for a low multipole range of the cleaned maps, power spectra contains some of the foregrounds residuals mainly from the W band.

We present a circular cross-correlation tests for the phases of the Internal Linear Combination Map (ILC) and {\it WMAP}'s foregrounds for all K--W frequency bands at the range of multipoles $\ell\le 50$. We have found significant deviations from the expected Poissonian statistics for the ILC and the foregrounds phases. Our analysis shows that the low multipole range of the ILC power spectrum contains some of the foregrounds residues.

The Large Scale Structure in the galaxy distribution is investigated using The First Data Release of the Sloan Digital Sky Survey. Using the Minimal Spanning Tree technique we have extracted sets of filaments, of wall--like structures, of galaxy groups, and of rich clusters from this unique sample. The physical properties of these structures were then measured and compared with the statistical expectations based on the Zel'dovich' theory. The measured characteristics of galaxy walls were found to be consistent with those for a spatially flat $\Lambda$CDM cosmological model with $\Omega_m\approx$ 0.3 and $\Omega_\Lambda \approx$ 0.7, and for Gaussian initial perturbations with a Harrison -- Zel'dovich power spectrum. Furthermore, we found that the mass functions of groups and of unrelaxed structure elements generally fit well with the expectations from Zel'dovich' theory. We also note that both groups and rich clusters tend to prefer the environments of walls, which tend to be of higher density, rather than the environments of filaments, which tend to be of lower density.

We present a new method of measuring the power spectrum of initial perturbations to an unprecedently small scale of $\sim$ 10$h^{-1}$ kpc. We apply this method to a sample of 4500 Ly-$\alpha$ absorbers and recover the cold dark matter (CDM) like power spectrum at scales $\geq 300h^{-1}$kpc with a precision of $\sim$ 10%. However at scales $\sim 10 - 300 h^{-1}$kpc the measured and CDM--like spectra are noticeable different. This result suggests a complex inflation with generation of excess power at small scales. The magnitude and reliability of these deviations depend also upon the possible incompleteness of our sample and poorly understood process of formation of weak absorbers. Confirmation of the CDM--like shape of the initial power spectrum or detection of its distortions at small scales are equally important for widely discussed problems of physics of the early Universe, galaxy formation, and reheating of the Universe. Our method links the observed mass function of absorbers with the correlation function of the initial velocity field and therefore it avoids the Nyquist restrictions limiting the investigations based on the smoothed flux or density fields. The physical model of absorbers adopted here asserts that they are formed in the course of both linear and nonlinear adiabatic or shock compression of dark matter (DM) and gaseous matter. At scales $\geq 1h^{-1}$Mpc all characteristics of the DM component and, in particular, their redshift distribution are found to be consistent with theoretical expectations for Gaussian initial perturbations with a CDM--like power spectrum.

The paper is devoted to the methods of determination of the cosmological parameters from recent CMB observations. We show that the more complex models of kinetics of recombination with a few "missing" parameters describing the recombination process provide better agreement between measured and expected characteristics of the CMB anisotropy. In particular, we consider the external sources of the Ly-{alpha} and Ly-{c} radiation and the model with the strong clustering of baryonic component. These factors can constrain the estimates of the cosmological parameters usually discussed. We demonstrate also that the measurements of polarization can improve these estimates and, for the precision expected for the PLANCK mission, allow to discriminate a wide class of models.

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 Large Scale Structure (LSS) in the galaxy distribution is investigated using the Sloan Digital Sky Survey Early Data Release (SDSS EDR). Using the Minimal Spanning Tree technique we have extracted sets of filaments, of wall-like structures, of galaxy groups, and of rich clusters from this unique sample. The physical properties of these structures were then measured and compared with the expectations from Zel'dovich' theory. The measured characteristics of galaxy walls were found to be consistent with those for a spatially flat $\Lambda$CDM cosmological model with $\Omega_m\approx$ 0.3 and $\Omega_\Lambda \approx$ 0.7, and for Gaussian initial perturbations with a Harrison -- Zel'dovich power spectrum. Furthermore, we found that the mass functions of groups and of unrelaxed structure elements generally fit well with the expectations from Zel'dovich' theory, although there was some discrepancy for lower mass groups which may be due to incompleteness in the selected sample of groups. We also note that both groups and rich clusters tend to prefer the environments of walls, which tend to be of higher density, rather than the environments of filaments, which tend to be of lower density. Finally, we note evidence of systematic differences in the properties of the LSS between the Northern Galactic Cap stripe and the Southern Galactic Cap stripe -- in particular, in the physical properties of the walls, their spatial distribution, and the relative numbers of clusters embedded in walls. Because the mean separation of walls is $\approx$ 60 -- 70$h^{-1}$ Mpc, each stripe only intersects a few tens of walls. Thus, small number statistics and cosmic variance are the likely drivers of these systematic differences.

In this paper we discuss the possible distortions of the ionization history of the universe caused by an injection of non-thermal energy due to decays of hypothetical Super Heavy Dark Matter (SHDM) particles. These particles are usually considered as a possible source of Ultra High Energy Cosmic Rays (UHECRs) in the framework of the Top-Down model. Estimates of fraction of energy of decays converted to the UV range show that, for suitable parameters of SHDM particles, the significant distortions of power spectra of CMB anisotropy appear. Comparison with the observed power spectrum allows to restrict some properties of the SHDM particles. These decays can also increase of about 5 -- 10 times the degree of ionization of hydrogen at redshifts $z\sim$ 10 -- 50 that essentially accelerates the formation of molecules $H_2$ and first stars during "dark ages".

The observational tests for the period of reionization of the universe are discussed. We show that this period can be observed as {\it negative} intensity patches of the CMB radiation with the amplitude $\delta T/T\sim 10^{-5}$ and the angular sizes $\theta_T\sim$10 angular seconds in range of the wavelength 0.1 cm$\leq\lambda\leq$1 cm. The expected number density and frequency dependence of the amplitude permit to recognize this effect and to discriminate it from the noise. This method applied to the small scale variations of CMB temperature complements well the traditional spectral approach. The number density and the amplitude of observed 'negative' intensity patches depend upon the redshift of reionization that allow to estimate roughly this redshift. The ionized bubbles formed just before the period of reionization could also be seen as the highest peaks. The expected results are sensitive to the Jeans scale at the period of reionization, to the small scale shape of the primordial power spectrum and to the mass of dark matter particles.

In this paper we briefly discuss the problem of the origin of Ultra High Energy Cosmic Rays in the framework of Top-Down models. We show that, for high energy of decays and in a wide range of spectra of injected protons, their extragalactic flux is consistent with the observed fluxes of cosmic rays in the energy range 0.1 E_{GZK}< E < 10E_{GZK}. For suitable energy and spectra of injected protons, the contribution of galactic sources is moderate, in this energy range, but it dominates at smaller and larger energies. In such models we can expect that at these energies the anisotropy of cosmic rays distribution over sky will be especially small. Some possible manifestations of decays of super massive particles such as, for example, primordial black holes with masses M_{pbh} ~ 10^{-5} g, are considered. In particular, we show that partial conversion of energy released during these decays at redshifts z ~ 1000 to Ly-alpha photons can delay the hydrogen recombination and distort the spectrum of fluctuations of the cosmic microwave background radiation.

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'.

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.

An approximate statistical description of the formation and evolution of structure of the universe based on the Zel'dovich theory of gravitational instability is proposed. It is found that the evolution of DM structure shows features of self-similarity and the main structure characteristics can be expressed through the parameters of initial power spectrum and cosmological model. For the CDM-like power spectrum and suitable parameters of the cosmological model the effective matter compression reaches the observed scales $R_{wall}\sim $20 -- 25$h^{-1}$Mpc with the typical mean separation of wall-like elements $D_{SLSS}\sim $50 -- 70$h^{-1}$Mpc. This description can be directly applied to the deep pencil beam galactic surveys and absorption spectra of quasars. For larger 3D catalogs and simulations it can be applied to results obtained with the core-sampling analysis. It is shown that the interaction of large and small scale perturbations modulates the creation rate of early Zel'dovich pancakes and generates bias on the SLSS scale. For suitable parameters of the cosmological model and reheating process this bias can essentially improve the characteristics of simulated structure of the universe. The models with $0.3\leq \Omega_m \leq 0.5$ give the best description of the observed structure parameters. The influence of low mass "warm" dark matter particles, such as a massive neutrino, will extend the acceptable range of $\Omega_m$ and $h$.

Geometrical and statistical properties of polarization of CMB are analyzed. Singular points of the vector field which describes CMB polarization are found and classified. Statistical distribution of the singularities is studied. A possible signature of tensor perturbations in CMB polarization is discussed. For a further analysis of CMB statistics Minkowski functionals are used, which present a technically simple method to search for deviations from a Gaussian distribution.

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.

Analysis of the singularities of the polarization field of CMB, where polarization is equal to zero, is presented. It is found that the classification of the singular points differs from the usual three types known in the ordinary differential equations. The new statistical properties of polarization field are discussed, and new methods to detect the presence of primordial tensor perturbations are indicated.

The large scale gravitational potential distribution and its influence on the large-scale matter clustering is considered on the basis of six simulations. It is found that the mean separation between zero levels of the potential along random straight lines coincides with the theoretical expectations, but it scatters largely. A strong link of the initial potential and the structure evolution is shown. It is found that the under-dense and over-dense regions correlate with regions of positive and negative gravitational potential at large redshifts. The over-dense regions arise due to a slow matter flow into the negative potential regions, where more pronounced non-linear structures appear. Such regions are related to the formation of huge super-large scale structures seen in the galaxy distribution.

The large-scale matter distribution represents a complex network of structure elements such as voids, clusters, filaments, and sheets. This network is spanned by a point distribution. The global properties of the point process can be measured by different statistical methods, which, however, do not describe directly the structure elements. The morphology of structure elements is an important property of the point distribution. Here we apply the core-sampling method to various Voronoi tessellations. Using the core-sampling method we identify one- and two-dimensional structure elements (filaments and sheets) in these Voronoi tessellations and reconstruct their mean separation along random straight lines. We compare the results of the core-sampling method with the a priori known structure elements of the Voronoi tessellations under consideration and find good agreement between the expected and found structure parameters, even in the presence of substantial noise. We conclude that the core-sampling method is a potentially powerful tool to investigate the distribution of such structure elements like filaments and walls of galaxies.

Geometrical model of structure of the universe is examined to obtain analytical expression for the two points nonlinear correlation function. According to the model the objects (galaxies) are concentrated into two types of structure elements - filaments and sheets. We considered the filaments ( similar to galaxy filaments ) simply as straight lines and the sheets ( similar to superclusters of galaxies ) simply as planes. The homogeneously distributed objects are also taken into consideration. The spatial distribution of lines, planes and points is uncorrelated. The nonlinear correlation function depends on four parameters and is similar to the observed and simulated ones for different samples. It describes quite well the correlation of galaxies, clusters of galaxies and dark matter distribution. Possible interpretation of the parameters of nonlinear correlation function is discussed.

The CMB anisotropy on arc second range is examined to test the power spectrum of perturbations in the small scale region and, in particular, to estimate the mass of dominant dark matter particles. It is shown that for the simplest evolutionary history with standard recombination, three and four beam observations could discriminate the mass of dark matter particles in the interval $0.5 KeV\leq M_{DM}\leq 4 KeV$ with an antenna beam (0.5 - 0.25) arc minute with amplitude $\approx 10^{-7}$.

A quantitative theory, based on the Zeldovich approximation, to provide an approximate description of the evolution of structure is presented. We give an expression for the characteristic scale of superlarge-scale structure which can also be applied to power spectra with a Harrison-Zeldovich asymptotic. The evolution of the characteristic scale of large-scale structure is described with an analytical expression. It corresponds to the mean separation of filaments or chains of galaxies. We apply the theory to some simulations and show that the theoretically defined characteristic scale for LSS does also correspond to the 'physical measurement' of the mean separation of dark matter LSS elements. This supports our claim for a new theory for the quantitative description of the formation and evolution of large-scale structure in the universe.