This short review is dedicated to academician Yakov Borisovich Zeldovich, the science of his epoch and the creation of modern accretion theory.

Approximate gravitational potentials are often used to describe analytically the motion of particles near black holes (BHs), as well as to study the structure of an accretion disk. Such 'pseudo-Newtonian' potentials are used with the flat-metric equations. Here we consider the motion of a free particle near a non-rotating BH in the context of an exact `logarithmic' gravitational potential. We show how the logarithmic potential gives an exact solution for a mechanical problem and present the relativistic Bernoulli equation for the fluid in the Schwarzschild metric.

The observed evolution of the broad-band spectral energy distribution (SED) in NS X-ray Nova Aql X-1 during the rise phase of a bright FRED-type outburst in 2013 can be understood in the framework of thermal emission from unstationary accretion disc with temperature radial distribution transforming from a single-temperature blackbody emitting ring into the multi-colour irradiated accretion disc. SED evolution during the hard to soft X-ray state transition looks curious, as it can not be repro- duced by the standard disc irradiation model with a single irradiation parameter for NUV, Optical and NIR spectral bands. NIR (NUV) band is correlated with soft (hard) X-ray flux changes during the state transition interval, respectively. In our interpreta- tion, at the moment of X-ray state transition UV-emitting parts of the accretion disc are screened from direct X-ray illumination from the central source and are heated primary by hard X-rays (E > 10 keV), scattered in the hot corona or wind possibly formed above the optically-thick outer accretion flow; the outer edge of multi-colour disc, which emits in Optical-NIR, can be heated primary by direct X-ray illumination. We point out that future simultaneous multi-wavelength observations of X-ray Nova systems during the fast X-ray state transition interval are of great importance, as it can serve as 'X-ray tomograph' to study physical conditions in outer regions of accretion flow. This can provide an effective tool to directly test the energy-dependent X-ray heating efficiency, vertical structure and accretion flow geometry in transient LMXBs.

We present the results of 3D-hydrodynamical simulations of accretion flow in the eclipsing dwarf nova V1239 Her in quiescence. The model includes the optical star filling its Roche lobe, a gas stream emanating from the inner Lagrangian point of the binary system, and the accretion disc structure. A cold hydrogen gas stream is initially emitted towards a point-like gravitational centre. A stationary accretion disc is formed in about 15 orbital periods after the beginning of accretion. The model takes into account partial ionization of hydrogen and uses realistic cooling function for hydrogen. The light curve of the system is calculated as the volume emission of optically thin layers along the line of sight up to the optical depth \tau =2/3 calculated using Planck-averaged opacities. The calculated eclipse light curves show good agreement with observations, with the changing shape of pre-eclipse and post-eclipse light curves being explained entirely due to ~ 50% variations in the mass accretion rate through the gas stream.

The theory of quasi-spherical subsonic accretion onto magnetized rotating neutron star is reviewed. Different regimes of quasi-spherical accretion onto a neutron star: supersonic (Bondi) accretion, which takes place when the captured matter cools down rapidly and falls supersonically towards the neutron-star magnetosphere, and subsonic (settling) accretion which occurs when the plasma remains hot until it meets the magnetospheric boundary. In subsonic accretion, which works at X-ray luminosities $\lesssim 4\times 10^{36}$~erg~s$^{-1}$, a hot quasi-spherical shell must form around the magnetosphere, and the actual accretion rate onto the neutron star is determined by the ability of the plasma to enter the magnetosphere due to the Rayleigh-Taylor instability. We show how the dimensionless parameters of the theory can be determined from observations of equilibrium X-ray pulsars (Vela X-1, GX 301-2). We also discuss how in the settling accretion theory bright X-ray flares ($\sim 10^{38}-10^{40}$~ergs) observed in supergiant fast X-ray transients (SFXT) may be produced by sporadic capture of magnetized stellar-wind plasma. At sufficiently low accretion rates, magnetic reconnection can enhance the magnetospheric plasma entry rate, resulting in copious production of X-ray photons, strong Compton cooling and ultimately in unstable accretion of the entire shell. A bright flare develops on the free-fall time scale in the shell, and the typical energy released in an SFXT bright flare corresponds to the mass of the shell.

The vertical structure of stationary thin accretion discs is calculated from the energy balance equation with heat generation due to microscopic ion viscosity {\eta} and electron heat conductivity {\kappa}, both depending on temperature. In the optically thin discs it is found that for the heat conductivity increasing with temperature, the vertical temperature gradient exceeds the adiabatic value at some height, suggesting convective instability in the upper disc layer. There is a critical Prandtl number, Pr = 4/9, above which a Keplerian disc become fully convective. The vertical density distribution of optically thin laminar accretion discs as found from the hydrostatic equilibrium equation cannot be generally described by a polytrope but in the case of constant viscosity and heat conductivity. In the optically thick discs with radiation heat transfer, the vertical disc structure is found to be convectively stable for both absorption dominated and scattering dominated opacities, unless a very steep dependence of the viscosity coefficient on temperature is assumed. A polytropic-like structure in this case is found for Thomson scattering dominated opacity.

The RXTE observed four outbursts of the accreting X-ray binary transient source, GX 304-1 in 2010 and 2011. We present results of detailed 3-100 keV spectral analysis of 69 separate observations, and report a positive correlation between cyclotron line parameters, as well as other spectral parameters, with power law flux. The cyclotron line energy, width and depth versus flux, and thus luminosity, correlations show a flattening of the relationships with increasing luminosity, which are well described by quasi-spherical or disk accretion that yield the surface magnetic field to be ~60 keV. Since HEXTE cluster A was fixed aligned with the PCA field of view and cluster B was fixed viewing a background region 1.5 degrees off of the source direction during these observations near the end of the RXTE mission, the cluster A background was estimated from cluster B events using HEXTEBACKEST. This made possible the detection of the ~55 keV cyclotron line and an accurate measurement of the continuum. Correlations of all spectral parameters with the primary 2-10 keV power law flux reveal it to be the primary driver of the spectral shape. The accretion is found to be in the collisionless shock braking regime.

Using the anelastic approximation of linearized hydrodynamic equations, we investigate the development of axially symmetric small perturbations in thin Keplerian discs. Dispersion relation is found as a solution of general Sturm-Liouville eigenvalue problem for different values of relevant physical parameters (viscosity, heat conductivity, disc semithickness). The analysis reveals the appearance of overstable mode for Prandtl parameter higher than some critical value. These modes have a viscous-convective nature and can serve as a seed for turbulence in astrophysical discs even in the absence of magnetic fields.

We present a model of the non-stationary $\alpha$-disk with account for the irradiation and the vertical convection in the outer accretion disk where hydrogen is partially ionized. We include the viscous energy generation in the mix-length convection equations in accretion disks. The optical and X-ray light curves of X-ray nova A0620-00 are investigated in terms of this model. The turbulent viscosity parameter of the accretion disk is estimated, $\alpha = 0.5 \div 0.6$, which is necessary to explain the luminosity decay rate on the descending branch of the X-ray light curve for the A0620-00 1975 outburst. The secondary luminosity maximum on the light curves is explained by assuming an additional injection of matter into the accretion disk from the optical companion.

Using RXTE/ASM archival data, we investigate the behaviour of the spectral hardness ratio as a function of X-ray luminosity in a sample of six transient X-ray pulsars (EXO 2030+375, GX 304-1, 4U 0115+63, V 0332+63, A 0535+26 and MXB 0656-072). In all sources we find that the spectral hardness ratio defined as $F_{5-12\mathrm{keV}}/ F_{1.33-3\mathrm{keV}}$ increases with the ASM flux (1.33--12 keV) at low luminosities and then saturates or even slightly decreases above some critical X-ray luminosity falling into the range $\sim(3-7)\times10^{37}$~erg~s$^{-1}$. Two-dimensional structure of accretion columns in the radiation-diffusion limit is calculated for two possible geometries (filled and hollow cylinder) for mass accretion rates $\dot M$ ranging from $10^{17}$ to 1.2$\times 10^{18}$~g s$^{-1}$. The observed spectral behaviour in the transient X-ray pulsars with increasing $\dot M$ can be reproduced by a Compton saturated sidewall emission from optically thick magnetized accretion columns with taking into account the emission reflected from the neutron star atmosphere. At $\dot M$ above some critical value $\dot M_{cr}\sim (6-8)\times 10^{17}$~g~s$^{-1}$, the hight of the column becomes such that the contribution of the reflected component to the total emission starts decreasing, which leads to the saturation and even slight decrease of the spectral hardness. Hollow-cylinder columns have a smaller height than the filled-cylinder ones, and the contribution of the reflected component in the total emission does not virtually change with $\dot M$ (and hence the hardness of the continuum monotonically increases) up to higher mass accretion rates than $\dot M_{cr}$ for the filled columns.

Using the anelastic approximation of linearised hydrodynamic equations, we investigate the development of axially symmetric small perturbations in thin Keplerian discs. The sixth-order dispersion equation is derived and numerically solved for different values of relevant physical parameters (viscosity, heat conductivity, disc semi-thickness and vertical structure). The analysis reveals the appearance of two overstable modes which split out from the classical Rayleigh inertial modes in a wide range of the parameters in both ionized and neutral gases. These modes have a viscous-convective nature and can serve as a seed for turbulence in astrophysical discs even in the absence of magnetic fields.

Conditions of Velikhov-Chandrasekhar magneto-rotational instability in ideal and non-ideal plasmas are examined. Linear WKB analysis of hydromagnetic axially symmetric flows shows that in the Rayleigh-unstable hydrodynamic case where the angular momentum decreases with radius, the MRI branch becomes stable, and the magnetic field suppresses the Rayleigh instability at small wavelengths. We investigate the limiting transition from hydromagnetic flows to hydrodynamic flows. The Rayleigh mode smoothly transits to the hydrodynamic case, while the Velikhov-Chandrasekhar MRI mode completely disappears without the magnetic field. The effects of viscosity and magnetic diffusivity in plasma on the MRI conditions in thin accretion discs are studied. We find the limits on the mean free-path of ions allowing MRI to operate in such discs.

A viscous instability in shearing laminar axisymmetric hydrodynamic flows around a gravitating center is described. In the linearized hydrodynamic equations written in the Boussinesq approximation with microscopic molecular transport coefficients, the instability arises when the viscous dissipation is taken into account in the energy equation. Using the local WKB approximation, we derive a third-order algebraic dispersion equation with two modes representing the modified Rayleigh modes R+ and R-, and the third X-mode. We show that in thin accretion flows the viscosity destabilizes one of the Rayleigh modes in a wide range of wavenumbers, while the X-mode always remains stable. In Keplerian flows, the instability increment is found to be a few Keplerian rotational periods at wavelengths with $kr\sim 10-50$. This instability may cause turbulence in astrophysical accretion discs even in the absence of magnetic field.

This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of high-mass X-ray binaries and ultra-luminous X-ray sources. For a summary, we refer to the paper.

We apply the model of subsonic settling accretion on to isolated neutron stars accreting from the interstellar medium (AINS). We show that in this regime the expected mean X-ray luminosity from AINS turns out to be 2-3 orders of magnitude as small as the maximum possible Bondi value, i.e. $10^{27}-10^{28}$~erg s$^{-1}$. The intrinsically unstable character of settling accretion due to long plasma cooling time leads to regular appearance of X-ray flares with a duration of about one hour and a maximum luminosity of about the Bondi value, $\sim10^{31}$~erg~s$^{-1}$. This feature can be used to distinguished AINS from other dim X-ray sources. With the sensitivity of the forthcoming all-sky X-ray surveys the expected number of the potentially detectable AINS can be from a few to ten.

In the settling accretion theory, which is applicable to quasi-spherical accreting slowly rotating magnetized neutron stars with X-ray luminosity $L_x\lesssim 4\times 10^{36}$~erg/s, bright X-ray flares ($\sim 10^{38}-10^{40}$~ergs) observed in supergiant fast X-ray transients (SFXT) may be produced by sporadic capture of magnetized stellar-wind plasma from the early-type supergiant. At sufficiently low steady accretion rates ($\lesssim 10^{15}$~g/s) through the shell around the neutron star magnetosphere at the settling accretion stage, magnetic reconnection can temporarily enhance the magnetospheric plasma entry rate, resulting in copious production of X-ray photons, strong Compton cooling, and ultimately in unstable accretion of the entire shell. A bright flare develops on the free-fall time scale in the shell, $R_B^{3/2}/\sqrt{GM}\sim 10^3-10^4$~s ($R_B$ is the classical Bondi capture radius), and the typical energy released in an SFXT bright flare corresponds to the mass of the shell.

Aims. We investigate the long-term evolution of the Cyclotron Resonance Scattering Feature (CRSF) in the spectrum of the binary X-ray pulsar Her X-1 and present evidence of a true long-term decrease in the centroid energy Ecyc of the cyclotron line in the pulse phase averaged spectra from 1996 to 2012. Methods. Our results are based on repeated observations of Her X-1 by those X-ray observatories capable of measuring clearly beyond the cyclotron line energy of about 40 keV. Results. The historical evolution of the pulse phase averaged CRSF centroid energy Ecyc since its discovery in 1976 is characterized by an initial value around 35 keV, an abrupt jump upwards to beyond about 40 keV between 1990 and 1994, and an apparent decay thereafter. Much of this decay, however, was found to be due to an artifact, namely a correlation between Ecyc and the X-ray luminosity Lx discovered in 2007. In observations after 2006, however, we now find a statistically significant true decrease in the cyclotron line energy. At the same time, the dependence of Ecyc on X-ray luminosity is still valid with an increase of about 5% in energy for a factor of two increase in luminosity. A decrease in Ecyc by 4.2 keV over the 16 years from 1996 to 2012 can either be modeled by a linear decay, or by a slow decay until 2006 followed by a more abrupt decrease thereafter. Conclusions. We speculate that the physical reason could be connected to a geometric displacement of the cyclotron resonant scattering region in the polar field or to a true physical change in the magnetic field configuration at the polar cap by the continued accretion.

We analyze spin-up/spin-down of the neutron star in Be X-ray binary system GX\,304-1 observed by \textit{Swift}/XRT and \textit{Fermi}/GBM instruments in the period of the source activity from April 2010 to January 2013 and discuss possible mechanisms of angular momentum transfer to/from the neutron star. We argue that the neutron star spin-down at quiescent states of the source with an X-ray luminosity of $L_x\sim 10^{34}$~erg s$^{-1}$ between a series of Type I outbursts and spin-up during the outbursts can be explained by quasi-spherical settling accretion onto the neutron star. The outbursts occur near the neutron star periastron passages where the density is enhanced due to the presence of an equatorial Be-disc tilted to the orbital plane. We also propose an explanation to the counterintuitive smaller spin-up rate observed at higher luminosity in a double-peak Type I outburst due to lower value of the specific angular momentum of matter captured from the quasi-spherical wind from the Be-star by the neutron star moving in an elliptical orbit with eccentricity $e\gtrsim 0.5$.

A review of wind accretion in HMXB is presented. We focus on different regimes of quasi-spherical accretion onto a NS: supersonic (Bondi) accretion, which takes place when the captured matter cools down rapidly and falls supersonically towards the NS magnetosphere, and subsonic (settling) accretion which occurs when the plasma remains hot until it meets the magnetospheric boundary. The two regimes are separated by a limit in X-ray luminosity at about 4 10^{36} erg/s. In subsonic accretion, which works a hot quasi-spherical shell must form around the magnetosphere, and the actual accretion rate onto the NS is determined by the ability of the plasma to enter the magnetosphere due to the Rayleigh-Taylor instability. Two regimes of subsonic accretion are possible, depending on the plasma cooling mechanism (Compton or radiative) near the magnetopshere. The transition from the high-luminosity regime with Compton cooling to the low-luminosity (L_x < 3\times 10^35 erg/s) regime with radiative cooling can be responsible for the onset of the 'off' states repeatedly observed in several X-ray pulsars, such as Vela X-1, GX 301-2 and 4U 1907+09. The triggering of the transition may be due to a switch in the X-ray beam pattern in response to a change in the optical depth in the accretion column with changing luminosity. We also show that in the settling accretion theory, bright X-ray flares (10^{38}-10^{40} ergs) observed in SFXT may be produced by sporadic capture of magnetized stellar-wind plasma. At sufficiently low accretion rates, magnetic reconnection can enhance the magnetospheric plasma entry rate, resulting in copious production of X-ray photons, strong Compton cooling and ultimately in unstable accretion of the entire shell. A bright flare develops on the free-fall time scale in the shell, and the typical energy released in an SFXT bright flare corresponds to the mass of the shell.

At steady low-luminosity states, Supergiant Fast X-ray Transients (SFXTs) can be at the stage of quasi-spherical settling accretion onto slowly rotating magnetized NS from the OB-companion winds. At this stage, a hot quasi-static shell is formed above the magnetosphere, the plasma entry rate into magnetosphere is controlled by (inefficient) radiative plasma cooling, and the accretion rate onto the NS is suppressed by a factor of \sim 30 relative to the Bondi-Hoyle-Littleton value. Changes in the local wind velocity and density can only slightly increase the mass accretion rate (a factor of \sim 10) bringing the system into the Compton cooling dominated regime and led to the production of moderately bright flares (L_x\lesssim 10^{36} erg/s). To interpret the brightest flares (L_x>10^{36}~erg/s) displayed by the SFXTs, we propose that a larger increase in the mass accretion rate can be produced by sporadic capture of magnetized stellar wind plasma. At sufficiently low accretion rates, magnetic reconnection can enhance the magnetospheric plasma entry rate, resulting in copious production of X-ray photons, strong Compton cooling and ultimately in unstable accretion of the entire shell. A bright flare develops on the free-fall time scale in the shell, and the typical energy released in an SFXT bright flare corresponds to the mass of the shell. This view is consistent with the energy released in SFXT bright flares (\sim 10^{38}-10^{40}~ergs), their typical dynamic range (\sim 100), and with the observed dependence of these characteristics on the average unflaring X-ray luminosity of SFXTs. Thus the flaring behavior of SFXTs, as opposed to steady HMXBs, may be primarily related to their low X-ray luminosity allowing sporadic magnetic reconnection to occur during magnetized plasma entry into the magnetosphere.

Accretion of matter onto the surface of a freely precessing NS with a complex non-dipole magnetic field can explain the change of X-ray pulse profiles of Her X-1 observed by RXTE with the phase of the 35-day cycle. We demonstrate this using all available measurements of X-ray pulse profiles in the 9-13 keV energy range obtained with the RXTE/PCA. The measured profiles guided the elaboration of a geometrical model and the definition of locations of emitting poles, arcs, and spots on the NS surface which satisfactorily reproduce the observed pulse profiles and their dependence on free precession phase. We have found that the observed trend of the times of the 35-day turn-ons on the O-C diagram, which can be approximated by a collection of consecutive linear segments around the mean value, can be described by our model by assuming a variable free precession period, with a fractional period change of about a few percent. We propose that the 2.5% changes in the free precession period that occur on time scales of several to tens of 35-day cycles can be related to wandering of the principal inertia axis of the NS body due to variations in the patterns of accretion onto the NS surface. The closeness of periods of the disk precession and the NS free precession can be explained by the presence of a synchronization mechanism in the system, which modulates the dynamical interaction of the gas streams and the accretion disk with the NS free precession period.

Strong magnetic field of accreting neutron stars ($10^{14}$ G) is hard to probe by X-ray spectroscopy but can be indirectly inferred from spin-up/spin-down measurement in X-ray pulsars. The existing observations of slowly rotating X-ray pulsars are discussed. It is shown that magnetic fields of neutron stars derived from these observations (or lower limits in some cases) fall within the standard $10^{12}$-$10^{13}$ G range. Claims about the evidence for accreting magnetars are critically discussed in the light of recent progress in understanding of accretion onto slowly rotating neutron stars in the subsonic regime.

A review of wind accretion in high-mass X-ray binaries is presented. We focus attention to different regimes of quasi-spherical accretion onto the neutron star: the supersonic (Bondi) accretion, which takes place when the captured matter cools down rapidly and falls supersonically toward NS magnetospghere, and subsonic (settling) accretion which occurs when plasma remains hot until it meets the magnetospheric boundary. Two regimes of accretion are separated by an X-ray luminosity of about $4\times10^{36}$ erg/s. In the subsonic case, which sets in at low luminosities, a hot quasi-spherical shell must be formed around the magnetosphere, and the actual accretion rate onto NS is determined by ability of the plasma to enter the magnetosphere due to Rayleigh-Taylor instability. We calculate the rate of plasma entry the magnetopshere and the angular momentum transfer in the shell due to turbulent viscosity appearing in the convective differentially rotating shell. We also discuss and calculate the structure of the magnetospheric boundary layer where the angular momentum between the rotating magnetosphere and the base of the differentially rotating quasi-spherical shell takes place. We show how observations of equilibrium X-ray pulsars Vela X-1 and GX 301-2 can be used to estimate dimensionless parameters of the subsonic settling accretion theory, and obtain the width of the magnetospheric boundary layer for these pulsars.

Quasi-spherical subsonic accretion onto slowly rotating magnetized NS is considered, when the accreting matter settles down subsonically onto the rotating magnetosphere, forming an extended quasi-static shell. The shell mediates the angular momentum transfer to/from the rotating NS magnetosphere by large-scale convective motions, which lead to an almost iso-angular-momentum rotation law inside the shell. The accretion rate through the shell is determined by the ability of the plasma to enter the magnetosphere due to Rayleigh-Taylor instability while taking cooling into account. The settling regime of accretion is possible for moderate X-ray luminosities L <4 10^36 erg/s. At higher luminosities a free-fall gap above the NS magnetosphere appears due to rapid Compton cooling, and accretion becomes highly non-stationary. From observations of spin-up/spin-down rates of wind accreting equilibrium XPSRs with known orbital periods (GX 301-2, Vela X-1), the main dimensionless parameters of the model and be determined and the NS magnetic field can be estimated. For equilibrium pulsars with independent measurements of the magnetic field the velocity of the stellar wind can be estimated without the use of complicated spectroscopic measurements. For non-equilibrium pulsars, a maximum possible value of the spin-down rate of the accreting neutron star exists. From observations of the spin-down rate and the X-ray luminosity in such pulsars (e. g. GX 1+4, SXP 1062 and 4U 2206+54) a lower limit on the neutron star magnetic field is obtained, which in all cases turns out to be close to the standard 10^12-10^13 G value, in agreement with cyclotron line measurements. The model explains both the spin-up/spin-down of the pulsar frequency on large time-scales and the irregular short-term frequency fluctuations, which may correlate or anti-correlate with the X-ray luminosity fluctuations.

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.

We performed a detailed pulse phase resolved spectroscopy of the accreting binary X-ray pulsar Her X-1 in the energy range 3.5-75 keV and have established pulse phase profiles for all spectral parameters. For the centroid of the cyclotron line, the photon index and the flux of the 6.4 keV iron line, we have studied the variation as a function of 35 d phase. We analyzed RXTE observations of the Main-On of November 2002. Four different time intervals of about 1 d duration were selected to provide a good coverage of a complete Main-On. The intervals are centered at 35 d phase 0.03, 0.10, 0.15, and 0.20, respectively. All spectral parameters show a strong modulation with pulse phase. While the centroid energy of the cyclotron line follows roughly the shape of the pulse profile, both the photon index and the iron line intensity exhibit distinct minima around the peak of the X-ray pulse. With respect to variations of the observed profiles with 35 d phase, we find that there is a clear evolution of the shape of the pulse profiles (flux versus pulse phase), a moderate increase of the maximum cyclotron line energy (found around pulse phase 0.7), but no significant evolution of the shape of the pulse phase profiles of the cyclotron line energy, the spectral power law index or the iron line intensity. The variation of spectral parameters as a function of the pulse phase provides important information about the system: 1. the disappearance of the Fe line flux near the highest continuum flux may be an indication of a hollow cone geometry of the accretion structure; ii. the apparent non-dependence of the cyclotron line energy profiles on 35 d phase provides a new possibility to test the model of free precession of the neutron star, proposed to be responsible for the systematic variations in the pulse profiles.

The accreting X-ray pulsar Her X-1 shows two types of long-term variations, both with periods of ~35 days: 1) Turn-on cycles, a modulation of the flux}, with a ten-day long Main-On and a five-day long Short-On, separated by two Off-states, and 2) a systematic variation in the shape of the 1.24 s pulse profile. While there is general consensus that the flux modulation is due to variable shading of the X-ray emitting regions on the surface of the neutron star by the precessing accretion disk, the physical reason for the variation in the pulse profiles has remained controversial. Following the suggestion that free precession of the neutron star may be responsible for the variation in the pulse profiles, we developed a physical model of strong feedback interaction between the neutron star and the accretion disk in order to explain the seemingly identical values for the periods of the two types of variations, which were found to be in basic synchronization. In a deep analysis of pulse profiles observed by several different satellites over the last three decades we now find that the clock behind the pulse profile variations shows exactly the same erratic behavior as the turn-on clock, even on short time scales (a few 35 d cycles), suggesting that there may in fact be only one 35 d clock in the system. If this is true, it raises serious questions with respect to the idea of free precession of the neutron star, namely how the neutron star can change its precessional period every few years by up to 2.5% and how the feedback can be so strong, such that these changes can be transmitted to the accretion disk on rather short time scales.

Quasi-spherical subsonic accretion can be realized in slowly rotating wind-fed X-ray pulsars (XPSRs) at X-ray luminosities <4 10^{36} erg/s. In this regime the accreting matter settles down subsonically onto the rotating magnetosphere, forming an extended quasi-static shell. The shell mediates the angular momentum removal from the rotating NS magnetosphere by shear turbulent viscosity in the boundary layer or via large-scale convective motions. In the last case the differential rotation law in the shell is close to iso-angular-momentum rotation. The accretion rate through the shell is determined by the ability of the plasma to enter the magnetosphere due to Rayleigh-Taylor instabilities while taking cooling into account. Measurements of spin-up/spin-down rates of quasi-spherically wind accreting XPSRs in equilibrium with known orbital periods (like e.g. GX 301-2 and Vela X-1) enable determination of the main dimensionless parameters of the model and the NS magnetic field. For equilibrium pulsars with independent measurements of the magnetic field, the stellar wind velocity from the companion can be estimated without the use of complicated spectroscopic measurements. For non-equilibrium pulsars, a maximum possible spin-down torque exerted on the accreting NS exists. From observations of the spin-down rate and X-ray luminosity in such pulsars (GX 1+4, SXP 1062, 4U 2206+54, etc.) a lower limit on the NS magnetic field is derived, which in all cases turns out to be close to the standard one and in agreement with cyclotron line measurements. The model explains the existence of super slowly rotating XPSRs without the need to hypothesize on additional accretion properties and magnetar-like magnetic fields in accreting neutron stars.

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.

We elaborate on a recently proposed model for subsonic quasi-spherical accretion onto slowly rotating pulsars, in which accretion is mediated through a hot quasi-static shell above the neutron star magnetosphere. We show that under the same external conditions, two regimes of subsonic accretion are possible, depending on if plasma cooling in the transition zone is dominated by Compton or radiative processes. We suggest that a transition from the higher luminosity Compton cooling regime to the lower luminosity radiative cooling regime can be responsible for the onset of the `off'-states repeatedly observed in several low luminosity slowly accreting pulsars, such as Vela X-1, GX 301-2 and 4U 1907+09. We further suggest that the triggering of the transition may be due to a switch in the X-ray beam pattern in response to a change in the optical depth in the accretion column with changing luminosity.

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.

Quasi-spherical accretion in wind-fed X-ray pulsars is discussed. At X-ray luminosities <4 10^{36} erg/s, a hot convective shell is formed around the neutron star magnetosphere, and subsonic settling accretion regime sets in. In this regime, accretion rate onto neutron star is determined by the ability of plasma to enter magnetosphere via Rayleigh-Taylor instability. A gas-dynamic theory of settling accretion is constructed taking into account anisotropic turbulence. The angular momentum can be transferred through the quasi-static shell via large-scale convective motions initiating turbulence cascade. The angular velocity distribution in the shell is found depending on the turbulent viscosity prescription. Comparison with observations of long-period X-ray wind-fed pulsars shows that an almost iso-angular-momentum distribution is most likely realized in their shells. The theory explains long-term spin-down in wind- fed accreting pulsars (e.g. GX 1+4) and properties of short-term torque-luminosity correlations. The theory can be applied to slowly rotating low-luminosity X-ray pulsars and non-stationary accretion phenomena observed in some SFXTs.

The accreting X-ray pulsar Her X-1 shows two types of long-term variations, both with a period of 35 days: 1) A modulation of the flux with a ten day long Main-On and a 5 d long Short-On, separated by two Off-states, and 2) A systematic variation of the shape of the 1.24 s pulse profile. While there is general consensus that the flux modulation is due to variable shading of the X-ray emitting regions on the surface of the neutron star by the precessing accretion disk, the physical reason for the variation of the pulse profiles had remained controversial. Following the suggestion by Truemper et al. (1986) that free precession of the neutron star may be responsible for the variation of the pulse profiles, we had developed physical models of strong feedback interaction between the neutron star and the accretion disk in order to explain the seemingly identical values for the periods of the two types of variations. In a deep analysis of pulse profiles observed by several different satellites over the last three decades we find now that the clock behind the pulse profile variations shows exactly the same erratic behavior as the turn-on clock, even on short time scales (100 d), suggesting that there may in fact be only one 35 d clock in the system. If this is true, it presents a serious challenge for the idea of free precession of the neutron star and calls for alternative physical models for the variation in pulse shape.

A theoretical model for quasi-spherical subsonic accretion onto slowly rotating magnetized neutron stars is constructed. In this model the accreting matter subsonically settles down onto the rotating magnetosphere forming an extended quasi-static shell. This shell mediates the angular momentum removal from the rotating neutron star magnetosphere during spin-down episodes by large-scale convective motions. The accretion rate through the shell is determined by the ability of the plasma to enter the magnetosphere. The settling regime of accretion can be realized for moderate accretion rates $\dot M< \dot M_*\simeq 4\times 10^{16}$ g/s. At higher accretion rates a free-fall gap above the neutron star magnetosphere appears due to rapid Compton cooling, and accretion becomes highly non-stationary. From observations of the spin-up/spin-down rates (the angular rotation frequency derivative $\dot \omega^*$, and $\partial\dot\omega^*/\partial\dot M$ near the torque reversal) of X-ray pulsars with known orbital periods, it is possible to determine the main dimensionless parameters of the model, as well as to estimate the magnetic field of the neutron star. We illustrate the model by determining these parameters for three wind-fed X-ray pulsars GX 301-2, Vela X-1, and GX 1+4. The model explains both the spin-up/spin-down of the pulsar frequency on large time-scales and the irregular short-term frequency fluctuations, which can correlate or anti-correlate with the X-ray flux fluctuations in different systems. It is shown that in real pulsars an almost iso-angular-momentum rotation law with $\omega \sim 1/R^2$, due to strongly anisotropic radial turbulent motions sustained by large-scale convection, is preferred.

The properties of wind accretion in symbiotic X-ray binaries (SyXBs) consisting of red-giant and magnetized neutron star (NS) are discussed. The spin-up/spin-down torques applied to NS are derived based on a hydrodynamic theory of quasi-spherical accretion onto magnetized NSs. In this model, a settling subsonic accretion proceeds through a hot shell formed around the NS magnetosphere. The accretion rate onto the NS is determined by the ability of the plasma to enter the magnetosphere.Due to large Reynolds numbers in the shell, the interaction of the rotating magnetosphere with plasma initiates a subsonic turbulence. The convective motions are capable of carrying the angular momentum through the shell. We carry out a population synthesis of SyXBs in the Galaxy with account for the spin evolution of magnetized NS. The Galactic number of SyXBs with bright (M_v<1) low-mass red-giant companion is found to be from \sim 40 to 120, and their birthrate is \sim 5\times 10^{-5}-10^{-4} per year. According to our model, among known SyXBs, Sct X-1 and IRXS J180431.1-273932 are wind-fed accretors. GX 1+4 lies in the transition from the wind-fed SyXBs to SyXBs in which the giants overflow their Roche lobe. The model successfully reproduces very long NS spins (such as in IGR J16358-4724 and 4U 1954+31) without the need to invoke very strong magnetic fields.

We have investigated the influence of X-ray irradiation on the vertical structure of the outer accretion disk in low-mass X-ray binaries by performing a self-consistent calculation of the vertical structure and X-ray radiation transfer in the disk. Penetrating deep into the disk, the field of scattered X-ray photons with energy $E\gtrsim10$\,keV exerts a significant influence on the vertical structure of the accretion disk at a distance $R\gtrsim10^{10}$\,cm from the neutron star. At a distance $R\sim10^{11}$\,cm, where the total surface density in the disk reaches $\Sigma_0\sim20$\,g\,cm$^{-2}$, X-ray heating affects all layers of an optically thick disk. The X-ray heating effect is enhanced significantly in the presence of an extended atmospheric layer with a temperature $T_{atm}\sim(2\div3)\times10^6$\,K above the accretion disk. We have derived simple analytic formulas for the disk heating by scattered X-ray photons using an approximate solution of the transfer equation by the Sobolev method. This approximation has a $\gtrsim10$\,% accuracy in the range of X-ray photon energies $E<20$\,keV.

Monitoring of pulse period variations in accreting binary pulsars is an important tool to study the interaction between the magnetosphere of the neutron star and the accretion disk. While the X-ray flux of the brightest X-ray pulsars have been successfully monitored over many years (e.g. with RXTE/ASM, CGRO/BATSE, Swift/BAT), the possibility to monitor their pulse timing properties continuously has so far been very limited. In our work we use Swift/BAT observations to study one of the most enigmatic X-ray pulsars, Hercules X-1. For the first time, a quasi-continuous monitoring of the pulse period and the pulse period derivative of Her X-1, is achieved over a long time (> 4 yrs). We argue that together with the long-term decrease of the orbital period in Her X-1 the measured pulse period behaviour requires the presence of mass ejection from the inner parts of the accretion disk along the open magnetic field lines. The mass ejection episodes probably take place during strong spin-down episodes which are associated with the low X-ray luminosity.

We present evidence for an identical behavior of the precession of the accretion disk and that of the neutron star (NS) in Her X-1, based on investigating the well established 35 day modulation in Her X-1 in two different ways: 1) following the turn-ons, thought to be due to the precession of the accretion disk, and 2) following the re-appearance of the shape of the pulse profiles, which we assume to be due to precession of the NS. The turn-on evolution and the evolution of the phase-zero values of the precessing NS (as determined from the pulse profiles) track each other very closely. Since the turn-on evolution is strongly correlated with the pulse period evolution, this means that there is also a strong correlation between the spin and the precession of the NS. There is a very strong physical coupling between the NS and the accretion disk, we suggest through physical feedback in the binary system. The apparent long-term stability of the 35 d clock may be due to the interior of the NS, the coupling of which to the observable surface effects is of general importance for the physics of super-dense, highly magnetized material.

We re-examine the infrared counterpart of the dipping low-mass x-ray binary 4U1323-619. New X-ray data available from the XMM and Chandra observatories combined with archival IR observations from the ESO 3.6m New Technology Telescope allow us to define a new possible counterpart. We present here its photometric properties and compare them with a simple analytical model of an accretion disc illuminated by the hot central corona known to be present in the binary system.

Context: Monitoring of pulse period variations in accreting binary pulsars is an important tool to study the interaction between the magnetosphere of the neutron star and the accretion disk. While the X-ray flux of the brightest X-ray pulsars have been successfully monitored over many years (e.g. with RXTE/ASM, CGRO/BATSE, Swift/BAT), the possibility to monitor their pulse timing properties continuously has so far been very limited. Aims: In our work we show that the Swift/BAT observations can be used to monitor coherent pulsations of bright X-ray sources and use the Swift archival data to study one of the most enigmatic X-ray pulsars, Hercules X-1. A quasi-continuous monitoring of the pulse period and the pulse period derivative of an X-ray pulsar, here Her X-1, is achieved over a long time (<~ 4 yrs). We compare our observational results with predictions of accretion theory and use them to test different aspects of the physical model of the system. Methods: In our analysis we use the data accumulated with Swift/BAT starting from the beginning of 2005 (shortly after launch) until the present time. To search for pulsations and for their subsequent analysis we used the count rate measured by the BAT detector in the entire field of view. Results: The slope of the correlation between the locally determined spin-up rate and the X-ray luminosity is measured for Her X-1 and found to be in agreement with predictions of basic accretion torque theory. The observed behaviour of the pulse period together with the previously measured secular decrease of the system's orbital period is discussed in the frame of a model assuming ejection of matter close to the inner boundary of the accretion disk.

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.

We fit outbursts of two X-ray novae (Nova Monocerotis 1975=A0620-00 and Nova Muscae GS 1991=1124-683) using a time-dependent accretion disk model. The model is based on a new solution for a diffusion-type equation for the non-steady-state accretion and describes the evolution of a viscous alpha-disk in a binary system after the peak of an outburst, when matter in the disk is totally ionized. The accretion rate in the disk decreases according to a power law. We derive formulas for the accretion rate and effective temperature of the disk. The model has three free input parameters: the mass of the central object M, the turbulence parameter alpha, and the normalization parameter delta t. Results of the modeling are compared with the observed X-ray and optical B and V light curves. The resulting estimates for the turbulence parameter $\alpha$ are similar: 0.2-0.4 for A 0620-00 and 0.45-0.65 for GS 1124-683, suggesting a similar nature for the viscosity in the accretion disks around the compact objects in these sources. We also derive the distances to these systems as functions of the masses of their compact objects.

We present evidence for the existence of two ~35 day clocks in the Her X-1/HZ Her binary system. ~35 day modulations are observed 1) in the Turn-On cycles with two on- and two off-states, and 2) in the changing shape of the pulse profiles which re-appears regularly. The two ways of counting the 35 day cycles are generally in synchronization. This synchronization did apparently break down temporarily during the long Anomalous Low (AL3) which Her X-1 experienced in 1999/2000, in the sense that there must have been one extra Turn-On cycle. Our working hypothesis is that there are two clocks in the system, both with a period of about ~35 days: precession of the accretion disk (the less stable "Turn-On clock") and free precession of the neutron star (the more stable "Pulse profile clock"). We suggest that free precession of the neutron star is the master clock, and that the precession of the accretion disk is basically synchronized to that of the neutron star through a feed-back mechanism in the binary system. However, the Turn-On clock can slip against its master when the accretion disk has a very low inclination, as is observed to be the case during AL3. We take the apparent correlation between the histories of the Turn-Ons, of the Anomalous Lows and of the pulse period evolution, with a 5 yr quasi-periodicity, as evidence for strong physical interaction and feed-back between the major components in the system. We speculate that the 5 yr (10 yr) period is either due to a corresponding activity cycle of HZ Her or a natural ringing period of the physical system of coupled components. The question whether free precession really exists in neutron stars is of great importance for the understanding of matter with supra-nuclear density.

The dynamics of small global perturbations in the form of linear combination of a finite number of non-axisymmetric eigenmodes is studied in two-dimensional approximation. The background flow is assumed to be an axisymmetric perfect fluid with the adiabatic index $\gamma=5/3$ rotating with power law angular velocity distribution $\Omega \propto r^{-q}$, $1.5<q<2.0$, confined by free boundaries in the radial direction. The substantial transient growth of acoustic energy of optimized perturbations is discovered. An optimal energy growth $G$ is calculated numerically for a variety of parameters. Its value depends essentially on the perturbation azimuthal wavenumber $m$ and increases for higher values of $m$. The closer the rotation profile to the Keplerian law, the larger growth factors can be obtained but over a longer time. The highest acoustic energy increase found numerically is of order $\sim 10^2$ over $\sim 6$ typical Keplerian periods. Slow neutral eigenmodes with corotation radius beyond the outer boundary mostly contribute to the transient growth. The revealed linear temporal behaviour of perturbations may play an important role in angular momentum transfer in toroidal flows near compact relativistic objects.

We address the task of modeling soft X-ray and optical light curves of X-ray novae in the high/soft state. The analytic model of viscous evolution of an externally truncated accretion \alpha-disk is used. Relativistic effects near a Kerr black hole and self-irradiation of an accretion disk are taken into account. The model is applied to the outbursts of X-ray nova Monocerotis 1975 (A 0620-00) and X-ray nova Muscae 1991 (GRS 1124-68). Comparison of observational data with the model yields constraints on the angular momentum (the Kerr parameter) of the black holes in A 0620-00 and GRS 1124-68: 0.3-0.6 and \leq 0.4, and on the viscosity parameter \alpha of the disks: 0.7-0.95 and 0.55-0.75. We also conclude that the accretion disks should have an effective geometrical thickness 1.5-2 times greater than the theoretical value of the distance between the photometric layers.

Aims: We investigate the X-ray spectral and timing properties of the accreting X-ray pulsar Her X-1 observed with the INTEGRAL satellite in July-August 2005. Methods: The data analyzed in this work cover a substantial part of one main-on state of the source. The short-time scale pulse period development is measured. X-ray pulse profiles for different energy ranges and time intervals are constructed. Pulse-averaged and pulse-phase resolved broad band X-ray spectra are studied. Spectral changes during X-ray dips are explored. Results: The X-ray pulse profiles are found to change significantly during the period of observations. For the first time a strong spinup is measured within one 35 d cycle. Spectral characteristics observed during the X-ray dips are consistent with their interpretaion as due to partial covering as has been reported by several authors. The fundamental cyclotron absorption line is firmly observed in both pulse-averaged and pulse-phase resolved X-ray spectra. The energy, width, and the depth of the line are found to vary significantly with pulse phase.

Analysis and interpretation of the Her X-1 X-ray light curve obtained with the ASM onbord RXTE over the period 1996 February to 2004 September is presented. We report that the features found previously in the averaged X-ray lightcurve are confirmed by the new RXTE/ASM data. In particular, anomalous dips and post-eclipse recoveries in two successive orbits in the short-on state are found to be prominent and stable details of the light curve. We argue that a change of the tilt of the accretion disk over the 35-day period is necessary to account for these features and show that our numerical model can explain such a behavior of the disk and reproduce the observed details of the light curve.

The instability of non-homoentropic axisymmetric flow of perfect fluid with respect to non-axisymmetric infinitesimal perturbations was investigated by numerical integration of hydrodynamical differential equations in two-dimensional approximation. The non-trivial influence of entropy gradient on unstable sound and surface gravity waves was revealed. In particular, both decrease and growth of entropy against the direction of effective gravitational acceleration $g_{eff}$ give rise to growing surface gravity modes which are stable with the same parameters in the case of homoentropic flow. At the same time increment of sound modes either grows monotonically while the rate of entropy decrease against $g_{eff}$ gets higher or vanishes at some values of positive and negative entropy gradient in the basic flow. The calculations have showed also that growing internal gravity modes appear only in the flow unstable to axisymmetric perturbations. At last, the analysis of boundary problem with free boundaries uncovered that's incorrect to set the entropy distribution according to polytropic law with polytropic index different from adiabatic value, since in this case perturbations don't satisfy the free boundary conditions.

First results of observations of the low mass X-ray binary Her X-1/HZ Her performed by the INTEGRAL satellite in July-August 2005 are presented. A significant part of one 35 day main-on state was covered. The cyclotron line in the X-ray spectrum is well observed and its position and shape, as well as its variability with time and phase of the 1.24 s pulsation are explored. X-ray pulse profiles for different energy bands are studied throughout the observation. The pulse period is found to vary on short time scales revealing a dynamical spin-up/spin-down behavior. Results of simultaneous optical observations of HZ Her are also discussed.

We have studied the long-term (1971-2005) behaviour of the 1.24 sec pulse period and the 35 day precession period of Her X-1 and show that both periods vary in a highly correlated way (see also Staubert et al. 1997 and 2000). When the spin-up rate decreases, the 35 day turn-on period shortens.This correlation is most evident on long time scales (~2000 days),e.g.around four extended spin-down episodes, but also on shorter time scales (a few 100 days) on which quasi-periodic variations are apparent. We argue that the likely common cause is variations of the mass accretion rate onto the neutron star.The data since 1991 allow a continuous sampling and indicate a lag between the turn-on behaviour and the spin behaviour, in the sense that changes are first seen in the spin, about one cycle later in the turn-on. Both the coronal wind model (Schandl & Meyer 1994) as well as the stream-disk model (Shakura et al.999) predict this kind of behaviour.

We present the results of ten years of repeated measurements of the Cyclotron Resonance Scattering Feature (CRSF) in the spectrum of the binary X-ray pulsar Her X-1 and report the discovery of a positive correlation of the centroid energy of this absorption feature in pulse phase averaged spectra with source luminosity.Our results are based on a uniform analysis of observations bythe RXTE satellite from 1996 to 2005, using sufficiently long observations of 12 individual 35-day Main-On states of the source. The mean centroid energy E_c of the CRSF in pulse phase averaged spectra of Her X-1 during this time is around 40 keV, with significant variations from one Main-On state to the next. We find that the centroid energy of the CRSF in Her X-1 changes by ~5% in energy for a factor of 2 in luminosity. The correlation is positive, contrary to what is observed in some high luminosity transient pulsars. Our finding is the first significant measurement of a positive correlation between E_c and luminosity in any X-ray pulsar. We suggest that this behaviour is expected in the case of sub-Eddington accretion and present a calculation of a quantitative estimate, which is very consistent with the effect observed in Her X-1.We urge that Her X-1 is regularly monitored further and that other X-ray pulsars are investigated for a similar behaviour.

We have performed an analysis and interpretation of the X-ray light curve of the accreting neutron star Her X-1 obtained with the ASM RXTE over the period 1996 February to 2004 September. The averaged X-ray light curves are constructed by means of adding up light curves corresponding to different 35 day cycles. A numerical model is introduced to explain the properties of the averaged light curves. We argue that a change of the tilt of the accretion disk over the 35 day period is necessary to account for the observed features and show that our numerical model can explain such a behavior of the disk and reproduce the details of the light curve.

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 present the analysis of Rossi X-ray Timing Explorer (RXTE) observations of the turn-on phase of a 35 day cycle of the X-ray binary Her X-1. During the early phases of the turn-on, the energy spectrum is composed of X-rays scattered into the line of sight plus heavily absorbed X-rays. The energy spectra in the 3-17 keV range can be described by a partial covering model, where one of the components is influenced by photoelectric absorption and Thomson scattering in cold material plus an iron emission line at 6.5 keV. In this paper we show the evolution of spectral parameters as well as the evolution of the pulse profile during the turn-on. We describe this evolution using Monte Carlo simulations which self-consistently describe the evolution of the X-ray pulse profile and of the energy spectrum.

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

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

Analytical investigation of time-dependent accretion in disks is carried out. We consider a time-dependent disk in a binary system at outburst which has a fixed tidally-truncated outer radius. The standard Shakura-Sunyaev model of the disk is considered. The vertical structure of the disk is accurately described in two regimes of opacity: Thomson and free-free. Fully analytical solutions are obtained, characterized by power-law variations of accretion rate with time. The solutions supply asymptotic description of disk evolution in flaring sources in the periods after outbursts while the disk is fully ionized. The X-ray flux of multicolor (black-body) alpha-disk is obtained as varying quasi-exponentially. Application to X-ray novae is briefly discussed concerning the observed faster-than-power decays of X-ray light curves. The case of time-dependent advective disk when the exponential variations of accretion rate can occur is discussed.

The X-ray binary pulsar Her X-1 shows a wide variety of long and short term variabilities in the X-ray light curve. The 35 d variability of the source is interpreted as the influence of a warped, inclined, and twisted accretion disk periodically covering the line of sight to the central neutron star. In 1997 September we observed the ``turn-on'' of a 35 d cycle with the Rossi X-ray Timing Explorer (RXTE). Spectral analysis reveals that during early phases of the turn-on the overall spectrum is composed of X-rays scattered into the line of sight plus heavily absorbed X-rays. This interpretation is consistent with the variation of the pulse profile observed at the same time. The overall shape of the pulse profile is not changing, but towards earlier phases of the turn-on the pulse signature is steadily ``washed out''. This behavior can be understood as an influence of scattering and absorption due to the presence of the accretion disk rim. Using a Monte Carlo code we simulate the influence of both processes on a time variable, beamed emission characteristic, similar to the pulse profile of Her X-1. By comparing the results of the simulation with the observed profiles we determine the amount of scattered radiation, absorbed radiation, and the size of the scattering region.

A careful analysis of X-ray light curves and pulse profiles of Her X-1 obtained over more than 20 years strongly evidences for free precession of a magnetized neutron star with rotational axis inclined to the orbital plane as a central clock underlying the observed 35-day period. Strong asymmetric X-ray illumination of the optical star atmosphere leads to the formation of gaseous streams coming out of the orbital plane and forming a tilted accretion disk around the neutron star. Such a disk precesses due to tidal forces and dynamical action of gaseous streams from the secondary companion. The locking of these torques with neutron star precession makes the net disk precession period to be very close to that of the neutron star free precession.

A failed main-on in the 35d cycle of Her X-1 was observed with the Rossi X-Ray Timing Explorer (RXTE) on 1999 April 26. Exceptions to the normal 35d cycle have been seen only twice before; in 1983 and again 1993. We present timing and spectral results of this latest Anomalous Low State (ALS) along with comparisons to the main-on and normal low states. Pulsations were observed in the 3-18 keV band with a fractional RMS variation of (0.037+-0.003). Spectral analysis indicates that the ALS spectrum has the same shape as the main-on but is modified by heavy absorption and scattering. We find that 70% of the observed emission has passed through a cold absorber (N_H=5.0x10^{23}cm^{-2}). This partially absorbing spectral fit can be applied to the normal low state with similar results. We find that the ALS observations may be interpreted as a decrease in inclination of the accretion disk causing the central X-Ray source to be obscured over the entire 35d cycle.

A strong X-ray illumination of the optical star atmosphere in Her X-1, asymmetric because of a partial shadowing by the tilted twisted accretion disk around central neutron star, leads to the formation of matter flows coming out of the orbital plane and crossing the line of sight before entering the disk. We suggest that the absorption of X-ray emission by this flow leads to the formation of pre-eclipse and anomalous dips of type I. These dips are observed during several orbits after turn-on both in the main-on and short-on state. Almost coherent action of tidal torques and matter streams enhances the disk wobbling which causes the disk edge to shield the X-ray source after the turn-on. Anomalous dips of type II and post-eclipse recovery appear due to this process only on the first orbit after turn-on.

An analysis of the publically available RXTE archive on Her X-1 including data on 23 34.85-day cycles is performed. The turn-on times for these cycles are determined. The number of cycles with a duration of 20.5 orbits has been found exceedingly larger than of shorter (20 orbits) or longer (21 orbits). A correlation between the duration of a cycle and its mean X-ray flux is noted. The mean X-ray light curve shows a very distinct short-on state. The anomalous X-ray absorption dip is found during the first orbit after the turn-on in the main-on state for the cycles starting near the binary phase 0.25, while is present during two successive orbits in the low-on state. The post-eclipse recovery feature have not been found in the main-on state but appears at least for two orbits during the low-on state. The pre-eclipse dips are present both in main-on and low-on state and demonstrate the behaviour as in early observations. The comparison of durations of the main-on and short-on states enabled us to constrain the accretion disk semi-thickness and its inclination to the orbital plane.

We show that the free precession of a triaxial body can naturally explain the anomalously rapid change of the X-ray pulse profile of Her X-1 observed by the HEAO-1 in September 1978 without requiring a large change in the moment of inertia.