See comments. We calculate analytically an effect of the thermal Comptonization on the cosmic microwave background (CMB) photons when they propagate through galaxy clusters. We estimate a deformation of the CMB spectrum due to multiple Compton scattering of the photons off electrons of an optically-thin hot spherical plasma cloud. This new approach allows us to disentangle and determine the optical depth and the electron temperature of a galaxy cluster using the CMB observations.

We investigate aperiodic X-ray flux variability in accreting highly magnetized neutron stars - X-ray pulsars (XRPs). The X-ray variability is largely determined by mass accretion rate fluctuations at the NS surface, which replicate accretion rate fluctuations at the inner radius of the accretion disc. The variability at the inner radius is due to fluctuations arising all over the disc and propagating inwards under the influence of viscous diffusion. The inner radius varies with mean mass accretion rate and can be estimated from the known magnetic field strength and accretion luminosity of XRPs. Observations of transient XRPs covering several orders of magnitude in luminosity give a unique opportunity to study effects arising due to the changes of the inner disc radius. We investigate the process of viscous diffusion in XRP accretion discs and construct new analytical solutions of the diffusion equation applicable for thin accretion discs truncated both from inside and outside. Our solutions are the most general ones derived in the approximation of Newtonian mechanics. We argue that the break observed at high frequencies in the power density spectra of XRPs corresponds to the minimal time scale of the dynamo process, which is responsible for the initial fluctuations. Comparing data from the bright X-ray transient A 0535+26 with our model, we conclude that the time scale of initial variability in the accretion disc is a few times longer than the local Keplerian time scale.

We present a model for a super-Eddington accretion disc around a magnetized neutron star taking into account advection of heat and the mass loss by the wind. The model is semi-analytical and predicts radial profiles of all basic physical characteristics of the accretion disc. The magnetospheric radius is found as an eigenvalue of the problem. When the inner disc is in radiation-pressure-dominated regime but does not reach its local Eddington limit, advection is mild, and the radius of the magnetosphere depends weakly on the accretion rate. Once approaching the local Eddington limit, the disc becomes advection-dominated, and the scaling for the magnetospheric radius with the mass accretion rate is similar to the classical Alfven relation. Allowing for the mass loss in a wind leads to an increase of the magnetospheric radius. Our model may be applied to a large variety of magnetized neutron stars accreting close to or above their Eddington limits: ultra-luminous X-ray pulsars, Be/X-ray binaries in outbursts, and other systems. In the context of our model we discuss the observational properties of NGC 5907~X-1, the brightest ultra-luminous pulsar known so far, and NGC 300~ULX-1 which is apparently a Be/X-ray binary experiencing a very bright super-Eddington outburst.

Approximate gravitational potentials are often used to describe analytically the motion of particles near black holes (BHs), as well as to study the structure of an accretion disk. Such 'pseudo-Newtonian' potentials are used with the flat-metric equations. Here we consider the motion of a free particle near a non-rotating BH in the context of an exact `logarithmic' gravitational potential. We show how the logarithmic potential gives an exact solution for a mechanical problem and present the relativistic Bernoulli equation for the fluid in the Schwarzschild metric.

The Multimessenger discovery of the merger of two neutron stars on August 17, 2017, GW170817 / GRB170817A, accompanied by a gamma-ray burst and an optical kilonova, is a triumph of the ideas about the evolution of the baryon component in the Universe. Despite the current uniqueness of this observation, the obtained variety of experimental data makes it possible right now to draw important theoretical conclusions about the origin of the double neutron star, their merger, and the subsequent flare-up of the electromagnetic radiation. We present that the discovery of the merger at a distance of 40 Mpc is in full agreement with the very first calculations of the Scenario Machine (Lipunov et al. 1987). In modern terms, the predicted rate is ~ 10 000 Gpc-3.

We investigate the viscous evolution of the accretion disk in 4U 1543-47, a black hole binary system, during the first 30 days after the peak of the 2002 burst by comparing the observed and theoretical accretion rate evolution $\dot M(t)$. The observed $\dot M(t)$ is obtained from spectral modelling of the archival RXTE/PCA data. Different scenarios of disk decay evolution are possible depending on a degree of self-irradiation of the disk by the emission from its centre. If the self-irradiation, which is parametrized by factor $C_\mathrm{irr}$, had been as high as $\sim 5\times10^{-3}$, then the disk would have been completely ionized up to the tidal radius and the short time of the decay would have required the turbulent parameter $\alpha\sim 3$. We find that the shape of the $\dot M(t)$ curve is much better explained in a model with a shrinking high-viscosity zone. If $C_\mathrm{irr}\approx(2-3)\times 10^{-4}$, the resulting $\alpha$ lie in the interval $0.5-1.5$ for the black hole masses in the range $6-10~\mathrm{M}_\odot$, while the radius of the ionized disk is variable and controlled by irradiation. For very weak irradiation, $C_\mathrm{irr} < 1.5 \times10^{-4}$, the burst decline develops as in normal outbursts of dwarf novae with $\alpha \sim 0.08-0.32$. The optical data indicate that $C_\mathrm{irr}$ in 4U 1543-47 (2002) was not greater than approximately $(3-6)\times10^{-4}$. Generally, modelling of an X-ray nova burst allows one to estimate $\alpha$ that depends on the black hole parameters. We present the public 1-D code Freddi to model the viscous evolution of an accretion disk. Analytic approximations are derived to estimate $\alpha$ in X-ray novae using $\dot M(t)$.

We report the early discovery of the optical afterglow of gamma-ray burst (GRB) 140801A in the 137 deg$^2$ 3-$\sigma$ error-box of the Fermi Gamma-ray Burst Monitor (GBM). MASTER is the only observatory that automatically react to all Fermi alerts. GRB 140801A is one of the few GRBs whose optical counterpart was discovered solely from its GBM localization. The optical afterglow of GRB 140801A was found by MASTER Global Robotic Net 53 sec after receiving the alert, making it the fastest optical detection of a GRB from a GBM error-box. Spectroscopy obtained with the 10.4-m Gran Telescopio Canarias and the 6-m BTA of SAO RAS reveals a redshift of $z=1.32$. We performed optical and near-infrared photometry of GRB 140801A using different telescopes with apertures ranging from 0.4-m to 10.4-m. GRB 140801A is a typical burst in many ways. The rest-frame bolometric isotropic energy release and peak energy of the burst is $E_\mathrm{iso} = 5.54_{-0.24}^{+0.26} \times 10^{52}$ erg and $E_\mathrm{p, rest}\simeq280$ keV, respectively, which is consistent with the Amati relation. The absence of a jet break in the optical light curve provides a lower limit on the half-opening angle of the jet $\theta=6.1$ deg. The observed $E_\mathrm{peak}$ is consistent with the limit derived from the Ghirlanda relation. The joint Fermi GBM and Konus-Wind analysis shows that GRB 140801A could belong to the class of intermediate duration. The rapid detection of the optical counterpart of GRB 140801A is especially important regarding the upcoming experiments with large coordinate error-box areas.

We find the Green's functions for the accretion disk with the fixed outer radius and time-independent viscosity. With the Green's functions, a viscous evolution of the disk with any initial conditions can be described. Two types of the inner boundary conditions are considered: the zero stress tensor and the zero accretion rate. The variable mass inflow at the outer radius can also be included. The well-known exponential decline of the accretion rate is a part of the solution with the inner zero stress tensor. The solution with the zero central accretion rate is applicable to the disks around stars with the magnetosphere's boundary exceeding the corotation radius. Using the solution, the viscous evolution of disks in some binary systems can be studied. We apply the solution with zero inner stress tensor to outbursts of short-period X-ray transients during the time around the peak. It is found that for the Kramers' regime of opacity and the initial surface density proportional to the radius, the rise time to the peak is t_rise ~ 0.15 r_out^2/nu_out and the e-folding time of the decay is t_exp ~ 0.45 r_out^2/nu_out. Comparison to non-stationary alpha-disks shows that both models with the same value of viscosity at the outer radius produce similar behaviour on the viscous time-scale. For six bursts in X-ray novae, which exhibit fast-rise-exponential-decay (FRED) and are fitted by the model, we find a way to restrict the turbulent parameter alpha.

We present results of optical polarization observations performed with the MASTER robotic net for three types of objects: gamma-ray bursts, supernovae, and blazars. For the Swift gamma-ray bursts GRB100906A, GRB110422A, GRB121011A, polarization observations were obtained during very early stages of optical emission. For GRB100906A it was the first prompt optical polarization observation in the world. Photometry in polarizers is presented for Type Ia Supernova 2012bh during 20 days, starting on March 27, 2012. We find that the linear polarization of SN 2012bh at the early stage of the envelope expansion was less than 3%. Polarization measurements for the blazars OC 457, 3C 454.3, QSO B1215+303, 87GB 165943.2+395846 at single nights are presented. We infer the degree of the linear polarization and polarization angle. The blazars OC 457 and 3C 454.3 were observed during their periods of activity. The results show that MASTER is able to measure substantially polarized light; at the same time it is not suitable for determining weak polarization (less than 5%) of dim objects (fainter than 16$^m$). Polarimetric observations of the optical emission from gamma-ray bursts and supernovae are necessary to investigate the nature of these transient objects.

We present results of the prompt, early, and afterglow optical observations of five gamma-ray bursts, GRBs 100901A, 100902A, 100905A, 100906A, and 101020A, made with the Mobile Astronomical System of TElescope-Robots in Russia (MASTER-II net), the 1.5-m telescope of Sierra-Nevada Observatory, and the 2.56-m Nordic Optical Telescope. For two sources, GRB 100901A and GRB 100906A, we detected optical counterparts and obtained light curves starting before cessation of gamma-ray emission, at 113 s and 48 s after the trigger, respectively. Observations of GRB 100906A were conducted with two polarizing filters. Observations of the other three bursts gave the upper limits on the optical flux; their properties are briefly discussed. More detailed analysis of GRB 100901A and GRB 100906A supplemented by Swift data provides the following results and indicates different origins of the prompt optical radiation in the two bursts. The light curves patterns and spectral distributions suggest a common production site of the prompt optical and high-energy emission in GRB 100901A. Results of spectral fits for GRB 100901A in the range from the optical to X-rays favor power-law energy distributions with similar values of the optical extinction in the host galaxy. GRB 100906A produced a smoothly peaking optical light curve suggesting that the prompt optical radiation in this GRB originated in a front shock. This is supported by a spectral analysis. We have found that the Amati and Ghirlanda relations are satisfied for GRB 100906A. An upper limit on the value of the optical extinction on the host of GRB 100906A is obtained.

We study statistical properties of long gamma-ray bursts (GRBs) produced by the collapsing cores of WR stars in binary systems. Fast rotation of the cores enables a two-stage collapse scenario, implying the formation of a spinar-like object. A burst produced by such a collapse consists of two pulses, whose energy budget is enough to explain observed GRBs. We calculate models of spinar evolution using results from a population synthesis of binary systems (done by the `Scenario Machine') as initial parameters for the rotating massive cores. Among the resulting bursts, events with the weaker first peak, namely, precursor, are identified, and the precursor-to-main-pulse time separations fully agree with the range of the observed values. The calculated fraction of long GRBs with precursor (about 10 per cent of the total number of long GRBs) and the durations of the main pulses are also consistent with observations. Precursors with lead times greater by up to one order of magnitude than those observed so far are expected to be about twice less numerous. Independently of a GRB model assumed, we predict the existence of precursors that arrive up to >~ 10^3 s in advance of the main events of GRBs.

We examine the jets and the disc of SS 433 at super-Eddington luminosities with 600 times Eddington critical accretion rate by time-dependent two-dimensional radiation hydrodynamical calculations, assuming alpha-model for the viscosity. One-dimensional supercritical accretion disc models with mass loss or advection are used as the initial configurations of the disc. As a result, from the initial advective disc models with alpha =0.001 and 0.1, we obtain the total luminosities 2.5x10^{40} and 2.0x10^{40} erg/s. The total mass-outflow rates are 4x10^{-5} and 10^{-4} solar-mass/yr and the rates of the relativistic axial outflows in a small half opening angle of 1 degree are about 10^{-6} solar-mass/yr: the values generally consistent with the corresponding observed rates of the wind and the jets, respectively. From the initial models with mass loss but without advection, we obtain the total mass-outflow and axial outflow rates smaller than or comparable to the observed rates of the wind and the jets respectively, depending on alpha. In the advective disc model with alpha=0.1, the initially radiation-pressure dominant, optically thick disc evolves to the gas-pressure dominated, optically thin state in the inner region of the disc, and the inner disc is unstable. Consequently, we find remarkable modulations of the disc luminosity and the accretion rate through the inner edge. These modulations manifest themselves as the recurrent hot blobs with high temperatures and low densities at the disc plane, which develop outward and upward and produce the QPOs-like variability of the total luminosity with an amplitude of a factor of 2 and quasi-periods of 10 -- 25 s. This may explain the massive jet ejection and the QPOs phenomena observed in SS 433.

We fit outbursts of two X-ray novae (Nova Monocerotis 1975=A0620-00 and Nova Muscae GS 1991=1124-683) using a time-dependent accretion disk model. The model is based on a new solution for a diffusion-type equation for the non-steady-state accretion and describes the evolution of a viscous alpha-disk in a binary system after the peak of an outburst, when matter in the disk is totally ionized. The accretion rate in the disk decreases according to a power law. We derive formulas for the accretion rate and effective temperature of the disk. The model has three free input parameters: the mass of the central object M, the turbulence parameter alpha, and the normalization parameter delta t. Results of the modeling are compared with the observed X-ray and optical B and V light curves. The resulting estimates for the turbulence parameter $\alpha$ are similar: 0.2-0.4 for A 0620-00 and 0.45-0.65 for GS 1124-683, suggesting a similar nature for the viscosity in the accretion disks around the compact objects in these sources. We also derive the distances to these systems as functions of the masses of their compact objects.

We address the task of modeling soft X-ray and optical light curves of X-ray novae in the high/soft state. The analytic model of viscous evolution of an externally truncated accretion \alpha-disk is used. Relativistic effects near a Kerr black hole and self-irradiation of an accretion disk are taken into account. The model is applied to the outbursts of X-ray nova Monocerotis 1975 (A 0620-00) and X-ray nova Muscae 1991 (GRS 1124-68). Comparison of observational data with the model yields constraints on the angular momentum (the Kerr parameter) of the black holes in A 0620-00 and GRS 1124-68: 0.3-0.6 and \leq 0.4, and on the viscosity parameter \alpha of the disks: 0.7-0.95 and 0.55-0.75. We also conclude that the accretion disks should have an effective geometrical thickness 1.5-2 times greater than the theoretical value of the distance between the photometric layers.

We derive the luminosity-temperature relation for the super-critically accreting black holes (BHs) and compare it to the data on ultraluminous X-ray sources (ULXs). At super-Eddington accretion rates, an outflow forms within the spherization radius. We construct the accretion disc model accounting for the advection and the outflow, and compute characteristic disc temperatures. The bolometric luminosity exceeds the Eddington luminosity L_Edd by a logarithmic factor 1+0.6 ln mdot (where mdot is the accretion rate in Eddington units) and the wind kinetic luminosity is close to L_Edd. The apparent luminosity for the face-on observer is 2-7 times higher because of geometrical beaming. Such an observer has a direct view of the inner hot accretion disc, which has a peak temperature T_max of a few keV in stellar-mass BHs. The emitted spectrum extends as a power-law F_E ~ E**{-1} down to the temperature at the spherization radius T_sp ~ mdot**(-1/2) keV. We associate T_max with a few keV spectral components and T_sp with the soft, 0.1-0.2 keV components observed in ULXs. An edge-on observer sees only the soft emission from the extended envelope, with the photosphere radius exceeding the spherization radius by orders of magnitude. The dependence of the photosphere temperature on luminosity is consistent with that observed in the super-Eddington accreting BHs SS 433 and V4641 Sgr. Strong outflows combined with the large intrinsic X-ray luminosity of the central BH explain naturally the presence of the photoionized nebulae around ULXs. An excellent agreement between the model and the observational data strongly argues in favour of ULXs being super-critically accreting, stellar-mass BHs similar to SS 433, but viewed close to the symmetry axis.

Analytical investigation of time-dependent accretion in disks is carried out. We consider a time-dependent disk in a binary system at outburst which has a fixed tidally-truncated outer radius. The standard Shakura-Sunyaev model of the disk is considered. The vertical structure of the disk is accurately described in two regimes of opacity: Thomson and free-free. Fully analytical solutions are obtained, characterized by power-law variations of accretion rate with time. The solutions supply asymptotic description of disk evolution in flaring sources in the periods after outbursts while the disk is fully ionized. The X-ray flux of multicolor (black-body) alpha-disk is obtained as varying quasi-exponentially. Application to X-ray novae is briefly discussed concerning the observed faster-than-power decays of X-ray light curves. The case of time-dependent advective disk when the exponential variations of accretion rate can occur is discussed.

The stages that follow the merging of two neutron stars are discussed. It is shown that if a rapidly rotating gravitationally bound object is formed after the merging (a spinar or a massive neutron star), then the characteristic time of its evolution is determined by a fundamental value t_{spin} = 3 k (m_p e^2 \hbar^{1/2}) / (m_e^3 c^{5/2} G^{1/2}) \approx 7e+5 s k, where the dimensionless value k depends on the exact equation of state of nuclear matter. The hypothesis is discussed as to whether the residual optical emission of the gamma-ray burst GRB 970228 is pulsar-like and its evolution driven by magnetodipole energy losses. It is shown that binary neutron star mergings can be accompanied by two gravitational wave burst separated either by the time interval t_{spin} or neutron star cooling time (~10 s), depending on the masses of neutron stars.

We consider supercritical accretion onto a compact object when the accretion rate exceeds its value derived from the Eddington limit: L_{Edd}/c^2. We use the scenario that was proposed by Shakura and Sunyaev (1973) for the supercritical regime, according to which matter flows out of the gas-disk surface under radiation pressure. An analytic solution of the accretion-disk structure is obtained for this scenario by using classical equations. An attempt is made to construct a model for nonconservative supercritical accretion with advection. To this end, we develop a numerical scheme for calculating the structure of an advective disk with mass outflow. The goal of solving the problem of a supercritical disk with mass loss is to explain accreting sources with outflowing envelopes, for example, SS 433 and galactic nuclei. The structure and observational manifestations of the envelope which forms around a supercritical disk are considered by taking into account bremsstrahlung absorption and Thomson scattering of photons.

The possible low-frequency radio emission from the progenitors of gamma ray bursts can experience a delay from tens of seconds to hours on the way to the observer due to the dispersion in galactic and extragalactic plasma, and thus reach the observer as a radio afterglow of the burst. This opens a unique possibility (peculiar "Time Machine") of seeing what happened in that place before the catastrophe.

The pattern of variations in the intensity of magnetodipole losses is studied with the relativistic effect of magnetic-field dissipation during collapse into a black hole taken into account. A burst-type solution can be obtained both for a direct collapse and for the formation of a rapidly-rotating, self-gravitating object - a spinar - using a simple model. Analytical dependences on radius describing an electromagnetic burst are derived. The time dependence of the burst shape for an infinitely distant observer and the maximum energy of relativistic particles accelerated by an electric field are numerically calculated. The objects under consideration are of particular interest because particles in their vicinity can be accelerated up to the Planck energies. Possible astrophysical applications to the theory of active galactic nuclei (AGNs) and QSOs are briefly discussed. It is shown for the first time that a spinar can be produced by a merger of neutron stars; this possibility is considered in and without connection with the formation of gamma-ray bursts.