It has been recently discovered that the transition of X-ray pulsars to the low luminosity state ($L\lesssim 10^{35}\,{\rm erg\,s^{-1}}$) is accompanied by a dramatic spectral changes. Namely, the typical power-law-like spectrum with high energy cutoff transforms into a two-component structure with a possible cyclotron absorption feature on top of it. It was proposed that these spectral characteristics can be explained qualitatively by the emission of cyclotron photons in the atmosphere of the neutron star caused by collisional excitation of electrons to upper Landau levels and further comptonization of the photons by electron gas. The latter is expected to be overheated in a thin top layer of the atmosphere. In this paper, we perform Monte Carlo simulations of the radiative transfer in the atmosphere of an accreting neutron star while accounting for a resonant scattering of polarized X-ray photons by thermally distributed electrons. The spectral shape is shown to be strongly polarization-dependent in soft X-rays ($\lesssim 10\,{\rm keV}$) and near the cyclotron scattering feature. The results of our numerical simulations are tested against the observational data of X-ray pulsar A 0535+262 in the low luminosity state. We show that the spectral shape of the pulsar can be reproduced by the proposed theoretical model. The applications of the discovery to the observational studies of accreting neutron stars are discussed.

We present results of the first dedicated observation of the young X-ray pulsar SXP 1062 in the broad X-ray energy band obtained during its 2019 outburst with the NuSTAR and XMM-Newton observatories. The analysis of the pulse-phase averaged and phase-resolved spectra in the energy band from 0.5 to 70 keV did not reveal any evidence for the presence of a cyclotron line. The spin period of the pulsar was found to have decreased to 979.48+/-0.06 s implying a ~10% reduction compared to the last measured period during the monitoring campaign conducted about five years ago, and is puzzling considering that the system apparently did not show major outbursts ever since. The switch of the pulsar to the spin-up regime supports the common assumption that torques acting on the accreting neutron star are nearly balanced and thus SXP 1062 likely also spins with a period close to the equilibrium value for this system. The current monitoring of the source revealed also a sharp drop of its soft X-ray flux right after the outburst, which is in drastic contrast to the behaviour during the previous outburst when the pulsar remained observable for years with only a minor flux decrease after the end of the outburst. This unexpected off state of the source lasted for at most 20 days after which SXP 1062 returned to the level observed during previous campaigns. We discuss this and other findings in context of the modern models of accretion onto strongly magnetized neutron stars.

We report on analysis of observations of the bright transient X-ray pulsar \src obtained during its 2017-2018 giant outburst with Insight-HXMT, \emph{NuSTAR}, and \textit{Swift} observatories. We focus on the discovery of a sharp state transition of the timing and spectral properties of the source at super-Eddington accretion rates, which we associate with the transition of the accretion disk to a radiation pressure dominated (RPD) state, the first ever directly observed for magnetized neutron star. This transition occurs at slightly higher luminosity compared to already reported transition of the source from sub- to super-critical accretion regime associate with onset of an accretion column. We argue that this scenario can only be realized for comparatively weakly magnetized neutron star, not dissimilar to other ultra-luminous X-ray pulsars (ULPs), which accrete at similar rates. Further evidence for this conclusion is provided by the non-detection of the transition to the propeller state in quiescence which strongly implies compact magnetosphere and thus rules out magnetar-like fields.

We report the discovery of absorption features in the X-ray spectrum of the transient X-ray pulsar GROJ2058+42. The features are detected around $\sim10$, $\sim20$ and $\sim30$ keV in both NuSTAR observations carried out during the source type II outburst in spring 2019. The most intriguing property is that the deficit of photons around these energies is registered only in the narrow phase interval covering around 10% of the pulsar spin period. We interpret these absorption lines as a cyclotron resonant scattering line (fundamental) and two higher harmonics. The measured energy allow us to estimate the magnetic field strength of the neutron star as $\sim10^{12}$ G.

The X-ray Polarization Probe (XPP) is a second generation X-ray polarimeter following up on the Imaging X-ray Polarimetry Explorer (IXPE). The XPP will offer true broadband polarimetery over the wide 0.2-60 keV bandpass in addition to imaging polarimetry from 2-8 keV. The extended energy bandpass and improvements in sensitivity will enable the simultaneous measurement of the polarization of several emission components. These measurements will give qualitatively new information about how compact objects work, and will probe fundamental physics, i.e. strong-field quantum electrodynamics and strong gravity.

Deep NuSTAR observation of X-ray pulsar A 0535+262, performed at a very low luminosity of $\sim7\times10^{34}$ erg s$^{-1}$, revealed the presence of two spectral components. We argue that the high-energy component is associated with cyclotron emission from recombination of electrons collisionally excited to the upper Landau levels. The cyclotron line energy of $E_{\rm cyc}=47.7\pm0.8$ keV was measured at the luminosity of almost an order of magnitude lower than what was achieved before. The data firmly exclude a positive correlation of the cyclotron energy with the mass accretion rate in this source.

The luminosity of X-ray pulsars powered by accretion onto magnetized neutron stars covers a wide range over a few orders of magnitude. The brightest X-ray pulsars recently discovered as pulsating ultraluminous X-ray sources reach accretion luminosity above $10^{40}\,{\rm erg\,s^{-1}}$ which exceeds the Eddington value more than by a factor of ten. Most of the energy is released within small regions in the vicinity of magnetic poles of accreting neutron star - in accretion columns. Because of the extreme energy release within a small volume accretion columns of bright X-ray pulsars are ones of the hottest places in the Universe, where the internal temperature can exceed 100 keV. Under these conditions, the processes of creation and annihilation of electron-positron pairs can be influential but have been largely neglected in theoretical models of accretion columns. In this letter, we investigate properties of a gas of electron-positron pairs under physical conditions typical for accretion columns. We argue that the process of pairs creation can crucially influence both the dynamics of the accretion process and internal structure of accretion column limiting its internal temperature, dropping the local Eddington flux and increasing the gas pressure.

The X-ray pulsar GRO J1744-28 is a unique source which shows both pulsations and type-II X-ray bursts, allowing studies of the interaction of the accretion disk with the magnetosphere at huge mass accretion rates exceeding $10^{19}$ g s$^{-1}$ during its super-Eddington outbursts. The magnetic field strength in the source, $B\approx 5\times 10^{11}$ G, is known from the cyclotron absorption feature discovered in the energy spectrum around 4.5 keV. Here, we explore the flux variability of the source in context of interaction of its magnetosphere with the radiation-pressure dominated accretion disk. Particularly, we present the results of the analysis of noise power density spectra (PDS) using the observations of the source in 1996-1997 by RXTE. Accreting compact objects commonly exhibit a broken power-law shape of the PDS with a break corresponding to the Keplerian orbital frequency of matter at the innermost disk radius. The observed frequency of the break can thus be used to estimate the size of the magnetosphere. We found, however, that the observed PDS of GRO J1744-28 differs dramatically from the canonical shape. Furthermore, the observed break frequency appears to be significantly higher than what is expected based on the magnetic field estimated from the cyclotron line energy. We argue that these observational facts can be attributed to the existence of the radiation-pressure dominated region in the accretion disk at luminosities above $\sim$2$\times 10^{37}$ erg s$^{-1}$. We discuss a qualitative model for the PDS formation in such disks, and show that its predictions are consistent with our observational findings. The presence of the radiation-pressure dominated region can also explain the observed weak luminosity-dependence of the inner radius, and we argue that the small inner radius can be explained by a quadrupole component dominating the magnetic field of the neutron star.

We investigate aperiodic X-ray flux variability in accreting highly magnetized neutron stars - X-ray pulsars (XRPs). The X-ray variability is largely determined by mass accretion rate fluctuations at the NS surface, which replicate accretion rate fluctuations at the inner radius of the accretion disc. The variability at the inner radius is due to fluctuations arising all over the disc and propagating inwards under the influence of viscous diffusion. The inner radius varies with mean mass accretion rate and can be estimated from the known magnetic field strength and accretion luminosity of XRPs. Observations of transient XRPs covering several orders of magnitude in luminosity give a unique opportunity to study effects arising due to the changes of the inner disc radius. We investigate the process of viscous diffusion in XRP accretion discs and construct new analytical solutions of the diffusion equation applicable for thin accretion discs truncated both from inside and outside. Our solutions are the most general ones derived in the approximation of Newtonian mechanics. We argue that the break observed at high frequencies in the power density spectra of XRPs corresponds to the minimal time scale of the dynamo process, which is responsible for the initial fluctuations. Comparing data from the bright X-ray transient A 0535+26 with our model, we conclude that the time scale of initial variability in the accretion disc is a few times longer than the local Keplerian time scale.

It has recently been discovered that a fraction of ultra-luminous X-ray sources (ULXs) exhibit X-ray pulsations, and are therefore powered by super-Eddington accretion onto magnetized neutron stars (NSs). For typical ULX mass accretion rates ($\gtrsim 10^{19}\,{\rm g\,s^{-1}}$), the inner parts of the accretion disc are expected to be in the supercritical regime, meaning that some material is lost in a wind launched from the disc surface, while the rest forms an optically thick envelope around the NS as it follows magnetic field lines from the inner disc radius to the magnetic poles of the star. The envelope hides the central object from a distant observer and defines key observational properties of ULX pulsars: their energy spectrum, polarization, and timing features. The optical thickness of the envelope is affected by the mass losses from the disc. We calculate the mass loss rate due to the wind in ULX pulsars, accounting for the NS magnetic field strength and advection processes in the disc. We argue that the detection of strong outflows from ULX pulsars can be considered evidence of a relatively weak dipole component of the NS magnetic field. We estimate the influence of mass losses on the optical thickness of the envelope and analyze how the envelope affects broadband aperiodic variability in ULXs. We show that brightness fluctuations at high Fourier frequencies can be strongly suppressed by multiple scatterings in the envelope and that the strength of suppression is determined by the mass accretion rate and geometrical size of the magnetosphere.

We report on the discovery of a dramatic change in the energy spectrum of the X-ray pulsar GX 304-1 appearing at low luminosity. Particularly, we found that the cutoff power-law spectrum typical for accreting pulsars, including GX 304-1 at higher luminosities of $L_{\rm X}\sim 10^{36} - 10^{37}$ erg s$^{-1}$, transformed at lower luminosity of $L_{\rm X}\sim 10^{34}$ erg s$^{-1}$ to a two-component spectrum peaking around 5 and 40 keV. We suggest that the observed transition corresponds to a change of the dominant mechanism responsible for the deceleration of the accretion flow. We argue that the accretion flow energy at low accretion rates is released in the atmosphere of the neutron star, and the low-energy component in the source spectrum corresponds to the thermal emission of the optically thick, heated atmospheric layers. The most plausible explanations for the high-energy component are either the cyclotron emission reprocessed by the magnetic Compton scattering or the thermal radiation of deep atmospheric layers partly Comptonized in the overheated upper layers. Alternative scenarios are also discussed.

In the work we present the results of two deep broad-band observations of the poorly studied X-ray pulsar IGR J19294+1816 obtained with the NuSTAR observatory. The source was observed during Type I outburst and in the quiescent state. In the bright state a cyclotron absorption line in the energy spectrum was discovered at $E_{\rm cyc}=42.8\pm0.7$ keV. Spectral and timing analysis prove the ongoing accretion also during the quiescent state of the source. Based on the long-term flux evolution, particularly on the transition of the source to the bright quiescent state with luminosity around $10^{35}$ erg s$^{-1}$, we concluded that IGR J19294+1816 switched to the accretion from the "cold" accretion disc between Type I outbursts. We also report the updated orbital period of the system.

We report on the first X-ray observation of the 0.28 s isolated radio pulsar PSR J1154--6250 obtained with the XMM-Newton observatory in February 2018. A point-like source is firmly detected at a position consistent with that of PSR J1154--6250. The two closest stars are outside the 3$\sigma$ confidence limits of the source position and thus unlikely to be responsible for the observed X-ray emission. The energy spectrum of the source can be fitted equally well either with an absorbed power-law with a steep photon index $\Gamma\approx 3.3$ or with an absorbed blackbody with temperature $kT=0.21\pm 0.04$~keV and emitting radius $R_\mathrm{BB} \approx 80$ m (assuming a distance of 1.36~kpc). The X-ray luminosity of $4.4\times 10^{30}$ erg s$^{-1}$ derived with the power-law fit corresponds to an efficiency of $\eta_X = L^\mathrm{unabs}_X/\dot E = 4.5\times 10^{-3}$, similar to those of other old pulsars. The X-ray properties of PSR J1154--6250 are consistent with an old age and suggest that the spatial coincidence of this pulsar with the OB association Cru OB1 is due to a chance alignment.

We report on the monitoring of the final stage of the outburst from the first Galactic ultraluminous X-ray pulsar Swift J0243.6+6124, which reached $\sim$40 Eddington luminosities at the peak of the outburst. The main aim of the monitoring program with the {\it Swift}/XRT telescope was to measure the magnetic field of the neutron star using the luminosity of transition to the "propeller" state. The visibility constraints, unfortunately, did not permit us to observe the source down to the fluxes low enough to detect such a transition. The tight upper limit on the propeller luminosity $L_{\rm prop}<6.8\times10^{35}$ erg s$^{-1}$ implies the dipole component of the magnetic field $B<10^{13}$ G, under plausible assumptions on the effective magnetosphere size. On the other hand, the observed evolution of the pulse profile and of the pulsed fraction with flux points to a change of the emission region geometry at the critical luminosity $L_{\rm crit}\sim3\times10^{38}$ erg s$^{-1}$ both in the rising and declining parts of the outburst. We associate the observed change with the onset of the accretion column, which allows us to get an independent estimate of the magnetic field strength close to the neutron stars surface of $B>10^{13}$ G. Given the existing uncertainty in the effective magnetosphere size, we conclude that both estimates are marginally compatible with each other.

The classical limit on the accretion luminosity of a neutron star is given by the Eddington luminosity. The advanced models of accretion onto magnetized neutron stars account for the appearance of magnetically confined accretion columns and allow the accretion luminosity to be higher than the Eddington value by a factor of tens. However, the recent discovery of pulsations from ultraluminous X-ray source (ULX) in NGC 5907 demonstrates that the accretion luminosity can exceed the Eddington value up to by a factor of 500. We propose a model explaining observational properties of ULX-1 in NGC 5907 without any ad hoc assumptions. We show that the accretion column at extreme luminosity becomes advective. Enormous energy release within a small geometrical volume and advection result in very high temperatures at the bottom of accretion column, which demand to account for the energy losses due to neutrino emission which can be even more effective than the radiation energy losses. We show that the total luminosity at the mass accretion rates above $10^{21}\,{\rm g\,s^{-1}}$ is dominated by the neutrino emission similarly to the case of core-collapse supernovae. We argue that the accretion rate measurements based on detected photon luminosity in case of bright ULXs powered by neutron stars can be largely underestimated due to intense neutrino emission. The recently discovered pulsating ULX-1 in galaxy NGC 5907 with photon luminosity of $\sim 10^{41}\,{\rm erg\,s^{-1}}$ is expected to be even brighter in neutrinos and is thus the first known Neutrino Pulsar.

The luminosity of accreting magnetised neutron stars can largely exceed the Eddington value due to appearance of accretion columns. The height of the columns can be comparable to the neutron star radius. The columns produce the X-rays detected by the observer directly and illuminate the stellar surface, which reprocesses the X-rays and causes additional component of the observed flux. The geometry of the column and the illuminated part of the surface determines the radiation beaming. Curved space-time affects the angular flux distribution. We construct a simple model of the beam patterns formed by direct and reflected flux from the column. We take into account the possibility of appearance of accretion columns, whose height is comparable to the neutron star radius. We argue that depending on the compactness of the star the flux from the column can be either strongly amplified due to gravitational lensing, or significantly reduced due to column eclipse by the star. The eclipses of high accretion columns result in specific features in pulse profiles. Their detection can put constraints on the neutron star radius. We speculate that column eclipses are observed in X-ray pulsar V 0332+53, leading us to the conclusion of large neutron star radius in this system ($\sim15\,{\rm km}$ if $M\sim 1.4M_\odot$). We point out that the beam pattern can be strongly affected by scattering in the accretion channel at high luminosity, which has to be taken into account in the models reproducing the pulse profiles.

Many statistical properties of X-ray aperiodic variability from accreting compact objects can be explained by the propagating fluctuations model applied to the accretion disc. The mass accretion rate fluctuations originate from variability of viscosity, which arises at every radius and causes local fluctuations of the density. The fluctuations diffuse through the disc and result in local variability of the mass accretion rate, which modulates the X-ray flux from the inner disc in the case of black holes, or from the surface in the case of neutron stars. A key role in the theoretical explanation of fast variability belongs to the description of the diffusion process. The propagation and evolution of the fluctuations is described by the diffusion equation, which can be solved by the method of Green functions. We implement Green functions in order to accurately describe the propagation of fluctuations in the disc. For the first time we consider both forward and backward propagation. We show that (i) viscous diffusion efficiently suppress variability at time scales shorter than the viscous time, (ii) local fluctuations of viscosity affect the mass accretion rate variability both in the inner and the outer parts of accretion disc, (iii) propagating fluctuations give rise not only to hard time lags as previously shown, but also produce soft lags at high frequency similar to those routinely attributed to reprocessing, (iv) deviation from the linear rms-flux relation is predicted for the case of very large initial perturbations. Our model naturally predicts bumpy power spectra.

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.

Magnetic field of accreting neutron stars determines their overall behaviour including the maximum possible luminosity. Some models require an above-average magnetic field strength (> 10^13 G) in order to explain super-Eddington mass accretion rate in the recently discovered class of pulsating ultraluminous X-ray sources (ULX). The peak luminosity of SMC X-3 during its major outburst in 2016-2017 reached ~2.5x10^39 erg/s comparable to that in ULXs thus making this source the nearest ULX-pulsar. SMC X-3 belongs to the class of transient X-ray pulsars with Be optical companions, and exhibited a giant outburst in July 2016 - February 2017. The source has been observed during the entire outburst with the Swift/XRT and Fermi/GBM telescopes, as well as the NuSTAR observatory. Collected data allowed us to estimate the magnetic field strength of the neutron star in SMC X-3 using several independent methods. Spin evolution of the source during and between the outbursts and the luminosity of the transition to so-called propeller regime in the range of (0.3 - 7)x10^35 erg/s imply relatively weak dipole field of (1 - 5)x10^12 G. On the other hand, there is also evidence for much stronger field in the immediate vicinity of the neutron star surface. In particular, transition from super- to sub-critical accretion regime associated with cease of the accretion column, absence of cyclotron absorption features in the broadband X-ray spectrum of the source obtained with NuSTAR and very high peak luminosity favor an order of magnitude stronger field. This discrepancy makes SMC X-3 a good candidate to posses significant non-dipolar components of the field, and an intermediate source between classical X-ray pulsars and accreting magnetars which may constitute an appreciable fraction of ULX population.

We report on the analysis of NuSTAR observations of the Be-transient X-ray pulsar V 0332+53 during the giant outburst in 2015 and another minor outburst in 2016. We confirm the cyclotron-line energy-luminosity correlation previously reported in the source and the line energy decrease during the giant outburst. Based on 2016 observations, we find that a year later the line energy has increased again essentially reaching the pre-outburst values. We discuss this behaviour and conclude that it is likely caused by a change of the emission region geometry rather than previously suggested accretion-induced decay of the neutron stars magnetic field. At lower luminosities, we find for the first time a hint of departure from the anticorrelation of line energy with flux, which we interpret as a transition from super- to sub-critical accretion associated with the disappearance of the accretion column. Finally, we confirm and briefly discuss the orbital modulation observed in the outburst light curve of the source.

Magnetized neutron stars power at least some ultra-luminous X-ray sources. The accretion flow in these cases is interrupted at the magnetospheric radius and then reaches the surface of a neutron star following magnetic field lines. Accreting matter moving along magnetic field lines forms the accretion envelope around the central object. We show that, in case of high mass accretion rates $\gtrsim 10^{19}\,{\rm g\,s^{-1}}$ the envelope becomes closed and optically thick, which influences the dynamics of the accretion flow and the observational manifestation of the neutron star hidden behind the envelope. Particularly, the optically thick accretion envelope results in a multi-color black-body spectrum originating from the magnetospheric surface. The spectrum and photon energy flux vary with the viewing angle, which gives rise to pulsations characterized by high pulsed fraction and typically smooth pulse profiles. The reprocessing of radiation due to interaction with the envelope leads to the disappearance of cyclotron scattering features from the spectrum. We speculate that the super-orbital variability of ultra-luminous X-ray sources powered by accreting neutron stars can be attributed to precession of the neutron star due to interaction of magnetic dipole with the accretion disc.

We present the results of the monitoring programmes performed with the Swift/XRT telescope and aimed specifically to detect an abrupt decrease of the observed flux associated with a transition to the propeller regime in two well known X-ray pulsars 4U 0115+63 and V 0332+53 during their giant outbursts in 2015. Such transitions were detected at the threshold luminosities of $(1.4\pm0.4)\times10^{36}$ erg s$^{-1}$ and $(2.0\pm0.4)\times10^{36}$ erg s$^{-1}$ for 4U 0115+63 and V 0332+53, respectively. Spectra of the sources are shown to be significantly softer during the low state. In both sources, the accretion at rates close to the aforementioned threshold values briefly resumes during the periastron passage following the transition into propeller regime. The strength of the dipole component of the magnetic field required to inhibit the accretion agrees well with estimates based on the position of the cyclotron lines in their spectra, thus excluding presence of a strong multipole component of the magnetic field in the vicinity of the neutron star.

Compton scattering of polarized radiation in a strong magnetic field is considered. The recipe for calculation of the scattering matrix elements, the differential and total cross sections based on quantum electrodynamic (QED) second order perturbation theory is presented for the case of arbitrary initial and final Landau level, electron momentum along the field and photon momentum. Photon polarization and electron spin state are taken into account. The correct dependence of natural Landau level width on the electron spin state is taken into account in general case of arbitrary initial photon momentum for the first time. A number of steps in calculations were simplified analytically making the presented recipe easy-to-use. The redistribution functions over the photon energy, momentum and polarization states are presented and discussed. The paper generalizes already known results and offers a basis for accurate calculation of radiation transfer in strong $B$-field, for example, in strongly magnetized neutron stars.

We present here the first convincing observational manifestation of a magnetar-like magnetic field in an accreting neutron star in binary system - the first pulsating ultra-luminous X-ray source X-2 in the galaxy M82. Using the Chandra X-ray observatory data we show that the source exhibit the bimodal distribution of the luminosity with two well-defined peaks separated by a factor of 40. This behaviour can be interpreted as the action of the "propeller regime" of accretion. The onset of the propeller in a 1.37 s pulsar at luminosity of ~$10^{40}$ erg/s implies the dipole component of the neutron star magnetic field of ~$10^{14}$ G.

We study properties of luminous X-ray pulsars using a simplified model of the accretion column. The maximal possible luminosity is calculated as a function of the neutron star (NS) magnetic field and spin period. It is shown that the luminosity can reach values of the order of $10^{40}\,{\rm erg/s}$ for the magnetar-like magnetic field ($B\gtrsim 10^{14}\,{\rm G}$) and long spin periods ($P\gtrsim 1.5\,{\rm s}$). The relative narrowness of an area of feasible NS parameters which are able to provide higher luminosities leads to the conclusion that $L\simeq 10^{40}\,\,{\rm erg/s}$ is a good estimate for the limiting accretion luminosity of a NS. Because this luminosity coincides with the cut-off observed in the high mass X-ray binaries luminosity function which otherwise does not show any features at lower luminosities, we can conclude that a substantial part of ultra-luminous X-ray sources are accreting neutron stars in binary systems.

Cyclotron resonance scattering features observed in the spectra of some X-ray pulsars show significant changes of the line centroid energy with the pulsar luminosity. Whereas for bright sources above the so called critical luminosity these variations are established to be connected with the appearance of the high accretion column above the neutron star surface, at low, sub-critical luminosities the nature of the variations (but with the opposite sign) has not been discussed widely. We argue here that the cyclotron line is formed when the radiation from a hotspot propagates through the plasma falling with a mildly relativistic velocity onto the neutron star surface. The position of the cyclotron resonance is determined by the Doppler effect. The change of the cyclotron line position in the spectrum with luminosity is caused by variations of the velocity profile in the line-forming region affected by the radiation pressure force. The presented model has several characteristic features: (i) the line centroid energy is positively correlated with the luminosity; (ii) the line width is positively correlated with the luminosity as well; (iii) the position and the width of the cyclotron absorption line are variable over the pulse phase; (iv) the line has a more complicated shape than widely used Lorentzian or Gaussian profiles; (v) the phase-resolved cyclotron line centroid energy and the width are negatively and positively correlated with the pulse intensity, respectively. The predictions of the proposed theory are compared with the variations of the cyclotron line parameters in the X-ray pulsar GX 304-1 over a wide range of sub-critical luminosities as seen by the INTEGRAL observatory.

We present the results of the pulse phase- and luminosity-resolved spectroscopy of the transient X-ray pulsar V0332+53, performed for the first time in a wide luminosity range (1-40)x10^{37} erg/s during a giant outburst observed by the RXTE observatory in Dec 2004 - Feb 2005. We characterize the spectra quantitatively and built the detailed "three-dimensional" picture of spectral variations with pulse phase and throughout the outburst. We show that all spectral parameters are strongly variable with the pulse phase, and the pattern of this variability significantly changes with luminosity directly reflecting the associated changes in the structure of emission regions and their beam patterns. Obtained results are qualitatively discussed in terms of the recently developed reflection model for the formation of cyclotron lines in the spectra of X-ray pulsars.

The accretion flow around X-ray pulsars with a strong magnetic field is funnelled by the field to relatively small regions close to the magnetic poles of the neutron star (NS), the hotspots. During strong outbursts regularly observed from some X-ray pulsars, the X-ray luminosity can be so high, that the emerging radiation is able to stop the accreting matter above the surface via radiation-dominated shock, and the accretion column begins to rise. This border luminosity is usually called the "critical luminosity". Here we calculate the critical luminosity as a function of the NS magnetic field strength $B$ using exact Compton scattering cross section in strong magnetic field. Influence of the resonant scattering and photon polarization is taken into account for the first time. We show that the critical luminosity is not a monotonic function of the B-field. It reaches a minimum of a few 10^{36} erg s^{-1} when the cyclotron energy is about 10 keV and a considerable amount of photons from a hotspot have energy close to the cyclotron resonance. For small B, this luminosity is about 10^{37} erg s^{-1}, nearly independent of the parameters. It grows for the B-field in excess of 10^{12} G because of the drop in the effective cross-section of interaction below the cyclotron energy. We investigate how different types of the accretion flow and geometries of the accretion channel affect the results and demonstrate that the general behaviour of the critical luminosity on B-field is very robust. The obtained results are shown to be in a good agreement with the available observational data and provide a necessary ground for the interpretation of upcoming high quality data from the currently operating and planned X-ray telescopes.

Cyclotron resonance scattering features observed in the spectra of some X-ray pulsars show significant changes of the line energy with the pulsar luminosity. At high luminosities, these variations are often associated with the onset and growth of the accretion column, which is believed to be the origin of the observed emission and of the cyclotron lines. However, this scenario inevitably implies large gradient of the magnetic field strength within the line-forming region, which makes the formation of the observed line-like features problematic. Moreover, the observed variation of the cyclotron line energy is much smaller than could be anticipated for the corresponding luminosity changes. We argue here that a more physically realistic situation is that the cyclotron line forms when the radiation emitted by the accretion column is reflected from the neutron star surface, where the gradient of the magnetic field strength is significantly smaller. We develop here the reflection model and apply it to explain the observed variations of the cyclotron line energy in a bright X-ray pulsar V 0332+53 over a wide range of luminosities.

Cyclotron resonance scattering features are observed in the spectra of some X-ray pulsars and show significant changes in the line energy with the pulsar luminosity. In a case of bright sources, the line centroid energy is anti-correlated with the luminosity. Such a behaviour is often associated with the onset and growth of the accretion column, which is believed to be the origin of the observed emission and the cyclotron lines. However, this scenario inevitably implies large gradient of the magnetic field strength within the line-forming region, and it makes the formation of the observed line-like features problematic. Moreover, the observed variation of the cyclotron line energy is much smaller than could be anticipated for the corresponding luminosity changes. We argue that a more physically realistic situation is that the cyclotron line forms when the radiation emitted by the accretion column is reflected from the neutron star surface. The idea is based on the facts that a substantial part of column luminosity is intercepted by the neutron star surface and the reflected radiation should contain absorption features. The reflection model is developed and applied to explain the observed variations of the cyclotron line energy in a bright X-ray pulsar V 0332+53 over a wide range of luminosities.

We derive the relativistic kinetic equation for Compton scattering of polarized radiation in strong magnetic field using the Bogolyubov method. The induced scattering and the Pauli exclusion principle are taken into account. The electron polarization is also considered in the general form of the kinetic equation. The special forms of the equation for the cases of the non-polarized electrons, the rarefied electron gas and the two polarization mode description of radiation are found. The derived equations are valid for any photon and electron energies and the magnetic field strength below about 10^{16} G. These equations provide the basis for formulation of the equation for polarized radiation transport in atmospheres and magnetospheres of strongly magnetized neutron stars.