We calculate models of extended stellar atmospheres with a temperature in the range of 12000-40000K and a mass loss rate of $10^{-6}-10^{-4} M_{\odot}$ yr$^{-1}$. A large number of objects with emission spectra, such as luminous blue variables (LBV), FeII - emission line stars, Of/late-WN stars, and even ultraluminous X-ray sources (ULXs) often have effective temperatures in this range. The paper presents the results of model grids calculating in the form of equivalent width diagrams for the selected lines of hydrogen, He, Si, and Fe, as well as the results of studies of some emission objects using the calculated models.

We present preliminary results of a study of the ultraluminous X-ray source UGC6456 ULX. The object is highly variable in both the X-ray (more than 100 times, 0.3-10 keV) and optical (2 times, V band) ranges. The peak X-ray luminosity of UGC6456 ULX exceeds $10^{40}$ erg/s, the absolute magnitude in the bright state exceeds M$_V=-7.6$, which makes it one of the brightest ULXs in the optical range. We found a correlation between the optical and X-ray fluxes, which may indicate that the optical emission is produced by re-processing of the X-rays in outer parts of the optically-thick wind coming from the supercritical accretion disk. Optical spectra of UGC6456 ULX show broad and variable hydrogen and helium emission lines, which also confirms the presence of the strong wind.

We analyzed the X-ray power density spectra of five ultraluminous X-ray sources (ULXs) NGC5408 X-1, NGC6946 X-1, M82 X-1, NGC1313 X-1 and IC342 X-1 that are the only ULXs which display both flat-topped noise (FTN) and quasi-periodic oscillations (QPO). We studied the QPO frequencies, fractional root-mean-square (rms) variability, X-ray luminosity and spectral hardness. We found that the level of FTN is anti-correlated with the QPO frequency. As the frequency of the QPO and brightness of the sources increase, their fractional variability decreases. We propose a simple interpretation using the spherizarion radius, viscosity time and $\alpha$-parameter as basic properties of these systems. The main physical driver of the observed variability is the mass accretion rate which varies >3 between different observations of the same source. As the accretion rate decreases the spherization radius reduces and the FTN plus the QPO move toward higher frequencies resulting in a decrease of the fractional rms variability. We also propose that in all ULXs when the accretion rate is low enough (but still super-Eddington) the QPO and FTN disappear. Assuming that the maximum X-ray luminosity depends only on the black hole (BH) mass and not on the accretion rate (not considering the effects of either the inclination of the super-Eddington disc nor geometrical beaming of radiation) we estimate that all the ULXs have about similar BH masses, with the exception of M82 X-1, which might be 10 times more massive.

We have found three new LBV candidates in the star-forming galaxy NGC 4736. They show typical well-known LBV spectra, broad and strong hydrogen lines, He I lines, many Fe II lines, and forbidden [Fe II] and [Fe III]. Using archival Hubble Space Telescope and ground-based telescope data, we have estimated the bolometric magnitudes of these objects from -8.4 to -11.5, temperatures, and reddening. Source NGC 4736_1 (Mv = -10.2 +/- 0.1 mag) demonstrated variability between 2005 and 2018 as Delta V = 1.1 mag and Delta B = 0.82 mag, the object belongs to LBV stars. NGC 4736_2 (Mv < -8.6 mag) shows P Cyg profiles and its spectrum has changed from 2015 to 2018. The brightness variability of NGC 4736_2 is Delta V = 0.5 mag and Delta B = 0.4 mag. In NGC 4736_3 (Mv = -8.2 +/- 0.2 mag), we found strong nebular lines, broad wings of hydrogen; the brightness variation is only 0.2 mag. Therefore, the last two objects may reside to LBV candidates.

Ultraluminous X-ray sources (ULXs) represent a class of binary systems that are more luminous than any black hole in our Galaxy. The nature of these objects remained unclear for a long time. The most popular models for the ULXs involve either intermediate mass black holes (IMBHs) or stellar-mass black holes accreting at super-Eddington rates. In the last few years our understating of these objects was significantly improved, which made the model of super-Eddington accretion much preferable. Both the X-ray and optical spectra provide evidence for the strong outflows coming from supercritical accretion disk. Another surprising result was discovery of pulsations in four ULXs, which claims that these systems must host neutron stars. Besides the presence of pulsations, there is no sharp difference between ultraluminous pulsars and normal ULXs. This fact implies that significant number of known ULXs might eventually be neutron stars.

We study the X-ray variability of SS433 based on data from the ASCA observatory and the MAXI and RXTE/ASM monitoring missions. Based on the ASCA data, we have constructed the power spectrum of SS433 in the frequency range from $10^{-6}$ to 0.1 Hz, which confirms the presence of a flat portion (flat-topped noise) in the spectrum at frequencies $3\times 10^{-5}$ - $10^{-3}$ Hz. The periodic variability (precession, nutation, eclipses) begins to dominate significantly over the stochastic variability at lower frequencies, which does not allow the stochastic variability to be studied reliably. The best agreement with the observations is reached by the model with the flat portion extending to $9.5\times10^{-6}$ Hz and a power-law spectrum with index of 2.6 below that frequency. The jet nutation with a period of about three days suggests that the time for the passage of material through the disk is less than this value. Therefore, at frequencies below $4\times10^{-6}$ Hz, the power spectrum probably does not reflect the disk structure. It may depend on other factors, for example, a variable mass accretion rate supplied by the donor. The flat portion may arise from a rapid decrease in the viscous time in the supercritical or radiative disk zones. It could be related to the variability of the X-ray jets that are formed in the supercritical region.

We study the optical counterparts of ultraluminous X-ray sources NGC4559 X-10 and NGC4395 ULX-1. Their absolute magnitudes, after taking the reddening into account, are $M_V \approx -5.3$ and $M_V \approx -6.2$, respectively. The spectral energy distribution of the NGC4559 X-10 counterpart is well fitted by a spectrum of an F-type star, whereas NGC4395 ULX-1 has a blue power-law spectrum. Optical spectroscopy of NGC4395 ULX-1 has shown a broad and variable HeII~$\lambda$4686 emission, which puts this object in line with all the other spectrally-studied ULXs. Using the Swift archival X-ray data for NGC4395 ULX-1, we have found a period of $62.8\pm 2.3$ days. The X-ray phase curve of the source is very similar to the precession curve of SS433. The optical variation of the counterpart (between two accurate measurements) amounts to 0.10 mag. Analyzing the absolute magnitudes of 16 well-studied ULX counterparts one may suggest that as the original accretion rate decreases (but nevertheless remains supercritical), the optical luminosity of the wind becomes dimmer and the donor star dominates. However, an observational bias may also influence the distribution.

The origin of Ultraluminous X-ray sources (ULXs) in external galaxies whose X-ray luminosities exceed those of the brightest black holes in our Galaxy by hundreds and thousands of times is mysterious. The most popular models for the ULXs involve either intermediate mass black holes (IMBHs) or stellar-mass black holes accreting at super-Eddington rates. Here we review the ULX properties, their X-ray spectra indicate a presence of hot winds in their accretion disks supposing the supercritical accretion. However, the strongest evidences come from optical spectroscopy. The spectra of the ULX counterparts are very similar to that of SS 433, the only known supercritical accretor in our Galaxy.

We identified two optical counterparts of brightest ultraluminous X-ray Sources (ULXs) in galaxies NGC5474 and NGC3627 (M66). The counterparts in Hubble Space Telescope images are very faint, their V magnitudes are 24.7 ($M_V \approx -4.5$) and 25.9 ($M_V \approx -4.2$), respectively. NGC5474 X-1 changes the X-ray flux more than two orders of magnitude, in its bright state it has $L_X \approx 1.6 \times 10^{40}$ erg/s, the spectrum is best fitted by an absorbed power-law model with a photon index $\Gamma \approx 0.94$. M66 X-1 varies in X-rays with a factor of ~2.5, its maximal luminosity being $2.0 \times 10^{40}$ erg/s with $\Gamma \approx 1.7$. Optical spectroscopy of the NGC5474 X-1 has shown a blue spectrum, which however was contaminated by a nearby star of 23 mag, but the counterpart has a redder spectrum. Among other objects captured by the slit are a background emission-line galaxy (z=0.359) and a new young cluster of NGC5474. We find that these two ULXs have largest X-ray-to-optical ratios of $L_X/L_{opt}$ ~ 7000 for NGC5474 X-1 (in its bright state) and 8000 for M66 X-1 both with the faintest optical counterparts ever measured. Probably their optical emission originates from the donor star. If they have super-Eddington accretion discs with stellar-mass black holes, they may also have the lowest mass accretion rates among ULXs such as in M81 X-6 and NGC1313 X-1.

We study a stochastic variability of SS433 in the $10^{-4} - 5\times 10^{-2}$ Hz frequency range based on RXTE data, and on simultaneous observations with RXTE and optical telescopes. We find that the cross-correlation functions and power spectra depend drastically on the precession phase of the supercritical accretion disc. When the wind funnel of the disc is maximally open to the observer, a flat part emerges in the power spectrum; a break is observed at the frequency $1.7\times10^{-3}$ Hz, with a power-law index $\beta \approx 1.67$ at higher frequencies. The soft emission forming mostly in the jets, lags behind the hard and optical emission. When the observer does not see the funnel and jets (the `edge-on' disc), the power spectrum is described by a single power-law with $\beta \approx 1.34$ and no correlations between X-ray ranges are detected. We investigated two mechanisms to explain the observed variability at the open disc phase, 1) reflection of radiation at the funnel wall (X-rays and optical) and 2) the gas cooling in the jets (X-rays only). The X-ray variability is determined by the contribution of both mechanisms, however the contribution of the jets is much higher. We found that the funnel size is $(2-2.5)\times10^{12}$ cm, and the opening angle is $\vartheta_f\sim 50^\circ$. X-ray jets may consist of three fractions with different densities: $8\times10^{13}$, $3\times10^{13}$ and $5\times10^{11}$ cm$^{-3}$, with most of the jet's mass falling within the latter fraction. We suppose that revealed flat part in the power spectrum may be related to an abrupt change in the disc structure and viscous time-scale at the spherization radius, because the accretion disc becomes thick at this radius, $h/r \sim 1$. The extent of the flat spectrum depends on the variation of viscosity at the spherization radius.

We describe a model of spectral energy distribution in supercritical accretion disks (SCAD) based on the conception by Shakura and Sunyaev. We apply this model to five ultra-luminous X-ray sources (ULXs). In this approach, the disk becomes thick at distances to the center less than the spherization radius, and the temperature dependence is T \propto r^{-1/2}. In this region the disk luminosity is L_bol ~ L_Edd ln(Mdot/Mdot_Edd), and strong wind arises forming a wind funnel above the disk. Outside the spherization radius, the disk is thin and its total luminosity is Eddington, L_Edd. The thin disk heats the wind from below. From the inner side of the funnel the wind is heated by the supercritical disk. In this paper we do not consider Comptonization in the inner hot winds which must cover the deep supercritical disk regions. Our model is technically similar to the DISKIR model of Gierlinski et al. The models differ in disk type (standard - supercritical) and irradiation (disk - wind). We propose to distinguish between these two models in the X-ray region, ~0.3 - 1 keV, where the SCAD model has a flat nu F_nu spectrum, and the DISKIR model never has a flat part, as it is based on the standard alpha-disk. An important difference between the models can be found in their resulting black hole masses. In application to the ULX spectra, the DISKIR model yields black hole masses of a few hundred solar masses, whereas the SCAD model produces stellar-mass black holes ~10 M_sun.