Review of the genuine black holes silhouettes. The isolated black hole is an invisible object in general relativity due to the infinite red-shift of photons moving from the black event horizon to a distant observer. Nevertheless, the dark shadow (silhouette) of the astrophysical black hole may be visible by a distant observer on the luminous background lensed in the black hole gravitational field. The black hole shadow is a projection on the celestial sphere of the photon capture cross-section. The classical black hole shadow has the maximal size if the luminous background is placed far behind the black hole or, more definitely, on the distance greatly exceeding the corresponding event horizon radius. The shadow of minimal size may be viewed if the same black hole is highlighted by the internal part of the luminous accretion disk, adjoining the black hole event horizon. In the case of a thin accretion disk this black hole shadow is the lensed silhouette of either the northern or southern hemisphere of its event horizon (depending on the orientation of the black hole rotation axis with respect to a distant observer). The dark silhouette of the southern hemisphere of the event horizon is seen on the first image of the supermassive black hole M87* obtained by the Event Horizon Telescope. The brightness of accretion matter near the black hole highly exceeds the corresponding one of the usual astrophysical distant background in the form of stars or hot gas clouds. For this reason, it is difficult to observe the classical shadow of the accreting black hole.

We propose the simple new method for extracting the value of the black hole spin from the direct high-resolution image of black hole by using a thin accretion disk model. In this model the observed dark region on the first image of the supermassive black hole in the galaxy M87, obtained by the Event Horizon Telescope, is a silhouette of the black hole event horizon. The outline of this silhouette is the equator of the event horizon sphere. The dark silhouette of the black hole event horizon is placed within the expected position of the black hole shadow, which is not revealed on the first image. We calculated numerically the relation between the observed position of the black hole silhouette and the brightest point in the thin accretion disk, depending on the black hole spin. From this relation we derive the spin of the supermassive black hole M87*, $a=0.75\pm0.15$.

We demonstrate that a dark silhouette of the black hole illuminated by a thin accretion disk and seen by a distant observer is, in fact, a silhouette of the event horizon hemisphere. The boundary of this silhouette is a contour of the event horizon equatorial circle if a thin accretion disk is placed in the black hole equatorial plane. A luminous matter plunging into black hole from different directions provides the observational opportunity for recovering a total silhouette of the invisible event horizon globe. The event horizon silhouette is projected on the celestial sphere within a position of the black hole shadow. A relative position of brightest point in the accretion disk with respect to the position of event horizon silhouette in the image of black hole in the galaxy M87, observed by the Event Horizon Telescope, corresponds to a rather high value of the black hole spin, $a\simeq0.75$.

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.

We argue that a genuine image of the black hole viewed by a distant observer is not its shadow, but a more compact event horizon image probed by the luminous matter plunging into black hole. The external border of the black hole shadow is washed out by radiation from matter plunging into black hole and approaching the event horizon. This effect will crucially influence the results of future observations by the Event Horizon Telescope. We show that gravitational lensing of the luminous matter plunging into black hole provides the possibility for visualization of the event horizon. The lensed image of the event horizon is formed by the highly red-shifted photons emitted by the plunging matter very near the black hole event horizon and detected by a distant observer. The resulting event horizon image is a gravitationally lensed projection on the celestial sphere of the whole black hole event horizon sphere. Seemingly, black holes are the unique objects in the Universe which may be viewed by distant observers at once from both the front and back sides.

How the supermassive black hole SgrA* in the Milky Way Center looks like for a distant observer? It depends on the black hole highlighting by the surrounding hot matter. The black hole shadow (the photon capture cross-section) would be viewed if there is a stationary luminous background. The black hole event horizon is invisible directly (per se). Nevertheless, a more compact (with respect to black hole shadow) projection of the black hole event horizon on the celestial sphere may be reconstructed by detecting the highly red-shifted photons emitted by the non-stationary luminous matter plunging into the black hole and approaching to the event horizon. It is appropriate to call this reconstructed projection of the event horizon on the celestial sphere for a distant observer as the ``lensed event horizon image'', or simply the ``event horizon image''. This event horizon image is placed on the celestial sphere within the position of black hole shadow. Amazingly, the event horizon image is a gravitationally lensed projection on the celestial sphere of the whole surface of the event horizon globe. In result, the black holes may be viewed at once from both the front and back sides. The lensed event horizon image may be considered as a genuine silhouette of the black hole. For example, a dark northern hemisphere of the event horizon image is the simplest model for a black hole silhouette in the presence of a thin accretion disk.

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.

The gravitational lensing of a finite star moving around a rotating Kerr black hole has been numerically calculated. Calculations for the direct image of the star and for the first and second light echoes have been performed for the star moving with an orbital period of 3.22 h around the supermassive black hole SgrA* at the Galactic Center. Time dependences for the observed star position on the celestial sphere, radiation flux from the star, frequency of detected radiation, major and minor semiaxes of the lensed star image have been calculated and plotted. The detailed observation of such lensing requires a space interferometer such as the Russian Millimetron project.

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.

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.

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.

Neutralino annihilation in the Galactic halo is the most definite observational signature proposed for indirect registration of the SUSY Dark Matter (DM) candidate particles. The corresponding annihilation signal (in the form of gamma-rays, positrons and antiprotons) may be boosted for one or three orders of magnitude due to the clustering of cold DM particles into the small-scale and very dense self-gravitating clumps. We discuss the formation of these clumps from the initial density perturbations and their successive fate in the Galactic halo. Only a small fraction of these clumps, $\sim0.1$%, in each logarithmic mass interval $\Delta\log M\sim1$ survives the stage of hierarchical clustering. We calculate the probability of surviving the remnants of dark matter clumps in the Galaxy by modelling the tidal destruction of the small-scale clumps by the Galactic disk and stars. It is demonstrated that a substantial fraction of clump remnants may survive through the tidal destruction during the lifetime of the Galaxy. 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. The survived dense remnants of tidally destructed clumps provide an amplification (boosting) of the annihilation signal with respect to the diffuse DM in the Galactic halo. We describe the anisotropy of clump distribution caused by the tidal destruction of clumps in the Galactic disk.

Quasi-Periodic Oscillations (QPOs) of the hot plasma spots or clumps orbiting an accreting black hole contain information on the black hole mass and spin. The promising observational signatures for the measurement of black hole mass and spin are the latitudinal oscillation frequency of the bright spots in the accretion flow and the frequency of black hole event horizon rotation. Both of these frequencies are independent of the accretion model and defined completely by the properties of the black hole gravitational field. Interpretation of the known QPO data by dint of a signal modulation from the hot spots in the accreting plasma reveals the Kerr metric rotation parameter, $a=0.65\pm0.05$, and mass, $M=(4.2\pm0.2)10^6M_\odot$, of the supermassive black hole in the Galactic center. At the same time, the observed 11.5 min QPO period is identified with a period of the black hole event horizon rotation, and, respectively, the 19 min period is identified with a latitudinal oscillation period of hot spots in the accretion flow. The described approach is applicable to black holes with a low accretion rate, when accreting plasma is transparent up to the event horizon region.

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 consider test planet and photon orbits of the third kind inside a black hole, which are stable, periodic and neither come out of the black hole nor terminate at the singularity. Interiors of supermassive black holes may be inhabited by advanced civilizations living on planets with the third-kind orbits. In principle, one can get information from the interiors of black holes by observing their white hole counterparts.

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.

We describe the formation of a seed massive black hole (MBH) inside a supermassive star (SMS) in a distant galactic nucleus. The short-lived SMS is naturally formed due to collision destructions of normal stars in the evolving galactic nucleus. The neutron stars (NSs) and stellar-mass black holes form a compact self-gravitating subsystem deep inside a SMS. This subsystem is short-lived in comparison with a host SMS and collapses finally into the MBH. Just before gravitational collapse of compact subsystem the frequent NS collisions are accompanied by the generation of numerous ultra-relativistic fireballs. A combined ram pressure of multiple coexisting fireballs produces a quasi-stationary rarefied cavity in the central part of SMS. The protons are accelerated in the fireballs and by relativistic shocks in the cavity. All secondary particles, produced in collisions, except the high-energy neutrinos are absorbed in the SMS interiors. An estimated high-energy neutrino signal from this hidden source can be detected by the neutrino telescope with an effective area 1 km^2 providing the evidence for MBH formation in a distant galactic nucleus. A corresponding lifetime of this high-energy hidden neutrino source is ~0.1-1 yr.

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 study the gravitational radiation from a rotating monopole-antimonopole pair connected by a string. While at not too high frequencies the emitted gravitational spectrum is described asymptotically by $P_n\propto n^{-1}$, the spectrum is exponentially suppressed in the high-frequency limit, $P_n\propto \exp(-n/n_{\rm cr})$. Below $n_{\rm cr}$, the emitted spectrum of gravitational waves is very similar to the case of an oscillating monopole pair connected by a string, and we argue therefore that the spectrum found holds approximately for any moving monopole-string system. As application, we discuss the stochastic gravitational wave background generated by monopole-antimonopole pairs connected by strings in the early Universe and gravitational wave bursts emitted at present by monopole-string networks. We confirm that advanced gravitational wave detectors have the potential to detect a signal for string tensions as small as $G\mu\sim 10^{-13}$.

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 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 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 calculate the gravitational radiation from a cusp of the chiral cosmic strings as a function of the current on the string in the limit of small values of the current. The smoothing of the cosmic string cusp due to the presence of the superconducting current on the string leads to the different behavior of the gravitational radiation from the cusp as compared with the non-superconducting (ordinary) string. The difference between the gravitational energy radiation from the chiral cusp and the ordinary one is proportional to the value of the current on the string. It is shown that there is a high-frequency cutoff in the spectrum of the gravitational radiation from the chiral string. The frequency cutoff depends on the value of the string current. This effect is crucial for the detection of gravitation wave bursts from cosmic strings: the rather large current on the string would lower the amplitude of the incoming signal.

We present the generalization of the Sedov-Taylor self-similar strong spherical shock solution for the case of a central energy source varying in time, $E=A t^k$, where $A$ and $k$ are constants. The known Sedov-Taylor solution corresponds to a particular adiabatic case of $k=0$ or \emph{instant shock} with an instant energy source of the shock, $E=A$. The self-similar hydrodynamic flow in the nonadiabatic $k\neq0$ case exists only under the appropriate local entropy (energy) input which must be supported by some radiative mechanism from the central engine. The specific case of $k=1$ corresponds to a permanent energy injection into the shock, or injection shock with a central source of constant luminosity, $L=A$, $E=A t$. The generalized self-similar shock solution may be applied to astrophysical objects in which the duration of central source activity is longer than the shock expansion time, e.g. the early phase of SN explosions, strong wind from stars and young pulsars, non-steady spherical outflow from black holes and collapsing dense stellar clusters with numerous neutron star collisions.

We derived expressions for energy, momentum and angular momentum losses due to gravitational and electromagnetic radiation from the closed superconducting chiral cosmic strings of arbitrary form. The expressions for corresponding radiation rates into the unit solid angle have the form of four-dimensional integrals. In the special case of piece-wise linear strings these formulas are reduced to sums over the kinks. We calculate numerically the total radiation rates for three examples of string loops in dependence of current along the string.

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.

We propose the model of a short-lived very powerful source of high energy neutrinos. It is formed as a result of the dynamical evolution of a galactic nucleus prior to its collapse into a massive black hole and formation of high-luminosity AGN. This stage can be referred to as ``pre-AGN''. A dense central stellar cluster in the galactic nucleus on the late stage of evolution consists of compact stars (neutron stars and stellar mass black holes). This cluster is sunk deep into massive gas envelope produced by destructive collisions of a primary stellar population. Frequent collisions of neutron stars in a central stellar cluster are accompanied by the generation of ultrarelativistic fireballs and shock waves. These repeating fireballs result in a formation of the expanding rarefied cavity inside the envelope. The charged particles are effectively accelerated in the cavity and, due to pp-collisions in the gas envelope, they produce high energy neutrinos. All high energy particles, except neutrinos, are absorbed in the thick envelope. Duration of this pre-AGN phase is about 10 yr, the number of the sources can be \~ 10 per cosmological horizon. High energy neutrino signal can be detected by underground neutrino telescope with effective area ~1 km^2.

The gravitational lensing of a relativistic fireball can produce the time delayed multiple images with quite different spectra and temporal patterns in contrast with a nonrelativistic source and hence can imitate the source recurrence. In particular, the enigmatic four multiple gamma-ray bursts detected during two days at October 1996 in the same sky region may be due to a single fireball event multiply imaged by the foreground galactic nucleus or cluster of galaxies.

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.