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.

Recent observational data on the type Ia supernova rates are in excellent agreement with the old prediction of the population synthesis of binary stars and confirm that the overwhelming majority of type Ia supernovas (~99%) in elliptical galaxies form via mergers of binary white dwarfs with the total mass exceeding the Chandrasekhar limit.

Hard X-ray INT observations of SS 433 carried out during 2003-2005 years with an analysis of precessional and orbital variability is presented. The width of X-ray eclipse in the 25-50 keV range at the precessional phase $\psi=0.1$ (accretion disk is open to observer) is higher than that in the Ginga 18.4-27.6 keV range. This fact suggests existance the presence of hot extended corona around the supercritical accretion disk. Spectrum of hard X-rays in the range 10-200 keV does not change with the precessional phase which also suggests that hard X-ray flux is generated in the hot extended corona around the accretion disk. The parameters of this hot corona are: kT=23-25 keV, \tau = 1.8-2.8. Mass ratio estimated from the analysis of the ingress part of the eclipse light curve is in the range q=m_x/m_v=0.3-0.5.

Results of simultaneous INTEGRAL and optical observations of galactic microquasar SS433 in May 2003 and INTEGRAL/RXTE observations in March 2004 are presented. Persistent precessional variability with a maximum to minimum uneclipsed hard X-ray flux ratio of \sim 4 is discovered. The 18-60 keV X-ray eclipse is found to be in phase with optical and near infrared eclipses. The orbital eclipse observed by INTEGRAL in May 2003 is at least two times deeper and apparently wider than in soft X-ray band. The broadband X-ray spectrum 2-100 keV simultaneously detected by RXTE/INTEGRAL in March 2004 can be described by bremsstrahlung emission from optically thin thermal plasma with kT\sim 30 keV. The optical spectroscopy with the 6-m SAO BTA telescope confirmed the optical companion to be an A5-A7 supergiant. For the first time, spectorscopic indications of a strong heating effect in the optical star atmosphere are found. The measurements of absorption lines which are presumably formed in the non-illuminated side of the supergiant yield its radial velocity semi-amplitude K_v=132\pm 9 km/s. The analysis of the observed hard X-ray light curve and the eclipse duration, combined with spectroscopically found optical star radial velocity corrected for the strong heating effect, allows us to model SS433 as a massive X-ray binary. Assuming that the hard X-ray source in SS433 is eclipsed by the donor star that exactly fills its Roche lobe, the masses of the optical and compact components in SS433 are suggested to be M_v\approx 30 M_\odot and M_x\approx 9M_\odot, respectively. This provides further evidence that SS433 is a massive binary system with supercritical accretion onto a black hole.

Results of simultaneous {\it INTEGRAL} and optical observations of galactic microquasar SS433 in May 2003 are presented. The analysis of the X-ray and optical eclipse duration and hard X-ray spectra obtained by {\it INTEGRAL} together with optical spectroscopy obtained on the 6-m telescope allows us to construct a model of SS433 as a massive X-ray binary. X-ray eclipse in hard X-rays has a depth of $\sim 80%$ and extended wings. The optical spectroscopy allows us to identify the optical companion as a A5-A7 supergiant and to measure its radial velocity semi-amplitude $K_v=132$ km/s. A strong heating effect in the optical star atmosphere is discovered spectroscopically. The observed broadband X-ray spectrum 2-100 keV can be described by emission from optically thin thermal plasma with $kT\sim 15-20 keV$

Observations of SS433 by INTEGRAL carried out in March -- May 2003 are presented. SS433 is evidently detected on the INTEGRAL images of the corresponding sky region in the energy bands 25-50 and 50-100 keV. The precessional variability of the hard X-ray flux is clearly seen. The X-ray eclipse caused by the binary orbital motion is also detected. A possible origin of the hard continuum is briefly discussed.

We estimate the relativistic binaries merger rate and redshift distribution population synthesis of an ensemble of evolving close binaries. Results of such simulations definitely depend on the cosmic star formation rate history, which is different in galaxies of different types. This leads to a difference in merger rate in spiral and elliptical galaxies, which in principle can be an observational test for GRB models. Also a fit of BATSE long bursts \lns is performed, showing that only a wide (2-3 orders of magnitude) luminosity function provides a good fit.

We argue that a nonthermally looking spectrum of a gamma-ray burst (GRB) can be formed as a superposition of a set of thermal blackbody spectra. This superposition may be done by time integration which is present even in `time resolved' GRB spectroscopy. A nonthermal spectrum can be obtained also by the space integration which should take place unless all the emission comes from a plane front moving exactly towards the observer. We propose a model of the gamma-ray burst spectrum formation based on this idea. This model allows the GRB radiation to be optically thick and to have higher values of baryon load. Thus the latter is limited by the energy considerations only, and not by the condition of a small optical depth.

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.

We examine and develop the model of gamma ray bursts origin proposed by Shaviv and Dar (astro-ph/9606032), according to which the strong gamma ray emission is produced by the interaction of the baryonic relativistic jet (outflow) with the soft photon field in a dense stellar region. It is shown by the simulations of the burst profiles that these profiles are very sensitive to the jet geometry, and the 'clarity' of the GRB profile is (in this model) somehow connected with the gamma rays beaming. Also some event-rate based restrictions are derived and discussed.

Binary neutron stars mergers that are expected to be the most powerful source of energy in the Universe definitely exist in nature, as is proven by the observed behavior of the Hulse-Taylor binary radio pulsar. Though most of energy in such events is radiated in gravitational waves, there probably exist several mechanisms giving also electromagnetic radiation. We propose a new one, involving a revival of the radio pulsar several orbital cycles before the merger.

We present a WWW-version of the {\it Scenario Machine} - a computer code designed to calculate the evolution of close binary stellar systems. The Internet users can directly access to the code and calculate binary evolutionary tracks with parameters at the user's will. The program is running on the {\it Pentium} server of the Division of the Relativistic Astrophysics of the Sternberg Astronimical Institute (http://xray.sai.msu.su/ ). The results are presented both in the form of tables and graphic diagrams. The work is always in progress. More possibilities for Internet users are intended to become available in the near future.

Evolution of the coalescence rate of double neutron stars (NS) and neutron star -- black hole (BH) binaries are computed for model galaxies with different star formation rates. Assuming gamma-ray bursts (GRB) to originate from NS+NS or NS+BH merging in distant galaxies, theoretical logN--logS distributions and <V/V_max> tests of gamma-ray bursts (GRB) are calculated for the first time taking the computed merging rates into account. We use a flat cosmological model (Omega=1) with different values of the cosmological constant Lambda and under various assumptions about the star formation history in galaxies. The calculated source evolution predicts a 5-10 times increase of the source statistics at count rates 3-10 times lower than the exising BATSE sensitivity limit. The most important parameter in fitting the 2nd BATSE catalogue is the initial redshift of star formation, which is found to be z_*=2-5 depending on a poorly determined average spectral index of GRB.