We consider the black hole (BH) shadow images which can be restored by data processing and image recovery procedures in future space Very Large Baseline Interferometry (VLBI) missions. For Kerr BHs with masses and coordinates of SgrA*, M87* and M31*, illuminated by light source behind them, we consider three kinds of observation: the ground-based interferometer (similar to the Event Horizon Telescope), space-ground interferometer with a satellite in geocentric orbit and space-ground interferometer with a satellite located in Lagrange point $L_2$. The significant difference between the images produced by the ground-based telescope alone and the one of the space VLBI with an added low-orbit satellite is caused by both the increased baseline and the improved of (u,v) coverage. The near-Earth configuration of the radio interferometer for the observation of BH shadows is the most preferable for the BH shadow observations among considered cases. With further increasing the orbit radius up to the Lagrange point $L_2$ the density of the (u,v) filling decreases and the results appear less reliable. Model images for all the cases are presented.

We investigate a possibility to find an accreting isolated black hole (IBH) with mass $1-100~\mathrm{M}_{\odot}$ within Central Galactic Molecular Zone (CMZ) in the submillimetre and IR spectral range with help of planned space observatories James Webb Space Telescope (JWST) and Millimetron (MM). We assume the spherical mode of accretion. We develop the simplest possible quantitative model of the formation of radiation spectrum in this range due to synchrotron emission and show that it fully agrees with the more complicated model of Ipser and Price 1982 for expected values of accretion rate. If a substantial fraction of LIGO events was caused by merger of primordial black holes, the JWST would be able to find them provided that there is a cusp in distribution of dark matter in the Galaxy and that the accretion efficiency parameter $\lambda $ defined as the ratio of accretion rate onto IBH to its Bondi-Hoyle-Lyttleton value is larger than $\sim 10^{-2}$. A comparable amount of IBHs is also predicted by recent models of their formation due to stellar evolution. MM capabilities are hampered by the effect of confusion due to distant submillimetre galaxies, but it can also be used for such purposes if the confusion effect is properly dealt with. In case of efficient accretion with $\lambda \sim 1$, both instruments could detect IBHs even when their number densities are as low as $10^{-6}~\mathrm{pc}^{-3}$.

Sensitivity of future far infrared space telescopes like Millimetron will be limited by a confusion noise created by distant galaxies. We construct a model of the Cosmic Infrared Background (CIB) aimed at exploration of methods of prediction and reducing the confusion noise. The model is based on a public available eGALICS simulation. For each simulated galaxy we construct a spectral energy distribution with the help of public GRASIL and CHE_EVO codes. In this paper, in order to put our model in the context of current CIB investigations, we compare the outputs of the model: luminosity and mass functions, source counts as a function of flux and redshift, spectrum of the CIB, prediction of confusion limit, with the available observational data and with three other models: one is a well known `backwards evolution' model of Bethermin et al. 2011 and two others are based on a simple mass-luminosity (M-L) relation applied to simulated dark matter halo catalogues. We conclude that our model reproduces the observational data reasonably well. All four models show significant differences in the predictions of the distribution of sources on the flux-redshift plane, especially at high redhifts. The predicted confusion noise on the wavelengths 70--350 microns is consistent between models, while for 650--2000 microns there are significant differences and M-L models should not be used.

The paper presents a catalog of supermassive black holes (SMBH) shapened for the interferometric observations in millimeter and submillimeter wavelength ranges and based on the open sources. The catalog includes the name of the object, coordinates, angular distance, the mass, the angular size of the gravitational radius of SMBH, the integral flux of radiosource, related with SMBH, in the range $20\div 900$~GHz which is considered to be used in Event Horizon Telescope, future space mission Millimetron and others. The catalog is intended for use during the planning of the interferometric observations of SMBH shadows.

The extragalactic background radiation produced by distant galaxies emitting in the far infrared limits the sensitivity of telescopes operating in this range due to confusion. We have constructed a model of the infrared background based on numerical simulations of the large-scale structure of the Universe and the evolution of dark matter halos. The predictions of this model agree well with the existing data on source counts. We have constructed maps of a sky field with an area of 1 deg$^2$ directly from our simulated observations and measured the confusion limit. At wavelengths $100-300$ $\mu$m the confusion limit for a 10-m telescope has been shown to be at least an order of magnitude lower than that for a 3.5-m one. A spectral analysis of the simulated infrared background maps clearly reveals the large-scale structure of the Universe. The two-dimensional power spectrum of these maps has turned out to be close to that measured by space observatories in the infrared. However, the fluctuations in the number of intensity peaks observed in the simulated field show no clear correlation with superclusters of galaxies; the large-scale structure has virtually no effect on the confusion limit.

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.

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.

We use the phenomenological approach to study properties of space-time in the vicinity of the Schwarzschild black-hole singularity. Requiring finiteness of the Schwarzschild-like metrics we come to the notion of integrable singularity that is, in a sense, weaker than the conventional singularity and allows the (effective) matter to pass to the white-hole region. This leads to a possibility of generating a new universe there. Thanks to the gravitational field of the singularity, this universe is already born highly inflated ('singularity-induced inflation') before the ordinary inflation starts.

We discuss status of the singularity problem in General Relativity and argue that the requirement that a physical solution must be completely free of singularities may be too strong. As an example, we consider properties of the integrable singularities and show that they represent light horizons separating T-regions of black and white holes. Connecting an astrophysical black hole to a white hole, they lead to a natural mechanism of generating new universes. Under favorable conditions the new universes will also contain black holes which, in their turn, will give rise to another generation of universes. In this case the cosmological evolutionary tree will continue to grow to form the "hyperverse". This scenario essentially differs from other known mechanisms, such as bounce, birth from "nothing", baby-universe scenario, etc.

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.

In this review, the formation, evolution, and decay of the large-scale structure of the Universe is discussed in the context of observational data, numerical simulations, and the Cosmological Standard Model (CSM). Problems concerning measuring and interpreting cosmological parameters, determining the composition of matter, and normalizing density perturbation spectra are especially highlighted.

We review different approaches to the problem of generating the observed Hubble flow, whose discussion was pioneered by A.D. Sakharov. Extrapolating the Cosmological Standard Model to the past makes it possible to determine physical properties of and conditions in the early universe. We discuss a new cosmogenesis paradigm based on studying geodesically complete geometries of black/white holes with integrable singularities.

We briefly review the problem of generating cosmological flows of matter in GR (the genesis of universes), analyze models' shortcomings and their basic assumptions yet to be justified in physical cosmology. We propose a paradigm of cosmogenesis based on the class of spherically symmetric solutions with {\it integrable} singularity $r=0$. They allow for geodesically complete geometries of black/white holes, which may comprise space-time regions with properties of cosmological flows.

Since the discovery of quasars in papers often appeared and appear the assertions that the redshift quasar distribution includes a periodic component with the period $\Delta z = 0.063$ or 0.11. A statement of such kind, if it is correct, may manifest the existence of a far order in quasar distribution in cosmological time, that might lead to a fundamental revision all the cosmological paradigm. In the present time there is a unique opportunity to check this statement with a high precision, using the rich statictics of 2dF and SDSS catalogues (about 85000 quasars). Our analysis indicates that the periodic component in distribution of quasar redshifts is absent at high confidence level.

We discuss the questions related to dark energy in the Universe. We note that in spite of the effect of dark energy, large-scale structure is still being generated in the Universe and this will continue for about ten billion years. We also comment on some statements in the paper ``Dark energy and universal antigravitation'' by A.D. Chernin [4].

Observational cosmology is on the verge of new discoveries that will change the essence of our world-view. The matter concerns origin of initial conditions and physics of dark matter.

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.

In spectra of many Seyfert galaxies there is a wide emission line of Fe $K_\alpha$. The line profile with two maxima supposes that the line emerges in innermost regions of an accretion disk around a black hole, hence, it is necessary to take into account General Relativity (GR) effects. In order to determine GR processes which occur in active galactic nuclei (AGN) an inverse problem of reconstructing the accreting system parameters from the line profile has to be solved quickly. In this paper we present a numerical approximation of the emission line Fe $K_\alpha$ with analytical functions. The approximation is accomplished for a range of the disk radial coordinate $r$ and the angle $\theta$ between line of sight and perpendicular to the disk and allows one to decrease computing time by $10^6$ times in certain astrophysical problems taking into account all GR effects. The approximation results are available in the Internet at http://www.iki.rssi.ru/people/repin/approx

An elementary derivation of the fundamental relation $T/S=4\gamma$ between the tensor and scalar modes of cosmological perturbations in the early universe is given. Statements by L.P.Grishchuk on this problem are commented on.

We modify the Press-Schechter formalism to calculate the mass function of gravitationally bounded objects in a local Universe with size $L> 10h^{-1}$Mpc, i.e. separately in voids ($\Omega_{m,L}<\Omega_m$) and superclusters ($\Omega_{m,L}>\Omega_m$). The dependence of the abundance of gravitationally bounded objects on local matter density is analyzed.

Observations of AGNs and microquasars by ASCA, RXTE, Chandra and XMM-Newton indicate the existence of wide X-ray emission lines of heavy ionized elements in their spectra. The emission can arise in the inner parts of accretion discs where the effects of General Relativity (GR) must be counted, moreover such effects can dominate. We describe a procedure to estimate an upper limit of the magnetic fields in the regions where X-ray photons are emitted. We simulate typical profiles of the iron $K_\alpha$ line in the presence of a magnetic field and compare them with observational data. As an illustration we find $H < 10^{10} - 10^{11}$ Gs for Seyfert galaxy MCG--6--30--15. Using the perspective facilities of measurement devices (e.g. Constellation-X mission) a better resolution of the blue peak structure of iron $K_\alpha$ line will allow to find the value of the magnetic fields if the latter are high enough.

The time evolution of the number density of galaxy clusters and their mass and temperature functions are used to constrain cosmological parameters in the spatially flat dark matter models containing a fraction of hot particles (massive neutrino) additional to cold and baryonic matter. We test the modified MDM models with cosmic gravitational waves and show that they neither pass the cluster evolution test nor reproduce the observed height of the first acoustic peak in $\Delta T/T$ spectrum, and therefore should be ruled out. The models with a non-zero cosmological constant are in better agreement with observations. We estimate the free cosmological parameters in $\Lambda$MDM with a negligible abundance of gravitational waves, and find that within the parameter ranges $h\in (0.6, 0.7)$, $n\in (0.9, 1.1)$, (i) the value of $\Omega_\Lambda$ is strongly affected by a small fraction of hot dark matter, $f_\nu\equiv\Omega_\nu /\Omega_m\in (0, 0.2)$: $0.45 <\Omega_\Lambda <0.7$ ($1\sigma$ CL), and (ii) the redshift evolution of galaxy clusters alone reveals the following explicit correlation between $\Omega_\Lambda$ and $f_\nu$: $\Omega_\Lambda +0.5f_\nu =0.65\pm 0.1$. The present accuracy of observational data allows only to bound the fraction of hot matter, $f_\nu\in (0, 0.2)$ (the number of massive neutrino species remains undelimited, $N_\nu =1, 2, 3$).

The galaxy cluster power spectrum and mass/temperature functions are currently the most precise observational tools for constraining the theory of the formation of large scale structure (LSS) in the Universe. Complementing these tests by the observational data at larger (cosmic microwave backgroud anisortopy (CMBA)) and smaller (distribution of $Ly_\alpha$ clouds) scales opens the way to a straightforward determination of the cosmological parameters in simplest dark matter (DM) models. We argue that such a 'minimal data set' is free from systematic effects and can indicate quite precisely the parameters of spatially flat mixed DM model with a positive cosmological constant and no cosmic gravitational waves.

The time evolution of galaxy cluster abundance is used to constrain cosmological parameters in dark matter models containing a fraction of hot particles (massive neutrino). We test the modified MDM models with cosmic gravitational waves which are in agreement with observational data at $z=0$, and show that they do not pass the cluster evolution test and therefore should be ruled out. The models with a non-zero cosmological constant are in better agreement with the evolution test. We estimate $\Omega_\Lambda$ and find that it is strongly affected by a small fraction of hot dark matter: $0.4 <\Omega_\Lambda <0.8$ for $\Omega_H /\Omega_M <0.2$.

The current status of the theoretical and observational cosmology is reviewed.

We argue that a $\Lambda$-inflation model can ensure large relative contribution of cosmic gravity waves into the $\Delta T/T$ at COBE scale preserving at the same time a near scale-invariant spectrum of cosmological density perturbations favored by observational data ($n_S\simeq 1$). High efficiency of these models to meet observational tests is discussed.

We propose a Lambda-inflation model which explains a significant part of the COBE signal by primordial cosmic gravitational waves. The primordial density perturbations fulfill both the constraints of large-scale microwave background and galaxy cluster normalization. The model is tested against the galaxy cluster power spectrum and the high-multipole angular CMB anisotropy.

An analysis of cosmological models in spatially flat Friedmann Universe with cosmic gravitational wave background and zero $\Lambda$-term is presented. The number of free parameters is equal to 5, they are $\sigma_8$, $n$, $\Omega_\nu$, $\Omega_b$, and $h$. The normalization of the spectrum of density perturbations on galaxy cluster abundance ($\sigma_8 = 0.52\pm 0.04$) has been used to calculate numerically the value of the large scale CMB anisotropy ($\ell\simeq 10$) and the relative contribution of cosmological gravitational waves T/S. Increasing $\Omega_\nu$ weaken the requirements to the value of T/S, however even for $\Omega_\nu\le 0.4$ the models with $h+n\ge 1.5$ suggest considerable abundance of gravitational waves: T/S${}^>_\sim 0.3$. In models with $\Omega_\nu\le 0.4$ and scale-invariant spectrum of density perturbations ($n=1$): T/S${}^>_\sim 10(h-0.47)$. Minimization of the value T/S is possible only in the range of the red spectra ($n<1$) and small $h$ ($<0.6$). It is shown that the models with T/S$\in [0, 3]$ admit both moderate red and blue spectra of density perturbations, $n\in[0.9,1.2]$, with rather high abundance hot dark matter, $\Omega_\nu\in [0.2,0.4]$. Any condition, $n<0.9$ or $\Omega_\nu<0.2$, decreases the relative amplitude of the first acoustic peak for more than 30% in comparison with its hight in the standard CDM normalized by COBE data.

We propose a $\Lambda$-inflation model which explains a large fraction of the COBE signal by cosmic gravitational waves. The primordial density perturbations fulfil both the contraints of large-scale microwave background and galaxy cluster normalization. The model is tested against the galaxy cluster power spectrum and the high-multipole angular CMB anisotropy.

The possibility of determining cosmological parameters on the basis of a wide set of observational data including the Abell-ACO cluster power spectrum and mass function, peculiar velocities of galaxies, the distribution of Ly-$\alpha$ clouds and CMB temperature fluctuations is analyzed. Using a $\chi^2$ minimization method, assuming $\Omega_{\Lambda}+\Omega_m =1$ and no contribution from gravity waves, we found that a tilted $\Lambda$MDM model with one sort of massive neutrinos and the parameters $n\approx 1.12$, $\Omega_m\approx 0.4$ ($\Omega_{\Lambda}\approx 0.6$), $\Omega_{cdm}\approx 0.3$, $\Omega_{\nu}\approx 0.06$, $\Omega_b\approx 0.04$ and $h\approx 0.7$ matches observational data best. $\Omega_{\nu}$ is higher for more species of massive neutrinos, $\sim 0.1$ for two and $\sim 0.13$ for three species. $\Omega_m$ raises by $\sim 0.08$ and $\sim 0.15$ respectively. Varying only a subset of parameters and fixing the others shows also that the observational data set used here rules out pure CDM models with $h\ge 0.5$, scale invariant primordial power spectrum, zero cosmological constant and spatial curvature at a very high confidence level, $>99.99%$. The corresponding class of MDM models are ruled out at $\sim 95%$ C.L. It is notable also that this data set determines the amplitude of scalar fluctuations approximately at the same level as COBE four-year data. It indicates that a possible tensor component in the COBE data cannot be very substantial.

The possibility of determining cosmological parameters on the basis of a wide set of observational data including the Abell-ACO cluster power spectrum and mass function, peculiar velocities of galaxies, the distribution of Ly-$\alpha$ clouds and CMB temperature fluctuations is analyzed. Using a $\chi^2$ minimization method, assuming $\Omega_{\Lambda}+\Omega_{\rm{matter}} =1 $ and no contribution from gravity waves, we found that a tilted $\Lambda$MDM model with one sort of massive neutrinos and the parameters $n=1.12\pm 0.10$, $\Omega_m=0.41\pm 0.11$ ($\Omega_{\Lambda}=0.59\pm0.11$), $\Omega_{cdm}=0.31\pm 0.15$, $\Omega_{\nu}=0.059\pm 0.028$, $\Omega_b=0.039\pm 0.014$ and $h=0.70\pm 0.12$ matches observational data best. The 1$\sigma$ (68.3%) confidence limits on each cosmological parameter, which are obtained by marginalizing over the other parameters, are $0.82\le n\le1.39$, $0.19\le\Omega_m\le 1$ ($0\le\Omega_{\Lambda}\le 0.81$), $0\le\Omega_{\nu}\le 0.17$, $0.021\le \Omega_b\le 0.13$ and $0.38\le h\le 0.85$ $1.5\le b_{cl}\le 3.5$. Here $b_{cl}$ is the cluster bias parameter. The best-fit parameters for 31 models which are inside of $1\sigma$ range of the best model are presented. It is shown also that observational data set used here rules out the class of CDM models with $h\ge 0.5$, scale invariant primordial power spectrum, zero cosmological constant and spatial curvature at very high confidence level, $>99.99%$. The corresponding class of MDM models are ruled out at $\sim 95%$ C.L.

In these lectures we dwell upon the cosmological corner-stones of the Very Early Universe (VEU) theory: Parametric Amplification Effect (PAE) responsible for the generation of Primordial Cosmological Perturbations (PCPs), Chaotic and Stochastic Inflation, Principal Tests of VEU, and others.

We explore a broad class of three-parameter inflationary models, called the $\Lambda$-inflation, and its observational predictions: high abundance of cosmic gravitational waves consistent with the Harrison-Zel'dovich spectrum of primordial cosmological perturbations, the non-power-law wing-like spectrum of matter density perturbations, high efficiency of these models to meet current observational tests, and others. We show that a parity contribution of the gravitational waves and adiabatic density perturbations into the large-scale temperature anisotropy, T/S $\sim 1$, is a common feature of $\Lambda$-inflation; the maximum values of T/S (basically not larger than 10) are reached in models where (i) the local spectrum shape of density perturbations is flat or slightly red ($n_S{}_\sim^< 1$), and (ii) the residual potential energy of the inflaton is near the GUT scale ($V_0^{{1/4}} \sim 10^{16} GeV$). The conditions to find large T/S in the paradigm of cosmic inflation and the relationship of T/S to the ratio of the power spectra, $r$, and to the inflationary $\gamma$ and Hubble parameters, are discussed. We argue that a simple estimate, T/S$\simeq 3r\simeq 12\gamma \simeq (\frac{H}{6\times 10^{13}{\rm GeV}})^2$, is true for most known inflationary solutions and allows to relate straightforwardly the important parameters of observational and physical cosmology.

An accurate analytic approximation of the transfer function for the power spectra of primordial density perturbations in mixed dark matter models is presented. The fitting formula in a matter-dominated Universe ($\Omega_0=\Omega_M=1$) is a function of wavenumber $k$, redshift $z$ and four cosmological parameters: the density of massive neutrinos, $\Omega_{\nu}$, the number of massive neutrino species, $N_{\nu}$, the baryon density, $\Omega_{b}$ and the dimensionless Hubble constant, $h$. Our formula is accurate in a broad range of parameters: $k\le 100 h/Mpc$, $z\le 30$, $\Omega_{\nu}\le 0.5$, $N_{\nu}\le 3$, $\Omega_{b}\le 0.3$, $0.3\le h\le 0.7$. The deviation of the variance of density fluctuations calculated with our formula from numerical results obtained with CMBfast is less than 6% for the entire range of parameters. It increases with $\Omega_bh^{2}$ and is less than $\le 3%$ for $\Omega_bh^{2}\le 0.05$. The performance of the analytic approximation of MDM power spectra proposed here is compared with other approximations found in the literature. Our approximation turns out to be closest to numerical results in the parameter space considered here.

We present an analysis of cosmological models with mixed dark matter (spectrum slope of density perturbations $n_S \in [0.6, 1.6]$, hot dark matter abundance $\Omega_\nu \in [0.0, 0.6]$) in spatially flat Universe, based on Press-Schechter formalism (P\&S). We found that the models with $n_S > 1$ and $\Omega_\nu < 0.5$ are preferable. Additionally, for all considered models we simulated $\Delta T/T\vert_{10^0}$ and compared it with the COBE measured amplitude of CMB anisotropy. We found that the models favourable for galaxy cluster abundance need large amount of cosmological gravitational waves (T/S $> 1$).

In this brief review I will touch upon the main points and thends of observational and theoretical cosmology, that influence and change today our understanding of the Universe.

We argue that gravitational wave contribution to the cosmic microwave background anisotropy at angular scale $\sim 10^0$ may exceed 50% for some models of hybrid inflation producing standard cosmology with the density perturbation slope $n \simeq 1$.

Recently, it was suggested that large concentrations or groups of quasars may trace sites of enhanced matter density at medium and high redshifts analogous to how galaxy clusters trace them in nearby space. We have checked existing quasar data for the presence of such groups. Large Quasar Groups (LQGs) were identified using a well-known cluster analysis technique. Our final list consists of 12 groups, including one reported previously. Most of the quasars in these groups come from deep homogeneous surveys. We have found that: (i) the probability that the detected groups are random is rather small (generally a few per cent); (ii) their sizes range from $\sim 70$ to $\sim 160h^{-1} Mpc$, which is comparable to sizes of nearby rich superclusters; (iii) the detected groups all have redshifts $0.5<z\leq 2$; (iv) the abundance of the LQGs is comparable with the abundance of large superclusters at $z\sim 0$, which is consistent with the idea that quasar groups and superclusters can be evolutionarily related. We argue that quasar groups could be a common feature of the spatial distribution of medium redshift quasars, and that the quasars in groups may belong to concentrations of young galaxy clusters and groups (distant superclusters) and hence be biased tracers of large-scale structure of matter distribution in the early Universe. Theoretical implications as well as other observations needed to test this point are discussed.

The majority of bright distant quasars $(z {}^{>}_{\sim} 1)$ may form in massive mergers appearing in compact galaxy groups in/and young clusters. The expected tests are (i) large correlation signal for medium-$z$ QSOs $(1 < z < 2)$ and (ii) direct search for quasar groups (QGs) indicating positions of distant pre-superclusters which later will evolve to the "systems" like the local Great Attractor or Shapley concentration. We discuss large QGs with more than ten members within regions $\le l_{LS} \sim 100-150 \, h^{-1} Mpc$, tracing the enhanced density regions at $z {}^{<}_{\sim} 2$. These early large scale structures (i) provide a natural way to "bias" the distribution of Abell clusters, and (ii) suggest that the spectrum of primordial density perturbations is nearly flat at scales encompassing both cluster and GAs, $l = \pi k^{-1} \in(10,100)h^{-1} Mpc: \Delta^{2}_{k} \sim k^{3}P(k) \sim k^{\gamma}, \gamma = 1^{+0.6}_{-0.4}$.

We report the investigation of spatial distribution of quasars using two different methods: the statistical study by means of 2-point correlation function and direct search for structures such as previously reported Large Quasar Groups. Having analyzed the combined sample of eight homogeneous quasar surveys (totalling near 2200 objects), we detected strong clustering at scales $r\leq 20h^{-1}$ Mpc ($H_0=100h$ and $q_0=0.5$) which seems, however, to be redshift dependent. Using reasonable selection criteria we have also sought for Large Quasar Groups (LQG) similar to those which were earlier found observationally. We discuss the properties of the detected groups and their implications for cosmology. We also discuss the suitability of the present data for getting the reliable statistical inference on quasar clustering.