In this paper we investigate the relation between the potential, and the geometrical time delays in the gravitational lensing. In the original paper of Shapiro (1964), it is stated that there is a time delay in the radar signals between Earth and Venus, passing near a massive object (the Sun), relative to a case of the absence of such mass. The reason is connected with the influence of the gravity on the speed of a light rays in a gravitational potential. The contribution from the change of the path length, which happens to be of a second order, is considered as negligible. Nevertheless, in the gravitational lens theory the geometrical delay, related to the change of path length, is routinely taken into account along with the potential term. Inspired by this apparent discrepancy, in this work we address the contribution of the geometric part of the time delay in different situations. We introduce a unified treatment with two limiting regimes of lensing. One of them corresponds to the time delay experiments near the Sun where geometrical delay is shown to be negligible. The second corresponds to the typical gravitational lens scenario with multiple imaging where the geometrical delay is shown to be significant. In particular, it is shown that the geometrical delay is non-negligible in the case of good alignment between source, lens and observer, because in such a case it becomes a first order quantity.

One of the problem revealed recently in cosmology is a so-called Hubble tension (HT), which is the difference between values of the present Hubble constant, measured by observation of the universe at redshift z smaller than 1, and by observations of a distant universe with CMB fluctuations originated at z about 1100. In this paper we suggest, that this discrepancy may be explained by deviation of the cosmological expansion from a simple Friedman model of a flat dusty universe during the period after recombination at z less than 1100, due to action of additional small component of a dark energy of a different origin.. We suppose, that a dark matter (DM) has a common origin with a small component of a dark energy (DEV). DE presently may have two components, one of which is the Einstein constant, and another, smaller component DEV comes from the remnants of a scalar fields responsible for inflation. Due to common origin and interconnections the ratio between DEV and DM densities is supposed to remain constant during, at least, the time after recombination. This part of the dark energy in not connected with the cosmological constant Lambda, but is defined by existence of scalar fields with a variable density. Taking into account the influence of DEV on the universe expansion we find the value of this ratio, which should remove the HT problem.

The formation of the chemical composition of neutron star envelopes, at high densities is considered. As hot matter is compressed in the process of collapse, which leads to the explosion of a core collapse supernova, the stage of nuclear equilibrium with free neutrino escape, kinetic equilibrium in \b{eta}-processes, and, as a result, the establishment of limited nuclear equilibrium with a fixed number of nuclei takes place. Cold matter is compressed at a fixed number of nuclei whose atomic weight initially does not change and subsequently decreases. A pycnonuclear reaction of the fusion of available nuclei and a decrease in their number begin at the end. The compression of cold matter is accompanied by an increase in the mass fraction of free neutrons. In this case, the chemical composition of the envelope differs significantly from the equilibrium one and contains a considerable store of nuclear energy. Nonequilibrium \b{eta}-reactions proceed at densities exceeding the upper bound for the non-equilibrium layer density, which lead to heating, nuclear energy release, and the possible attainment of a state of complete thermodynamic equilibrium. The thermodynamics of nonequilibrium \b{eta}-processes, which lead to the heating of matter as neutrinos escape freely, is considered.

In this paper, we investigate the influence of the plasma surrounding the gravitational lens on the effect of microlensing. In presence of plasma around the lens, the deflection angle is determined by both the gravitational field of the lens and the chromatic refraction in the inhomogeneous plasma. We calculate microlensing light curves numerically for point-mass lens surrounded by power-law density distribution of plasma. A variety of possible curves is revealed, depending on the plasma density and frequency of observations. In the case of significant influence of plasma, the shape of microlensing light curve is strongly deformed in comparison with vacuum case. If the refractive deflection is large enough to compensate or to overcome the gravitational deflection, microlensing images can completely disappear for the observer. In this case, the remarkable effect occurs: formation of a 'hole' instead of a 'hill' in the center of microlensing light curve. Observational prospects of 'hill-hole' effect in different microlensing scenarios are discussed.

Periodic changes in a thermal soft X-ray flux of a rotating neutron star indicate a non-uniform distribution of the surface temperature. A possible cause of this phenomenon is a suppression of the heat flux across the magnetic field lines in a crust and an envelope of magnetized neutron stars. In this paper we study three-dimensional effects, associated with non-axisymmetric magnetic fields in neutron stars. We calculate the surface temperature distribution by solving numerically a three dimensional heat transfer equation in a magnetized neutron star crust. We adopt an anisotropic (tensorial) electron thermal conductivity coefficient, which is derived as an analytical solution of the Boltzmann equation with a Chapman-Enskog method. To calculate the surface temperature distribution, we construct a local one-dimensional plane-parallel model ("Ts-Tb"-relationship) of a magnetized neutron star envelope. We then use it as an outer boundary condition for the three-dimensional problem in the crust to find the self-consistent solution. To study possible observational manifestations from anisotropic temperature distributions we calculate light curves with a composite black-body model. Our calculations show, that a non-axisymmetric magnetic field distribution can lead to the irregular non-sinusoidal shape of a pulse profile as well as in some cases a significant amplification of pulsations of the thermal flux in comparison to the pure-dipolar magnetic field configurations.

X-ray binaries are cauldrons of fundamental physical processes which appear along practically the whole electromagnetic spectrum. The sub-class of X-ray transient sources show multifrequency behaviour which deserve particular attention in order to understand the causing physics. These binary systems consist of a compact star and an optical star, therefore there is a mutual influence between these two stars that drive the low energy (LE) (i.e. radio, IR, optical) and high energy (HE) (i.e. UV, X-ray, $\gamma$-ray) processes. The LE processes are produced mostly on the optical star and the HE processes mostly on the compact star, typically a neutron star. Thus it appears evident that through the study of LE processes it is possible to understand also the HE processes and vice versa. In this paper we will discuss this problem starting from the experimental evidence of a delay between LE and HE processes detected for the first time in the X-ray/Be system A0535+26/HDE245770 (e.g. Giovannelli \& Sabau-Graziati, 2011; Giovannelli, Bisnovatyi-Kogan \& Klepnev, 2013 (here after GBK13); Giovannelli et al., 2015b). This delay is common in cataclysmic variables (CVs) and other binary systems with either a neutron star or a black hole. Since a delay between LE processes and HE processes has been experimentally observed in several active galactic nuclei (AGNs), we will discuss also the tidal disruption of stars by massive BHs, following the original idea of Rees (1998): stars in galactic nuclei can be captured or tidally disrupted by a central black hole. Some debris would be ejected at high speed, the remainder would be swallowed by the hole, causing a bright flare lasting at most a few years.

The exact time-dependent solution is obtained for a magnetic field growth during a spherically symmetric accretion into a black hole (BH) with a Schwarzschild metric. Magnetic field is increasing with time, changing from the initially uniform into a quasi-radial field. Equipartition between magnetic and kinetic energies in the falling gas is established in the developed stages of the flow. Estimates of the synchrotron radiation intensity are presented for the stationary flow. The main part of the radiation is formed in the region $r \leq 7 r_g$, here $r_g$ is a BH gravitational radius. The two-dimensional stationary self-similar magnetohydrodynamic solution is obtained for the matter accretion into BH, in a presence of a large-scale magnetic field, when the magnetic field far from the BH is homogeneous and does not influence the flow. At the symmetry plane perpendicular to the direction of the distant magnetic field, the quasi-stationary disk is formed around BH, which structure is determined by dissipation processes. Parameters of the shock forming due to matter infall onto the disk are obtained. The radiation spectrum of the disk and the shock are obtained for the $10\,\, M_\odot$ BH. The luminosity of such object is about the solar one, for a characteristic galactic gas density, with possibility of observation at distances less than 1 kpc. The spectra of a laminar and a turbulent disk structure around BH are very different. The turbulent disk emits a large part of its flux in the infrared. It may occur that some of the galactic infrared star-like sources are a single BH in the turbulent accretion state. The radiative efficiency of the magnetized disk is very high, reaching $\sim 0.5\,\dot M\,c^2$ so it was called recently as a magnetically arrested disk (MAD). Numerical simulations of MAD, and its appearance during accretion into neutron stars are considered and discussed.

In this article, we present an overview of the new developments in problems of the plasma influence on the effects of gravitational lensing, complemented by pieces of new material and relevant discussions. Deflection of light in the presence of gravity and plasma is determined by a complex combination of various physical phenomena: gravity, dispersion, refraction. In particular, the gravitational deflection itself, in a homogeneous plasma without refraction, differs from the vacuum one and depends on the frequency of the photon. In an inhomogeneous plasma, chromatic refraction also takes place. We describe chromatic effects in strong lens systems including a shift of angular position of image and a change in magnification. We also investigate high-order images that arise when lensing on a black hole surrounded by homogeneous plasma. The recent results of analytical studies of the effect of plasma on the shadow of the Schwarzschild and Kerr black holes are presented.

Jets are observed in young stellar objects, X-ray sources, active galactic nuclei (AGN). The mechanisms of jet formation may be divided in regular, acting continuously for a long time, and explosive ones. Continuous mechanisms are related with electrodynamics and radiation pressure acceleration, hydrodynamical acceleration in the nozzle inside a thick disk, acceleration by relativistic beam of particles. Explosive jet formation is connected with supernovae, gamma ray bursts and explosive events in galactic nuclei. Mechanisms of jet collimation may be connected with magnetic confinement, or a pressure of external gas. Explosive formation of jets in the laboratory is modeled in the experiments with powerful laser beam, and plasma focus.

The idea of the magnetorotational explosion mechanism is that the energy of rotation of the neutron star formed in the course of a collapse is transformed into the energy of an expanding shock wave by means of a magnetic field. In the two-dimensional case, the time of this transformation depends weakly on the initial strength of the poloidal magnetic field because of the development of a magnetorotational instability. Differential rotation leads to the twisting and growth of the toroidal magnetic-field component, which becomes much stronger than the poloidal component. In the case where the initial configuration of the magnetic field is close to a dipole configuration, the ejection of matter has a jet character, whereas, in the case of a quadrupole configuration, there arises an equatorial ejection. In either case, the energy release is sufficient for explaining the observed average energy of supernova explosion. Neutrinos are emitted as the collapse and the formation of a rapidly rotating neutron star proceeds. In addition, neutrino radiation arises in the process of magnetorotational explosion owing to additional rotational-energy losses. In order to explain an interval of 4.5 hours between the two observed neutrino signals from SN 1987A, it is necessary to assume a weakening of the magnetorotional instability and a small initial magnetic field (10(9)-10(10)G) in the newly formed rotating neutron star. The existence of a black hole in the SN 1987A remnant could explain the absence of any visible pointlike source at the center of the explosion.

The results of MHD numerical simulations of the formation and development of magnetized jets are presented. Similarity criteria for comparisons of the results of laboratory laser experiments and numerical simulations of astrophysical jets are discussed. The results of laboratory simulations of jets generated in experiments at the Neodim laser installation are presented.

A solution to the Boltzmann equation is obtained for a magnetized plasma with strongly degenerate nonrelativistic electrons and nondegenerate nuclei. The components of the diffusion, thermal diffusion and diffusion thermoeffect tensors in a nonquantizing magnetic field are calculated in the Lorentz approximation without allowance for electron-electron collisions, which is asymptotically accurate for plasma with strongly degenerate electrons. Asymptotically accurate analytical expressions for the electron diffusion, thermal diffusion and diffusion thermoeffect tensors in the presence of a magnetic field are obtained for the first time. The expressions reveal a considerably more complicated dependence on magnetic field than analogous dependences derived in the previous publications on this subject.

If the large scale structure of the Universe was created, even partially, via Zeldovich pancakes, than the fluctuations of the CMB radiation should be formed due to bulk comptonization of black body spectrum on the contracting pancake. Approximate formulas for the CMB energy spectrum after bulk comptonization are obtained. The difference between comptonized energy spectra of the CMB due to thermal and bulk comptonization may be estimated by comparison of the plots for the spectra in these two cases.

A mechanism of formation of gravitational waves in the Universe is considered for a nonspherical collapse of matter. Nonspherical collapse results are presented for a uniform spheroid of dust and a finite-entropy spheroid. Numerical simulation results on core-collapse supernova explosions are presented for the neutrino and magnetorotational models. These results are used to estimate the dimensionless amplitude of the gravitational wave with a frequency \nu ~1300 Hz, radiated during the collapse of the rotating core of a pre-supernova with a mass of 1:2M(sun) (calculated by the authors in 2D). This estimate agrees well with many other calculations (presented in this paper) that have been done in 2D and 3D settings and which rely on more exact and sophisticated calculations of the gravitational wave amplitude. The formation of the large-scale structure of the Universe in the Zel'dovich pancake model involves the emission of very long-wavelength gravitational waves. The average amplitude of these waves is calculated from the simulation, in the uniform spheroid approximation, of the nonspherical collapse of noncollisional dust matter, which imitates dark matter. It is noted that a gravitational wave radiated during a core-collapse supernova explosion in our Galaxy has a sufficient amplitude to be detected by existing gravitational wave telescopes.

Self-similar solution is obtained for propagation of a strong shock, in a flat expanding dusty Friedman universe. Approximate analytic solution was obtained earlier, using relation between self-similar variables, equivalent to the exact energy conservation integral, which was obtained by L.I. Sedov for the strong explosion in the static uniform medium. Numerical integration of self-similar equation is made here, giving an exact solution of the problem, which is rather close to the approximate analytic one. The differences between these solutions are most apparent in the vicinity of the shock. For polytropic equation of state, self-similar solutions exist in more narrow interval of the adiabatic power than in the static case.

It is shown that investigation of ionization zones left after supernova explosions can yield information on the mechanism active in those explosions.

Propagation of the strong shock in the flat expanding Friedman universe is investigated using methods of dimension and similarity. Exact analytic solution of self-similar equations is obtained, determining dependences of the radius and velocity of the shock wave on time and radius. It is obtained, that in the expanding medium the velocity of shock decreases as $\sim t^{-1/5}$, what is slower than the shock velocity in the static uniform medium $\sim t^{-3/5}$. The radius of the shock wave in the expanding self-gravitating medium increases $\sim t^{4/5}$, more rapidly than the shock wave radius in the uniform non-gravitating medium $\sim t^{2/5}$. So, the shock propagates in the direction of decreasing density with larger speed, that in the static medium, due to accelerating action of the decreasing density, even in the presence of a self-gravitation.

The observational properties of Soft Gamma Repeaters and Ano\-malous X-ray Pulsars (SGR/AXP) indicate to necessity of the energy source different from a rotational energy of a neutron star. The model, where the source of the energy is connected with a magnetic field dissipation in a highly magnetized neutron star (magnetar) is analyzed. Some observational inconsistencies are indicated for this interpretation. The alternative energy source, connected with the nuclear energy of superheavy nuclei stored in the nonequilibrium layer of low mass neutron star is discussed.

Electron-positron annihilation line from a Galactic center direction was discovered by the balloon-borne germanium gamma-ray telescope, and confirmed by OSSE experiment of CGRO mission. Extensive observations by INTEGRAL observatory permit to determine properties of the annihilation line from Galactic bulge. Possible sources of $e^{+}$ already discussed are supernovae explosions, microquasars, gamma-ray bursts, tidal disruption events, activity near black hole of Sgr A*, $e^{+}$ generation by subrelativistic cosmic rays, or even dark matter. One remarkable feature of the line emission is an absence of resolved point like sources. Any model should take into account this feature. We consider flares of low-mass stars as a possible cumulative $e^{+}$ source of the observed annihilation line from the bulge. Our estimations show that $e^{+}$ production by numerous flaring stars in Galaxy bulge can be significant, or even the main source of $e^{+}$ responsible for the annihilation line observed from the central part of our Galaxy.

We develop models for time lag between maxima of the source brightness in different wavelengths during a transient flash of luminosity connected with a short period of increase of the mass flux onto the central compact object. We derive a simple formula for finding the time delay among events in different wavelengths, valid in general for all disk accreting cosmic sources, and discuss quantitatively a model for time lag formation in AGNs. In close binaries with accretion disks the time lag is connected with effects of viscosity defining a radial motion of matter in the accretion disk. In AGN flashes, the falling matter has a low angular momentum, and the time lag is defined by the free fall time to the gravitating center. We show the validity of these models by means of several examples of galactic and extragalactic accreting sources.

The solution of Boltzmann equation for plasma in magnetic field, with arbitrarily degenerate electrons and non-degenerate nuclei, is obtained by Chapman-Enskog method. Functions, generalizing Sonin polynomials are used for obtaining an approximate solution. Fully ionized plasma is considered. The tensor of the heat conductivity coefficients in non-quantized magnetic field is calculated. For non-degenerate and strongly degenerate plasma the asymptotic analytic formulas are obtained, which are compared with results of previous authors. The Lorentz approximation, with neglecting of electron-electron encounters, is asymptotically exact for strongly degenerate plasma. We obtain, for the first time, in three polynomial approximation, with account of electron-electron collisions, analytical expressions for the heat conductivity tensor for non-degenerate electrons, in presence of a magnetic field. Account of the third polynomial improved substantially the precision of results. In two polynomial approximation our solution coincides with the published results. For strongly degenerate electrons we obtain, for the first time, an asymptotically exact analytical solution for the heat conductivity tensor in presence of a magnetic field. This solution has considerably more complicated dependence on the magnetic field than those in previous publications, and gives several times smaller relative value of a thermal conductivity across the magnetic field at $\omega\tau\gtrsim 0.8$.

Strong discontinuities in solutions of the gas dynamic equations under isentropic conditions, i.e., with continuity of entropy at the discontinuity, are examined. Solutions for a standard shock wave with continuity of energy at the discontinuity are compared with those for an isentropic 'shock wave'. It is shown that numerical simulation of astrophysical problems in which high-amplitude shock waves are encountered (supernova explosions, modelling of jets) with conservation of entropy, rather than of energy, leads to large errors in the shock calculations. The isentropic equations of gas dynamics can be used only when there are no strong discontinuities in the solution or when the intensity of the shocks is not high and they do not significantly affect the flow.

We compare polarization properties of the cyclotron, and relativistic dipole radiation of electrons moving in the magnetic field on a helix with ultra-relativistic longitudinal and non-relativistic transverse velocity components. The applicability of these models in the case of accretion onto a neutron star is discussed. The test, based on polarization observations is suggested, to distinguish between the cyclotron, and relativistic dipole origin of features, observed in X-ray spectra of some X-ray sources, among which the Her X-1 is the most famous.

Estimations of magnetic fields of neutron stars, observed as radio and X-ray pulsars, are discussed. It is shown, that theoretical and observational values for different types of radiopulsars are in good correspondence. Radiopulsars in close binaries and millisecond pulsars, which have passed the stage of disk accretion (recycled radiopulsars), have magnetic fields 2-4 orders of magnitude smaller than ordinary single pulsars. Most probably, the magnetic field of the neutron star was screened by the infalling material. Several screening models are considered. Formation of single recycled pulsars loosing its companion is discussed. Magnetic fields of some X-ray pulsars are estimated from the cyclotron line energy. In the case of Her X-1 this estimation exceeds considerably the value of its magnetic field obtained from long term observational data related to the beam structure evolution. Another interpretation of the cyclotron feature, based on the relativistic dipole radiation mechanism, could remove this discrepancy. Observational data about soft gamma repeaters and their interpretation as magnetars are critically analyzed.

We analytically calculate the influence of a plasma on the shadow of a black hole (or of another compact object). We restrict to spherically symmetric and static situations, where the shadow is circular. The plasma is assumed to be non-magnetized and pressure-less. We derive the general formulas for a spherically symmetric plasma density on an unspecified spherically symmetric and static spacetime. The formalism applies not only to black holes but also, e.g., to wormholes. As examples for the underlying spacetime model, we consider the Schwarzschild spacetime and the Ellis wormhole. In particular, we treat the case that the plasma is in radial free fall from infinity onto a Schwarzschild black hole. The perspectives of actually observing the influence of a plasma on the shadows of supermassive black holes are discussed.

The influence of plasma on different effects of gravitational lensing is reviewed. Using the Hamiltonian approach for geometrical optics in a medium in the presence of gravity, an exact formula for the photon deflection angle by a black hole (or another body with a Schwarzschild metric) embedded in plasma with a spherically symmetric density distribution is derived. The deflection angle in this case is determined by the mutual combination of different factors: gravity, dispersion, and refraction. While the effects of deflection by the gravity in vacuum and the refractive deflection in a nonhomogeneous medium are well known, the new effect is that, in the case of a homogeneous plasma, in the absence of refractive deflection, the gravitational deflection differs from the vacuum deflection and depends on the photon frequency. In the presence of a plasma nonhomogeneity, the chromatic refractive deflection also occurs, so the presence of plasma always makes gravitational lensing chromatic. In particular, the presence of plasma leads to different angular positions of the same image if it is observed at different wavelengths. It is discussed in detail how to apply the presented formulas for the calculation of the deflection angle in different situations. Gravitational lensing in plasma beyond the weak deflection approximation is also considered.

In this paper we review a recently developed approximate method for investigation of dynamics of compressible ellipsoidal figures. Collapse and subsequent behaviour are described by a system of ordinary differential equations for time evolution of semi-axes of a uniformly rotating, three-axis, uniform-density ellipsoid. First, we apply this approach to investigate dynamic stability of non-spherical bodies. We solve the equations that describe, in a simplified way, the Newtonian dynamics of a self-gravitating non-rotating spheroidal body. We find that, after loss of stability, a contraction to a singularity occurs only in a pure spherical collapse, and deviations from spherical symmetry prevent the contraction to the singularity through a stabilizing action of nonlinear non-spherical oscillations. The development of instability leads to the formation of a regularly or chaotically oscillating body, in which dynamical motion prevents the formation of the singularity. We find regions of chaotic and regular pulsations by constructing a Poincare diagram. A real collapse occurs after damping of the oscillations because of energy losses, shock wave formation or viscosity. We use our approach to investigate approximately the first stages of collapse during the large scale structure formation. The theory of this process started from ideas of Ya. B. Zeldovich, concerning the formation of strongly non-spherical structures during nonlinear stages of the development of gravitational instability, known as 'Zeldovich's pancakes'. In this paper the collapse of non-collisional dark matter and the formation of pancake structures are investigated approximately. We estimate an emission of very long gravitational waves during the collapse, and discuss the possibility of gravitational lensing and polarization of the cosmic microwave background by these waves.

We consider the motion of classical spinning test particles in Schwarzschild and Kerr metrics and investigate innermost stable circular orbits (ISCO). The main goal of this work is to find analytically the small-spin corrections for the parameters of ISCO (radius, total angular momentum, energy, orbital angular frequency) of spinning test particles in the case of vectors of black hole spin, particle spin and orbital angular momentum being collinear to each other. We analytically derive the small-spin linear corrections for arbitrary Kerr parameter $a$. The cases of Schwarzschild, slowly rotating and extreme Kerr black hole are considered in details. For a slowly rotating black hole the ISCO parameters are obtained up to quadratic in $a$ and particle's spin $s$ terms. From the formulae obtained it is seen that the spin-orbital coupling has attractive character when spin and angular momentum are parallel and repulsive when they are antiparallel. For the case of the extreme Kerr black hole with co-rotating particle we succeed to find the exact analytical solution for the limiting ISCO parameters for arbitrary spin. It has been shown that the limiting values of ISCO radius and frequency do not depend on the particle's spin while values of energy and total angular momentum depend on it. We have also considered circular orbits of arbitrary radius and have found small-spin linear corrections for the total angular momentum and energy at given radius. System of equations for numerical calculation of ISCO parameters for arbitrary $a$ and $s$ is also explicitly written.

In this paper we address the question of why the observed periods of the Anomalous X-ray Pulsars (AXPs) and Soft Gamma-ray Repeaters (SGRs) are clustered in the range 2-12 s. We explore a possibility to answer this question assuming that AXPs and SGRs are the descendants of High Mass X-ray Binaries (HMXBs) which have been disintegrated in the core-collapse supernova explosion. The spin period of neutron stars in HMXBs evolves towards the equilibrium period, averaging around a few seconds. After the explosion of its massive companion, the neutron star turns out to be embedded into a dense gaseous envelope, the accretion from which leads to the formation of a residual magnetically-levitating (ML) disk. We show that the expected mass of a disk in this case is 10^-7 - 10^-8 M_sun which is sufficient to maintain the process of accretion at the rate 10^14 - 10^15 g/s over a time span of a few thousand years. During this period the star manifests itself as an isolated X-ray pulsar with a number of parameters resembling those of AXPs and SGRs. Period clustering of such pulsars can be provided if the lifetime of the residual disk does not exceed the spin-down timescale of the neutron star.

The optical behaviour of the Be star in the high-mass X-ray transient A0535+26/HDE245770 shows that at periastron the luminosity is typically enhanced by 0.02 to a few tenths magnitude, and the X-ray outburst occurs eight days after the periastron. Indeed, at the periastron an increase of the mass flux occurs. This sort of flush reaches the external part of the temporary accretion disk around the neutron star and moves to the hot central parts of the accretion disk and the neutron star surface. The time necessary for this way is dependent on the turbulent viscosity in the accretion disk, as discussed by Giovannelli, Bisnovatyi-Kogan, and Klepnev (2013) (GBK13). In this paper we will show the behaviour of this system in optical band around the predicted periastron passage on 21st February 2014, by using the GBK13 ephemeris that we used to schedule our spectroscopic and photometric optical observations. Spectroscopic unusual activity detected in the Balmer lines and the enhancement in the emission in B, V, and R bands around the periastron passage, and the subsequent X-ray event definitively demonstrate the existence of about 8 day delay between optical and X-ray flares.

Core-collapse supernovae are accompanied by formation of neutron stars. The gravitation energy is transformed into the energy of the explosion, observed as SN II, SN Ib,c type supernovae. We present results of 2-D MHD simulations, where the source of energy is rotation, and magnetic field serves as a "transition belt" for the transformation of the rotation energy into the energy of the explosion. The toroidal part of the magnetic energy initially grows linearly with time due to differential rotation. When the twisted toroidal component strongly exceeds the poloidal field, magneto-rotational instability develops, leading to a drastic acceleration in the growth of magnetic energy. Finally, a fast MHD shock is formed, producing a supernova explosion. Mildly collimated jet is produced for dipole-like type of the initial field. At very high initial magnetic field no MRI development was found.

The optical behaviour of the Be star in the high mass X-ray transient A0535+26/HDE245770 shows that at the periastron typically there is an enhancement in the luminosity of order 0.05 to few tenths mag, and the X- ray outburst happens about 8 days after the periastron. We construct a quantitative model of this event, basing on a nonstationary accretion disk behavior, connected with a high ellipticity of the orbital motion. We explain the observed time delay between the peaks of the optical and X-ray outbursts in this system by the time of radial motion of a matter in the accretion disk, after increase of the mass flux in the vicinity of a periastral point in the binary. This time is determined by the turbulent viscosity, with the parameter $ \alpha=0.1-0.3$, estimated from the comparison of the model with observational data

We explore a possibility to explain the phenomenon of the Anomalous X-ray Pulsars (AXP) and Soft Gamma-ray Repeaters (SGR) within the scenario of fall-back magnetic accretion onto a young isolated neutron star. The X-ray emission of the pulsar in this case is originated due to accretion of matter onto the surface of the neutron star from the magnetic slab surrounding its magnetosphere. The expected spin-down rate of the neutron star within this approach is close to the observed value. We show that these neutron stars are relatively young and are going through a transition from the propeller state to the accretor state. The pulsars activity in the gamma-rays is connected with their relative youth and is provided by the energy stored in the non-equilibrium layer located in the crust of low-mass neutron stars. This energy can be released due to mixing of matter in the neutron star crust with super heavy nuclei approaching its surface and getting unstable. The nuclei fission in the low-density region initiates chain reactions leading to the nuclear explosion. The outbursts are likely to be triggered by an instability developing in the region where the matter accreted by the neutron star is accumulated at the magnetic pole regions.

Context. The external regions of galaxy clusters may be under strong influence of the dark energy, discovered by observations of the SN Ia at redshift z < 1. Aims. The presence of the dark energy in the gravitational equilibrium equation, with the Einstein $\Lambda$ term, contrasts the gravity, making the equilibrium configuration more extended in radius. Methods. In this paper we derive and solve the kinetic equation for an equilibrium configuration in presence of the dark energy, by considering the Newtonian regime, being the observed velocities of the galaxies inside a cluster largely smaller than the light velocity. Results. The presence of the dark energy in the gravitational equilibrium equation leads to wide regions in the $W_0-\rho_{\Lambda}$ diagram where the equilibrium solutions are not permitted, due to the prevalence of the effects of the dark energy on the gravity.

We investigate the influence of plasma presence on relativistic images formed by Schwarzschild black hole lensing. When a gravitating body is surrounded by a plasma, the lensing angle depends on a frequency of the electromagnetic wave due to refraction properties, and the dispersion properties of the light propagation in gravitational field in plasma. The last effect leads to difference, even in uniform plasma, of gravitational deflection angle in plasma from vacuum case. This angle depends on the photon frequency, what resembles the properties of the refractive prism spectrometer. Here we consider the case of a strong deflection angle for the light, traveling near the Schwarzschild black hole, surrounded by a uniform plasma. Asymptotic formulae are obtained for the case of a very large deflection angle, exceeding $2\pi$. We apply these formulae for calculation of position and magnification of relativistic images in a homogeneous plasma, which are formed by the photons performing one or several revolutions around the central object. We conclude that the presence of the uniform plasma increases the angular size of relativistic rings or the angular separation of point images from the gravitating center. The presence of the uniform plasma increases also a magnification of relativistic images. The angular separation and the magnification become significantly larger than in the vacuum case, when the photon frequency goes to a plasma frequency.

The optical behaviour of the Be star in the high mass X-ray transient A0535+26/HDE245770 shows that at the periastron typically there is an enhancement in the luminosity of order 0.02 to few tenths mag, and the X-ray outburst happens about 8 days after the periastron. We construct a quantitative model of this event, basing on the a nonstationary accretion disk behavior, connected with a high ellipticity of the orbital motion. The ephemeris used in this paper -- JD$_{\rm opt-outb}$ = JD$_0$(2,444,944) $\pm$ n(111.0 $\pm$ 0.4) days are derived from the orbital period of the system P$_{\rm orb} = 111.0 \pm 0.4$ days, determined by Priedhorsky & Terrell (1983), and from the optical flare of December 5, 1981 (Giovannelli et al., 1985) (here after 811205-E; E stands for the Event occurred at that date) that triggered the subsequent X-ray outburst of December 13, 1981 (Nagase et al., 1982) (here after 811213-E). We explain the observed time delay between the peaks of the optical and X-ray outbursts in this system by the time of radial motion of the matter in the accretion disk, after an increase of the mass flux in the vicinity of a periastral point in the binary. This time is determined by the turbulent viscosity, with the parameter $\alpha=0.1-0.3$. The increase of the mass flux is a sort of flush that reaches the external part of the accretion disk around the neutron star, producing an enhancement in the optical luminosity. The consequent X-ray flare happens when the matter reaches the hot central parts of the accretion disk, and the neutron star surface.

The coefficients that determine the electron heat transfer and diffusion in the crust of neutron stars are calculated on the basis of a solution of the Boltzmann equation with allowance for degeneracy.

We obtain a solution for the hydrodynamic outflow of the polytropic gas from the gravitating center, in presence of the uniform Dark Energy (DE). The antigravity of DE is enlightening the outflow and make the outflow possible at smaller initial temperature, at the same density. The main property of the wind in presence of DE is its unlimited acceleration after passing the critical point. In application of this solution to the winds from galaxy clusters we suggest that collision of the strongly accelerated wind with another galaxy cluster, or with another galactic cluster wind could lead to the formation of a highest energy cosmic rays.

{We consider the Coma cluster of galaxies as a gravitationally bound physical system embedded in the perfectly uniform static dark energy background as implied by the $\Lambda$CDM cosmology.} {We ask if the density of dark energy is high enough to affect the structure of a large rich cluster of galaxies?} {We use recent observational data on the cluster together with our theory of local dynamical effects of dark energy, including the zero-gravity radius $R_{\rm ZG}$ of the local force field as the key parameter.} {1) Three masses are defined which characterize the structure of a regular cluster: the matter mass $M_{\rm M}$, the dark-energy effective mass $M_{\rm DE}$ ($<0$) and the gravitating mass $M_{\rm G}$ ($= M_{\rm M} + M_{\rm DE}$). 2) A new matter density profile is suggested which reproduces well the observational data for the Coma cluster in the radius range from 1.4 Mpc to 14 Mpc and takes into account the dark energy background. 3) Using this profile, we calculate upper limits for the total size of the Coma cluster, $R \le R_{\rm ZG} \approx 20$ Mpc, and its total matter mass, $M_{\rm M} \la M_{\rm M}(R_{\rm ZG}) = 6.2 \times 10^{15} M_{\odot}$.} {The dark energy antigravity affects strongly the structure of the Coma cluster at large radii $R \ga 14$ Mpc and should be taken into account when its total mass is derived.}

The problem of steady-state accretion to nonrotating black holes is examined. Advection is included and generalized formulas for the radiation pressure in both the optically thick and thin cases are used. Special attention is devoted to models with a high accretion rate. Global solutions for accretion disks are studied which describe a continuous transition between an optically thick outer region and an optically thin inner region. It is shown that there is a maximum disk temperature for the model with a viscosity parameter {\alpha} = 0.5. For the model with {\alpha} = 0.1, no optically thin regions are found to exist for any accretion rate.

We discuss the problem of the formation of a large-scale magnetic field in the accretion disks around black holes, taking into account the non-uniform vertical structure of the disk. The high electrical conductivity of the outer layers of the disk prevents the outward diffusion of the magnetic field. This implies a stationary state with a strong magnetic field in the inner parts of the accretion disk close to the black hole, and zero radial velocity at the surface of the disk. Structure of advective accretion disks is investigated, and conditions for formation of optically thin regions in central parts of the accretion disk are found. The problem of jet collimation by magneto-torsion oscillations is considered.

We study physics of clusters of galaxies embedded in the cosmic dark energy background. Under the assumption that dark energy is described by the cosmological constant, we show that the dynamical effects of dark energy are strong in clusters like the Virgo cluster. Specifically, the key physical parameters of the dark mater halos in clusters are determined by dark energy: 1) the halo cut-off radius is practically, if not exactly, equal to the zero-gravity radius at which the dark matter gravity is balanced by the dark energy antigravity; 2) the halo averaged density is equal to two densities of dark energy; 3) the halo edge (cut-off) density is the dark energy density with a numerical factor of the unity order slightly depending on the halo profile. The cluster gravitational potential well in which the particles of the dark halo (as well as galaxies and intracluster plasma) move is strongly affected by dark energy: the maximum of the potential is located at the zero-gravity radius of the cluster.

The analytic solution is obtained describing kinetic equilibrium of the $\beta$-processes in the nucleonic plasma with relativistic pairs. The nucleons $(n,p)$ are supposed to be non-relativistic and non-degenerate, while the electrons and positrons are ultra-relativistic due to high temperature $(T>6\cdot 10^9$K), or high density $(\rho>\mu 10^6$g/cm$^3$), or both, where $\mu$ is a number of nucleons per one electron. The consideration is simplified because of the analytic connection of the density with the electron chemical potential in the ultra-relativistic plasma, and Gauss representation of Fermi functions. Electron chemical potential and number of nucleons per one initial electron are calculated as functions of $\rho$ and $T$.

We explore a detailed model in which the active galactic nucleus (AGN) obscuration results from the extinction of AGN radiation in a global flow driven by the pressure of infrared radiation on dust grains. We assume that external illumination by UV and soft X-rays of the dusty gas located at approximately 1pc away from the supermassive black hole is followed by a conversion of such radiation into IR. Using 2.5D, time-dependent radiation hydrodynamics simulations in a flux-limited diffusion approximation we find that the external illumination can support a geometrically thick obscuration via outflows driven by infrared radiation pressure in AGN with luminosities greater than $0.05\, L_{\rm edd}$ and Compton optical depth, $\tau_{\rm T}\gtrsim 1$.

Thermal balance of the jet in the source SS433 is considered with account of radiative and adiabatic cooling, and different heating mechanisms. We consider jet heating by the inverse Compton effect of coronal hard X-ray quanta on jet electrons, the influence of shock wave propagation along the jet, and jet kinetic energy transformation into heat via Coulomb collisions of jet and corona protons. The most important heating mechanism for the source SS433 turns out to be Coulomb collisions of jet particles with the surrounding medium.

We construct a radiation-hydrodynamics model for the obscuring toroidal structure in active galactic nuclei. In this model the obscuration is produced at parsec scale by a dense, dusty wind which is supported by infrared radiation pressure on dust grains. To find the distribution of radiation pressure, we numerically solve the 2D radiation transfer problem in a flux limited diffusion approximation. We iteratively couple the solution with calculations of stationary 1D models for the wind, and obtain the z-component of the velocity. Our results demonstrate that for AGN luminosities greater than 0.1 Ledd external illumination can support a geometrically thick obscuration via outflows driven by infrared radiation pressure. The terminal velocity of marginally Compton-thin models (0.2<tau_T<0.6), is comparable to or greater than the escape velocity. In Compton thick models the maximum value of the vertical component of the velocity is lower than the escape velocity, suggesting that a significant part of our torus is in the form of failed wind. The results demonstrate that obscuration via normal or failed infrared-driven winds is a viable option for the AGN torus problem and AGN unification models. Such winds can also provide an important channel for AGN feedback.

We solve the equation of the equilibrium of the gravitating body, with a polytropic equation of state of the matter $P=K\rho^{\gamma}$, with $\gamma=1+1/n$, in the frame of the Newtonian gravity, with non-zero cosmological constant $\Lambda$. We consider the cases with $n=1,\,\,1.5,\,\,3$ and construct series of solutions with a fixed value of $\Lambda$. For each value of $n$, the non-dimensional equation of the static equilibrium has a family of solutions, instead of the unique solution of the Lane-Emden equation at $\Lambda=0$. The equilibrium state exists only for central densities $\rho_0$ larger than the critical value $\rho_c$. There are no static solutions at $\rho_0 < \rho_c$. We find the values of $\rho_c$ for each value of $n$ and show that the presence of dark energy decrease the dynamic stability of the configuration. We apply our results for analyzing the possibility of existence of equilibrium states for cluster of galaxies in the present universe with non-zero $\Lambda$.

In earlier works we pointed out that the disk's surface layers are non-turbulent and thus highly conducting (or non-diffusive) because the hydrodynamic and/or magnetorotational (MRI) instabilities are suppressed high in the disk where the magnetic and radiation pressures are larger than the plasma thermal pressure. Here, we calculate the vertical profiles of the {\it stationary} accretion flows (with radial and azimuthal components), and the profiles of the large-scale, magnetic field taking into account the turbulent viscosity and diffusivity and the fact that the turbulence vanishes at the surface of the disk. Also, here we require that the radial accretion speed be zero at the disk's surface and we assume that the ratio of the turbulent viscosity to the turbulent magnetic diffusivity is of order unity. Thus at the disk's surface there are three boundary conditions. As a result, for a fixed dimensionless viscosity $\alpha$-value, we find that there is a definite relation between the ratio ${\cal R}$ of the accretion power going into magnetic disk winds to the viscous power dissipation and the midplane plasma-$\beta$, which is the ratio of the plasma to magnetic pressure in the disk. For a specific disk model with ${\cal R}$ of order unity we find that the critical value required for a stationary solution is $\beta_c \approx 2.4r/(\alpha h)$, where $h$ the disk's half thickness. For weaker magnetic fields, $\beta > \beta_c$, we argue that the poloidal field will advect outward while for $\beta< \beta_c$ it will advect inward. Alternatively, if the disk wind is negligible (${\cal R} \ll 1$), there are stationary solutions with $\beta \gg \beta_c$.

A hydrodynamic model for steady state, spherically-symmetric winds driven by young stellar clusters with an exponential stellar density distribution is presented. Unlike in most previous calculations, the position of the singular point R_sp, which separates the inner subsonic zone from the outer supersonic flow, is not associated with the star cluster edge, but calculated self-consistently. When the radiative losses of energy are negligible, the transition from the subsonic to the supersonic flow occurs always at R_sp ~ 4 R_c, where R_c is the characteristic scale for the stellar density distribution, irrespective of other star cluster parameters. This is not the case in the catastrophic cooling regime, when the temperature drops abruptly at a short distance from the star cluster center and the transition from the subsonic to the supersonic regime occurs at a much smaller distance from the star cluster center. The impact from the major star cluster parameters to the wind inner structure is thoroughly discussed. Particular attention is paid to the effects which radiative cooling provides to the flow. The results of the calculations for a set of input parameters, which lead to different hydrodynamic regimes, are presented and compared to the results from non-radiative 1D numerical simulations and to those from calculations with a homogeneous stellar mass distribution.

We investigate dynamics of a jet collimated by magneto-torsional oscillations. The problem is reduced to an ordinary differential equation containing a singularity and depending on a parameter. We find a parameter range for which this system has stable periodic solutions and study bifurcations of these solutions. We use Poincar\'e sections to demonstrate existence of domains of regular and chaotic motions. We investigate transition from periodic to chaotic solutions through a sequence of period doublings.

We discuss the problem of the formation of a large-scale magnetic field in the accretion disks around black holes, taking into account the nonuniform vertical structure of the disk. The high electrical conductivity of the outer layers of the disk prevents the outward diffusion of the magnetic field. This implies a stationary state with a strong magnetic field in the inner parts of the accretion disk close to the black hole, and zero radial velocity at the surface of the disk. Magnetic jet collimation is considered, when the jet radius is hold due to magneto-torsional oscillations. The range of parameters is found where jet radius is oscillating, regularly or chaotically, within restricted values.

We analyze uncertainties in the cosmic X-ray background measurements performed by the INTEGRAL observatory. We find that the most important effect limiting the accuracy of the measurements is related to the intrinsic background variation in detectors. Taking into account all of the uncertainties arising during the measurements we conclude that the X-ray background intensity obtained in the INTEGRAL observations is compatible with the historic X-ray background observations performed by the HEAO-1 satellite.

We develop a model of gravitational lensing in a non-uniform plasma. When a gravitating body is surrounded by a plasma, the lensing angle depends on the frequency of the electromagnetic wave, due to dispersion properties of plasma, in presence of a plasma inhomogeneity, and of a gravity. The second effect leads, even in a uniform plasma, to a difference of the gravitational photon deflection angle from the vacuum case, and to its dependence on the photon frequency. We take into account both effects, and derive the expression for the lensing angle in the case of a strongly nonuniform plasma in presence of the gravitation. Dependence of the lensing angle on the photon frequency in a homogeneous plasma resembles the properties of a refractive prism spectrometer, which strongest action is for very long radiowaves. We discuss the observational appearances of this effect for the gravitational lens with a Schwarzschild metric, surrounded by a uniform plasma. We obtain formulae for the lensing angle and the magnification factors in this case and discuss a possibility of observation of this effect by the planned VLBI space project Radioastron. We also consider models with a nonuniform plasma distribution. For different gravitational lens models we compare the corrections to the vacuum lensing due to the gravitational effect in plasma, and due to the plasma inhomogeneity. We have shown that the gravitational effect could be detected in the case of a hot gas in the gravitational field of a galaxy cluster.

Activity of the nuclei of galaxies and stellar mass systems involving disk accretion to black holes is thought to be due to (1) a small-scale turbulent magnetic field in the disk (due to the magneto-rotational instability or MRI) which gives a large viscosity enhancing accretion, and (2) a large-scale magnetic field which gives rise to matter outflows and/or electromagnetic jets from the disk which also enhances accretion. An important problem with this picture is that the enhanced viscosity is accompanied by an enhanced magnetic diffusivity which acts to prevent the build up of a significant large-scale field. Recent work has pointed out that the disk's surface layers are non-turbulent and thus highly conducting (or non-diffusive) because the MRI is suppressed high in the disk where the magnetic and radiation pressures are larger than the thermal pressure. Here, we calculate the vertical ($z$) profiles of the stationary accretion flows (with radial and azimuthal components), and the profiles of the large-scale, magnetic field taking into account the turbulent viscosity and diffusivity due to the MRI and the fact that the turbulence vanishes at the surface of the disk. The stationary solutions we find indicate that a weak ($\beta > 1$) large-scale field does not diffuse away as suggested by earlier work.

We develop a Monte-Carlo technique based on L.B. Lucy's indivisible photon packets method to calculate X-ray continuum spectra of comptonized thermal plasma in arbitrary geometry and apply it to describe the broadband X-ray continuum of the galactic superaccreting microquasar SS433 observed by INTEGRAL. A physical model of the X-ray emitting region is proposed that includes thermal emission from the accretion disk, jets and hot corona where the photons of different origin are comptonized. From comparison with INTEGRAL observations, we estimate physical parameters of the complex X-ray emitting region in SS433 and present model spectra for different viewing angles of the object.

Gravitational lensing is predicted by general relativity and is found in observations. When a gravitating body is surrounded by a plasma, the lensing angle depends on a frequency of the electromagnetic wave due to refraction properties, and the dispersion properties of the light propagation in plasma. The last effect leads to dependence, even in the uniform plasma, of the lensing angle on the frequency, what resembles the properties of the refractive prism spectrometer. The strongest action of this spectrometer is for the frequencies slightly exceeding the plasma frequency, what corresponds to very long radiowaves.

Gravitational lensing by gravitational wave is considered. We notice that although final and initial direction of photons coincide, displacement between final and initial trajectories occurs. This displacement is calculated analytically for the plane gravitational wave pulse. Estimations for observations are discussed.

The motion with intersections of relativistic gravitating shells in the Schwarzschild gravitational field of a central body is considered. Formulas are derived for calculating parameters of the shells after intersection via their parameters before intersection. Such special cases as the Newtonian approximation, intersections of light shells, and intersections of a test shell with a gravitating shell are also considered. The ejection of one of the shells to infinity in the relativistic region is described. The equations of motion for the shells are analyzed numerically.

The formation of a neutrino pulse emitted during the relativistic collapse of a spherical supermassive star is considered. The free collapse of a body with uniform density in the absence of rotation and with the free escape of the emitted neutrinos can be solved analytically by quadrature. The light curve of the collapsing star and the spectrum of the emitted neutrinos at various times are calculated.

The properties of the relativistic rings which show up in images of a source when a black hole lies between the source and observer are examined. The impact parameters are calculated, along with the distances of closest approach of the rays which form a relativistic ring, their angular sizes, and their "magnification" factors, which are much less than unity.

The evolution of the primordial low mass black hole (PBH) in hot universe is considered. Increase of mass and decrease of PBH spin due to the accretion of radiation dominated matter are estimated with using of results of numerical simulation of PBH formation and approximate relations for accretion to a rotating black hole.

The following problems are discussed. 1. Pulsars and close binaries. 2. Hulse-Taylor pulsar. 3. Disrupted pulsar pairs. 4. RP statistics. 5. Enhanced evaporation: formation of single RP. 6. General relativity effects: NS+NS. 7. A Double pulsar system. 8. Checking general relativity. 9. Variability of the gravitational constant. 10. Space Watch.

We discuss results of 2D simulations of magnetorotational(MR) mechanism of core collapse supernova explosions. Due to the nonuniform collapse the collapsed core rotates differentially. In the presence of initial poloidal magnetic field its toroidal component appears and grows with time. Increased magnetic pressure leads to foramtion of compression wave which moves outwards. It transforms into the fast MHD shock wave (supernova shock wave). The shape of the MR supernova explosion qualitatively depends on the configuration of the initial magnetic field. For the dipole-like initial magnetic field the supernova explosion develops mainly along rotational axis forming mildly collimated jet. Quadrupole-like initial magnetic field leads to the explosion developing mainly along equatorial plane. Magnetorotational instability was found in our simulations. The supernova explosion energy is growing with increase of the initial core mass and rotational energy of the core, and corresponds to the observational data.

We solve equations, describing in a simplified way the newtonian dynamics of a selfgravitating nonrotating spheroidal body after loss of stability. We find that contraction to a singularity happens only in a pure spherical collapse, and deviations from the spherical symmetry stop the contraction by the stabilising action of nonlinear nonspherical oscillations. A real collapse happens after damping of the oscillations due to energy losses, shock wave formation or viscosity. Detailed analysis of the nonlinear oscillations is performed using a Poincar\'{e} map construction. Regions of regular and chaotic oscillations are localized on this map.

We reconsider the problem of the formation of a large-scale magnetic field in the accretion disks around black holes. In contrast with previous work we take into account the nonuniform vertical structure of the disk. The high electrical conductivity of the outer layers of the disk prevents the outward diffusion of the magnetic field. This implies a stationary state with a strong magnetic field in the inner parts of the accretion disk close to the black hole.

Many quasars and active galactic nuclei (AGN) appear in radio, optical, and X-ray maps, as a bright nuclear sources from which emerge single or double long, thin jets. When observed with high angular resolution these jets show structure with bright knots separated by relatively dark regions. Nonthermal nature of a jet radiation is well explained as the synchrotron radiation of the relativistic electrons in an ordered magnetic field. We consider magnetic collimation, connected with torsional oscillations of a cylinder with elongated magnetic field, and periodically distributed initial rotation around the cylinder axis. The stabilizing azimuthal magnetic field is created here by torsional oscillations, where charge separation is not necessary. Approximate simplified model is developed. Ordinary differential equation is derived, and solved numerically, what gives a possibility to estimate quantitatively the range of parameters where jets may be stabilized by torsional oscillations.

The present common view about GRB origin is related to cosmology, what is based on statistical analysis, and on measurements of the redshifts in the GRB optical afterglows of long GRB. Models of radiation, and models of the central machine for GRB are critically analyzed. Most of these models cannot produce the amount of energy required for the cosmological GRB. Collimation is needed to reduce GRB energy output. No correlation is found between redshifts, GRB spectrum, and total GRB fluence. Comparison of KONUS and BATSE data about statistics and hard X-ray lines is done, and some differences are noted. Hard gamma-ray afterglows, prompt optical spectra, hard X-ray lines measurements could be important for farther insight into GRB origin. Two bright optical afterglows had been detected, indicating that initial bright optical emission is directly connected with the GRB itself, and has a different origin from the subsequent faint optical emission. Results of these observations are represented. Observations indicate that GRB originate in star forming regions of galaxies with a high matter density. Interaction of the mighty gamma ray flux from the cosmological GRB with the dense surrounding molecular cloud lead to long optical afterglow, lasting up to 10 years, due to reradiation of the cloud. Results of numerical simulations of such reradiation for different forms of clouds are represented. Possible connection of short GRB with soft gamma repeaters is discussed.

Many quasars and active galactic nuclei (AGN) appear in radio, optical, and X-ray maps, as a bright nuclear sources from which emerge single or double long, thin jets (Thomson et al., 1993). When observed with high angular resolution these jets show structure with bright knots separated by relatively dark regions. High percentages of polarization, sometimes more then 50% in some objects, indicates the nonthermal nature of the radiation which is well explained as the synchrotron radiation of the relativistic electrons in an ordered magnetic field. A strong collimation of jets is most probably connected with ordered magnetic fields. The mechanism of magnetic collimation, first suggested by Bisnovatyi-Kogan et al. (1969), was based on the initial charge separation, leading to creation of oscillating electrical current, which produces azimuthal magnetic field, preventing jet expansion and disappearance. Here we consider magnetic collimation, connected with torsional oscillations of a cylinder with elongated magnetic field. Instead of initial blobs with charge separation, we consider a cylinder with a periodically distributed initial rotation around the cylinder axis. The stabilizing azimuthal magnetic field is created here by torsional oscillations, where charge separation is not necessary. Approximate simplified model is developed. Ordinary differential equation is derived, and solved numerically, what gives a possibility to estimate quantitatively the range of parameters where jets may be stabilized by torsional oscillations.

Binary radio pulsars, first discovered by Hulse and Taylor in 1974 [1], are a unique tool for experimentally testing general relativity (GR), whose validity has been confirmed with a precision unavailable in laboratory experiments. In particular, indirect evidence of the existence of gravitational waves has been obtained. Radio pulsars in binary systems (which have come to be known as recycled) have completed the accretion stage, during which neutron star spins reach millisecond periods and their magnetic fields decay 2 to 4 orders of magnitude more weakly than ordinary radio pulsars. Among about a hundred known recycled pulsars, many have turned out to be single neutron stars. The high concentration of single recycled pulsars in globular clusters suggests that close stellar encounters are highly instrumental in the loss of the companion. A system of one recycled pulsar and one 'normal' one discovered in 2004 is the most compact among binaries containing recycled pulsars [2]. Together with the presence of two pulsars in one system, this suggests new prospects for further essential improvements in testing GR. This paper considers theoretical predictions of binary pulsars, their evolutionary formation, and mechanisms by which their companions may be lost. The use of recycled pulsars in testing GR is discussed and their possible relation to the most intriguing objects in the universe, cosmic gamma-ray bursts, is examined.

The typical spectra of gamma-ray bursts (GRBs) are discussed in the context of the compactness problem for GRB sources and how it is resolved in the popular fireball model. In particular, observational (model-independent) constraints on the collimation of the gamma-rays and the dependence of the collimation angle on the photon energy are considered. The fact that the threshold for the creation of $e^{-}e^{+}$ pairs depends on the angle between the momenta of the annihilating photons in the GRB source provides an alternative solution to the compactness problem. A new approach to explaining GRBs, taking into account the angular dependence for pair creation, is proposed, and the main features of a scenario describing a GRB source with a total (photon) energy smaller or of the order of $10^{49}$ erg are laid out. Thus, we are dealing with an alternative to an ultra-relativistic fireball, if it turns out (as follows from observations) that all "long" GRBs are associated with normal (not peculiar) core-collapse supernovae. The effects of radiation pressure and the formation of jets as a consequence of even a small amount of anisotropy in the total radiation field in a (compact) GRB source are examined in this alternative model. Possible energy release mechanisms acting in regions smaller or of the order of $10^{8}$ cm in size (a compact model for a GRB) are discussed. New observational evidence for such compact energy release in the burst source is considered.

We present results of 2D simulation of magnetorotational (MR) supernova accompanied by jet formation in the core collapse supernova explosion. Initial magnetic field used in the simulations has dipole-like symmetry. Contrary to the simulations of MR supernova with initial quadrupole-like magnetic field, where the matter was ejected mainly near the equatorial plane, in presence of the dipole-like initial magnetic field the supernova explosion is developing preferably along a rotational axis, and leads to formation of a protojet. We expect that protojet propagation through the envelope of the star will be accompanied by its collimation. The magnetorotational instability (MRI) was found in simulations, similar to the earlier considered case of the quadrupole-like initial magnetic field. Our estimations show that the characteristic time for the reconnection of the magnetic field is much larger than the MRI development time. The supernova explosion energy for the dipole-like field is about $0.61\cdot 10^{51}$erg, and about $0.13M_\odot$ of mass was ejected during the explosion.

Typical observational gamma-ray burst (GRB) spectra are discussed and, in this connection, what is the origin of the compactness problem and how it was solved at first. If the threshold for $e^{-}e^{+}$ pair production depends on an angle between photon momenta, then another solution of the compactness problem is possible. We discuss a possibility of the $\gamma$-rays collimation and the dependence of photon beaming on photon energies. The list of basic assumptions of the scenario describing the GRB source with energy $< 10^{49}$ ergs are adduced: the matter is about an alternative to the ultrarelativistic fireball if {\it all} long-duration GRBs are related or physically connected with {\it normal/unpeculiar} core-collapse supernovae (SNe). Namely, we consider the questions about radiation pressure and how the jet arises on account of even small asymmetry of the radiation field in a compact GRB source. The possibility of a new approach to explanation of GRB phenomenon is shown. Possible mechanisms of their generation in regions of size $< 10^8$cm are discussed (a compact model of GRBs). Observational consequences of the compact GRB energy release are considered.

We derive a variational principle for the dynamical stability of a cluster as a gas sphere in a box. Newtonian clusters are always dynamically stable and, for relativistic clusters, the relation between dynamical and thermodynamical instabilities is analyzed. The boundaries between dynamically and thermodynamically stable and unstable models are found numerically for relativistic stellar systems with different cut off parameters. A criterion based on binding energy curve is used for determination of the boundary of dynamical stability.

The most precise measurements are done at present by timing of radiopulsars in binary systems with two neutron stars. The timing measurements of the Taylor-Hulse pulsar B1913+16 gave the most precise results on testing of general relativity (GR), finding implicit proof of existence of gravitational waves. We show that available results of existing measurements, obtained to the year 1993, in combination with the results of the Mariner 10 in (1992), give the boundaries for the variation of the gravitational constant ${\dot G}/{G}$ inside the limits $(-0.6 \div +2)\cdot 10^{-12}$ year{$^{-1}$}.

The collapse of non-collisional dark matter and the formation of pancake structures in the Universe are investigated approximately. Collapse is described by a system of ordinary differential equations, in the model of a uniformly rotating, three-axis, uniform density ellipsoid. Violent relaxation, mass and angular momentum losses are taken into account phenomenologically. The formation of the equilibrium configuration, secular instability and the transition from a spheroid to a three-axis ellipsoid are investigated numerically and analytically in this dynamical model.

Core-collapse supernovae are connected with formation of neutron stars. Part of the gravitation energy is transformed into the energy of the explosion, observed in SN II, SN Ib,c type supernovae. The mechanism of transformation is not simple, because the overwhelming majority of the energy is going into weakly interacting neutrino. The attempts to use this energy for the explosion were not successful during about 40 years of investigation. We consider the explosion mechanism in which the source of energy is the rotation, and magnetic field serves for the transformation of the rotation energy into the energy of explosion. 2-D MHD simulations of this mechanism were performed. After the collapse the core consists of a rapidly rotating proto-neutron star with a differentially rotating envelope. The toroidal part of the magnetic energy generated by the differential rotation grows as quadratic function with time at the initial stage of the evolution of the magnetic field. The linear growth of the toroidal magnetic field is terminated by the development of magnetohydrodynamic instability, when the twisted toroidal component strongly exceeds the poloidal field, leading to a drastic acceleration in the growth of magnetic energy. At the moment when the magnetic pressure becomes comparable to the gas pressure at the periphery of the proto-neutron star the MHD compression wave appears and goes through the envelope of the collapsed core. It transforms into the fast MHD shock and produces a supernova explosion. Our simulations give the energy of the explosion $0.6\cdot 10^{51}$ ergs. The amount of the mass ejected by the explosion is $\sim 0.14M_\odot$. The implicit numerical method, based on the Lagrangian triangular grid of variable structure, was used for the simulations.

An analytical model is developed for the screening of the external magnetic field of a rotating, axisymmetric neutron star due to the accretion of plasma from a disk. The decrease of the field occurs due to the electric current in the infalling plasma. The deposition of this current carrying plasma on the star's surface creates an induced magnetic moment with a sign opposite to that of the original magnetic dipole. The field decreases independent of whether the star spins-up or spins-down. The time-scale for an appreciable decrease (factor of $>100$) of the field is found to be $\sim 1.6 \times 10^7$ yr, for a mass accretion rate $\dot{M}=10^{-9} M_\odot/$yr and an initial magnetic moment $\mu_i = 10^{30}{\rm G cm}^3$ which corresponds to a surface field of $10^{12} $G if the star's radius is $10^6$ cm. The time-scale varies approximately as $\mu_i^{3.8}/\dot{M}^{1.9}$. The decrease of the magnetic field does not have a simple relation to the accreted mass. Once the accretion stops the field leaks out on an Ohmic diffusion time scale which is estimated to be $ > 10^9$ yr.

Observations of gamma ray bursts (GRB) afterglows in different spectral bands provide a most valuable information about their nature, as well as about properties of surrounding medium. Powerful infrared afterglow was observed from the strong GRB041219. Here we explain the observed IR afterglow in the model of a dust reradiation of the main GRB signal in the envelope surrounding the GRB source. In this model we do not expect appearance of the prompt optical emission which should be absorbed in the dust envelope. We estimate the collimation angle of the gamma ray emission, and obtain restrictions on the redshift (distance to GRB source), by fitting the model parameters to the observational data.

We present the results of 2D simulations of the magnetorotational model of a supernova explosion. After the core collapse the core consists of rapidly a rotating proto-neutron star and a differentially rotating envelope. The toroidal part of the magnetic energy generated by the differential rotation grows linearly with time at the initial stage of the evolution of the magnetic field. The linear growth of the toroidal magnetic field is terminated by the development of magnetohydrodynamic instability, leading to drastic acceleration in the growth of magnetic energy. At the moment when the magnetic pressure becomes comparable with the gas pressure at the periphery of the proto-neutron star $\sim 10-15$km from the star centre the MHD compression wave appears and goes through the envelope of the collapsed iron core. It transforms soon to the fast MHD shock and produces a supernova explosion. Our simulations give the energy of the explosion $0.6\cdot 10^{51}$ ergs. The amount of the mass ejected by the explosion is $\sim 0.14M_\odot$. The implicit numerical method, based on the Lagrangian triangular grid of variable structure, was used for the simulations.

A qualitative study is made of the evolutionary tracks of various dense clusters, with allowance for the evaporation of stars, and for head-on collisions, until relativistic collapse occures or the cluster breaks up.

The collapse of a rotating 3-axis ellipsoid is approximated by a system of ordinary differential equations. Violent relaxation, mass and angular momentum losses are taken into account phenomenologically. The formation of the equilibrium configuration and different types of instability are investigated.

We present 2D results of simulations of the magnetorotational core collapsed supernova. For the first time we obtain strong explosion for the core collapsed supernova. In 2D approximation we show that amplification of the toroidal magnetic field due to the differential rotation leads to the formation of MHD shockwave, which produces supernova explosion. The amounts of the ejected mass $0.1M_\odot$ and energy $\sim 0.5\div0.6 \cdot 10^{51}$ergs can explain the energy output for supernova type II or type Ib/c explosions. The shape of the explosion is qualitatively depends on the initial configuration of the magnetic field, and may form strong ejection neat the equatorial plane, or produce mildly collimated jets. Our simulation show that during the evolution of the magnetic field the magnetorotational instability appears and leads to exponential growth of the magnetic field strength.

We consider the effects of advection and radial gradients of pressure and radial drift velocity on the structure of accretion disks around black holes with proper description of optically thick/thin transitions. We concentrated our efforts on the models with large accretion rate. Contrary to disk models neglecting advection, we find that continuous solutions extending from the outer disk regions to the inner edge exist for all accretion rates we have considered. We show that the sonic point moves outward with increasing accretion rate, and that in the innermost disk region advection acts as a heating process that may even dominate over dissipative heating. Despite the importance of advection on it's structure, the disk remains geometrically thin. Global solutions of advective accretion disks, which describe continuously the transition between optically thick outer region and optically thin inner region are constructed and analyzed.

It is shown that during the motion of two initially gravitationally bound spherical shells, consisting of point particles moving along ballistic trajectories, one of the shell may be expelled to infinity at subrelativistic speed $v_{exp}\leq 0.25 c$. The problem is solved in Newtonian gravity. Motion of two intersecting shells in the case when they do not runaway shows a chaotic behaviour. We hope that this toy and oversimplyfied model can give nevertheless a qualitative idea on the nature of the mechanism of matter outbursts from the dense stellar clusters.

Interaction of cosmological gamma ray burst radiation with the dense interstellar medium of host galaxy is considered. Gas dynamical motion of interstellar medium driven by gamma ray burst is investigated in 2D approximation for different initial density distributions of host galaxy matter and different total energy of gamma ray burst. The maximum velocity of motion of interstellar medium is $1.8\cdot10^4$ km/s. Light curves of gamma ray burst afterglow are calculated for set of non homogeneous density, distribution gamma ray burst total energy, and different viewing angles. Spectra of gamma ray burst afterglow are modeled taking into account conversion of hard photons (soft X-ray, hard UV) to soft UV and optics photons.

Exact solution is obtained for electromagnetic field around a conducting cylinder of infinite length and finite radius, with a periodical axial current, when the wave length is much larger than the radius of the cylinder. The solution describes simultaneously the fields in the near zone close to the cylinder, and transition to the wave zone. Proper long-wave oscillations of such cylinder are studied. The electromagnetic energy flux from the cylinder is calculated. These solutions could be applied for description of the electromagnetic field around relativistic jets from active galactic nuclei and quasars and particle acceleration inside jets.

Astrophysical sources of high frequency gravitational radiation are considered in association with a new interest to very sensitive HFGW receivers required for the laboratory GW Hertz experiment. A special attention is paid to the phenomenon of primordial black holes evaporation. They act like black body to all kinds of radiation, including gravitons, and, therefore, emit an equilibrium spectrum of gravitons during its evaporation. Limit on the density of high frequency gravitons in the Universe is obtained, and possibilities of their detection are briefly discussed.

Accretion disks are observed in young stars, cataclysmic variables, binary X-ray sources et al. Accretion disk theory was first developed as a theory with the local heat balance, where the whole energy produced by a viscous heating was emitted to the sides of the disk. Important part of this theory was the phenomenological treatment of the turbulent viscosity, known the `` alpha'' prescription, where the $(r \phi)$ component of the stress tensor was connected with the pressure as $\alpha P$. Sources of turbulence in the accretion disk are discussed, including hydrodynamic turbulence, convection and magnetic field role. Optically thin solution and advective disks are considered. Related problems of mass ejection from magnetized accretion disks and jet formation are discussed.

The present common view about GRB origin is related to cosmology, what is based on statistical analysis, and on measurements of the redshifts in the GRB optical afterglows of long GRB. No correlation is found between redshifts, GRB spectrum, and total GRB fluence. Comparison of KONUS and BATSE data about statistics and hard X-ray lines is done, and some differences are noted. Hard gamma-ray afterglows, prompt optical spectra, hard X-ray lines measurements could be important for farther insight into GRB origin. Possible possible connection of short GRB with soft gamma repeaters is discussed.

After some historical remarks we discuss different criteria of dynamical stability of stars, and properties of the critical states where dynamical stability is lost, leading to collapse with formation of the neutron star or a black hole. At the end some observational and theoretical problems related to quark stars are discussed.

The present common view about GRB origin is related to cosmology. There are two evidences in favor of this interpretation. The first is connected with statistics, the second is based on measurements of the redshifts in the GRB optical afterglows. Redshifts in optical afterglows had been observed only in long GRB. Comparison of KONUS and BATSE data about statistics and hard X-ray lines is done, and some differences are noted. Hard gamma-ray afterglows, prompt optical spectra and polarization measurements could be very important for farther insight into GRB origin. Possibility of galactic origin of short GRB is discussed as well as their possible connection with soft gamma repeaters.

Results of 2D simulations of the magnetorotational mechanism of supernova type II are presented. Amplification of toroidal magnetic field of the star due to differential rotation of the star leads to the transformation of the rotational (gravitational) energy to the energy of the supernova explosion. In our simulation the energy of the explosion is $1.12 \cdot 10^{51}erg$. The explosion ejects about $0.11 M_\odot$.

Many phenomena such as stellar variability, stellar explosions, behavior of different kind of X-ray and gamma-ray sources, processes in active galactic nuclei and other astrophysical objects show stochastic features. Brief description of these phenomena is given in this review.

Collapse of the rotating spheroid is approximated by a system of ordinary differential equations describing its dynamics. The gravitational potential is approximated by the one of the iniform Maclaurin spheroid. Developement of gravitational instability and collapse in the dark matter medium do not lead to any shock formation or radiation, but is characterized by non-collisional relaxation, which is accompanied by the mass and angular momentum losses. Phenomenological account of these processes is done in this model. Formation of the equilibrium configuration dynamics of collapse is investigated for several parameters, characterizing the configuration. A very long gravitational wave emission during the collapse is estimated, and their possible connection with the observed gravitational lenses is discussed.

Slow magnetic field variations in stars and planets create a quasistationary electrical field which may be observed. It is supposed that the electrical field near the Earth surface may be partially connected with variation of the Earth magnetic field. Two examples of the electrical field distribution around the infinite cylinder, and the circular loop with a lineary growing with time electrical currents are given.

Accretion disk theory was first developed as a theory with the local heat balance, where the whole energy produced by a viscous heating was emitted to the sides of the disk. One of the most important new invention of this theory was the phenomenological treatment of the turbulent viscosity, known the `` alpha'' prescription, where the (r phi) component of the stress tensor was approximated by (alpha p) with a unknown constant alpha. This prescription played the role in the accretion disk theory as well important as the mixing-length theory of convection for stellar evolution. Sources of turbulence in the accretion disk are discussed, including hydrodynamical turbulence, convection and magnetic field role. In parallel to the optically thick geometrically thin accretion disk models, a new branch of the optically thin accretion disk models was discovered, with a larger thickness for the same total luminosity. The choice between these solutions should be done of the base of a stability analysis. The ideas underlying the necessity to include advection into the accretion disk theory are presented and first models with advection are reviewed. The present status of the solution for a low-luminous optically thin accretion disk model with advection is discussed and the limits for an advection dominated accretion flows (ADAF) imposed by the presence of magnetic field are analysed. Related problems of mass ejection from accretion disks and jet formation are discussed.

We revisit the problem of magnetic field generation in accretion flows onto black holes owing to the excess radiation force on electrons. This excess force may arise from the Poynting-Robertson effect. Instead of a recent claim of the generation of dynamically important magnetic fields, we establish the validity of earlier results from 1977 which show only small magnetic fields are generated. The radiative force causes the magnetic field to initially grow linearly with time. However, this linear growth holds for only a {\it restricted} time interval which is of the order of the accretion time of the matter. The large magnetic fields recently found result from the fact that the linear growth is unrestricted. A model of the Poynting-Robertson magnetic field generation close to the horizon of a Schwarzschild black hole is solved exactly using General Relativity, and the field is also found to be dynamically insignificant. These weak magnetic fields may however be important as seed fields for dynamos.

Jet formation is connected most probably with matter acceleration from the vicinity of rotating magnetized bodies. It is usually related to the mass outflows and ejection from accretion disks around black holes. Problem of jet collimation is discussed. Collapse of a rotating magnetized body during star formation or supernovae explosion may lead to a jet-like mass ejection for certain angular velocity and magnetic field distributions at the beginning of the collapse. Jet formation during magnetorotational explosion is discussed basing on the numerical simulation of collapse of magnetized bodied with quasi-dipole field.

The present common view about GRB origin is related to cosmology. There are two evidences in favour of this interpretation. The first is connected with statistics, the second is based on measurements of the redshifts in the GRB optical afterglows. Red shifts in optical afterglows had been observed only in long GRB. Statistical errors, and possibility of galactic origin of short GRB is discussed; their connection with Soft Gamma Repeaters (SGR) is analyzed.

It is shown that in presence of URCA shell pulsational energy losses due to neutrino emission and nonequilibrium beta heating are much less than energy losses by excitation of short-wavelength acoustic waves. Convective motion in presence of URCA shell is considered, and equations generalizing the mean free path model of convection are derived.

Recently, there appeared many papers devoted to the modeling of X-ray properties of Cygnus X-1 and other black hole accretion disk candidates: e.g. J.Poutanen, J.Krolik, F.Ryde, MNRAS 292 (1997) L21; E.Agol, J.Krolik, ApJ 507 (1998) 304; A.Beloborodov, in ``High Energy Processes in Accreting Black Holes'', ASP Conf. Series, 161, (1999), p.295. The goal of this electronic publication is to draw attention to our old papers where many ideas of recent discussions were anticipated (hot coronae, Comptonization, photon damping of waves, particle acceleration and matter ejection from accretion disks with large scale poloidal magnetic fields, etc.)

Global solutions of optically thick advective accretion disks around black holes are constructed. The solutions are obtained by solving numerically a set of ordinary differential equations corresponding to a steady axisymmetric geometrically thin disk. We pay special attention to consistently satisfy the regularity conditions at singular points of the equations. For this reason we analytically expand a solution at the singular point, and use coefficients of the expansion in our iterative numerical procedure. We obtain consistent transonic solutions in a wide range of values of the viscosity parameter alpha and mass acretion rate. We compare two different form of viscosity: one takes the shear stress to be proportional to the pressure, while the other uses the angular velocity gradient-dependent stress. We find that there are two singular points in solutions corresponding to the pressure-proportional shear stress. The inner singular point locates close to the last stable orbit around black hole. This point changes its type from a saddle to node depending on values of alpha and accretion rate. The outer singular point locates at larger radius and is always of a saddle-type. We argue that, contrary to the previous investigations, a nodal-type inner singular point does not introduce multiple solutions. Only one integral curve, which corresponds to the unique global solution, passes simultaneously the inner and outer singular points independently of the type of inner singular point. Solutions with the angular velocity gradient-dependent shear stress have one singular point which is always of a saddle-type and corresponds to the unique global solution. The structure of accretion disks corresponding to both viscosities are similar.

We generalize a formula for the Goldreich-Julian charge density (\rho), originally derived for rotating neutron star, for arbitrary oscillations of a neutron star with arbitrary magnetic field configuration under assumption of low current density in the inner parts of the magnetosphere. As an application we consider toroidal oscillation of a neutron star with dipole magnetic field and calculate energy losses. For some oscillation modes the longitudinal electric field can not be canceled by putting charged particles in the magnetosphere without a presence of strong electric current j = (c/(\omega r))\rho c. It is shown that the energy losses are strongly affected by plasma in the magnetosphere, and cannot be described by vacuum formulas.

We perform 2D numerical simulations of a magnetorotational explosion of a rotating magnetized gas cloud. We found that amplification of a toroidal magnetic field due to the differential rotation leads to a transformation of the part of the rotational energy of the cloud to the radial kinetic energy. Simulations have been made for 3 initial values of $\xi$ (the relation of magnetic energy to the gravitational energy of the cloud): $\xi =10^{-2},10^{-4},10^{-6}$. Part of the matter - $\sim 7%$ of the mass of the cloud ($\sim 3.3%$ of the final gravitational energy of the cloud) - gets radial kinetic energy which is larger than its potential energy and can be thrown away to the infinity. It carries about 30% of the initial angular momentum of the cloud. This effect is important for angular momentum loss in the processes of stellar formation, and for the magnetorotational mechanism of explosion suggested for supernovae. Simulations have been made on the basis of the Lagrangian 2D numerical implicit scheme on a triangular grid with grid reconstruction.

Estimations of magnetic fields of neutron stars, observed as radio and X-ray pulsars, are discussed. It is shown, that theoretical and observational values for different types of radiopulsars are in good correspondence. Magnetic fields of X-ray pulsars are estimated from the cyclotron line energy. In the case of Her X-1 this estimation exceeds considerably the value of its magnetic field obtained from long term observational data related to the beam structure evolution. Another interpretation of the cyclotron feature, based on the relativistic dipole radiation mechanism, could remove this discrepancy. Observational data about soft gamma repeators and their interpretation as magnetars are analyzed.

Standard accretion disk theory is formulated which is based on the local heat balance. The energy produced by a turbulent viscous heating is supposed to be emitted to the sides of the disc. Sources of turbulence in the accretion disc are connected with nonlinear hydrodynamic instability, convection, and magnetic field. In standard theory there are two branches of solution, optically thick, and optically thin. Advection in accretion disks is described by the differential equations what makes the theory nonlocal. Low-luminous optically thin accretion disc model with advection at some suggestions may become advectively dominated, carrying almost all the energy inside the black hole. The proper account of magnetic filed in the process of accretion limits the energy advected into a black hole, efficiency of accretion should exceed $\sim 1/4$ of the standard accretion disk model efficiency.

The spectrum of cyclotron radiation is calculated for anisotropically distributed relativistic electrons with a nonrelativistic velocity scattering across the magnetic field. It is shown that if such electrons are responsible for a formation of the "cyclotron" line in the spectrum of Her X-1, then the value of its magnetic field $(3-6)\cdot 10^{10}$ Gs resulting from this interpretation is in a good agreement with some other observations and theoretical estimations. Observations of a time dependence of the energy of this "cyclotron" line in the spectra of several X-ray pulsars is explained by a variability of the average longitude energy of the electrons, decreasing with increasing liminocity due to radiational braking of the accretion flow.

We investigated influence of cosmological GRB on the surrounding interstellar medium. It was shown that \gamma-radiation from the burst heats interstellar gas to the temperatures > 10^4 K up to the distance \sim 10 pc. For high density ISM optical and UV radiation of the heated gas can be observed on the Eath several years as a GRB`s counterpart

The launch of the Beppo-Sax satellite gave a unique opportunity to investigate gamma ray bursts (GRB) in different spectral regions. The large diversity of the afterglow behavior creates additional problems for the cosmological model with a fireball. Formation of the afterglow giving the observed diversity of properties is suggested in the Galactic neutron star model of GRBs. It is based on the transient accretion disc formation around the neutron star with a low-mass brown companion irradiated by the neutron star.

Recent papers discussing advection dominated accretion flows (ADAF) as a solution for astrophysical accretion problems should be treated with some caution because of their uncertain physical basis. The suggestions underlying ADAF involve ignoring the magnetic field reconnection in heating of the plasma flow, assuming electron heating due only to binary Coulomb collisions with ions. Here, we analyze the physical processes in optically thin accretion flows at low accretion rates including the influence of an equipartition turbulent magnetic field. For these conditions there is continuous destruction of magnetic flux by reconnection. The reconnection is expected to significantly heat the electrons which can efficiently emit magnetobremstrahlung radiation. Because of this electron emission, the radiative efficiency of the ADAF is not small. We suggest that the small luminosities of nearby galactic black holes is due to outflows rather than ADAF accretion.

A model is developed for magnetic `propeller'-driven outflows which cause a rapidly rotating magnetized star accreting from a disk to spin-down. Energy and angular momentum lost by the star goes into expelling most of the accreting disk matter. The theory gives an expression for the effective Alfven radius $R_A$ (where the inflowing matter is effectively stopped) which depends on the mass accretion rate, the star's mass and magnetic moment, and the star's rotation rate. The model points to a mechanism for `jumps' between spin-down and spin-up evolution and for the reverse transition, which are changes between two possible equilibrium configurations of the system. In for example the transistion from spin-down to spin-up states the Alfven radius $R_A$ decreases from a value larger than the corotation radius to one which is smaller. In this transistion the `propeller' goes from being ``on'' to being ``off.'' The ratio of the spin-down to spin-up torque (or the ratio for the reverse change) in a jump is shown to be of order unity.

We study axisymmetric accretion of slowly rotating matter onto a gravitation center. The process of the flow restructuring is presented as the angular velocity of accreting gas approaches the Keplerian velocity. For spherically-symmetric accretion we found conditions of existence of steady-state supersonic solutions for various external boundary conditions. These conditions are used for choosing proper initial and boundary conditions for accretion of the rotating matter. The study is performed on the basis of the numerical solution of the Euler gas dynamic equations using the second order extension of the Lax-Friedrichs scheme with the implicit approximation of the source term.

Accretion disc theory was first developed as a theory with the local heat balance, where the whole energy produced by a viscous heating was emitted to the sides of the disc. One of the most important new invention of this theory was a phenomenological treatment of the turbulent viscosity, known as ''alpha'' prescription, when the (r$\phi$) component of the stress tensor was approximated by ($\alpha$ P) with a unknown constant $\alpha$. This prescription played the role in the accretion disc theory as well important as the mixing-length theory of convection for stellar evolution. Sources of turbulence in the accretion disc are discussed, including nonlinear hydrodynamical turbulence, convection and magnetic field role. In parallel to the optically thick geometrically thin accretion disc models, a new branch of the optically thin accretion disc models was discovered, with a larger thickness for the same total luminosity. The choice between these solutions should be done of the base of a stability analysis. The ideas underlying the necessity to include advection into the accretion disc theory are presented and first models with advection are reviewed. The present status of the solution for a low-luminous optically thin accretion disc model with advection is discussed and the limits for an advection dominated accretion flows (ADAF) imposed by the presence of magnetic field are analysed.

Approximate approach based on the Galerkin method is suggested for the investigation of equilibrium stellar models, a relativistic collapse problem and black hole formation. Some results of its simplified version - energetic method- are reviewed, and equations for general Gelerkin method are presented.

We have considered a stationary outflowing envelope from the star accelerated by the radiative force in arbitrary optical depth case. Introduced approximations provide satisfactory description of the behavior of the matter flux with partially separated radiation at arbitrary optical depths. The obtained system of differential equations provides a continuous transition of the solution between optically thin and optically thick regions. We analytically derived approximate representation of the solution at the vicinity of the sonic point. Using this representation we numerically integrate the system of equations from the critical point to the infinity. Matching the boundary conditions we obtain the solution of the describing the problem system of differential equations. We have discuss in details the problem of the outer boundary condition. The situation of the photosphere was clarified for the obtained expanding solution. Theoretical approach submitted in this work could be useful for self - consistent simulations of the massive stars evolution with mass loss.

Stability of supermassive stars embedded in hot dark matter is investigated on the base of the energetic method. Stability effect of dark mater is compared with rotational stabilization and preference of the last factor is advocated.

Results of 2D numerical simulation of magnetorotational mechanism of supernova explosion are presented. It is shown that due to the differential rotation of the star toroidal component of magnetic field appears and grows with time. Angular momentum transfers outwards by the toroidal magnetic field. With the evolution of the process part of the envelope of the star is throwing away. The amount of thrown away mass and energy are estimated. The results of the simulation are qualitatively correspond to supernova explosion picture.

Selection effects, connected with stochastic errors in source flux and threshold value determination are analyzed. Normal and normal logarithmic distributions of stochastic deviations are considered. These two kind of distributions produce different effects on the source statistics. Applications to Gamma Ray Burst statistics are discussed. A physical test for checking a close neutron star model of GRBs is suggested.

Timing of Geminga gamma-ray pulsar is done using data of COS B and EGRET. It is shown, that errors in angular coordinates of sources similar to Geminga strongly influence a determination of $\ddot \nu$, so that at angular precision less then $10^{-3}$ arc sec determination of the value of $\ddot \nu$ by means of criterion gives error more then $100 %$. Attempt have been done to improve coordinates of a gamma-ray pulsar using timing analysis. In addition to search of $\nu$, $\dot\nu$ and $\ddot\nu$, a technique is first developed permitting search of two angular coordinates, absolute speed value and direction of a proper motion. In that way timing of gamma-ray pulsars gives amount of information compatible with radiopulsars, but using of real data gives much poorer precision. In gamma-ray sources with rare pulses the periodicity criteria are quite different from the ones in radio region. Data on the coordinates and proper motion of Geminga, obtained from timing studies, do not contradict inside the errors to its identification with $G''$ star and its proper motion. These errors are larger then ones in optical measurements, but are smaller then corresponding errors in X-ray and $\gamma$-ray data. Estimations of the gamma-ray pulsar coordinates and its proper motion could be obtained independently on its optical or radio component, and are available in their absence.

Equations are derived which describe a propagation of strong shocks in the interstellar matter, without any demands for a symmetry, in a thin layer approximation (2.5 dimensions). Using these equations permits to calculate a propagation of shock waves from nonsymmetric supernovae explosions in the medium with arbitrary density distribuion, formation of superbubbles in galaxies.

Formation of a neutron star is accompanied by neutrino emission carring about 10% of the rest energy of the star. Toroidal field produced by twisting of a dipole field in differentially rotating star is antisymmetric. Its summation with antisymmetric toroidal field results in braking of mirror symmetry of the magnetic field . For large magnetic field the neutron decay rate depends on its strength. Neutrino is emitted more in one direction leading to flux asymmetry and recoil of the neutron star. Estimations show that the neutron star can reach velocities $\sim 1000$ km/s for 3% asymmetry of the neutrino flux.

The neutron star formation during the collapse with the strong magnetic field may lead to a mirror symmetry violation and formation of an asymmetric magnetic field. Dependence of the week interaction cross-section on the magnetic field strength lead to the asymmetric neutrino flux and formation of rapidly mooving pulsars due to the recoil action as well as rapidly moving black holes.

Advection-dominated, high-temperature, quasi-spherical accretion flow onto a compact object, recently considered by a number of authors, assume that the dissipation of turbulent energy of the flow heats the ions and that the dissipated energy is advected inward. It is suggested that the efficiency of conversion of accretion energy to radiation can be very much smaller than unity. However, it is likely that the flows have an equipartition magnetic field with the result that dissipation of magnetic energy at a rate comparable to that for the turbulence must occur by Ohmic heating. We argue that this heating occurs as a result of plasma instabilities and that the relevant instabilities are current driven in response to the strong electric fields parallel to the magnetic field. We argue further that these instabilities heat predominantly the electrons. We conclude that the efficiency of conversion of accretion energy to radiation can be much smaller than unity only for the unlikely condition that the Ohmic heating of the electrons is negligible.

We consider the effects of advection and radial gradients of pressure and radial drift velocity on the structure of optically thick accretion disks. We concentrate our efforts on highly viscous disks, $\alpha=1.0$, with large accretion rates. Contrary to disk models neglecting advection, we find that continuous solutions extending from the outer disk regions to the inner edge exist for all accretion rates we have considered. We show that the sonic point moves outward with increasing accretion rate, and that in the innermost disk region advection acts as a heating process that may even dominate over dissipative heating. Despite the importance of advection on it's structure, the disk remains geometrically thin.

We present solutions to the accretion disk structure equations in which advective cooling is accounted for in a self-consistent way. It is shown that for high rates of accretion, above a critical value, when without advection there are no solutions extending continuously from large to small radii, there exist global solutions with self-consistent advective cooling. These solutions are quite different fromthose solutions obtained earlier by other authors, where advective cooling was not treated self-consistently.

An investigation is made of disk accretion of matter onto a rotating star with an aligned dipole magnetic field. A new aspect of this work is that when the angular velocity of the star and disk differ substantially we argue that the $\bf B$ field linking the star and disk rapidly inflates to give regions of open field lines extending from the polar caps of the star and from the disk. The open field line region of the disk leads to the possibility of magnetically driven outflows. An analysis is made of the outflows and their back affect on the disk structure assuming an ``$\ap$" turbulent viscosity model for the disk and a magnetic diffusivity comparable to this viscosity. The outflows are found to extend over a range of radial distances inward to a distance close to $r_{to}$, which is the distance of the maximum of the angular rotation rate of the disk. We find that $r_{to}$ depends on the star's magnetic moment, the accretion rate, and the disk's magnetic diffusivity. The outflow regime is accompanied in general by a spin-up of the rotation rate of the star. When $r_{to}$ exceeds the star's corotation radius $r_{cr} = (GM/\om_*^2)^{1\ov 3}$, we argue that outflow solutions do not occur, but instead that ``magnetic braking" of the star by the disk due to field-line twisting occurs in the vicinity of $r_{cr}$. The magnetic braking solutions can give spin-up or spin-down (or no spin change) of the star depending mainly on the star's magnetic moment and the mass accretion rate. For a system with $r_{to}$ comparable to $r_{cr}$, bimodal behavior is possible where extraneous perturbations cause the system to flip between spin-up and spin-down.