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

We report on \textit{NuSTAR} observations of the well-known wind-accreting X-ray pulsar \source\ during a strong spin-up episode that took place in January-March 2019. A high luminosity of the source in a most recent observation allowed us to detect a positive correlation of the cyclotron line energy with luminosity. Beyond that, only minor differences in spectral and temporal properties of the source during the spin-up, presumably associated with the formation of a transient accretion disk, and the normal wind-fed state could be detected. We finally discuss conditions for the formation of the disk and possible reasons for lack of any appreciable variations in most of the observed source properties induced by the change of the accretion mechanism, and conclude that the bulk of the observed X-ray emission is still likely powered by direct accretion from the wind.

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.

We present results of investigation of the poorly studied X-ray pulsar Swift J1816.7--1613 during its transition from the type I outburst to the quiescent state. Our studies are based on the data obtained from X-ray observatories \textit{Swift}, \textit{NuSTAR} and \textit{Chandra} alongside with the latest IR data from UKIDSS/GPS and \textit{Spitzer}/GLIMPSE surveys. The aim of the work is to determine parameters of the system: the strength of the neutron star magnetic field and the distance to the source, which are required for the interpretation of the source behaviour in the framework of physically motivated models. No cyclotron absorption line was detected in the broad-band energy spectrum. However, the timing analysis hints at the typical for the X-ray pulsars magnetic field from a few $\times 10^{11}$ to a few $\times 10^{12}$ G. We also estimated type of the IR-companion as a B0-2e star located at distance of 7--13~kpc.