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.

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.

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.

We aim both to complement the existing data on the spin history of the peculiar accreting X-ray pulsar GX 1+4 with more past and current data from BeppoSAX, INTEGRAL, and Fermi and to interpret the evolution in the framework of accretion theory. We used source light curves obtained from BeppoSAX/WFC and INTEGRAL/ISGRI to derive pulse periods using an epoch-folding analysis. Fermi/GBM data were analyzed by fitting a constant plus a Fourier expansion to background-subtracted rates, and maximizing the Y2 statistic. We completed the sample with hard X-ray light curves from Swift/BAT. The data were checked for correlations between flux and changes of the pulsar spin on different timescales. The spin-down of the pulsar continues with a constant change in frequency, i.e., an apparently accelerating change in the period. Over the past three decades, the pulse period has increased by about ~50%. Short-term fluctuations on top of this long-term trend do show anti-correlation with the source flux. Possible explanations of the observed long-term frequency and its dependence on flux are discussed.

We report on the long-term evolution of the spin period of the symbiotic X-ray pulsar GX 1+4 and a possible interpretation within a model of quasi-spherical accretion. New period measurements from BeppoSAX, INTEGRAL and Fermi observations have been combined with previously published data from four decades of observations. During the 1970s GX 1+4 was spinning up with the fastest rate among the known X-ray pulsars at the time. In the mid 1980s it underwent a change during a period of low X-ray ux and started to spin down with a rate similar in magnitude to the previous spin up rate. The spin period has changed from ~110 s to ~160 s within the last three decades. Our results demonstrate that the overall spin down trend continues and is stronger than ever. We compare the observations with predictions from a model assuming quasi-spherical accretion from the slow wind of the M giant companion.

We report on the long-term evolution of the spin period of the symbiotic X-ray pulsar GX 1+4 and a possible interpretation within a model of quasi-spherical accretion. New period measurements from BeppoSAX/WFC, INTEGRAL/ISGRI and Fermi/GBM observations have been combined with previously published data from four decades of observations. During the 1970's GX 1+4 was spinning up with the fastest rate among the known X-ray pulsars at the time. In the mid 1980's it underwent a change during a period of low X-ray flux and started to spin down with a rate similar in magnitude to the previous spin up rate. The spin period has changed from ~110 s to ~160 s within the last three decades. Our results demonstrate that the overall spin down trend continues and is stronger than ever. We compare the observations with predictions from a model assuming quasi-spherical accretion from the slow wind of the M giant companion.

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.

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.