Magnetic fields of strongly magnetized stars can trap conducting matter due to frozen-in condition. In the force-free regime, the motion of the matter along the field lines may be considered in the "bead on a wire" approximation. Such a motion, if gravity and centrifugal forces are taken into account, has equilibrium points, some of which are stable. In most cases, stability is possible in about several per cent of the possible locations. Corresponding oscillation frequencies span the range from zero to $\sqrt{3}$ of the spin frequency. We suggest that this variability mode may be excited in some X-ray pulsars during the outbursts and create the peaked broad-band noise component near the break frequency in the power density spectrum, as well as produce some of the quasi-periodic oscillation features in this frequency range. Existence of this variability does not require any changes in mass accretion rate and involves only a small amount of matter infiltrating from the disc and magnetic flow due to interchange instabilities.

When the accretion disc around a weakly magnetised neutron star (NS) meets the stellar surface, it should brake down to match the rotation of the NS, forming a boundary layer. As the mechanisms potentially responsible for this braking are apparently inefficient, it is reasonable to consider this layer as a spreading layer (SL) with negligible radial extent and structure. We perform hydrodynamical 2D spectral simulations of an SL, considering the disc as a source of matter and angular momentum. Interaction of new, rapidly rotating matter with the pre-existing, relatively slow material co-rotating with the star leads to instabilities capable of transferring angular momentum and creating variability on dynamical timescales. For small accretion rates, we find that the SL is unstable for heating instability that disrupts the initial latitudinal symmetry and produces large deviations between the two hemispheres. This instability also results in breaking of the axial symmetry as coherent flow structures are formed and escape from the SL intermittently. At enhanced accretion rates, the SL is prone to shearing instability and acts as a source of oblique waves that propagate towards the poles, leading to patterns that again break the axial symmetry. We compute artificial light curves of an SL viewed at different inclination angles. Most of the simulated light curves show oscillations at frequencies close to 1kHz. We interpret these oscillations as inertial modes excited by shear instabilities near the boundary of the SL. Their frequencies, dependence on flux, and amplitude variations can explain the high-frequency pair quasi-periodic oscillations observed in many low-mass X-ray binaries.

We report on the analysis of the spin evolution of a slow X-ray pulsar GX 301$-$2 along the orbit using long-term monitoring by Fermi/GBM. Based on the observationally confirmed accretion scenario and an analytical model for the accretion of angular momentum we demonstrate that in this system, the neutron star spins retrogradely, that is, in a direction opposite to the orbital motion. This first-of-a-kind discovery of such a system proves the principal possibility of retrograde rotation in accreting systems with suitable accretion torque, and might have profound consequences for our understanding of the spin evolution of X-ray pulsars, estimates of their initial spin periods, and the ultimate result of their evolution.

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

In standard accretion discs, outward angular momentum transfer by viscous forces is compensated by the inward motion of the accreting matter. However, the vertical structure of real accretion discs leads to meridional circulation with comparable amplitudes of poloidal velocities. Using thin-disc approximation, we consider different regimes of disc accretion with different vertical viscosity scalings. We show that, while gas-pressure-dominated discs can easily have a midplane outflow, standard thin radiation-pressure-dominated disc is normally moving inwards at all the heights. However, quasi-spherical scaling for pressure ($p\propto \varpi^{-5/2}$) leads to a midplane outflow for a very broad range of parameters. It particular, this may lead to a reversed, outward heat advection in geometrically thick discs when the temperature decreases rapidly enough with height. While the overall direction of heat advection depends on the unknown details of vertical structure and viscosity mechanisms, existence of the midplane counterflow in quasi-spherical flows is a robust result weakly dependent on the parameters and the assumptions of the model. Future models of thick radiatively inefficient flows should take meridional circulation into account.

The aim of this paper is to investigate the transition of a strongly magnetized neutron star into the accretion regime with very low accretion rate. For this purpose we monitored the Be-transient X-ray pulsar GRO J1008-57 throughout a full orbital cycle. The current observational campaign was performed with the Swift/XRT telescope in the soft X-ray band (0.5-10 keV) between two subsequent Type I outbursts in January and September 2016. The expected transition to the propeller regime was not observed. However, the transitions between different regimes of accretion were detected. In particular, after an outburst the source entered a stable accretion state characterised by the accretion rate of ~10^14-10^15 g/s. We associate this state with accretion from a cold (low-ionised) disc of temperature below ~6500 K. We argue that a transition to such accretion regime should be observed in all X-ray pulsars with certain combination of the rotation frequency and magnetic field strength. The proposed model of accretion from a cold disc is able to explain several puzzling observational properties of X-ray pulsars.

Most of ultraluminous X-ray sources are thought to be objects accreting above their Eddington limits. In the recently identified class of ultraluminous X-ray pulsars, accretor is a neutron star and thus has a fairly small mass with a small Eddington limit. The accretion disc structure around such an object affects important observables such as equilibrium period, period derivative and the size of the magnetosphere. We propose a model of a nearly-standard accretion disc interacting with the magnetosphere only in a thin layer near the inner disc rim. Our calculations show that the size of the magnetosphere may be represented as the classical Alfv\'en radius times a dimensionless factor $\xi$ which depends on the disc thickness only. In the case of radiation-pressure-dominated disc, the size of the magnetosphere does not depend on the mass accretion rate. In general, increasing the disc thickness leads to a larger magnetosphere size in units of the Alfv\'en radius. For large enough mass accretion rates and magnetic moments, it is important to take into account not only the pressure of the magnetic field and the radiation pressure inside the disc, but also the pressure of the radiation produced close to the surface of the neutron star in accretion column. The magnetospheric size may increase by up to factor of two as a result of the effects related to the disc thickness and the irradiation from the central source. Accounting for these effects reduces the estimate of the neutron star magnetic moment by a factor of several.

Quasar microlensing effects make it possible to measure the accretion disc sizes around distant supermassive black holes that are still well beyond the spatial resolution of contemporary instrumentation. The sizes measured with this technique appear inconsistent with the standard accretion disc model. Not only are the measured accretion disc sizes larger, but their dependence on wavelength is in most cases completely different from the predictions of the standard model. We suggest that these discrepancies may arise not from non-standard accretion disc structure or systematic errors, as it was proposed before, but rather from scattering and reprocession of the radiation of the disc. In particular, the matter falling from the gaseous torus and presumably feeding the accretion disc may at certain distances become ionized and produce an extended halo that is free from colour gradients. A simple analytical model is proposed assuming that a geometrically thick translucent inflow acts as a scattering mirror changing the apparent spatial properties of the disc. This inflow may be also identified with the broad line region or its inner parts. Such a model is able to explain the basic properties of the apparent disc sizes, primarily their large values and their shallow dependence on wavelength. The only condition required is to scatter significant portion of the luminosity of the disc. This can easily be fulfilled if the scattering inflow has large geometrical thickness and clumpy structure.

The GeV-range spectra of blazars are shaped not only by non-thermal emission processes internal to the relativistic jet but also by external pair-production absorption on the thermal emission of the accretion disc and the broad-line region (BLR). For the first time, we compute here the pair-production opacities in the GeV range produced by a realistic BLR accounting for the radial stratification and radiation anisotropy. Using photoionization modelling with the CLOUDY code, we calculate a series of BLR models of different sizes, geometries, cloud densities, column densities and metallicities. The strongest emission features in the model BLR are Ly$\alpha$ and HeII Ly$\alpha$. Contribution of recombination continua is smaller, especially for hydrogen, because Ly continuum is efficiently trapped inside the large optical depth BLR clouds and converted to Lyman emission lines and higher-order recombination continua. The largest effects on the gamma-ray opacity are produced by the BLR geometry and localization of the gamma-ray source. We show that when the gamma-ray source moves further from the central source, all the absorption details move to higher energies and the overall level of absorption drops because of decreasing incidence angles between the gamma-rays and BLR photons. The observed positions of the spectral breaks can be used to measure the geometry and the location of the gamma-ray emitting region relative to the BLR. Strong dependence on geometry means that the soft photons dominating the pair-production opacity may be actually produced by a different population of BLR clouds than the bulk of the observed broad line emission.

Standard accretion disc model relies upon several assumptions, the most important of which is geometrical thinness. Whenever this condition is violated, new physical effects become important such as radial energy advection and mass loss from the disc. These effects are important, for instance, for large mass accretion rates when the disc approaches its local Eddington limit. In this work, we study the upper limits for standard accretion disc approximation and find the corrections to the standard model that should be considered in any model aiming on reproducing the transition to super-Eddington accretion regime. First, we find that for thin accretion disc, taking into account relativistic corrections allows to increase the local Eddington limit by about a factor of two due to stronger gravity in General Relativity (GR). However, violation of the local Eddington limit also means large disc thickness. To consider consequently the disc thickness effects, one should make assumptions upon the two-dimensional rotation law of the disc. For rotation frequency constant on cylinders $r\sin\theta=const$, vertical gravity becomes stronger with height on spheres of constant radius. On the other hand, effects of radial flux advection increase the flux density in the inner parts of the disc and lower the Eddington limit. In general, the effects connected to disc thickness tend to increase the local Eddington limit even more. The efficiency of accretion is however decreased by advection effects by about a factor of several.

Black holes are spun up by accreting matter and possibly spun-down by magnetic fields. In our work we consider the effect on black hole rotation of the two electromagnetic processes, Blandford-Znajek and Direct Magnetic Link, that differ in their magnetic field configuration. The efficiency of these processes varies with mass accretion rate and accretion regime and generally result in an equilibrium spin parameter in the range from 0.35 to ~0.98. Magnetic field loses its energy while being accreted that may lead to an increase in equilibrium Kerr parameter for the case of advection-dominated disc. We find magnetic field decay decay can decrease electromagnetic term significantly thus increasing the Kerr parameter. We have performed Monte-Carlo simulations for a supermassive black hole population. Our simulations show broad distributions in Kerr parameter (0.1< a <0.98) with a peak at a~0.6. To explain the high observational Kerr parameter values of a>0.9, episodes of supercritical accretion are required. This implication does not however take into account black hole mergers (that play an important role for supermassive black hole evolution).

If accretion disc contains weak frozen-in entangled magnetic fields, their dynamical effect may be important inside the last stable orbit because of the decompression near the sonic point. Here, I consider the radial and vertical structure of a nearly free-falling flow inside the last stable orbit of a thin disc around a Kerr black hole. The thickness of such a flow is determined primarily by the vertical stress created by radial and azimuthal magnetic fields. The thickness is predicted to oscillate vertically around its equilibrium value determined by the magnetic field balance with gravity. For thin discs, this thickness is much larger than that of the accretion disc itself. Numerical simulations with HARM2d show the vertical structure is more complicated. In particular, magnetically supported disc seems to be unstable to segregation of matter into thinner streams with the vertical scale determined by thermal pressure or other processes.

Previously unremarkable, the extragalactic radio source GB 1310+487 showed a gamma-ray flare on 2009 November 18, reaching a daily flux of ~10^-6 photons/cm^2/s at energies E>100 MeV and becoming one of the brightest GeV sources for about two weeks. Its optical spectrum is not typical for a blazar, instead, it resembles those of narrow emission-line galaxies. We investigate changes of the object's radio-to-GeV spectral energy distribution (SED) during and after the prominent GeV flare with the aim to determine the nature of the object and constrain the origin of the variable high-energy emission. The data collected by the Fermi and AGILE satellites at gamma-ray energies, Swift at X-ray and ultraviolet, Kanata, NOT, and Keck telescopes at optical, OAGH and WISE at infrared, and IRAM 30m, OVRO 40m, Effelsberg 100m, RATAN-600, and VLBA at radio, are analysed together to trace the SED evolution on timescales of months. The gamma-ray/radio-loud narrow-line active galactic nucleus (AGN) is located at redshift z=0.638. It is shining through an unrelated foreground galaxy at z=0.500. The AGN light is likely amplified by a factor of a few because of gravitational lensing. The AGN SED shows a two-humped structure typical of blazars and gamma-ray-loud NLSy1 galaxies, with the high-energy (inverse-Compton) emission dominating by more than an order of magnitude over the low-energy (synchrotron) emission during gamma-ray flares. The difference between the two SED humps is smaller during the low-activity state. Fermi observations reveal a strong correlation between the gamma-ray flux and spectral index, with the hardest spectrum observed during the brightest gamma-ray state. If the gamma-ray flux is a mixture of synchrotron self-Compton (SSC) and external Compton (EC) emission, the observed GeV spectral variability may result from varying relative contributions of these two emission components.

The influence of disc radiation capture upon black hole rotational evolution is negligible for radiatively inefficient discs. For the standard thin disc model it is a slight but potentially important effect leading to the equilibrium spin parameter value of about 0.998. For optically thin discs, the fraction of disc radiation captured by the black hole is however about two times larger. In some disc radiation models, inner parts of the accretion flow are optically thin, advection-dominated flows, and the thin disc ends at some transition radius R_{tr}. The thermal energy of the disc stored in trapped radiation is released at this radius. Angular distribution of the radiation released at this radial photosphere facilitates its capture by the black hole. For accretion rates close to critical and disc truncation radius of (2..4) GM/c^2, radiation capture is most efficient in spinning the black hole down that may lead to a_{eq} ~ 0.996..0.997 or less depending on the mass accretion rate. For an accretion flow radiating some constant fraction \epsilon of dissipated energy, the equilibrium Kerr parameter is shown to obey the relation 1-a_{eq} \propto \epsilon^{3/2} as long as 1-a_{eq} << 1. Deviations from Keplerian law near the last stable orbit dominate over the radiation capture effect if they exceed 1..2%.

Microlensing by the stellar population of lensing galaxies provides an important opportunity to spatially resolve the accretion disc structure in strongly lensed quasars. Disc sizes estimated this way are on average larger than the predictions of the standard Shakura-Sunyaev accretion disk model. Analysing the observational data on microlensing variability allows to suggest that some fraction of lensed quasars (primarily, smaller-mass objects) are accreting in super-Eddington regime. Super-Eddington accretion leads to formation of an optically-thick envelope scattering the radiation formed in the disc. This makes the apparent disc size larger and practically independent of wavelength. In the framework of our model, it is possible to make self-consistent estimates of mass accretion rates and black hole masses for the cases when both amplification-corrected fluxes and radii are available.

Gravitational microlensing by the stellar population of lensing galaxies provides an important opportunity to spatially resolve the accretion disk structure in strongly lensed quasars. Some of the objects (like Einstein's cross) are reasonably consistent with the predictions of the standard accretion disk model. In other cases, the size of the emitting region is larger than predicted by the standard thin disk theory and practically independent on wavelength. This may be interpreted as an observational manifestation of an optically-thick scattering envelope possibly related to super-Eddington accretion with outflows.

Though the bulk of the observed optical flux from the discs of intermediate-redshift lensed quasars is formed well outside the region of strong relativistic boosting and light-bending, relativistic effects have important influence on microlensing curves. The reason is in the divergent nature of amplification factors near fold caustics increasingly sensitive to small spatial size details. Higher-order disc images produced by strong light bending around the black hole may affect the amplification curves, making a contribution of up to several percent near maximum amplification. In accordance with theoretical predictions, some of the observed high-amplification events possess fine structure. Here we consider three putative caustic crossing events, one by SBS1520+530 and two events for individual images of the Einstein's cross (QSO J2237+0305). Using relativistic disc models allows to improve the fits, but the required inclinations are high, about 70deg or larger. Such high inclinations apparently contradict the absence of any strong absorption that is likely to arise if a disc is observed edge-on through a dust torus. Still, the high inclinations are required only for the central parts of the disc, that allows the disc itself to be initially tilted by 60-90deg with respect to the black hole and aligned toward the black hole equatorial plane near the last stable orbit radius. For SBS1520+530, an alternative explanation for the observed amplification curve is a superposition of two subsequent fold caustic crossings. While relativistic disc models favour black hole masses ~10^10 solar (several times higher than the virial estimates) or small Eddington ratios, this model is consistent with the observed distribution of galaxies over peculiar velocities only if the black hole mass is about 3 10^8 solar.

We consider the luminous blue variable (LBV) star V532 in M33, also known as Romano's star, in two different spectral states: in the optical minimum of 2007/2008 and during a local brightening in 2005. Optical spectra of low and moderate resolution are modelled using the non-LTE model atmosphere code CMFGEN. All the observed properties of the object in the minimum are well described by a late nitrogen-sequence Wolf-Rayet (WN) star model with relatively high hydrogen abundance (H/He=1.9), while the spectrum during the outburst corresponds to the spectral class WN11 and is similar to the spectrum of P Cyg. The atmosphere is enriched in nitrogen by about a factor of 6 in both states. Most of the heavy element abundances are consistent with the chemical composition of M33. Bolometric luminosity is shown to vary between the two states by a factor of ~1.5. This makes V532 another example of an LBV that shows variations in its bolometric luminosity during an outburst.

We study and analyse low-resolution spectra of the unordinary late WN star FSZ35 in M33. We classify the object as a hydrogen-rich WN8 star. Using the radiative transfer code CMFGEN, we determine the physical parameters of this object and compare them to the parameters of other WN8 stars including the LBV star Romano's star during the minimum of brightness. Unlike Romano's star, the object is fairly stable both spectrally and photometrically, that may be attributed to its more advanced evolutionary stage or lower luminosity. FSZ35 is shown to possess a compact nebula producing faint but detectable [OIII]emission. Location of this object at a large distance (about 100 pc) from the nearest association suggests the object may be one more example of a massive runaway star.

Using the archives of the American Association of Variable Stars Observers and our own data, we analyse the long-term variability of several well-studied Luminous Blue Variables (LBVs) aiming on a general picture of stochastic variability of these objects. The power density spectra of all the selected objects may be generally described by a single power law contaminated by observational noise at higher frequencies. The slopes of the power-law component are close to p=2 (where PDS ~ f^{-p}, and f is frequency) for strongly variable flaring objects like AG Car and significantly smaller (p~1.3) for P Cyg where brightness variation amplitude is <~ 1 mag and dominated by slow low-amplitude variability. The slope holds for about two orders of magnitude in the frequency domain, though peaks and curvatures are present at f ~ 10^{-2}.. 10^{-3} d^{-1}. We show that pseudo-photosphere approach to variability may explain the power-law shape of the variability spectrum at higher frequencies. However, the observed spectra are actually rather "red" than "brown": flux variations are correlated up to tens of years that is much longer than the characteristic refreshment time scales of the pseudo-photosphere. We propose that several stochastic noise components produce the power spectra of LBV stars.

For a long time, the well-known supercritically accreting binary SS433 is being proposed as a prototype for a class of hypothetical bright X-ray sources that may be identified with the so-called Ultraluminous X-ray sources (ULXs) in nearby galaxies or at least with part of them. Like SS433, these objects should be associated with optical nebulae, powered by both radiation of the central source and its wind or jet activity. Indeed, around many ULXs, bright optical nebulae (ULX Nebulae, ULXNe) are found. Here, we use SS433 as a prototype for the power source creating the nebulae around ULXs. Though many factors are important such as the structure of the host star-forming region and the possible supernova remnant formed together with the accreting compact object, we show that most of the properties of ULXNe may be explained by an SS433-like system evolving for up to about one million years in a constant density environment. The basic stages of evolution of a ULX Nebula include a non-spherical HII-region with a central cavity created by non-radiative shock waves, an elongated or bipolar shock-powered nebula created by jet activity and a large-scale quasi-spherical bubble.

We combine archival spectral observations of the LBV star V532 (Romano's star) together with the existing photometric data in the B band. Spectroscopic data cover 15 years of observations (from 1992 to 2007). We show that the object in maximum of brightness behaves as an emission line supergiant while in minimum V532 moves along the sequence of late WN stars. In this sense, the object behaves similarly to the well-known Luminous Blue Variable (LBV) stars AG Car and R127, but is somewhat hotter in the minima. We identify about 100 spectral lines in the 3700..7300 Angstrom wavelength range. For today, our spectroscopy is the most comprehensive for this object. The velocity of the wind is derived using HeI triplet lines (360+/-30 km/s). Physical parameters of the nebula around V532 are estimated.

We present new spectral (FPI and long-slit) data on the Eastern optical filament of the well known radionebula W50 associated with SS433. We find that on sub-parsec scales different emission lines are emitted by different regions with evidently different physical conditions. Kinematical properties of the ionized gas show evidence for moderately high (V ~ 100 km/s) supersonic motions. [OIII]5007 emission is found to be multi-component and differs from lower-excitation [SII]6717 line both in spatial and kinematical properties. Indirect evidence for very low characteristic densities of the gas (n ~ 0.1cm^{-3}) is found. We propose radiative (possibly incomplete) shock waves in low-density, moderately high metallicity gas as the most probable candidate for the power source of the optical filament. Apparent nitrogen over-abundance is better understood if the location of W50 in the Galaxy is taken into account.

In this short paper we show that making turbulence two- rather than three-dimensional may increase effective turbulent viscosity by about 40%. Dimensionless hydrodynamical viscosity parameter up to \alpha\_{max} = 0.25 M_t^2 may be obtained in this approach, that is in better agreement with the observational data on non-stationary accretion than the values obtained in numerical simulations. However, the \alpha-parameter values known from observations are still several times higher.

Highly supercritical accretion discs are probable sources of dense optically thick axisymmetric winds. We introduce a new approach based on diffusion approximation radiative transfer in a funnel geometry and obtain an analytical solution for the energy density distribution inside the wind assuming that all the mass, momentum and energy are injected well inside the spherization radius. This allows to derive the spectrum of emergent emission for various inclination angles. We show that self-irradiation effects play an important role altering the temperature of the outcoming radiation by about 20% and the apparent X-ray luminosity by a factor of 2-3. The model has been successfully applied to two ULXs. The basic properties of the high ionization HII-regions found around some ULXs are also easily reproduced in our assumptions.

We present the results of optical panoramic and long-slit spectroscopy of the nebula MF16 associated with the Ultraluminous X-ray Source NGC6946 ULX-1. More than 20 new emission lines are identified in the spectra. Using characteristic line ratios we find the electron density n_e ~ 600cm^{-3}, electron temperature in the range from ~9000K to ~20 000K (for different diagnostic lines) and the total emitting gas mass M ~ 900 Msolar. We also estimate the interstellar extinction towards the nebula as A_V = 1.m54 somewhat higher than the Galactic absorption. Observed line luminosities and ratios appear to be inconsistent with excitation and ionization by shock waves so we propose the central object responsible for powering the nebula. We estimate the parameters of the ionizing source using photon number estimates and Cloudy modelling. Required EUV luminosity ($\sim 10^{40}$\ergl) is high even if compared with the X-ray luminosity. We argue that independently of their physical nature ULXs are likely to be bright UV and EUV sources. It is shown that the UV flux expected in the GALEX spectral range (1000-3000Angstroms) is quite reachable for UV photometry. Measuring the luminosities and spectral slopes in the UV range may help to distinguish between the two most popular ULX models.

One of the interesting features of Ultraluminous X-ray sources is that many of them are surrounded by luminous nebulae exhibiting diverse observational properties. In different cases the nebulae are photoionized or shock-powered. Generally, the two energy sources appear to coexist. ULX bubble nebulae may be considered a new class of shock-powered nebulae similar to upscaled versions of stellar wind bubbles. Their expansion rates support constant energy influx rather than single powerful events like Hypernova explosions.

We report the results of our observations of the nebular complex MH9/10/11, associated with the ULX HoIX X-1, with scanning Fabry-Perot Interferometer. Two regions differing by their kinematics and line ratios may be distinguished, roughly corresponding to the bubble nebula MH9/10 and fainter HII-region MH11. For MH9/10 we find the expansion rate of 20-70km/s that is different for the approaching and receding parts. MH11 is characterised by very low velocity dispersion (less than or about 15km/s) and nearly constant line-of-sight velocities. Properties of MH11 may be explained by photoionization of gas with hydrogen density of about 0.2cm^-3. Luminosity required for that should be of the order of 10^39erg/s. Similar power source is required to explain the expansion rate of MH9/10. Modelling results also indicate that oxygen abundance in MH11 is about solar.

Optical photometric and spectroscopic data are presented that show an association of an ultraluminous X-ray source in NGC 7331 with a young star cluster of mass 1.1e5 solar masses and age 4.25 Myr. If the ULX is part of the bright stellar cluster, then the mass of the progenitor of the compact accretor must have been greater than about 40-50 solar masses in order to already have evolved through the supernova stage to a compact object. The companion star is also likely an evolved massive star. The emission line spectrum of the nebula surrounding the cluster can be interpreted as a result of photoionization by the cluster OB stars with an additional source of shock excitation producing strong [SII], [OI] and NII lines. This additional source appears to be as much as five times more powerful than the supernovae and stellar winds in the cluster can provide. Additional mechanical energy input associated with the ULX itself can help explain the residual shock excited line luminosities of the emission region.

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

Here we present the results of panoramic and long-slit observations of eight ULX nebular counterparts held with the 6m SAO telescope. In two ULXNe we detected for the first time signatures of high excitation ([OIII]5007 / H\beta > 5). Two of the ULXs were identified with young (T ~ 5-10 Myr) massive star clusters. Four of the eight ULX Nebulae (ULXNe) show bright high-excitation lines. This requires existence of luminous (~ 10^{38} .. 10^{40} erg/s) UV/EUV sources coinciding with the X-ray sources. Other 4 ULXNe require shock excitation of the gas with shock velocities of 20-100km/s. However, all the studied ULXN spectra show signatures of shock excitation, but even those ULXNe where the shocks are prevailing show presence of a hard ionizing source with the luminosity at least ~10^{38} erg/s. Most likely shock waves, X-ray and EUV ionization act simultaneously in all the ULXNe, but they may be roughly separated in two groups, shock-dominated and photoionization-dominated ULXNe. The ULXs have to produce strong winds and/or jets powering their nebulae with \~10^{39} erg/s. Both the wind/jet activity and the EUV source needed are consistent with the suggestion that ULXs are high-mass X-ray binaries with the supercritical accretion disks of the SS433 type.

A radio-IR-optical-X-ray observation campaign of SS 433 has been performed in April 2006, when the jet axis is almost perpendicular to the line of sight. Five flares have been detected during the campaign by radio monitoring observation with RATAN-600. The X-ray astronomical satellite Suzaku observed the source in and out of eclipse. In the X-ray data out of eclipse, the flux shows a significant variation with a time scale of hours. The source seems to be in the active state during the campaign. The observation logs and preliminary results are presented.

SS433 is the only known persistent supercritical accretor, it may be very important for understanding ultraluminous X-ray sources (ULXs) located in external galaxies. We describe main properties of the SS433 supercritical accretion disk and jets. Basing on observational data of SS433 and published 2D simulations of supercritical accretion disks we estimate parameters of the funnel in the disk/wind of SS 433. We argue that the UV radiation of the SS433 disk (~50000 K, ~10^{40}erg/s) is roughly isotropic, but X-ray radiation (~10^7 K, ~10^{40}erg/s) of the funnel is midly anisotropic. A face-on SS433 object has to be ultraluminous in X-rays (10^{40-41}erg/s). Typical time-scales of the funnel flux variability are estimated. Shallow and very broad (0.1-0.3c) and blue-shifted absorption lines are expected in the funnel X-ray spectrum.

We present a study of a peculiar nebula MF16 associated with an Ultraluminous X-ray Source NGC6946 ULX-1. We use integral-field and long-slit spectral data obtained with the 6-m telescope (Russia). The nebula was for a long time considered powered by strong shocks enhancing both high-excitation and low-excitation lines. However, kinematical properties point to rather moderate expansion rates (V ~ 100-200 km/s). The total power of the emission-line source exceeds by one or two orders of magnitude the power observed expansion rate can provide, that points towards the existence of an additional source of excitation and ionization. Using CLOUDY96.01 photoionization code we derive the properties of the photoionizing source. Its total UV/EUV luminosity must be about 10^{40} erg/s.

Five well-known LBV stars in M33 were observed with the Multi-Pupil Fiber Spectrograph (MPFS) on the 6-m Russian telescope. We observed LBVs varA, varB, varC, var2 and var83. In three of them, var2, var83, varB, large-scale nebulae were found with sizes from 15 pc and larger. The nebula shapes are complex, like one-side tails or conical nebulae. They all are related to their LBV stars. In var2 and var83 stars we found radial velocity gradients 15-30 km/s across their nebulae. The stars varA and varC do not show extended nebulae, but nebular lines are certainty present in their spectra.

The ultraluminous X-ray sources (ULXs) were isolated in external galaxies for the last 5 years. Their X-ray luminosities exceed 100-10000 times those of brightest Milky Way black hole binaries and they are extremely variable. There are two models for the ULXs, the best black hole candidates. 1. They are supercritical accretion disks around a stellar mass black hole like that in SS433, observed close to the disk axes. 2. They are Intermediate Mass Black Holes (of 100-10000 solar masses). Critical observations which may throw light upon the ULXs nature come from observations of nebulae around the ULXs. We present results of 3D-spectroscopy of nebulae around several ULXs located in galaxies at 3-6 Mpc distances. We found that the nebulae to be powered by their central black holes. The nebulae are shocked and dynamically perturbed probably by jets. The nebulae are compared with SS433 nebula (W50).

We present a study of a peculiar Supernova Remnant MF16, associated with the Ultraluminous X-ray Source (ULX) NGC6946 ULX-1. Observations were taken with the MultiPupil Fiber Spectrograph (MPFS) with 6-m telescope on January 2005. The nebula is found to be highly asymmetric, one of the parts being much denser and colder. The two-component structure of the emission lines and radial velocity gradient argue for a non-spherical nebula, expanding with a velocity of about 100 km/s. Neither shock models nor the X-ray emission can adequately explain the actual emission line spectrum of MF16, so we suggest an additional ultraviolet source with a luminosity of about 10E40 erg/s. We confirm coincidence of the ULX with the central star, and identify radio emission observed by VLA with the most dense part of the nebula.

We suggest that the ultraluminous X-ray sources located in external galaxies (ULXs) are supercritical accretion disks like that in SS433, observed close to the disk axis. We estimate parameters of the SS433 funnel, where the relativistic jets are formed. Emergent X-ray spectrum in the proposed model of the multicolor funnel (MCF) is calculated. The spectrum can be compared with those of ULXs. We predict a complex absorption-line spectrum with broad and shallow K$\alpha$/Kc blends of the most abundant heavy elements and particular temporal variability. Another critical idea comes from observations of nebulae around the ULXs. We present results of 3D-spectroscopy of nebulae of two ULXs located in Holmberg II and NGC6946. In both cases the nebula is found to be powered by the central black hole. The nebulae are compared with SS433 nebula (W50).

The photometric and spectroscopic observations of nearby type IIP supernova 2004dj are presented. The $^{56}$Ni mass estimated from the light curve is $\approx0.02 M_{\odot}$. This estimate is found to be consistent with the H$\alpha$ luminosity. SN2004dj reveals a strong asymmetry of the H$\alpha$ emission line at the nebular epoch with the shift of the maximum of -1600 km s$^{-1}$. A similar asymmetric component is detected in H$\beta$, [O I] and [Ca II] lines. The line asymmetry is interpreted as a result of the asymmetry of $^{56}$Ni ejecta. The H$\alpha$ profile and its evolution are reproduced in the model of the asymmetric bipolar $^{56}$Ni and spherical hydrogen distributions. The mass of the front $^{56}$Ni jet is comparable to the central component and twice as larger compared to the rear $^{56}$Ni jet. We note that the asymmetric bipolar structure of of $^{56}$Ni ejecta is revealed also by SN1999em, another type IIP supernova.

We present broad band power spectra of variations of SS433 in radio, optical and X-ray spectral bands. We show that at frequencies lower than 1e-5 Hz the source demonstrates the same variability pattern in all these bands. The broad band power spectrum can be fitted by one power law down to frequencies ~1e-7 Hz with flattening afterwards. Such a flattening means that on time scales longer than ~1e7 sec the source variability becomes uncorrelated. This naturally leads to the appearance of quasi-poissonian flares in the source light curve, which have been regularly observed in radio and optical spectral bands. The radio flux power spectrum appears to have a second break at Fourier frequencies ~1e-5 Hz which can be caused by the smearing of the intrinsic radio variability on timescale of the light-crossing time of the radio emitting region. We find a correlation of the radio and optical fluxes of SS433 and the radio flux is delayed by about ~2 days with respect to the optical one. Power spectra of optical and X-ray variabilities continue with the same power law from 1e-7 Hz up to ~0.01-0.05 Hz. The broad band power spectrum of SS433 can be interpreted in terms of self-similar accretion rate modulations in the accretion disk proposed by Lyubarskii (1997) and elaborated by Churazov et al. (2001). We discuss a viscous time-scale in the accretion disk of SS433 in implication to the observed broad band power spectrum.