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.

INTEGRAL IBIS/ISGRI 18-60 keV observations of SS433 performed in 2003-2011 enabled the hard X-ray phase-resolved orbital and precessional light curves and spectra to be constructed. The spectra can be fitted by a power-law with photon index $\simeq 3.8$ and remain almost constant while the X-ray flux varies by a factor of a few. This suggests that the hard X-ray emission is produced in an extended quasi-isothermal hot 'corona' surrounding central parts of a supercritical accretion disc. A joint analysis of the broadband 18-60 keV orbital and precessional light curves was performed in the model assuming a significant Roche lobe overfilling by the optical star, up to its filling the outer Lagrangian surface enabling mass loss through the outer Lagrangian L$_2$ point. From this modeling, the relativistic-to-optical component mass ratio $q=M_x/M_v\gtrsim0.4\div 0.8$ is estimated. An analysis of the observed long-term stability of the orbital period of SS433 with an account of the recent observations of SS433 by the VLTI GRAVITY interferometer enabled an independent mass ratio estimate $q>0.6$. This estimate in combination with the radial velocity semi-amplitude for stationary He II emission, $K_x=168\pm 18$ km/s (Hillwig et al 2004) suggests the optical component mass in SS433 $M_v>12 M_\odot$. Thus, the mass of the relativistic component in SS433 is $M_x>7 M_\odot$, which is close to the mean mass of black holes in X-ray binaries ($\sim 8 M_\odot$). The large binary mass ratio in SS433 allows us to understand why there is no common envelope in this binary at the secondary mass transfer evolutionary stage and the system remains semi-detached (van den Heuvel et al. 2017). We also discuss unsolved issues and outline prospects for further study of SS433.

We revisit the determination of binary mass ratio in the Galactic microquasar SS433 based on recent GRAVITY VLTI measurements of mass and angular momentum outflow through a circumbinary disc. The new observations combined with the constancy of the binary orbital period over $\sim 30$ yrs confirm that the mass ratio in SS433 is $q=M_\mathrm{x}/M_\mathrm{v}\gtrsim 0.6$. For the assumed optical star mass $M_\mathrm{v}$ ranging from $\sim 8$ to 15 $M_\odot$ such a mass ratio suggests a low limit of the compact object mass of $M_\mathrm{x}\sim 5-9 M_\odot$, placing the compact object in SS433 as a stellar-mass black hole.

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.

A huge optical luminosity of the supercritical accretion disc and powerful stellar wind in the high-mass X-ray binary SS433 make it difficult to reliably estimate the mass ratio of the binary components from spectroscopic observations. We analyze different indirect methods of the mass ratio estimate. We show that with an account of the possible Roche lobe overflow by the optical star, the analysis of X-ray eclipses in the standard and hard X-ray bands suggests the estimate $q=M_\mathrm{x}/M_\mathrm{v}\gtrsim 0.3$. We argue that the double-peak hydrogen Brackett lines in SS433 should form not in the accretion disc but in a circumbinary envelope, suggesting a total mass of $M_\mathrm{v}+M_\mathrm{x}\gtrsim 40 M_\odot$. The observed long-term stability of the orbital period in SS433 $|\dot P_b/P_b|\le 1.793\times 10^{-14}$~s$^{-1}$ over $\sim 28$ year period is used to place an independent constraint of $q\gtrsim 0.6$ in SS433, confirming its being a Galactic microquasar hosting a superaccreting black hole.

By using available archival X-ray data, we significantly extended the list of times of X-ray minima. The new list includes 65 data points obtained by critically re-analyzing RXTE ASM data, 88 data points based on observations by MAXI, and two data points based on observations by SUZAKU and AstroSat. Analyzing the data along with times of X-ray minima available from the literature, we provide the most accurate estimate of the rate of period change to date. We do not confirm existence of a second derivative of the orbital period suggested by some authors earlier. Instead, we find that the changes in the period can be fit by a sum of quadratic and sinusoidal functions. The period of sinusoidal variations is 15.8 yr. They can be related either to apsidal motion in the close binary with eccentricity e=0.03 or to a presence of a third body with the mass of about 0.7 solar masses located at a distance about 16 au from the close binary. We also detect irregular and abrupt changes in the residuals between the best fit ephemeris and the data. While we discuss possible reasons for the changes, their origin remains unclear. A tentative period of about 188 days in X-ray flux variations was found. Such a period could be attributed to a small precessing disk around the compact object.

Photometric observations of the low-mass X-ray binary system A0620-00=V616 Mon are performed in the optical (unfiltered light, lambda_eff~6400A) and the near-infrared J and K-bands. The mean system flux, the orbital light curve shape and the flickering amplitude dependences on wavelength are examined for two activity stages of the system remaining in quiescence. In 2015-16 A0620-00 was in passive stage (as by Cantrell et al., 2008) exhibiting the regular orbital light curves and low flickering. In less than 230 days in 2016-17 the system switched into active stage: the brightness increased by ~0.2-0.3 mag, the orbital light curve changed while the flickering amplitude increased more than twice. The object regular orbital light curves were fitted by models with "cold" spots on the optical star surface and without those. These models reproduce the observed orbital light curves both in passive and in active stages. The dependence of the mean square flickering amplitude (in fluxes, extinction corrected) on wavelength is computed in the lambda 6400-22000AA range. In active stage, the observed flickering amplitude decreases over the whole studied range and may be represented as Delta F_fl~lambda^{-2} which corresponds to the free-free emission of optically thin high-temperature plasma. In passive stage, flickering obeys Delta F_fl~lambda^{-4} law in the range 6400-12500AA that corresponds to the thermal radiation of optically thick high-temperature plasma. At longer wavelengths the flickering amplitude dependence is flat which may imply existence of a synchrotron component of the relativistic jets emission. These flickering features let us propose that the flickering mechanism includes at least two components: thermal and, apparently, synchrotron, that agrees with the recent discovery of the variable linear polarization of the IR system emission (Russell et al., 2016).

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.

ASTRONIRCAM is a cryogenic-cooled slit spectrograph for the spectral range 1-2.5 mkm installed at the Nasmyth focus of the 2.5-meter telescope of the Caucasian observatory of Sternberg Astronomical Institute of Lomonosov Moscow State University. The instrument is equipped with the HAWAII-2RG 2048x2048 HgCdTe array. Grisms are used as dispersive elements. In the photometric mode ASTRONIRCAM allows for extended astronomical object imaging in the field of view of 4.6x4.6 arc minutes with the 0.269 arcsec/pixel scale in standard photometric bands J, H, K and Ks as well as in narrow-band filters CH_4, [Fe II], H_2 v=1-0 S(1), Br_gamma and CO. In the spectroscopic mode, ASTRONIRCAM takes spectra of extended or point-like sources with spectral resolution R=lambda/Delta lambda <= 1200. The general design, optical system, detector electronics and readout, amplification and digitization scheme are considered. The conversion factor GAIN measurement results are described as well as its dependence on the accumulated signal (non-linearity). The full transmission of the atmosphere-to-detector train ranges from 40 to 50\% in the wide-band photometry mode. The ASTRONIRCAM sensitivity at the 2.5-m telescope is characterized by the limiting J=20, K=19 star magnitudes measured with the 10% precision and 15 minutes integration at the 1 arcsec atmospheric seeing conditions. The references to first results published on the base of ASTRONIRCAM observations are given.

An exact calculation of CaI $\lambda 6439.075$A absorption profiles in the spectra of optical stars in low-mass X-ray binary systems with black holes (BH LMXBs) is carried out. We show that the approximation of a real Roche lobe filling star as disk with uniform local line profile and linear limb darkening law leads to overestimation of projected equatorial rotational velocity $V_{rot} \sin i$ and accordingly, underestimation of mass ratio $q=M_x/M_v$. Refined value of $q$ does not affect the mass of a black hole, but the mass of an optical star has shrunk $\sim 1.5$ times. We present refined components masses in low-mass X-ray binaries with black holes. Companion masses in BH LMXBs are found in the mass range 0.1 - 1.6 $M_{\odot}$ with the peak at $M_v \simeq 0.35 M_{\odot}$. This finding poses additional problem for the standard evolutionary scenarios of BH LMXBs formation with a common envelope phase (CE). We also discuss the implications of these masses into the evolutionary history of the binary.

There is a mismatch between modelled and observed distributions of optical stars masses in BH LMXB. Companion masses in BH LMXB are found in the mass range 0.1 - 1.6 $M_{\odot}$ with the peak at 0.6 $M_{\odot}$. The standard evolutionary scenarios require the donor mass distribution peaks to be at least $\sim 1\, M_{\odot}$ to eject a massive envelope of the black hole progenitor. Imperfect of the methods of optical stars masses determination may cause this difference. Two common used approximations of real Roche lobe filling star as distortion sources have been tested. On the one hand, there is the a approximation of real Roche lobe filling optical stars as sphere. We tested rotational broadening of absorption lines based on an exact calculation of CaI $\lambda 6439.075 \, \AA$ absorption profiles in the spectra of Roche lobe filling optical stars. There is overestimation of projected equatorial rotational velocity $V_{rot} \sin i$ and, accordingly, underestimation of mass ratio $q=M_x/M_v$ in the spherical star model. On the other hand, the approximation of a real Roche lobe filling star as disk with uniform local profile and linear limb darkening law is more rough. In this case overestimation of projected equatorial rotational velocity $V_{rot} \sin i$ is $ \sim 20\%$. Such overestimation of $V_{rot} \sin i$ can result in significant underestimation of mass ratio $q=M_x/M_v$ at hight value of $q=M_x/M_v$. Refined value of $q$ does not affect the mass of a black hole, but the mass of an optical star has shrunk 1.5 times. Hence, the masses of optical stars in BH LMXB must be corrected downward, despite the contradictions to the standard evolutionary scenarios.

This paper describes outstanding issues in astrophysics and cosmology that can be solved by astronomical observations in a broad spectral range from far infrared to millimeter wavelengths. The discussed problems related to the formation of stars and planets, galaxies and the interstellar medium, studies of black holes and the development of the cosmological model can be addressed by the planned space observatory Millimetron (the "Spectr-M" project) equipped with a cooled 10-m mirror. Millimetron can operate both as a single-dish telescope and as a part of a space-ground interferometer with very long baseline.

Resolved UBVRI photometry of RW Aur binary was performed on November 13/14, 2014 during the deep dimming of RW Aur with a newly installed 2.5 meter telescope of the Caucasus observatory of Lomonosov Moscow State University at the mount Shatzhatmaz. At that moment RW Aur A was $\simeq 3^m$ fainter than in November 1994 in all spectral bands. We explain the current RW Aur A dimming as a result of eclipse of the star by dust particles with size $>1 \mu m.$ We found that RW Aur B is also a variable star: it was brighter than 20 years ago at $0.7^m$ in each of UBVRI band (gray brightening).

We consider the relationship between the masses of the compact nuclear objects in the centers of disky galaxies -- supermassive black holes (SMBHs) or nuclear star clusters (NCs) -- and such parameters as the maximal velocity of rotation $V_{\textrm{max}}$, obtained from the rotation curves, indicative dynamical mass $M_{25}$, and the color index ($B{-}V$) of their parent galaxies. It was found that the mass of nuclear clusters $M_{\rm nc}$ correlates more closely with the velocity of rotation and total mass of galaxies than the mass of supermassive black holes $M_{\rm bh}$. The dependence of masses of the central objects on the color index is bimodal: galaxies of the red group (red-sequence), which have ($B{-}V) > 0.6{-}0.7$, differ from bluer galaxies, by higher values of $M_{\rm bh}$ for similar host-galaxy parameters. In contrast, in the diagrams for nuclear clusters the "blue" and "red" galaxies form unified sequences. It agrees with scenarios in which most red-group galaxies form as a result of loss of interstellar gas in a stage of high nuclear activity in galaxies whose central black-hole masses are high, exceeding $10^6 {-} 10^7 M_{\odot}$ (depending on the total mass of the galaxy). The active growth of nuclear star clusters possibly begun after the violent AGN activity.

We consider the relationship between the masses of the central objects in disky galaxies - supermassive black holes (SMBHs) and nuclear star clusters (NCs) - and various parameters of parent galaxies: velocity of rotation $V_{(2)}$ at $R= 2$ kpc, maximal velocity of rotation $V_{\textrm{max}}$, the indicative dynamical mass $M_{25}$, the total mass of the stellar population $M_{*}$, and the total color index ($B{-}V$). The mass of nuclear clusters $M_{\rm nc}$ correlates more closely with the kinematic parameters and total mass of galaxies than the mass of central black holes $M_{\rm bh}$, including correlation with the velocity $V_{\textrm{max}}$, which characterizes the virial mass of the dark halo. On average, lenticular galaxies are distinguished by higher masses $M_{\rm bh}$ compared to other types of galaxies with similar characteristics. The dependence of masses of central objects on the color index is bimodal: galaxies of the red group (red-sequence), which have ($B{-}V) > 0.6{-}0.7$, being mostly early-type galaxies, differ from bluer galaxies, by higher values of $M_{\rm nc}$ and $M_{\rm bh}$. The red-group galaxies have systematically higher $M_{\rm bh}$ values, even when the host-galaxy parameters are similar. In contrast, in the case of nuclear stellar clusters, the "blue" and "red" galaxies form unified sequences at the diagrams. The results agree with scenarios in which most red-group galaxies form as a result of the partial or complete loss of interstellar gas in a stage of high nuclear activity in galaxies whose central black-hole mass exceeds $10^6 {-} 10^7 M_{\odot}$ (depending on the total mass of a galaxy). Most of galaxies with $M_{\rm bh}> 10^7 M_{\odot}$ are lenticular galaxies (types S0, E/S0), whose disks are usually devoid of gas.

The radial velocity curves of optical components in X-ray binary systems can differ from the radial velocity curves of their barycenters due to tidal distortion, gravitational darkening, X-ray heating, etc. This motivated us to investigate how the semiamplitudes of the radial velocity curves of these optical components can depend on the binary-system parameters in a Roche model. The K-correction is taken to be the ratio of the radial velocity semiamplitude for a star in the Roche model to the corresponding value for the stellar barycenter. K-corrections are tabulated for the optical stars in the massive X-ray binaries Cen X-3, LMC X-4, SMC X-1, Vela X-1, and 4U 1538-52.

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.

We suggest a simple synthesis model of an eclipsing binary system which includes one component with strong stellar wind. Numerical simulations show that the shape of the light curve (and in particularly the widths of the minima) strongly depends on wind parameters. Wind effects are crucial in modelling light curves of binaries including e.g., WR stars.

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 have obtained the first estimates of the masses of the components of the Her X-1/HZ Her X-ray binary system taking into account non-LTE effects in the formation of the H_gamma absorption line: mx=1.8Msun and mv=2.5Msun. These mass estimates were made in a Roche model based on the observed radial-velocity curve of the optical star, HZ Her. The masses for the X-ray pulsar and optical star obtained for an LTE model lie are mx=0.85\pm0.15Msun and mv=1.87\pm0.13Msun. These mass estimates for the components of Her X-1/HZ Her derived from the radial-velocity curve should be considered tentative. Further mass estimates from high-precision observations of the orbital variability of the absorption profiles in a non-LTE model for the atmosphere of the optical component should be made.

We interpret the observed radial-velocity curve of the optical star in the low-mass X-ray binary 2S 0921-630 using a Roche model, taking into account the X-ray heating of the optical star and screening of X-rays coming from the relativistic object by the accretion disk. Consequences of possible anisotropy of the X-ray radiation are considered.We obtain relations between the masses of the optical and compact (X-ray) components, mv and mx, for orbital inclinations i=60, 75, 90 degrees. Including X-ray heating enabled us to reduce the compact object's mass by near 0.5-1Msun, compared to the case with no heating. Based on the K0III spectral type of the optical component (with a probable mass of mv=2.9Msun, we concluded that mx=2.45-2.55Msun (for i=75-90 degrees). If the K0III star has lost a substantial part of its mass as a result of mass exchange, as in the V404 Cyg and GRS 1905+105 systems, and its mass is $m_v=0.65-0.75Msun, the compact object's mass is close to the standard mass of a neutron star, mx=1.4Msun (for i=75-90 degrees). Thus, it is probable that the X-ray source in the 2S 0921-630 binary is an accreting neutron star.

High-accuracy multicolor light curves of the binary system HD 189733, which contains an exoplanet, are analyzed. We have determined the radii of the star and the planet in the binary system, as well as the orbital inclination. The limb-darkening coefficients of the stellar disk were obtained in 10 filters in the wavelength interval 5500-10500 AA. The uncertainties of the fitted parameters were estimated using the differential-correction and confidence-area methods. The wavelength dependence of the limb-darkening coefficients is compared to the corresponding theoretical function for a model thin stellar atmosphere. We confirm the wavelength dependence of the exoplanet's radius found by Pont et al. The exoplanet radius increases with decreasing wavelength, which seems to argue for the presence of an atmosphere around the planet.

The statistical relations between the masses of supermassive black holes (SMBHs) in disk galaxies and the kinematic properties of their host galaxies are analyzed. We use the radial velocity profiles for several galaxies obtained earlier at the 6-m telescope of the Special Astrophysical Observatory of the Russian Academy of Sciences parallel with the data for other galaxies taken from the literature. We demonstrate that the SMBH masses correlate well with the velocities of rotation of disks at a fixed distance R \approx 1 kpc (V1), which characterize the mean density of the central region of the galaxy. The SMBH masses correlate appreciably weaker with the asymptotic velocity at large distances from the center and with the angular velocity at the optical radius R_{25}. We suggest that the growth of the SMBH occurs inside of the forming "classical" bulge during a monolithic collapse of gas in the central kpc-size region of the protogalaxy. We have also found a correlation between the SMBH mass and the total (indicative) mass of the galaxy M_{25} within the optical radius R_{25}, which includes both baryonic and "dark" mass. The masses of the nuclear star clusters in early-type disk galaxies (based on the catalog of Seth et al.) are also scaled with the dynamical mass M_{25}, whereas the correlations with the luminosity and velocity of rotation of galaxies are practically absent for them. For a given M_{25} the masses of the nuclear clusters are, on average, nearly order of magnitude higher in S0-Sbc galaxies than in late-type galaxies.

We report the discovery of a new Galactic candidate Luminous Blue Variable (cLBV) via detection of an infrared circular nebula and follow-up spectroscopy of its central star. The nebula, MN112, is one of many dozens of circular nebulae detected at $24 \mu$m in the {\it Spitzer Space Telescope} archival data, whose morphology is similar to that of nebulae associated with known (c)LBVs and related evolved massive stars. Specifically, the core-halo morphology of MN112 bears a striking resemblance to the circumstellar nebula associated with the Galactic cLBV GAL 079.29+00.46, which suggests that both nebulae might have a similar origin and that the central star of MN112 is a LBV. The spectroscopy of the central star showed that its spectrum is almost identical to that of the bona fide LBV P Cygni, which also supports the LBV classification of the object. To further constrain the nature of MN112, we searched for signatures of possible high-amplitude ($\ga 1$ mag) photometric variability of the central star using archival and newly obtained photometric data covering a 45 year period. We found that the B magnitude of the star was constant ($\simeq$ 17.1$\pm$0.3 mag) over this period, while in the I band the star brightened by $\simeq 0.4$ mag during the last 17 years. Although the non-detection of large photometric variability leads us to use the prefix `candidate' in the classification of MN112, we remind that the long-term photometric stability is not unusual for genuine LBVs and that the brightness of P Cygni remains relatively stable during the last three centuries.

We report the serendipitous discovery of a ring nebula around a candidate Wolf-Rayet (WR) star, HBHA 4202-22, in Cygnus using the Spitzer Space Telescope archival data. Our spectroscopic follow-up observations confirmed the WR nature of this star (we named it WR138a) and showed that it belongs to the WN8-9h subtype. We thereby add a new example to the known sample of late WN stars with circumstellar nebulae. We analyzed the spectrum of WR138a by using the Potsdam Wolf-Rayet (PoWR) model atmospheres, obtaining a stellar temperature of 40 kK. The stellar wind composition is dominated by helium with 20 per cent of hydrogen. The stellar spectrum is highly reddened and absorbed (E_{B-V} = 2.4 mag, A_V = 7.4 mag). Adopting a stellar luminosity of log L/Lsun = 5.3, the star has a mass-loss rate of 10^{-4.7} Msun yr^{-1}, and resides in a distance of 4.2 kpc. We measured the proper motion for WR138a and found that it is a runaway star with a peculiar velocity of ~50 km/s. Implications of the runaway nature of WR138a for constraining the mass of its progenitor star and understanding the origin of its ring nebula are discussed.

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.

The analysis of hard X-ray INTEGRAL observations (2003-2008) of superaccreting galactic microquasar SS433 at precessional phases of the source with the maximum disk opening angle is carried out. It is found that the shape and width of the primary X-ray eclipse is strongly variable suggesting additional absorption in dense stellar wind and gas outflows from the optical A7I-component and the wind-wind collision region. The independence of the observed hard X-ray spectrum on the accretion disk precessional phase suggests that hard X-ray emission (20-100 keV) is formed in an extended, hot, quasi-isothermal corona, probably heated by interaction of relativistic jet with inhomogeneous wind outflow from the precessing supercritical accretion disk. A joint modeling of X-ray eclipsing and precessional hard X-ray variability of SS433 revealed by INTEGRAL by a geometrical model suggests the binary mass ratio $q=m_x/m_v\simeq 0.25\div 0.5$. The absolute minimum of joint orbital and precessional $\chi^2$ residuals is reached at $q\simeq 0.3$. The found binary mass ratio range allows us to explain the substantial precessional variability of the minimum brightness at the middle of the primary optical eclipse. For the mass function of the optical star $f_v=0.268 M_\odot$ as derived from Hillwig & Gies data, the obtained value of $q\simeq 0.3$ yields the masses of the components $m_x\simeq 5.3 M_\odot$, $m_v\simeq 17.7 M_\odot$, confirming the black hole nature of the compact object in SS433.

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.

The analysis of INTEGRAL observations (2003-2008) of superaccreting galactic microquasar SS433 at precessional phases with the maximum disk opening angle is carried out. The shape and width of the primary X-ray eclipse is found to be strongly variable suggesting additional absorption in dense stellar wind and gas outflows from the optical A7I-component. The joint modeling of X-ray eclipse and precessional X-ray variability by a geometrical model suggests the binary mass ratio q=m_x/m_v=0.3, allowing an explnation of peculiarities of the optical variability of SS433, in particular, the substantial precessional variability at the primary optical eclipse minimum. For the mass function of the optical star f_v=0.268 M_\odot as derived from Hillwig and Gies (2008) data, the obtained q yields the masses of the components m_x=5 M_\odot, m_v=15 M_\odot, confirming the black hole nature of the relativistic object in SS433. The independence of the observed hard X-ray spectrum on the precession phase suggests that hard X-ray emission is formed in an extended hot corona. The Monte-Carlo simulations of the broadband X-ray spectrum of SS433 at the maximum disk opening precessional phases allowed us to determine physical parameters of the corona (temperature T_{cor}=20 keV, Thomson optical depth \tau=0.2), and to estimate the jet mass outflow rate \dot M_j=3\times 10^{19} g/s yielding the kinetic power of the jets \sim 10^{39} erg/s.

We consider small-scale spheroidal clusters of weakly interacting massive particles in our Galaxy as non-compact gravitational microlenses and predict the appearance of caustics in the plane of a lensed source. The crossing of these caustics by a lensed star can produce a large variety of light curves, including some observed in actual microlensing events that have been interpreted as manifestations of binary gravitational lenses. We consider also observable effects during the gravitational microlensing of stars of non-zero angular size with a given brightness distribution across their disks by such an exotic objects as natural wormholes and objects whose space-time environment is described with the NUT metric. We demonstrate that, under certain conditions, the microlensing light curves, chromatic and polarizational effects due to the properties of the lens and the star disk brightness distributions can differ considerably from those observed for a Schwarzschild gravitational lens, so that their analysis can facilitate the identification of such objects.

BAT99-129 is a massive eclipsing binary system in the Large Magellanic Cloud (LMC) which consists of WN3(h) and O5V components. A broad-band MACHO light curve of the system is studied in the present paper. A dense and extended atmosphere of the Wolf-Rayet (WR) star does not allow one to analyze the light curve in terms of standard parametric models of Wilson-Devinney type. Distributions of brightness and absorption across the WR star disk are restored using direct solution of integral equations describing eclipses in the system. As a result, reliable estimates of the orbital inclination and component parameters have been obtained. The orbital inclination is 78 deg, the orbital separation is 28.5 solar radii, the radius of the O component is 7.1 solar radii. The size of the WR component core opaque in optical continuum is 3.4 solar radii. The brightness temperature in the center of the WR disk is about 45000K. Probable errors of the parameters are discussed. Evolutionary status of the system is discussed.

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.

We carried out numerical experiments on the evaluation of the possibilities of obtaining the information about brightness distributions for the components of eclipsing variables from the data of high-precision photometry expected for planned satellites COROT and Kepler. We examined a simple model of the eclipsing binary with the spherical components on circular orbits and the linear law of the limb darkening. The solutions of light curves have been obtained as by fitting of the nonlinear model, into the number of parameters of which included the limb darkening coefficients, so also by the solution of the ill-posed inverse problem of restoration of brightness distributions across the disks of stars without rigid model constraints on the form of these functions. The obtained estimations show that if the observational accuracy amounts to 0.0001 then the limb darkening coefficients can be found with the relative error approximately 0.01 . The brightness distributions across the disks of components can be restored also nearly with the same accuracy.

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.

We present a technique for the successive restoration of the branches of the one-dimensional strip brightness distribution across a quasar's accretion disk via the analysis of observations of high magnification events in measured fluxes from the multiple quasar images produced by a gravitational lens. Hypothesizing these events to by associated with microlensing by a fold caustic, the branches of brightness distribution are searched for on compact sets of non-negative, monotonically non-increasing, convex downward functions. The results of numerical simulations show that the solution obtained is stable against random noise. Analysis of the light curves of high magnification events in the fluxes from components C and A of the gravitational lens QSO 2237+0305, observed by the OGLE and GLITP groups, has yielded the forms of the strip brightness distributions across the accretion disk of the lensed quasar. The resulting sizes of the accretion disk are in agreement with results obtained earlier via model-fitting. The form of the brightness distribution is consistent with the expected appearance of an accretion disk rotating around supermassive black hole.

The relations between masses of Super-Massive Black Holes, Mbh, in galactic nuclei, maximal rotational velocities, Vm, and indicative masses, Mi, of galaxies are studied for galaxies with the available rotation curves. Mbh correlates with both Vm and Mi, although much weaker than with the central velocity dispersion S0. Masses Mbh for early-type galaxies (S0-Sab), which usually possess more luminous bulges, are, on average, larger than those for late-type galaxies with similar velocities of rotation. It enables to conclude that masses of black holes are determined mainly by the bulge properties, and - much weaker - by dark haloes of galaxies.

We present optical and NIR (near infrared) follow up observations of the GRB 001109 from 1 to 300 days after the burst. No transient emission was found at these wavelengths within this GRB's (Gamma Ray Burst) 50" radius BeppoSAX error box. Strong limits (3 sigma) are set with: R >~ 21, 10.2 hr after the GRB; I >~ 23, 11.4 hr after the GRB; H >~ 20.7, 9.9 hr after the GRB; and Ks >~ 20, 9.6 hours after the GRB. We discuss whether the radio source found in the GRB's error box (Taylor et al. 2000) might be related to the afterglow. We also present a multiwavelength study of a reddened starburst galaxy, found coincident with the potential radio and the X ray afterglow. We show that our strong I band upper limit makes of the GRB 001109 the darkest one localised by the BeppoSAX's NFI (Narrow Field Instrument), and it is one of the most constraining upper limits on GRB afterglows to date. Further to it, the implications of these observations in the context of dark GRBs are considered.

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$

Analysis of 18-120 keV images of the Sagittarius Arm Tangent region (SATR) obtained by IBIS telescope onboard INTEGTRAL observatory during the spring of 2003 is performed. In the 18-60 keV energy range, 28 sources have been detected with a flux level above 1.4 mCrab. Of these sources, 16 were identified earlier as galactic X-ray binary systems, 3 as extragalactic objects, 2 as pulsars inside supernova remnants, and 7 has unknown nature. The analysis revealed the presence of three previously unknown sources. Fourteen sources show significant flux in the 60-120 keV energy range

During the period Aug.23-Sept.24 2003, the INTEGRAL observatory performed an ultra deep survey of the Galactic Center region with a record sensitivity at energies higher than 20 keV. We have analized images of the Galactic Center region obtained with the ISGRI detector of the IBIS telescope (15-200 keV) and present here a catalog of detected sources. In total, 60 sources with a flux higher than 1.5 mCrab have been detected. 44 of them were earlier identified as Galactic binary systems, 3 are extragalactic objects. 2 new sources are discovered.

We analyze the mass distribution of stellar black holes derived from the light and radial velocity curves of optical stars in close binary systems using dynamical methods. The systematic errors inherent in this approach are discussed. These are associated primarily with uncertainties in models for the contribution from gaseous structures to the optical brightness of the systems under consideration. The mass distribution is nearly flat in the range 4-15M_sun. This is compared with the mass distribution for black holes in massive close binaries, which can be manifest as ultraluminous X-ray sources (L_x > 10^39 erg/s) observed in other galaxies. If the X-ray luminosities of these objects correspond to the Eddington limit, the black-hole mass distribution should be described by a power law, which is incompatible with the flat shape derived dynamically from observations of close binaries in our Galaxy. One possible explanation of this discrepancy is the rapid evaporation of stellar-mass black holes predicted in recent multi-dimensional models of gravity. This hypothesis can be verifed by measuring the stellar black-hole mass spectrum or finding isolated or binary black holes with masses below 3M_sun.

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 present synthetic Doppler maps of gaseous flows in binary IP Peg based on the results of 3D gasdynamical simulations. Using of gasdynamical calculations alongside with the Doppler tomography technique permits us to identify the main features of the flow on the Doppler maps without solution of an ill-posed inverse problem. Comparison of synthetic tomograms with observations shows that in quiescence there are two zones of high emission: a shock wave on the edge of the stream from L_1 caused by the interaction of the gas of circumbinary envelope and the stream, and the dense region in a apoastron of quasi-elliptical accretion disk. A single arm of the spiral shock and the stream itself give a minor input to the total brightness. During outburst the accretion disk dominates, and the most emitive regions are the two arms of the spiral shock.

Observed properties of GRBs, WR stars and their CO-cores in the end of evolution are analyzed. A possible bimodality of the observed GRB energy distribution ($10^{48}$ erg for GRB9809425; $3\times 10^{51}\div 2\times 10^{54}$ erg for others) is in accord with the bimodal mass distribution of relativistic objects ($M_{NS}=(1.35\pm 0.15) M_\odot$; $M_{BH}=(4\div 15) M_\odot$). The peculiarity of SN1998bw can be related to the rotation of the collapsing CO-core. The expected galactic collapse rate of CO-cores of most compact WR stars of type WO is $\sim 10^{-5}$ per year, only by one and a half order of magnitude higher than the GRB rate. The allowance for a gamma-ray beaming or random outcome of the CO-core collapse due to some instabilities brings this rate in accordance GRB rate. We argue that WR stars (most probably, of type WO) can be considered as progenitors of cosmic gamma-ray bursts. Two types of GRBs are predicted in correspondence with the bimodal mass distribution of the relativistic objects. Three types of optical afterglows should appear depending on which CO-core is collapsing: of a single WR star, of a WR star in a WR+O or a hypothetic WR+(A-M) binary system. We briefly discuss a model of GRB as a transient phenomenon occurring at early stages of galactic evolution ($z>1$), when very massive ($M>100 M_\odot$) low-metallicity stars could form. WR-galaxies can be most probable candidates for GRB host galaxies.

The results of three-dimensional numerical simulations of the flow of matter in non-magnetic semidetached binary systems are presented. Self-consistent solutions indicate the absence of a shock interaction between the stream of matter flowing from the inner Lagrange point and the accretion disk (a "hot spot"). At the same time, the interaction between the stream and the common envelope of the system forms an extended shock wave along the edge of the stream, whose observational properties are roughly equivalent to those of a hot spot in the disk. Comparison of synethesized and observed light curves confirm the adequacy of the proposed model for the flow of matter without formation of a hot spot. The model makes it possible for the first time to reproduce the variety of the observed light curves of cataclysmic variables in the framework of a single model, without invoking additional assumptions.