Optical, J, and K photometric observations of the KV UMa black hole X-ray nova in its quiescent state obtained in 2017-2018 are presented. A significant flickering within light curves was not detected, although the average brightness of the system faded by $\approx 0.1^m$ during 350 days. Changes in the average brightness were not accompanied with the increase or the decrease of the flickering. From the modelling of five light curves the inclination of the KV UMa orbit and the black hole mass were obtained: $i=74^{\circ}\pm 4^{\circ}$, $M_{BH}=(7.06\div 7.24)M_{\odot}$ dependently on the used mass ratio. The non-stellar component of the spectrum in the range $\lambda=6400\div 22000$\r{A} can be fitted by a power law $F_{\lambda}\sim \lambda^{\alpha}$, $\alpha\approx -1.8$. The accretion disk orientation angle changed from one epoch to another. The model with spots on the star was inadequate. Evolutionary calculations using the "Scenario Machine" code were performed for low mass X-ray binaries, a recently discovered anomalously rapid decrease of the orbital period was taken into account. We showed that the observed decrease can be consistent with the magnetic stellar wind of the optical companion which magnetic field was increased during the common envelope stage. Several constraints on evolutionary scenario parameters were done.

We investigated the time-variability and spectral properties of the eclipsing X-ray source Circinus Galaxy X-1 (GG X-1), using Chandra, XMM-Newton and ROSAT. We phase-connected the lightcurves observed over 20 years, and obtained a best-fitting period $P = (25,970.0 \pm 0.1)$ s $\approx$7.2 hr, and a period derivative $\dot{P}/P = ( 10.2\pm4.6) \times 10^{-7}$ yr$^{-1}$. The X-ray lightcurve shows asymmetric eclipses, with sharp ingresses and slow, irregular egresses. The eclipse profile and duration vary substantially from cycle to cycle. We show that the X-ray spectra are consistent with a power-law-like component, absorbed by neutral and ionized Compton-thin material, and by a Compton-thick, partial-covering medium, responsible for the irregular dips. The high X-ray/optical flux ratio rules out the possibility that CG X-1 is a foreground Cataclysmic Variable; in agreement with previous studies, we conclude that it is the first example of a compact ultraluminous X-ray source fed by a Wolf-Rayet star or stripped Helium star. Its unocculted luminosity varies between $\approx$4 $\times 10^{39}$ erg s$^{-1}$ and $\approx$3 $\times 10^{40}$ erg s$^{-1}$. Both the donor star and the super-Eddington compact object drive powerful outflows: we suggest that the occulting clouds are produced in the wind-wind collision region and in the bow shock in front of the compact object. Among the rare sample of Wolf-Rayet X-ray binaries, CG X-1 is an exceptional target for studies of super-critical accretion and close binary evolution; it is also a likely progenitor of gravitational wave events.

Recently discovered bursts of gravitational waves provide a good opportunity to verify the current view on the evolution of close binary stars. Modern population synthesis computer programs help to study this evolution from two main sequence stars up to the formation of compact remnants. To calculate the evolution of predecessors of black hole (BH) mergers we used the `Scenario Machine' code. The scenario modelling allowed to describe the evolution of systems for which the final stage is a BH+BH merger and showed the plausibility of modern evolutionary scenarios for binary stars and population synthesis models based on it. We discuss possible values of spins in BH mergers, and give arguments in favour of different values of spins of BH components in BH mergers (low spin + low spin, low spin + high spin, high spin + high spin). We predict that $\geq 50\%$ BH mergers originated through isolated binary evolution should possess spins of both BHs close to zero (low spin + low spin). Up to $\approx$10\% of BH mergers are able to consist of two BHs with spins close to one (high spin + high spin), predecessors of such binaries could be sources of two gamma ray bursts. The BH with the spin close to one could be formed as the result of the merger of two BHs formed via the collapse of a fast rotating star accompanied with a long gamma ray burst and/or a gravitational wave burst with smaller power in comparison with the merger of two BH remnants of binary components.

In 2002, 2004, and 2017 we conducted high precision CCD photometry observations of the eclipsing binary system AS~Cam. By the analysis of the light curves from 1967 to 2017 (our data + data from the literature) we obtained photometric elements of the system and found the change of the system's orbital eccentricity by $\Delta e=0.03 \pm 0.01$. This change can indicate that there is a third companion in the system in a highly inclined orbit with respect to the orbital plane of the central binary, and its gravitational influence may cause the discrepancy between the observed and theoretical apsidal motion rates of AS~Cam.

This study aims timing the eclipses of the short period low mass binary star AB And. The times of minima are taken from the literature and from our observations in October 2013 (22 times of minima) and in August 2014 (23 times of minima). We find and discuss an inaccuracy in the determination of the types of minima in the previous investigation by \citet{li2014}. We study the secular evolution of the central binary's orbital period and the possibility of the existence of third and fourth companions in the system.

The recently discovered burst of gravitational waves GW150914 provides a good new chance to verify the current view on the evolution of close binary stars. Modern population synthesis codes help to study this evolution from two main sequence stars up to the formation of two final remnant degenerate dwarfs, neutron stars or black holes [Massevich 1988]. To study the evolution of the GW150914 predecessor we use the "Scenario Machine" code presented by [Lipunov 1996]. The scenario modelling conducted in this study allowed to describe the evolution of systems for which the final stage is a massive BH+BH merger. We find that the initial mass of the primary component can be $100\div 140 M_{\odot}$ and the initial separation of the components can be $50\div 350 R_{\odot}$. Our calculations show the plausibility of modern evolutionary scenarios for binary stars and the population synthesis modelling based on it.

A short period eclipsing binary star CV Boo is tested for the possible existence of additional bodies in the system with a help of the light equation method. We use data on the moments of minima from the literature as well as from our observations during 2014 May--July. A variation of the CV Boo's orbital period is found with a period of $\approx 75$ d. This variation can be explained by the influence of a third star with a mass of $\approx 0.4M_{\odot}$ in an eccentric orbit with $e\approx 0.9$. A possibility that the orbital period changes on long time scales is discussed. The suggested tertiary companion is near the chaotic zone around the central binary, so CV Boo represents an interesting example to test its dynamical evolution. A list of 14 minima moments of the binary obtained from our observations is presented.

Evolutionary tracks for the X-ray binaries Cyg X-3, IC 10 X-1, NGC 300 X-1, SS 433, and M33 X-7 are computed using the Scenario Machine code. The compact objects in IC 10 X-1, NGC 300 X-1, and M33 X-7 are the most massive stellar mass black hole candidates. Cyg X-3, IC 10 X-1, and NGC 300 X-1 are the only currently known Wolf-Rayet stars with degenerate companions. SS 433 is the only known superaccretor in the Milky Way. Therefore, the stars studied provide excellent laboratories for testing scenarios for the evolution of binaries under extreme conditions. The classical evolutionary scenario is consistent with modern observational data. During the evolution of these binaries, hypernova explosions accompanied by the collapse of stellar cores with large angular momenta can occur, leading to long gamma-ray bursts. At the end of their evolution, Cyg X-3, IC 10 X-1, NGC 300 X-1, and SS 433 may form binary relativistic objects, which will subsequently merge due to the radiation of gravitational waves. The gravitational waves emitted during mergers of relativistic stars should be detectable by existing and future gravitational-wave antennas. In the course of its future evolution, M33 X-7 will pass through a Thorne-Zytkow stage. The formation of a Thorne-Zytkow object can also be accompanied by gravitational-wave radiation. Figure 2 contains two merging black holes with masses almost identical to the second LIGO event GW151226.

Light curves of the eclipsing binary FL Lyr acquired by the Kepler space telescope are analyzed. Eclipse timing measurements for FL Lyr testify to the presence of a third body in the system. Preliminary estimates of its mass and orbital period are > 4M_Jupiter and > 7 yrs. The times of primary minimum in the light curve of FL Lyr during the operation of the Kepler mission are presented.

We analyse a set of moments of minima of eclipsing variable V0873 Per. V0873 Per is a short period low mass binary star. Data about moments of minima of V0873 Per were taken from literature and our observations during 2013-2014. Our aim is to test the system on existence of new bodies using timing of minima of eclipses. We found the periodical variation of orbital period of V0873 Per. This variation can be explained by the gravitational influence of a third companion on the central binary star. The mass of third body candidate is $\approx 0.2 M_{\odot}$, its orbital period is $\approx 300$ days. The paper also includes a table with moments of minima calculated from our observations which can be used in future investigations of V0873 Per.

Based on multiyear INTEGRAL observations of SS433 in 2003-2011, a composite IBIS/ISGRI 18-60 keV orbital light curve is constructed around zero precessional phases psi_{pr}= 0 at the maximim accretion disk opening angle. It shows a peculiar shape with significant excess near the orbital phase phi_orb= 0.25, which is not seen in the softer 2-10 keV energy band. The 40-60 keV orbital light curve demonstrates two almost equal humps at phases \sim 0.25 and \sim 0.75, most likely due to nutation effects of the accretion disk. The nutational variability of SS433 in 15-50 keV with a period of 6.290 days is independently found from timing analysis of Swift/BAT data. The change of the off-eclipse 18-60 keV X-ray flux with the precessional phase shows a double-wave form with strong primary maximum at psi_{pr}= 0 and weak but significant secondary maximum at psi_{pr}= 0.6. A weak variability of the 18-60 keV flux in the middle of the orbital eclipse correlated with the disk precessional phase is also observed. The joint analysis of the broadband 18-60 keV orbital and precessional light curves confirms the presence of a hot extended corona in the central parts of the supercritical accretion disk and constrains the binary mass ratio in SS433 in the range 0.5>q>0.3, suggesting the black hole nature of the compact object.

Based on multiyear INTEGRAL observations of SS433, a composite IBIS/ISGRI 18-60 keV orbital light curve is constructed around zero precessional phase $\psi_{pr}= 0$. It shows a peculiar shape characterized by a significant excess near the orbital phase $\phi_{orb}= 0.25$, which is not seen in the softer 2-10 keV energy band. Such a shape is likely to be due to a complex asymmetric structure of the funnel in a supercritical accretion disk in SS433. The orbital light curve at 40-60 keV demonstrates two almost equal bumps at phases $\sim 0.25$ and $\sim 0.75$, most likely due to nutation effects of the accretion disk. The change of the off-eclipse 18-60 keV X-ray flux with the precessional phase shows a double-wave form with strong primary maximum at $\psi_{pr}= 0$ and weak but significant secondary maximum at $\psi_{pr}= 0.6$. A weak variability of the 18-60 keV flux in the middle of the orbital eclipse correlated with the disk precessional phase is also observed. The joint analysis of the broadband (18-60 keV) orbital and precessional light curves obtained by INTEGRAL confirms the presence of a hot extended corona in the central parts of the supercritical accretion disk and constrain the binary mass ratio in SS433 in the range $0.5\gtrsim q\gtrsim 0.3$, confirming the black hole nature of the compact object. Orbital and precessional light curves in the hardest X-ray band 40-60 keV, which is free from emission from thermal X-ray jets, are also best fitted by the same geometrical model with hot extended corona at $q\sim 0.3$, stressing the conclusions of the modeling of the broad-band X-ray orbital and precessional light curves.

We analyze the influence of the evolution of light absorption by gray dust in the host galaxies of type Ia supernovae (SNe Ia) and the evolution of the mean combined mass of close-binary carbon-oxygen white dwarfs merging due to gravitational waves (SNe Ia precursors) on the interpretation of Hubble diagrams for SNe Ia. A significant increase in the mean SNe Ia energy due to the higher combined masses of merging dwarfs should be observable at redshifts z > 2. The observed relation between the distance moduli and redshifts of SNe Ia can be interpreted not only as evidence for accelerated expansion of the Universe, but also as indicating time variations of the gray-dust absorption of light from these supernovae in various types of host galaxies, observational selection effects, and the decrease in mean combined masses of merging degenerate dwarfs.

We analyze the influence of the evolution of light absorbtion by grey dust in SNe Ia host galaxies and the influence of the evolution of average total mass of coalescing double carbon-oxygen white dwarfs (progenitors of SNe Ia) under the influence of gravitational radiation on the interpretation of Hubble diagrams of SNe Ia. Significant increase in the average energy of SNe Ia due to increase in the total mass of merging dwarfs can be observed at red shift z> 2. The observed dependence of the distance modulus from the red shift in observations of SNe Ia can be explained not only by the assumption about accelerated expansion of the Universe, but also by the evolution of the absorbtion of light by grey dust in various types of host galaxies of SNe Ia, by the effects of observational selection and by the decrease in the average mass of coalescing degenerate dwarfs.

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 have analyzed the observed radial-velocity curve for the X-ray binary M33 X-7 in a Roche model. We have analyzed the dependence between the component masses and the degree of filling of the optical star's Roche lobe to obtain the ratio of the masses of the optical star and compact object. For the most probable mass of the optical star, $m_v=70M_{\odot}$, the mass of the compact object is $m_x=15.55\pm3.20 M_{\odot}$. It has been shown that black holes with masses of $m_x=15M_{\odot}$ and even higher can form in binaries. We present characteristic evolutionary tracks for binary systems passing through an evolutionary stage with properties similar to M33 X-7 - type objects. According to population-synthesis analyses, such binaries should be present in galaxies with masses of at least $10^{11}M_{\odot}$. The present number of such systems in M33 should be of the order of unity. We have also studied the evolutionary status of the X-ray binary IC 10 X-1 with a Wolf-Rayet component, which may contain a massive black hole. The final stages of the evolution of the M33 X-7 and IC 10 X-1 systems should be accompanied by the radiation of gravitational waves.

We consider the possible existence of a common channel of evolution of binary systems, which results in a gamma-ray burst during the formation of a black hole or the birth of a magnetar during the formation of a neutron star. We assume that the rapid rotation of the core of a collapsing star can be explained by tidal synchronization in a very close binary. The calculated rate of formation of rapidly rotating neutron stars is qualitatively consistent with estimates of the formation rate of magnetars. However, our analysis of the binarity of newly-born compact objects with short rotational periods indicates that the fraction of binaries among them substantially exceeds the observational estimates. To bring this fraction into agreement with the statistics for magnetars, the additional velocity acquired by a magnetar during its formation must be primarily perpendicular to the orbital plane before the supernova explosion, and be large.

The "Scenario Machine" (a computer code designed for studies of the evolution of close binaries) was used to carry out a population synthesis for a wide range of merging astrophysical objects: main-sequence stars with main-sequence stars; white dwarfs with white dwarfs, neutron stars, and black holes; neutron stars with neutron stars and black holes; and black holes with black holes.We calculate the rates of such events, and plot the mass distributions for merging white dwarfs and main-sequence stars. It is shown that Type Ia supernovae can be used as standard candles only after approximately one billion years of evolution of galaxies. In the course of this evolution, the average energy of Type Ia supernovae should decrease by roughly 10%; the maximum and minimum energies of Type Ia supernovae may differ by no less than by a factor of 1.5. This circumstance should be taken into account in estimations of parameters of acceleration of the Universe. According to theoretical estimates, the most massive - as a rule, magnetic - white dwarfs probably originate from mergers of white dwarfs of lower mass. At least some magnetic Ap and Bp stars may form in mergers of low-mass main sequence stars (<1.5 mass of the Sun) with convective envelopes.

We present the results of population syntheses for binary stars carried out using the ``Scenario Machine'' code with the aim of analyzing events that may result in long gamma-ray bursts. We show that the observed distribution of morphological types of the host galaxies of long gamma-ray bursts can be explained in a model in which long gamma-ray bursts result from the core collapse of massive Wolf-Rayet stars in close binaries. The dependence of the burst rate on galaxy type is associated with an increase in the rate of stellar-wind mass-loss with increasing stellar metallicity. The separation of binary components at the end of their evolution increases with the stellar-wind rate, resulting in a reduction of the number of binaries that produce gamma-bursts.

Observations of eclipses of the radio pulsar B1259-63 by the disk of its Be-star companion SS 2883 provide an excellent opportunity to study the winds of stars of this type. The eclipses lead to variations in the radio flux (due to variations in the free-free absorption), dispersion measure, rotation measure, and linear polarization of the pulsar. We have carried out numerical modeling of the parameters of the Be-star wind and compared the results with observations.

We present here an updated description of the "Scenario Machine" code. This tool is used to carry out a population synthesis of binary stars. Previous version of the description can be found at http://xray.sai.msu.ru/~mystery//articles/review/contents.html

Using the ``Scenario Machine'' we have carried out a population synthesis of X-ray binaries for the purpose of modelling of X-ray luminosity functions (XLFs) in different types of galaxies: star burst, spiral, and elliptical. This computer code allows to calculate, by using Monte Carlo simulations, the evolution of a large ensemble of binary systems, with proper accounting for the spin evolution of magnetized neutron stars. We show that the XLF has no universal type. It depends on the star formation rate in the galaxy. Also it is of importance to take into account the evolution of binary systems and life times of X-ray stages in theoretical models of such functions. We have calculated cumulative and differential XLFs for the galaxy with the constant star formation rate. Also we have calculated cumulative luminosity functions for different intervals of time after the star formation burst in the galaxy and curves depicting the evolution of the X-ray luminosity after the star formation burst in the galaxy.

We analyze here late evolutionary stages of massive (with initial mass higher than 8 masses of the Sun) close binary stars. Our purposes are to study possible mechanisms of gamma ray bursts (GRBs) origin. We suppose in this paper that GRB phenomenon require formation of massive (approx. 1 M_sun) compact (approx. 10 km) accretion disks around Kerr black holes and neutron stars. Such Kerr black holes are products of collapse of Wolf-Rayet stars in extremely close binaries and merging of neutron stars with black holes and neutron stars with neutron stars in close binary systems. Required accretion disks also can be formed around neutron stars which were formed during collapse of accreting oxygen-neon white dwarfs. We have estimated frequencies of events which lead to a rotational collapse concerned with formation of rapidly rotating relativistic objects in the Galaxy. We made our calculations using the "Scenario Machine".

Using "Scenario Machine" we have carried out population synthesis of radio pulsar with black hole binaries (BH+Psr) in context of the most wide assumptions about star mass loss during evolution, binary stars mass ratio distribution, kick velocity and envelope mass lost during collapse. Our purpose is to display that under any suppositional parameters of evolution scenario BH+Psr population have to be abundant in Galaxy. It is shown that in the all models including models evolved by Heger et al. (2002), Woosley et al. (2002), Heger et al. (2003) expected number of the black holes paired with radio pulsars is sufficient enough to discover such systems within the next few years.