The accumulation of dark matter particles near the primordial black holes starts at the radiation-dominated cosmological stage and produces the central density spikes. The spikes can be the bright gamma-ray sources due to dark matter annihilation. We present the self-consistent derivation of the equation of motion of particle in the metrics of primordial black hole immersed into cosmological background. By numerical solution of this equation we find the central dark matter density profile. The density growth is suppressed in the central part of the profile compared with previous calculations.

A change in gravity law in some regimes is predicted in the modified gravity models that are actively discussed at present. In this paper, we consider a possibility that the signal recorded by the Geograv resonant gravitational-wave detector in 1987 during the explosion of SN 1987A was produced by an abrupt change in the metric during the passage of a strong neutrino flux through the detector. Such an impact on the detector is possible, in particular, in extended scalar-tensor theories in which the local matter density gradient affects the gravitational force. The first short neutrino pulse emitted at the initial stage of stellar core collapse before the onset of neutrino opacity could exert a major influence on the detector by exiting the detector response at the main resonance frequency. In contrast, the influence of the subsequent broad pulse (with a duration of several seconds) in the resonant detector is exponentially suppressed, despite the fact that the second pulse carries an order-of-magnitude more neutrino energy, and it could generate a signal in the LSD neutrino detector. This explains the time delay of 1.4s between the Geograv and LSD signals. The consequences of this effect of modified gravity for LIGO/Virgo observations are discussed.

The Primordial Black Holes (PBHs) are gradually involved into consideration as the phenomenon having reliable basis. We discuss here the possibility of their agglomeration into clusters that may have several prominent observable features. The clusters can form due to closed domain walls appearance in the natural and the hybrid inflation with subsequent evolution and gravitational collapse. Early dustlike stages of dominance of heavy metastable dissipative particles, at which star-like objects are formed, can also naturally lead to formation of black hole clusters, remaining in the Universe after decay of particles, from which they have originated. The dynamical evolution of such clusters discussed here is of the crucial importance. Such a model inherits all the advantages of the single PBHs like possible explanation of existence of supermassive black holes (origin of the early quasars), binary BH merges registered by LIGO/Virgo through gravitational waves, contribution to reionization of the Universe, but also has additional benefits. The cluster could alleviate or completely avoid existing constraints on the single PBH abundance making PBHs a real dark matter candidate. The most of existing constraints on (single) PBH density should be re-considered as applied to the clusters. Also unidentified cosmic gamma-ray point-like sources could be (partially) accounted for by them. One can conclude, that it seems really to be much more viable model with respect to the single PBHs.

We study a black hole in an expanding Universe during the radiation-dominated stage. In the case when the black hole radius is much smaller than the cosmological horizon, we present a solution of the Einstein equations for the metric, the matter density and velocity distributions. At distances much smaller than the cosmological horizon the solution features quasi-stationary accretion, while at large distances the solution asymptotes the homogeneous radiation dominated FLRW. We discuss how our solution is related to the McVittie solution. The obtained results can be applied, in particular, for the formation of dark matter density spikes around primordial black holes, and for the evolution of dark matter clumps during the radiation-dominated stage.

Previously, it has been established that axion dark matter (DM) is clustered to form clumps (axion miniclusters) with masses $M\sim10^{-12}M_\odot$. The passages of such clumps through the Earth are very rare events occurring once in $10^5$ years. It has also been shown that the Earth's passage through DM streams, which are the remnants of clumps destroyed by tidal gravitational forces from Galactic stars, is a much more probable event occurring once in several years. In this paper we have performed details calculations of the destruction of miniclusters by taking into account their distribution in orbits in the Galactic halo. We have investigated two DM halo models, the Navarro-Frenk-White and isothermal density profiles. Apart from the Galactic disk stars, we have also taken into account the halo and bulge stars. We show that about 2-5% of the axion miniclusters are destroyed when passing near stars and transform into axion streams, while the clump destruction efficiency depends on the DM halo model. The expected detection rate of streams with an overdensity exceeding an order of magnitude is 1-2 in 20 years. The possibility of detecting streams by their tidal gravitational effect on gravitational-wave interferometers is also considered.

It was shown by (Nakamura et al. 1997), (Ioka et al. 1998), and (Sasaki et al. 2016) that primordial black holes (PBHs) binaries can form effectively at the cosmological stage of radiation dominance, and the merge of the PBHs in pairs can explain the gravitational wave burst GW150914. In this paper, the model is re-examined by considering the effect of inflationary dark matter density perturbations which produce additional tidal forces. As a result, the merge rate of PBHs binaries and the corresponding rate of the gravitational bursts are suppressed by the factor $\sim1.5-2$ in comparison with previous calculations. This rate matches the LIGO data if the PBHs constitute the $f\sim5\times10^{-4}-5\times10^{-3}$ fraction of dark matter.

We propose the model describing the observed multiple fast radio bursts due to the close encounters and collisions of neutron stars in the central clusters of the evolved galactic nuclei. The subsystem of neutron star cluster may originate in the dense galactic nucleus evolutionary in the combined processes of stellar and dynamical evolution. The neutron stars in the compact cluster can produce the short living binaries with the highly eccentric orbits, and finally collide after several orbital revolutions. Fast radio bursts may be produced during the close periastron approach and at the process of the final binary merging. In the sufficiently dense star cluster the neutron stars collisions can be very frequent. Therefore, this model can explain in principle the observed recurrent (multiple) fast radio bursts, analogous to the observed ones from the source FRB 121102. Among the possible observational signatures of the proposed model may be the registration of the gravitational wave bursts by the laser interferometers LIGO/VIRGO or by the next generation of gravitational wave detectors.

We show that density spikes begin to form from dark matter particles around primordial black holes immediately after their formation at the radiation-dominated cosmological stage. This follows from the fact that in the thermal velocity distribution of particles there are particles with low velocities that remain in finite orbits around black holes and are not involved in the cosmological expansion. The accumulation of such particles near black holes gives rise to density spikes. These spikes are considerably denser than those that are formed later by the mechanism of secondary accretion. The density spikes must be bright gamma-ray sources. Comparison of the calculated signal from particle annihilation with the Fermi-LAT data constrains the present-day cosmological density parameter for primordial black holes with masses $M_{\rm BH}\geq10^{-8}M_\odot$ from above by values from $\Omega_{\rm BH}\leq1$ to $\Omega_{\rm BH}\leq10^{-8}$, depending on $M_{\rm BH}$. These constraints are several orders of magnitude more stringent than other known constraints.

The classical Vaidya metric is transformed to the special diagonal coordinates in the case of the linear mass function allowing rather easy treatment. We find the exact analytical expressions for metric functions in these diagonal coordinates. Using these coordinates, we elaborate the maximum analytic extension of the Vaidya metric with a linear growth of the black hole mass and construct the corresponding Carter-Penrose diagrams for different specific cases. The derived global geometry seemingly is valid also for a more general behavior of the black hole mass in the Vaidya metric.

We review the physical processes that occur at the center of the Galaxy and that are related to the supermassive black hole Sgr A* residing there. The discovery of high-velocity S0 stars orbiting Sgr A* for the first time allowed measuring the mass of this supermassive black hole, the closest one to us, with a 10\% accuracy, with the result $M_h=(4.1\pm0.4)\times 10^6M_\odot$. Further monitoring can potentially discover the Newtonian precession of the S0 star orbits in the gravitational field of the black hole due to invisible distributed matter. This will yield the "weight" of the elusive dark matter concentrated there and provide new information for the identification of dark matter particles. The weak accretion activity of the "dormant quasar" at the Galactic center occasionally shows up as quasiperiodic X-ray and near-IR oscillations with mean periods of $11$ and $19$ min. These oscillations can possibly be interpreted as related to the rotation frequency of the Sgr A* event horizon and to the latitude oscillations of hot plasma spots in the accretion disk. Both these frequencies depend only on the black hole gravitational field and not on the accretion model. Using this interpretation yields quite the accurate values for both the mass $M_h$ and the spin $a$ (Kerr rotation parameter) of Sgr A*: $M_h=(4.2\pm0.2)\times 10^6M_\odot$ and $a=0.65\pm0.05$.

A promising method for measuring the total mass of the dark matter near a supermassive black hole at the center of the Galaxy based on observations of nonrelativistic precession of the orbits of fast S0 stars together with constraints on the annihilation signal from the dark matter particles has been discussed. An analytical expression for the precession angle has been obtained under the assumption of a power-law profile of the dark matter density. In the near future, modern telescopes will be able to measure the precession of the orbits of S0 stars or to obtain a strong bound on it. The mass of the dark matter necessary for the explanation of the observed excess of gamma radiation owing to the annihilation of the dark matter particles has been calculated with allowance for the Sommerfeld effect.

We examine the compatibility of the mirror matter concept with the non-orientable wormholes. If any particle (or classical object) is traversing through the non-orientable wormhole, it turns into a corresponding mirror particle and vice versa. The critical sonic point in the hydrodynamical flow of the perfect fluid is located exactly at the throat of the wormhole. This means that wormholes are traversable for the fluid. The astrophysical signatures of the presence of non-orientable wormholes in the Universe are discussed. Some non-trivial aspects of the non-orientable wormholes electrodynamics are considered.

The new, rapidly developing field of theoretical research --- studies of dark energy interacting with black holes (and, in particular, accreting onto black holes) --- is reviewed. The term `dark energy' is meant to cover a wide range of field theory models, as well as perfect fluids with various equations of state, including cosmological dark energy. Various accretion models are analyzed in terms of the simplest test field approximation or by allowing back reaction on the black-hole metric. The behavior of various types of dark energy in the vicinity of Schwarzschild and electrically charged black holes is examined. Nontrivial effects due to the presence of dark energy in the black hole vicinity are discussed. In particular, a physical explanation is given of why the black hole mass decreases when phantom energy is being accreted, a process in which the basic energy conditions of the famous theorem of nondecreasing horizon area in classical black holes are violated. The theoretical possibility of a signal escaping from beneath the black hole event horizon is discussed for a number of dark energy models. Finally, the violation of the laws of thermodynamics by black holes in the presence of noncanonical fields is considered.

Small-scale clumps of dark matter are gravitationally bounded structures that have masses comparable to or lower than stellar masses and consist of noninteracting or weakly interacting dark matter particles. In this paper, the current knowledge about the formation and evolution of such structures is reviewed, various types of spectra of primordial cosmological perturbations are considered, and various dark matter models are discussed. Depending on the particular spectrum type, dark matter clumps may differ considerably in their formation processes and ultimate characteristics. The role of clumps in experiments on indirect detection of dark matter particles via their annihilation products is discussed. A number of astrophysical problems and phenomena that are related to dark matter clumps are examined.

We propose the new dark matter particle candidate - the "black hole atom", which is an atom with the charged black hole as an atomic nucleus and electrons in the bound internal quantum states. As a simplified model we consider the the central Reissner-Nordstrom black hole with the electric charge neutralized by the internal electrons in bound quantum states. For the external observers these objects would look like the electrically neutral Schwarzschild black holes. We suppose the prolific production of black hole atoms under specific conditions in the early universe.

We discuss the formation mechanisms and structure of the superdense dark matter clumps (SDMC) and ultracompact minihaloes (UCMH), outlining the differences between these types of DM objects. We define as SDMC the gravitationally bounded DM objects which have come into virial equilibrium at the radiation-dominated (RD) stage of the universe evolution. Such objects can be formed from the isocurvature (entropy) density perturbations or from the peaks in the spectrum of curvature (adiabatic) perturbation. The axion miniclusters (Kolb and Tkachev 1994) are the example of the former model. The system of central compact mass (e. g. in the form of SDMC or primordial black hole (PBH)) with the outer DM envelope formed in the process of secondary accretion we refer to as UCMH. Therefore, the SDMC can serve as the seed for the UCMH in some scenarios. Recently, the SDMC and UCMH were considered in the many works, and we try to systematize them here. We consider also the effect of asphericity of the initial density perturbation in the gravitational evolution, which decreases the SDMC amount and, as the result, suppresses the gamma-ray signal from DM annihilation.

We consider observational signatures of the collisions and partial destructions of giant exoplanets at the chaotic stage of the planetary systems evolution. The rate of these collisions in the Galaxy is estimated to be ~0.01-1 per yr. In the inelastic collision with a small relative velocity and small impact parameter, the planets may sufficiently lose their kinetic energy and merge together. Otherwise, if the planet is experienced a shallow tangential collision, it flews safety away. At the same time, the planets loss some part of their gas envelopes during mutual collisions. Collisions of the giant planets must be accompanied also by the radiation of energy from the radio to optical bands. The optical and near UV flashes result from the collision heating of the planet surface layers. Additionally, the collision compression and collapse of the powerful magnetosphere of giant planets are possible, with a successive generation of the radio bursts. According to our estimations, the corresponding rate the radio bursts is ~0.01-1 bursts per year with the maximum spectral flux ~30mJy at the frequencies ~3GHz and with a duration of the each burst ~1.5 hours. These signals are available for registration by the existing radio telescopes.

We discuss the multiband energy release in a model of a collapsing galactic nucleus, and we try to interpret the unique super-long cosmic gamma-ray event Swift J164449.3+573451 (GRB 110328A by early classification) in this scenario. Neutron stars and stellar-mass black holes can form evolutionary a compact self-gravitating subsystem in the galactic center. Collisions and merges of these stellar remnants during an avalanche contraction and collapse of the cluster core can produce powerful events in different bands due to several mechanisms. Collisions of neutron stars and stellar-mass black holes can generate gamma-ray bursts (GRBs) similar to the ordinary models of short GRB origin. The bright peaks during the first two days may also be a consequence of multiple matter supply (due to matter release in the collisions) and accretion onto the forming supermassive black hole. Numerous smaller peaks and later quasi-steady radiation can arise from gravitational lensing, late accretion of gas onto the supermassive black hole, and from particle acceleration by shock waves. Even if this model will not reproduce exactly all the Swift J164449.3+573451 properties in future observations, such collapses of galactic nuclei can be available for detection in other events.

We develop a model of production of the very dense clumps of DM in RD epoch due to the accretion of DM on the loops of cosmic strings as the seeds. At some time the loops disappear, for example due to the gravitational radiation, and the remaining dense clumps produce the enhancement of the annihilation signal. We take into account the velocity distribution of the strings, and consider the two extreme regimes of DM annihilation: fast decay and continuous evaporation. The produced annihilation flux of gamma radiation is detectable, and for some parameters of DM particles and the strings can exceed the extragalactic flux of the gamma-radiation observed by Fermi. For the fixed parameters of DM particles (e.g. neutralino with fixed masses and cross-section of annihilation) one can obtain the limits on the basic string parameter, tension \mu, which is stronger than (more general) limits obtained from WMAP observations, cosmological nucleosynthesis and gravitational lensing. In particular for the neutralino with 100 GeV mass we exclude the interval 5\times10^{-10}<G{\mu}/c^2<5.1\times10^{-9}.

We consider the possibility of observing supernova (SN) - gamma-ray burst (GRB) pairs separated by several years. These energetic events correspond to the second SN explosion in a double star system followed by the formation of the very short-lived pair of compact objects and finally by the GRB generation after their merge. Possible connection of SN 1997ds and GRB 990425 is discussed.

We describe spherically symmetric steady-state accretion of perfect fluid in the Reissner-Nordstrom metric. We present analytic solutions for accretion of a fluid with the linear equations of state and of the Chaplygin gas. It is also shown that, under reasonable physical conditions, there is no steady-state accretion of a perfect fluid onto a Reissner-Nordstrom naked singularity. Instead, a static atmosphere of fluid is formed. We discuss a possibility of violation of the third law of black hole thermodynamics for a phantom fluid accretion.

The probability of observing "supernova - gamma-ray burst" (GRB) pair events and recurrent GRBs from one galaxy in a time interval of several years has been estimated. Supernova explosions in binary systems accompanied by the formation of a short-lived pair of compact objects can be the sources of such events. If a short GRB is generated during the collision of a pair, then approximately each of ~300 short GRBs with redshift z must have an optical precursor - a supernova in the observer's time interval <2(1+z)yr. If the supernova explosion has the pattern of a hypernova, then a successive observation of long and short GRBs is possible. The scenario for the generation of multiple GRBs in collapsing galactic nuclei is also discussed.

A model for the formation of supermassive primordial black holes in galactic nuclei with the simultaneous suppression of the formation of intermediate-mass black holes is presented. A bimodal mass function for black holes formed through phase transitions in a model with a "Mexican hat" potential has been found. The classical motion of the phase of a complex scalar field during inflation has been taken into account. Possible observational manifestations of primordial black holes in galaxies and constraints on their number are discussed.

Superheavy dark matter (SHDM) exchanges energy with its environment much slower than particles with masses close to the electroweak (EW) scale and has therefore different small-scale clustering properties. Using the neutralino as candidate for the SHDM, we find that free-streaming allows the formation of DM clumps of all masses down to $\sim 260 m_\chi$ in the case of bino. If small-scale clumps evolve from a non-standard, spiky spectrum of perturbations, DM clumps may form during the radiation dominated era. These clumps are not destroyed by tidal interactions and can be extremely dense. In the case of a bino, a "gravithermal catastrophe" can develop in the central part of the most dense clumps, increasing further the central density and thus the annihilation signal. In the case of a higgsino, the annihilation signal is enhanced by the Sommerfeld effect. As a result annihilations of superheavy neutralinos in dense clumps may lead to observable fluxes of annihilation products in the form of ultrahigh energy particles, for both cases, higgsinos and binos, as lightest supersymmetric particles.

The formation and evolution of superdense clumps (or subhalos) is studied. Such clumps of dark matter (DM) can be produced by many mechanisms, most notably by spiky features in the spectrum of inflationary perturbations and by cosmological phase transitions. Being produced very early during the radiation dominated epoch, superdense clumps evolve as isolated objects. They do not belong to hierarchical structures for a long time after production, and therefore they are not destroyed by tidal interactions during the formation of larger structures. For DM particles with masses close to the electroweak (EW) mass scale, superdense clumps evolve towards a power-law density profile $\rho(r) \propto r^{-1.8}$ with a central core. Superdense clumps cannot be composed of standard neutralinos, since their annihilations would overproduce the diffuse gamma radiation. If the clumps are constituted of superheavy DM particles and develop a sufficiently large central density, the evolution of their central part can lead to a 'gravithermal catastrophe.' In such a case, the initial density profile turns into an isothermal profile with $\rho \propto r^{-2}$ and a new, much smaller core in the center. Superdense clumps can be bserved by gamma radiation from DM annihilations and by gravitational wave detectors, while the production of primordial black holes and cascade nucleosynthesis constrain this scenario.

The rate of gravitational wave bursts from the mergers of massive primordial black holes in clusters is calculated. Such clusters of black holes can be formed through phase transitions in the early Universe. The central black holes in clusters can serve as the seeds of supermassive black holes in galactic nuclei. The expected burst detection rate by the LISA gravitational wave detector is estimated.

An analytic solution for the accretion of ultra-hard perfect fluid onto a moving Kerr-Newman black hole is found. This solution is a generalization of the previously known solution by Petrich, Shapiro and Teukolsky for a Kerr black hole. Our solution is not applicable for an extreme black hole due to violation of the test fluid approximation. We also present a stationary solution for a massless scalar field in the metric of a Kerr-Newman naked singularity.

What happened to the central cores of tidally destructed dark matter clumps in the Galactic halo? We calculate the probability of surviving of the remnants of dark matter clumps in the Galaxy by modelling the tidal destruction of the small-scale clumps. It is demonstrated that a substantial fraction of clump remnants may survive through the tidal destruction during the lifetime of the Galaxy if the radius of a core is rather small. The resulting mass spectrum of survived clumps is extended down to the mass of the core of the cosmologically produced clumps with a minimal mass. Since the annihilation signal is dominated by the dense part of the core, destruction of the outer part of the clump affects the annihilation rate relatively weakly and the survived dense remnants of tidally destructed clumps provide a large contribution to the annihilation signal in the Galaxy. The uncertainties in minimal clump mass resulting from the uncertainties in neutralino models are discussed.

Model of supermassive black holes formation inside the clusters of primordial black holes is developed. Namely, it is supposed, that some mass fraction of the universe ~10^-3 is composed of the compact clusters of primordial (relic) black holes, produced during phase transitions in the early universe. These clusters are the centers of dark matter condensation. We model the formation of protogalaxies with masses about 2*10^8M_sun at the redshift z=15. These induced protogalaxies contain central black holes with mass ~10^5M_sun and look like dwarf spheroidal galaxies with central density spike. The subsequent merging of induced protogalaxies and ordinary dark matter haloes corresponds to the standard hierarchical clustering scenario of large-scale structure formation. The coalescence of primordial black holes results in formation of supermassive black holes in the galactic centers. As a result, the observed correlation between the masses of central black holes and velocity dispersion in the galactic bulges is reproduced.

The origin of supermassive black holes in the galactic nuclei is quite uncertain in spite of extensive set of observational data. We review the known scenarios of galactic and cosmological formation of supermassive black holes. The common drawback of galactic scenarios is a lack of time and shortage of matter supply for building the supermassive black holes in all galaxies by means of accretion and merging. The cosmological scenarios are only fragmentarily developed but propose and pretend to an universal formation mechanism for supermassive black holes.

We describe the anisotropy of dark matter clump distribution caused by tidal destruction of clumps in the Galactic disk. A tidal destruction of clumps with orbit planes near the disk plane occurs more efficiently as compared with destruction of clumps at near-polar orbits. A corresponding annihilation of dark matter particles in small-scale clumps produces the anisotropic gamma-ray signal with respect to the Galactic disk. This anisotropy is rather small, 9%, and superimposed on that due to off-centering position of the Sun in the Galaxy. The anisotropy of annihilation signal with respect to the Galactic disk provides a possibility to discriminate dark matter annihilation from the diffuse gamma-ray backgrounds of other origin.

A mass function of small-scale dark matter clumps is calculated in the standard cosmological scenario with an inflationary-produced primordial fluctuation spectrum and with a hierarchical clustering. We take into account the tidal destruction of clumps at early stages of structure formation starting from a time of clump detachment from the Universe expansion. Only a small fraction of these clumps, ~0.1-0.5%, in each logarithmic mass interval $\Delta\log M\sim1$ survives the stage of hierarchical clustering. The surviving clumps can be disrupted further in the galaxies by tidal interactions with stars. We performed the detailed calculations of the tidal destruction of clumps by stars in the Galactic bulge and halo and by the Galactic disk itself. It is shown that the Galactic disc provides the dominant contribution to the tidal destruction of small-scale clumps outside the bulge. The results obtained are crucial for calculations of the dark matter annihilation signal in the Galaxy.

The stationary, spherically symmetric accretion of dark energy onto a Schwarzschild black hole is considered in terms of relativistic hydrodynamics. The approximation of an ideal fluid is used to model the dark energy. General expressions are derived for the accretion rate of an ideal fluid with an arbitrary equation of state p=p(\rho) onto a black hole. The black hole mass was found to decrease for the accretion of phantom energy. The accretion process is studied in detail for two dark energy models that admit an analytical solution: a model with a linear equation of state, p=\alpha(\rho-\rho_0), and a Chaplygin gas. For one of the special cases of a linear equation of state, an analytical expression is derived for the accretion rate of dark energy onto a moving and rotating black hole. The masses of all black holes are shown to approach zero in cosmological models with phantom energy in which the Big Rip scenario is realized.

We describe the model of protogalaxy formation around the cluster of primordial black holes with a minimum extension of standard cosmological model. Namely, it is supposed, that a mass fraction of the universe ~10^-3 is composed of the compact clusters of primordial (relict) black holes produced during the phase transitions in the early universe. These clusters are the centers of the dark matter (DM) condensations. As a result the protogalaxies with a mass 2x10^8Msun form at the redshift z=15. These induced protogalaxies contain the central black holes of mass ~10^5Msun and look like the dwarf spheroidal galaxies with a central density spike. Subsequent merging of the induced protogalaxies and ordinary DM haloes leads to the standard scenario of the large scale structure formation. Black holes merging gives the nowadays supermassive black holes and reproduces the observed correlations between their masses and velocity dispersions in the bulges.

We propose the model of first quasars formation around the cluster of rimordial black holes (PBHs). It is supposed, that mass fraction of the universe ~10^-3 is composed of the compact clusters of PBHs, produced during the phase transitions in the early universe. The clusters of PBHs become the centers of dark matter condensation. As a result, the galaxies with massive central black holes are formed. In the process of galaxies formation the central black holes are growing due to accretion. This accretion is acompaned by the early quasar activity.

We study the cosmological origin of small-scale DM clumps in the hierarchical scenario with the most conservative assumption of adiabatic Gaussian fluctuations. The mass spectrum of small-scale clumps with M<10^3Msun is calculated with tidal destruction of the clumps taken into account within the hierarchical model of clump structure. Only 0.1-0.5% of small clumps survive the stage of tidal destruction in each logarithmic mass interval. The mass distribution of clumps has a cutoff at Mmin due to diffusion of DM particles out of a fluctuation and free streaming at later stage. Mmin is a model dependent quantity. In the case the neutralino DM, considered as a pure bino, Mmin~10^-8 Msun. The evolution of density profile in a DM clump does not result in the singularity because of formation of the core under influence of tidal interaction. The radius of the core is ~0.1R, where R is radius of the clump. The applications for annihilation of DM particles in the Galactic halo are studied. The number density of clumps as a function of their mass, radius and distance to the Galactic center is presented. The enhancement of annihilation signal due to clumpiness, valid for arbitrary DM particles, is calculated. In spite of small survival probability, the global annihilation signal in most cases is dominated by clumps, with major contribution given by small clumps. The enhancement due to large clumps with M>10^6 Msun is very small.

Dark energy with the usually used equation of state $p=w\rho$, where $w=const<0$ is hydrodynamically unstable. To overcome this drawback we consider the cosmology of a perfect fluid with a linear equation of state of a more general form $p=\alpha(\rho-\rho_0)$, where the constants $\alpha$ and $\rho_0$ are free parameters. This non-homogeneous linear equation of state provides the description of both hydrodynamically stable ($\alpha>0$) and unstable ($\alpha<0$) fluids. In particular, the considered cosmological model describes the hydrodynamically stable dark (and phantom) energy. The possible types of cosmological scenarios in this model are determined and classified in terms of attractors and unstable points by the using of phase trajectories analysis. For the dark energy case there are possible some distinctive types of cosmological scenarios: (i) the universe with the de Sitter attractor at late times, (ii) the bouncing universe, (iii) the universe with the Big Rip and with the anti-Big Rip. In the framework of a linear equation of state the universe filled with an phantom energy, $w<-1$, may have either the de Sitter attractor or the Big Rip.

Solution for a stationary spherically symmetric accretion of the relativistic perfect fluid with an equation of state $p(\rho)$ onto the Schwarzschild black hole is presented. This solution is a generalization of Michel solution and applicable to the problem of dark energy accretion. It is shown that accretion of phantom energy is accompanied with the gradual decrease of the black hole mass. Masses of all black holes tend to zero in the phantom energy universe approaching to the Big Rip.

Production of small-scale DM clumps is studied in the standard cosmological scenario with an inflation-produced primeval fluctuation spectrum. Special attention is given to three following problems: (i) The mass spectrum of small-scale clumps with $M \lesssim 10^3 M_{\odot}$ is calculated with tidal destruction of the clumps taken into account within the hierarchical model of clump structure. Only 0.1 - 0.5% of small clumps survive the stage of tidal destruction in each logarithmic mass interval $\Delta\ln M\sim1$. (ii) The mass distribution of clumps has a cutoff at $M_{\rm min}$ due to diffusion of DM particles out of a fluctuation and free streaming at later stage. $M_{\rm min}$ is a model dependent quantity. In the case the neutralino, considered as a pure bino, is a DM particle, $M_{\rm min} \sim 10^{-8} M_{\odot}$. (iii) The evolution of density profile in a DM clump does not result in the singularity because of formation of the core under influence of tidal interaction. The radius of the core is $R_c \sim 0.1 R$, where $R$ is radius of the clump. The applications for annihilation of DM particles in the Galactic halo are studied. The number density of clumps as a function of their mass, radius and distance to the Galactic center is presented. The enhancement of annihilation signal due to clumpiness, valid for arbitrary DM particles, is calculated. In spite of small survival probability, the annihilation signal in most cases is dominated by clumps. For observationally preferable value of index or primeval fluctuation spectrum $n_p \approx 1$, the enhancement of annihilation signal is described by factor 2 - 5 for different density profiles in a clump.

We argue that the observed correlations between central black holes masses M_{BH} and galactic bulge velocity dispersions \sigma_e in the form M_{BH}\propto\sigma_e^4 may witness on the pregalactic origin of massive black holes. Primordial black holes would be the centers for growing protogalaxies which experienced multiple mergers with ordinary galaxies. This process is accompanied by the merging of black holes in the galactic nuclei.

To account for the microlensing events observed in the Galactic halo, Gurevich, Zybin, and Sirota have proposed a model of gravitationally bound, noncompact objects with masses of 0.01-1M_\odot. These objects are formed in the expanding Universe from adiabatic density perturbations and consist of weakly interacting particles of dark matter, for example, neutralinos. They assumed the perturbation spectrum on some small scale to have a distinct peak. We show that the existence of this peak would inevitably give rise to a large number of primordial black holes (PBHs) with masses of \sim10^5M_\odot at the radiation-dominated evolutionary stage of the Universe. Constraints on the coefficient of nonlinear contraction and on the compactness parameter of noncompact objects were derived from constraints on the PBH number density. We show that noncompact objects can serve as gravitational lenses only at a large PBH formation threshold, \delta_c > 0.5, or if noncompact objects are formed from entropic density perturbations.

We consider the cosmological model in which a part of the Universe \Omega_h\sim 10^-5 is in the form of primordial black holes with mass \sim 10^5M_\odot. These primordial black holes would be centers for growing protogalaxies which experienced multiple mergers with ordinary galaxies. This process of galaxies formation is accompanied by the merging of central black holes in the galactic nuclei. It is shown that recently discovered correlations between the central black holes and bulges of galaxies are naturally reproduced in this scenario.

Non-steady and eruptive phenomena in quasars are thought to be associated with the Eddington or super-Eddington luminous stage. Although there is no lack in hypotheses about the total duration of such a stage, the latter remains essentially unknown. We calculate the duration of quasar luminous phase in dependence upon the initial mass of a newborn massive black hole (MBH) by comparing the observed luminosity- and redshift distributions of quasars with mass distribution of the central MBHs in normal galactic nuclei. It is assumed that, at the quasar stage, each MBH goes through a single (or recurrent) phase(s) of accretion with, or close to, the Eddington luminosity. The mass distributions of quasars is found to be connected with that of MBHs residing in normal galaxies by a one-to-one corrrespondence through the entire mass range of the inferred MBHs if the accretion efficiency of mass-to-energy transformation $\eta \sim 0.1$.

Recently we have elaborated a new cosmological model of gamma-ray burst (GRB) origin (1998, ApJ 502, 192), which employs the dynamical evolution of central dense stellar clusters in the galactic nuclei. Those clusters inevitably contain a large fraction of compact stellar remnants (CSRs), such as neutron stars (NSs) and stellar mass black holes (BHs), and close encounters between them result in radiative captures into short-living binaries, with subsequent merging of the components, thereby producing GRBs (typically at large distances from the nucleus). In the present paper, we calculate the redshift distribution of the rate of GRBs produced by close encounters of NSs in distant galactic nuclei. To this end, the following steps are undertaken: (i) we establish a connection between the parameters of the fast evolving central stellar clusters (i.e. those for which the time of dynamical evolution exceeds the age of the Universe) with masses of the forming central supermassive black holes (SMBHs) using a dynamical evolution model; (ii) we connect these masses with the inferred mass distributions of SMBHs in the galactic nuclei and the redshift distribution of quasars by assuming a certain `Eddington luminosity phase' in their activity; (iii) we incorporate available observational data on the redshift distribution of quasars as well as a recently found correlation between the masses of galaxies and their central SMBHs. The resulting redshift distribution of the GRB rate, which accounts for both fast and slowly evolving galactic nuclei is consistent with that inferred from the BATSE data if the fraction of fast evolving galactic nuclei is in the range $0.016-0.16$.

A new cosmological scenario for the origin of gamma ray bursts (GRBs) is proposed. In our scenario, a highly evolved central core in the dense galactic nucleus is formed containing a subsystem of compact stellar remnants (CSRs), such as neutron stars and black holes. Those subsystems result from the dynamical evolution of dense central stellar clusters in the galactic nuclei through merging of stars, thereby forming (as has been realized by many authors) the short-living massive stars and then CSRs. We estimate the rate of random CSR collisions in the evolved galactic nuclei by taking into account, similar to Quinlan & Shapiro (1987), the dissipative encounters of CSRs, mainly due to radiative losses of gravitational waves, which results in the formation of intermediate short-living binaries, with further coalescence of the companions to produce GRBs. We also consider how the possible presence of a central supermassive black hole, formed in a highly evolved galactic nucleus, influences the CSR binary formation. This scenario does not postulate ad hoc a required number of tight binary neutron stars in the galaxies. Instead, it gives, for the most realistic parameters of the evolved nuclei, the expected rate of GRBs consistent with the observed one, thereby explaining the GRB appearance in a natural way of the dynamical evolution of galactic nuclei. In addition, this scenario provides an opportunity for a cosmological GRB recurrence, previously considered to be a distinctive feature of GRBs of a local origin only. We also discuss some other observational tests of the proposed scenario.