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 long-term multi-wavelength study of optical, UV and X-ray variability of the nearby changing-look (CL) Seyfert NGC 1566 observed with the Swift Observatory and the MASTER Global Robotic Network from 2007 to 2019. We started spectral observations with South African Astronomical Observatory 1.9-m telescope soon after the brightening was discovered in July 2018 and present here the data for the interval between Aug. 2018 to Sep. 2019. This paper concentrates on the remarkable post-maximum behaviour after July 2018 when all bands decreased with some fluctuations. We observed three significant re-brightenings in the post-maximum period during 17 Nov. 2018 - 10 Jan. 2019, 29 Apr. - 19 Jun. 2019 and 27 Jul.- 6 Aug. 2019. An X-ray flux minimum occurred in Mar. 2019. The UV minimum occurred about 3 months later. It was accompanied by a decrease of the Luv/Lx ratio. New post-maximum spectra covering (31 Nov. 2018 - 23 Sep. 2019) show dramatic changes compared to 2 Aug. 2018, with fading of the broad lines and [Fe X]6374 until Mar. 2019. These lines became somewhat brighter in Aug.-Sep. 2019. Effectively, two CL states were observed for this object: changing to type 1.2 and then returning to the low state as a type 1.8 Sy. We suggest that the changes are due mostly to fluctuations in the energy generation. The estimated Eddington ratios are about 0.055% for minimum in 2014 and 2.8% for maximum in 2018.

We present the discovery with Keck of the extremely infrared (IR) luminous transient AT 2017gbl, coincident with the Northern nucleus of the luminous infrared galaxy (LIRG) IRAS 23436+5257. Our extensive multi-wavelength follow-up spans ~900 days, including photometry and spectroscopy in the optical and IR, and (very long baseline interferometry) radio and X-ray observations. Radiative transfer modelling of the host galaxy spectral energy distribution shows the presence of a hitherto undetected dust obscured active galactic nucleus (AGN). The optical and near-IR spectra show broad 2000 km/s hydrogen, He I and O I emission features that decrease in flux over time. Radio imaging shows a fast evolving compact source of synchrotron emission spatially coincident with AT 2017gbl. We infer a lower limit for the radiated energy of 7.3 x 10^50 erg from the IR photometry. An extremely energetic supernova would satisfy this budget, but is ruled out by the radio counterpart evolution. Instead, we propose AT 2017gbl is related to an accretion event by the central supermassive black hole, where the spectral signatures originate in the AGN broad line region and the IR photometry is consistent with re-radiation by polar dust. Given the fast evolution of AT 2017gbl, we deem a tidal disruption event (TDE) of a star a more plausible scenario than a dramatic change in the AGN accretion rate. This makes AT 2017gbl the third TDE candidate to be hosted by a LIRG, in contrast to the so far considered TDE population discovered at optical wavelengths and hosted preferably by post-starburst galaxies.

We report on detailed spectral analysis of emission from X-ray pulsar 4U 1901+03 using the data obtained by the Insight-HXMT and NuSTAR observatories during the 2019 outburst of the source. Thanks to the extensive coverage of the outburst by Insight-HXMT, we were able to investigate spectral evolution of the source as a function of flux, and compare these results to the previous reports with focus on the properties of a putative absorption feature at around 10 keV. In particular, we demonstrate that the broadband X-ray continuum of 4U 1901+03 can be well described with a two-component continuum model without absorption line at 10 keV, which casts doubt on its interpretation as a cyclotron line. The high quality of the data has also allowed us to perform both phase-averaged and phase-resolved spectral analysis as function of luminosity. Finally, we performed a detailed investigation of another absorption feature in the spectrum of the source around 30 keV recently reported in the NuSTAR data. In line with the previous findings, we show that this feature appears to be significantly detected, albeit only in the pulse-phase resolved spectra at relatively high luminosities. The non-detection of the line in the averaged spectrum is possibly associated with the observed strong variations of its parameters with the pulse phase.

We use the In data collected during the 2019 outburst from X-ray pulsar 4U 1901+03 to complement the orbital parameters reported by Fermi/GBM. Using the Insight-HXMT, we examine the correlation between the derivative of the intrinsic spin frequency and bolometric flux based on accretion torque models. It was found that the pulse profiles significantly evolve during the outburst. The existence of two types of the profile's pattern discovered in the Insight-HXMT data indicates that this source experienced transition between a super-critical and a sub-critical accretion regime during its 2019 outburst. Based on the evolution of the pulse profiles and the torque model, we derive the distance to 4U 1901+03 as 12.4+-0.2 kpc.

PSRJ 2032+4127 is only the second known gamma-ray binary where it is confirmed that a young radio pulsar is in orbit around a Be-star. The interaction of the pulsar wind with the mass outflow from the companion leads to broad band emission from radio up to TeV energies. In the current paper we present results of optical monitoring of the 2017 periastron passage with the Nordic Optical Telescope. These observations are complemented by X-ray (Swift/XRT, NuSTAR) and GeV (Fermi/LAT) monitoring. Joint analysis of the evolution of the parameters of the Halpha line and the broadband (X-ray to TeV) spectral shape allows us to propose a model linking the observed emission to the interaction of the pulsar and Be-star winds under the assumption of the inclined disc geometry. Our model allows the observed flux and spectral evolution of the system to be explained in a self-consistent way.

We present observations from the short-term intensive optical campaign (from Sep2019 to Jan2020) of the changing-look Seyfert NGC 3516. This active galactic nucleus is known to have strong optical variability and has changed its type in the past. It has been in the low-activity state in the optical since 2013, with some rebrightening from the end of 2015 to the beginning of 2016, after which it remained dormant. We aim to study the photometric and spectral variability of NGC 3516 from the new observations in U- and B-bands and examine the profiles of the optical broad emission lines in order to demonstrate that this object may be entering a new state of activity. NGC 3516 has been monitored intensively for the past 4 months with an automated telescope in U and B filters, enabling accurate photometry of 0.01 precision. Spectral observations were triggered when an increase in brightness was spotted. We support our analysis of past-episodes of violent variability with the UV and X-ray long-term light curves constructed from the archival Swift data. An increase of the photometric magnitude is seen in both U and B filters to a maximum amplitude of 0.25mag and 0.1 mag, respectively. During the flare, we observe stronger forbidden high-ionization iron lines than reported before, as well as the complex broad Ha and Hb lines. This is especially seen in Ha which appears to be double-peaked. It seems that a very broad component of ~10,000 km/s in width in the Balmer lines is appearing. The trends in the optical, UV, and X-ray light curves are similar, with the amplitudes of variability being significantly larger in the case of UV and X-ray bands. The increase of the continuum emission, the variability of the coronal lines, and the very broad component in the Balmer lines may indicate that the AGN of NGC 3516 is finally leaving the low-activity state in which it has been for the last ~3 years.

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

We present 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.

IGR J17591-2342 is a recently INTEGRAL discovered accreting millisecond X-ray pulsar that went into outburst around July 21, 2018. To better understand the physics acting in these systems during the outburst episode we performed detailed temporal-, timing- and spectral analyses across the 0.3-300 keV band using data from NICER, XMM-Newton, NuSTAR and INTEGRAL. The hard X-ray 20-60 keV outburst profile is composed of four flares. During the maximum of the last flare we discovered a type-I thermonuclear burst in INTEGRAL JEM-X data. We derived a distance of 7.6+/-0.7 kpc, adopting Eddington luminosity limited photospheric radius expansion burst emission and assuming anisotropic emission. In the timing analysis using all NICER 1-10 keV monitoring data we observed a rather complex behaviour starting with a spin-up period, followed by a frequency drop, a episode of constant frequency and concluding with irregular behaviour till the end of the outburst. The 1-50 keV phase distributions of the pulsed emission, detected up to $\sim$ 120 keV using INTEGRAL ISGRI data, was decomposed in three Fourier harmonics showing that the pulsed fraction of the fundamental increases from ~10% to ~17% going from ~1.5 to ~4 keV, while the harder photons arrive earlier than the soft photons for energies <10 keV. The total emission spectrum of IGR J17591-2342 across the 0.3-150 keV band could adequately be fitted in terms of an absorbed compPS model yielding as best fit parameters a column density of N_H=(2.09+/-0.05) x 10^{22} /cm2, a blackbody seed photon temperature kT_bb,seed of 0.64+/- 0.02 keV, electron temperature kT_e=38.8+/-1.2 keV and Thomson optical depth Tau_T=1.59+/-0.04. The fit normalisation results in an emission area radius of 11.3+/-0.5 km adopting a distance of 7.6 kpc. Finally, the results are discussed within the framework of accretion physics- and X-ray thermonuclear burst theory.

We report on our analysis of the 2019 outburst of the X-ray accreting pulsar 4U 1901+03 observed with Insight-HXMT and NICER. Both spectra and pulse profiles evolve significantly in the decaying phase of the outburst. Dozens of flares are observed throughout the outburst. They are more frequent and brighter at the outburst peak. We find that the flares, which have a duration from tens to hundreds of seconds, are generally brighter than the persistent emission by a factor of $\sim$ 1.5. The pulse profile shape during the flares can be significantly different than that of the persistent emission. In particular, a phase shift is clearly observed in many cases. We interpret these findings as direct evidence of changes of the pulsed beam pattern, due to transitions between the sub- and super-critical accretion regimes on a short time scale. We also observe that at comparable luminosities the flares' pulse profiles are rather similar to those of the persistent emission. This indicates that the accretion on the polar cap of the neutron star is mainly determined by the luminosity, i.e., the mass accretion rate.

The hard X-ray transient source Swift J1845.7-0037 was discovered in 2012 by Swift/BAT. However, at that time no dedicated observations of the source were performed. On Oct 2019 the source became active again, and X-ray pulsations with a period of ~199s were detected with Swift/XRT. This triggered follow-up observations with NuSTAR. Here we report on the timing and spectral analysis of the source properties using NuSTAR and Swift/XRT. The main goal was to confirm pulsations and search for possible cyclotron lines in the broadband spectrum of the source to probe its magnetic field. Despite highly significant pulsations with period of 207.379(2) were detected, no evidence for a cyclotron line was found in the spectrum of the source. We therefore discuss the strength of the magnetic field based on the source flux and the detection of the transition to the "cold-disc" accretion regime during the 2012 outburst. Our conclusion is that, most likely, the source is a highly magnetized neutron star with B 1e13G at a large distance of d~10 kpc. The latter one consistent with the non-detection of a cyclotron line in the NuSTAR energy band.

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.

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 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.

We report results of the spectral and timing analysis of the Be/X-ray pulsar SXP 4.78 using the data obtained during its recent outburst with NuSTAR, Swift, Chandra and NICER observatories. Using an overall evolution of the system luminosity, spectral analysis and variability power spectrum we obtain constraints on the neutron star magnetic field strength. We found a rapid evolution of the variability power spectrum during the rise of the outburst, and absence of the significant changes during the flux decay. Several low frequency quasi-periodic oscillation features are found to emerge on the different stages of the outburst, but no clear clues on their origin were found in the energy spectrum and overall flux behaviour. We use several indirect methods to estimate the magnetic field strength on the neutron star surface and found that most of them suggest magnetic field $B \lesssim 2 \times10^{12}$ G. The strictest upper limit comes from the absence of the cyclotron absorption features in the energy spectra and suggests relatively weak magnetic field $B < 6 \times 10^{11}$ G.

We report a 72 ks XMM-Newton observation of the Be/X-ray pulsar (BeXRP) RX J0812.4-3114 in quiescence ($L_X \approx 1.6 \times 10^{33}~\mathrm{erg~s^{-1}}$). Intriguingly, we find a two component spectrum, with a hard power-law ($\Gamma \approx 1.5$) and a soft blackbody-like excess below $\approx 1~\mathrm{keV}$. The blackbody component is consistent in $kT$ with a prior quiescent Chandra observation reported by Tsygankov et al. and has an inferred blackbody radius of $\approx 10~\mathrm{km}$, consistent with emission from the entire neutron star (NS) surface. There is also mild evidence for an absorption line at $\approx 1~\mathrm{keV}$ and/or $\approx 1.4~\mathrm{keV}$. The hard component shows pulsations at $P \approx 31.908~\mathrm{s}$ (pulsed fraction $0.84 \pm 0.10$), agreeing with the pulse period seen previously in outbursts, but no pulsations were found in the soft excess (pulsed fraction $\lesssim 31\%$). We conclude that the pulsed hard component suggests low-level accretion onto the neutron star poles, while the soft excess seems to originate from the entire NS surface. We speculate that, in quiescence, the source switches between a soft thermal-dominated state (when the propeller effect is at work) and a relatively hard state with low-level accretion, and use the propeller cutoff to estimate the magnetic field of the system to be $\lesssim 8.4 \times 10^{11}~\mathrm{G}$. We compare the quiescent thermal $L_X$ predicted by the standard deep crustal heating model to our observations and find that RX J0812.4-3114 has a high thermal $L_X$, at or above the prediction for minimum cooling mechanisms. This suggests that RX J0812.4-3114 either contains a relatively low-mass NS with minimum cooling, or that the system may be young enough that the NS has not fully cooled from the supernova explosion.

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 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 report on analysis of properties of the X-ray binary pulsar GROJ1750-27 based on X-ray (Chandra, Swift, and Fermi/GBM), and near-infrared (VVV and UKIDSS surveys) observations. An accurate position of the source is determined for the first time and used to identify its infrared counterpart. Based on the VVV data we investigate the spectral energy distribution (SED) of the companion, taking into account a non-standard absorption law in the source direction. A comparison of this SED with those of known Be/X-ray binaries and early type stars has allowed us to estimate a lower distance limit to the source at $>12$ kpc. An analysis of the observed spin-up torque during a giant outburst in 2015 provides an independent distance estimate of $14-22$ kpc, and also allows to estimate the magnetic field on the surface of the neutron star at $B\simeq(3.5-4.5)\times10^{12}$ G. The latter value is in agreement with the possible transition to the propeller regime, a strong hint for which was revealed by Swift/XRT and Chandra. We conclude that GROJ1750-27 is located far behind the Galactic Center, which makes it one of the furthest Galactic X-ray binaries known.

Aims. The optical emission of black hole transients increases by several magnitudes during the X-ray outbursts. Whether the extra light arises from the X-ray heated outer disc, from the inner hot accretion flow, or from the jet is currently debated. Optical polarisation measurements are able to distinguish the relative contributions of these components. Methods. We present the results of BVR polarisation measurements of the black hole X-ray binary MAXI J1820+070 during the period of March-April 2018. Results. We detect small, $\sim$0.7%, but statistically significant polarisation, part of which is of interstellar origin. Depending on the interstellar polarisation estimate, the intrinsic polarisation degree of the source is between $\sim$0.3% and 0.7%, and the polarisation position angle is between $\sim10\deg-30\deg$. We show that the polarisation increases after MJD 58222 (2018 April 14). The change is of the order of 0.1% and is most pronounced in the R band. The change of the source Stokes parameters occurs simultaneously with the drop of the observed V-band flux and a slow softening of the X-ray spectrum. The Stokes vectors of intrinsic polarisation before and after the drop are parallel, at least in the V and R filters. Conclusions. We suggest that the increased polarisation is due to the decreasing contribution of the non-polarized component, which we associate with the the hot flow or jet emission. The low polarisation can result from the tangled geometry of the magnetic field or from the Faraday rotation in the dense, ionised, and magnetised medium close to the black hole. The polarized optical emission is likely produced by the irradiated disc or by scattering of its radiation in the optically thin outflow.

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

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.

We present a study of optical, UV and X-ray light curves of the nearby changing look active galactic nucleus in the galaxy NGC 1566 obtained with the Neil Gehrels Swift Observatory and the MASTER Global Robotic Network over the period 2007 - 2018. We also report on our optical spectroscopy at the South African Astronomical Observatory with the 1.9-m telescope on the night 2018 August 2-3. A substantial increase in X-ray flux by 1.5 orders of magnitude was observed following the brightening in the UV and optical bands during the last year. After a maximum was reached at the beginning of 2018 July the fluxes in all bands decreased with some fluctuations. The amplitude of the flux variability is strongest in the X-ray band and decreases with increasing wavelength. Low-resolution spectra reveal a dramatic strengthening of the broad emission as well as high-ionization [FeX]6374 A lines. These lines were not detected so strongly in the past published spectra. The change in the type of the optical spectrum was accompanied by a significant change in the X-ray spectrum. All these facts confirm NGC 1566 to be a changing look Seyfert galaxy.

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 results of a spectral and timing analysis of the X-ray pulsar XTE J1829-098 using data obtained with the NuSTAR observatory during an outburst in August 2018. A strong absorption feature was detected at the energy of $E_{cyc}\simeq 15$ keV in the source spectrum. This feature was interpreted as a cyclotron resonance scattering line corresponding to the magnetic field strength of the neutron star surface $B\simeq1.7\times10^{12}$ G. The pulse phase-resolved spectroscopy shows that the cyclotron line is significantly detected at all phases of the pulse and its energy and other parameters are variable over the pulse period. The timing analysis of the source emission revealed strong pulsations with a period of $P = 7.84480(2)$ s. The pulsed fraction is changed with the energy, including its local increase in the vicinity of the cyclotron line. Using the archival data of the RXTE observatory the presence of the cyclotron line in the spectrum of XTE J1829-098 was independently confirmed.

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.

The closest known ultraluminous X-ray source (ULX), M33 X-8, has been recently observed with NuSTAR during its Extragalactic Legacy program, which includes a hard X-ray survey of the M33 galaxy. We present results of two long observations of M33 taken in 2017 March and July, with M33 X-8 in the field of view. The source demonstrates a nearly constant flux during the observations, and its 3-20 keV spectrum can be well described by two distinct components: a standard accretion disc with a temperature of ~1 keV at the inner radius and a power law with a photon index ~3, which is significantly detected up to 20 keV. There is also an indication of a high-energy cutoff in the spectrum, corresponding to a temperature of the Comptonizing medium of >10 keV. The broad-band spectral properties of M33 X-8 resemble black hole X-ray binaries in their very high states, suggesting that M33 X-8 is a black hole accreting at a nearly Eddington rate, in contrast to super-Eddington accretion believed to take place in more luminous ULXs.

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.

Tidal disruption events (TDEs) are transient flares produced when a star is ripped apart by the gravitational field of a supermassive black hole (SMBH). We have observed a transient source in the western nucleus of the merging galaxy pair Arp 299 that radiated >1.5x10^52 erg in the infrared and radio, but was not luminous at optical or X-ray wavelengths. We interpret this as a TDE with much of its emission re-radiated at infrared wavelengths by dust. Efficient reprocessing by dense gas and dust may explain the difference between theoretical predictions and observed luminosities of TDEs. The radio observations resolve an expanding and decelerating jet, probing the jet formation and evolution around a SMBH.

We take advantage of a long (with a total exposure time of 120 ks) X-ray observation of the unique Galactic microquasar SS 433, carried out with the XMM-Newton space observatory, to search for a fluorescent line of neutral (or weakly ionized) nickel at the energy 7.5 keV. We consider two models of the formation of fluorescent lines in the spectrum of SS 433: 1) due to reflection of hard X-ray radiation from a putative central source on the optically thick walls of the accretion disk "funnel"; and 2) due to scattering of the radiation coming from the hottest parts of the jets in the optically thin wind of the system. It is shown, that for these cases, the photon flux of Ni I K$_{\alpha}$ fluorescent line is expected to be 0.45 of the flux of Fe I K$_{\alpha}$ fluorescent line at 6.4 keV, for the relative nickel overabundance $Z_{Ni}/Z = 10$, as observed in the jets of SS 433. For the continuum model without the absorption edge of neutral iron, we set a 90 per cent upper limit on the flux of the narrow Ni I K$_{\alpha}$ line at the level of $0.9 \times 10^{-5}$ ph s$^{-1}$ cm$^{-2}$. For the continuum model with the absorption edge, the corresponding upper limit is $2.5 \times 10^{-5}$ ph s$^{-1}$ cm$^{-2}$. At the same time, for the Fe I K$_{\alpha}$ line, we measure the flux of $9.9_{8.4}^{11.2} \times 10^{-5}$ ph s$^{-1}$ cm$^{-2}$. Taken at the face value, the results imply that the relative overabundance of nickel in the wind of the accretion disc should be at least 1.5 times less than the corresponding excess of nickel observed in the jets of SS 433.

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.

We present the orbital solution for the newly discovered transient Be X-ray binary Swift J0243.6+6124 based on the data from gamma-ray burst monitor onboard Fermi obtained during the Oct 2017 outburst. We model the Doppler induced and intrinsic spin variations of the neutron star assuming that the later is driven by accretion torque and discuss the implications of the observed spin variations for the parameters of the neutron star and the binary. In particular we conclude that the neutron star must be strongly magnetized, and estimate the distance to the source at $\sim$5 kpc.

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.

Optical and near-infrared photometry, optical spectroscopy, and soft X-ray and UV monitoring of the changing-look active galactic nucleus NGC 2617 show that it continues to have the appearance of a type-1 Seyfert galaxy. An optical light curve for 2010$-$2017 indicates that the change of type probably occurred between 2010 October and 2012 February and was not related to the brightening in 2013. In 2016 and 2017 NGC 2617 brightened again to a level of activity close to that in 2013 April. However, in 2017 from the end of the March to end of July 2017 it was in very low level and starting to change back to a Seyfert 1.8. We find variations in all passbands and in both the intensities and profiles of the broad Balmer lines. A new displaced emission peak has appeared in H$\beta$. X-ray variations are well correlated with UV$-$optical variability and possibly lead by $\sim$ 2$-$3 d. The $K$ band lags the $J$ band by about 21.5 $\pm$ 2.5 d and lags the combined $B + J$ bands by $\sim$ 25 d. $J$ lags $B$ by $\sim$ 3 d. This could be because $J$-band variability arises predominantly from the outer part of the accretion disc, while $K$-band variability is dominated by thermal re-emission by dust. We propose that spectral-type changes are a result of increasing central luminosity causing sublimation of the innermost dust in the hollow bi-conical outflow. We briefly discuss various other possible reasons that might explain the dramatic changes in NGC 2617.

The dependence of the spin frequency derivative $\dot{\nu}$ of accreting neutron stars with a strong magnetic field (X-ray pulsars) on the mass accretion rate (bolometric luminosity, $L_{bol}$) has been investigated for eight transient pulsars in binary systems with Be stars. Using data from the Fermi/GBM and Swift/BAT telescopes, we have shown that for seven of the eight systems the dependence $\dot{\nu}$ can be fitted by the model of angular momentum transfer through an accretion disk, which predicts the relation $\dot{\nu}\sim L^{6/7}_{bol}$. Hysteresis in the dependence $\dot{\nu}(L_{bol})$ has been confirmed in the system V 0332+53 and has been detected for the first time in the systems KS 1947+300, GRO J1008-57, and 1A 0535+26. The radius of the neutron star magnetosphere in all of the investigated systems have been estimated. We show that this quantity varies from pulsar to pulsar and depends strongly on the analytical model and the estimates for the neutron star and binary system parameters.

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.

Observations of accreting neutron stars (NS) with strong magnetic fields can be used not only for studying the accretion flow interaction with NS magnetospheres, but also for understanding the physical processes inside NSs and for estimating their fundamental parameters. Of particular interest are (i) the interaction of a rotating neutron star (magnetosphere) with the in-falling matter at different accretion rates, and (ii) the theory of deep crustal heating and the influence of a strong magnetic field on this process. Here, we present results of the first systematic investigation of 16 X-ray pulsars with Be optical companions during their quiescent states, based on data from the Chandra, XMM-Newton and Swift observatories. The whole sample of sources can be roughly divided into two distinct groups: i) relatively bright objects with a luminosity around ~10^34 erg/s and (hard) power-law spectra, and ii) fainter ones showing thermal spectra. X-ray pulsations were detected from five objects in group i) with quite a large pulse fraction of 50-70 per cent. The obtained results are discussed within the framework of the models describing the interaction of the in-falling matter with the neutron star magnetic field and those describing heating and cooling in accreting NSs.

LSI +61 303 is one of the few X-ray binaries with a Be star companion from which radio, X-rays and high-energy gamma-ray (GeV and TeV) emission have been observed. The nature of the high energy activity of the system is not yet fully understood, but it is widely believed that it is generated due to the interaction of the relativistic electrons leaving the compact object with the photons and non-relativistic wind of the Be star. The superorbital variability of the system has been observed in the radio, optical and X-ray domains and could be due to the cyclic change of the Be star disk size. In this paper we systematically review all publicly available data from Suzaku, XMM-Newton, Chandra and Swift observatories in order to measure the absorption profile of the circumstellar Be disk as a function of orbital and superorbital phases. We also discuss short-term variability of the system, found during the analysis and its implications for the understanding of the physical processes in this system.

The International Gamma-Ray Astrophysics Laboratory (INTEGRAL) continues to successfully work in orbit after its launch in 2002. The mission provides the deepest ever survey of hard X-ray sources throughout the Galaxy at energies above 20 keV. We report on a catalogue of new hard X-ray source candidates based on the latest sky maps comprising 14 years of data acquired with the IBIS telescope onboard INTEGRAL in the Galactic Plane (|b|<17.5 deg). The current catalogue includes in total 72 hard X-ray sources detected at S/N>4.7 sigma and not known to previous INTEGRAL surveys. Among them, 31 objects have also been detected in the on-going all-sky survey by the BAT telescope of the Swift observatory. For 26 sources on the list, we suggest possible identifications: 21 active galactic nuclei, two cataclysmic variables, two isolated pulsars or pulsar wind nebulae, and one supernova remnant; 46 sources from the catalogue remain unclassified.

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.

Optical and near-infrared photometry, optical spectroscopy, and soft X-ray and UV monitoring of the changing look active galactic nucleus NGC 2617 show that it continues to have the appearance of a type-1 Seyfert galaxy. An optical light curve for 2010-2016 indicates that the change of type probably occurred between 2010 October and 2012 February and was not related to the brightening in 2013. In 2016 NGC 2617 brightened again to a level of activity close to that in 2013 April. We find variations in all passbands and in both the intensities and profiles of the broad Balmer lines. A new displaced emission peak has appeared in H$\beta$. X-ray variations are well correlated with UV-optical variability and possibly lead by $\sim$ 2-3 d. The $K$ band lags the $J$ band by about 21.5 $\pm$ 2.5 d. and lags the combined $B+J$ filters by $\sim$ 25 d. $J$ lags $B$ by about 3 d. This could be because $J$-band variability arises from the outer part of the accretion disc, while $K$-band variability comes from thermal re-emission by dust. We propose that spectral-type changes are a result of increasing central luminosity causing sublimation of the innermost dust in the hollow biconical outflow. We briefly discuss various other possible reasons that might explain the dramatic changes in NGC 2617.

Our simultaneous three-colour ($BVR$) polarimetric observations of the low-mass black hole X-ray binary V404 Cyg show a small but statistically significant change of polarization degree ($\Delta p \sim 1$ per cent) between the outburst in June 2015 and the quiescence. The polarization of V404 Cyg in the quiescent state agrees within the errors with that of the visually close (1\farcs4) companion ($p_{R} = 7.3\pm 0.1$ per cent), indicating that it is predominantly of interstellar origin. The polarization pattern of the surrounding field stars supports this conclusion. From the observed variable polarization during the outburst we show that polarization degree of the intrinsic component peaks in the $V$-band, $p_{V} = 1.1\pm 0.1$ per cent, at the polarization position angle of $\theta_{V}= -7\deg\pm 2\deg$, which is consistent in all three passbands. We detect significant variations in the position angle of the intrinsic polarization in $R$ band from $-30\deg$ to $\sim 0\deg$ during the outburst peak. The observed wavelength dependence of the intrinsic polarization does not support non-thermal synchrotron emission from a jet as a plausible mechanism, but is in better agreement with the combined effect of electron (Thomson) scattering and absorption in a flattened plasma envelope or outflow surrounding the illuminating source. Alternatively, the polarization signal can be produced by scattering of the disc radiation in a mildly relativistic polar outflow. The position angle of the intrinsic polarization, nearly parallel to the jet direction (i.e. perpendicular to the accretion disc plane), is in agreement with these interpretations.

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 report the results of the monitoring campaign of the transient X-ray pulsar SMC X-2 performed with the Swift/XRT telescope over the period of 2015 September - 2016 January during the Type II outburst. During this event, the bolometric luminosity of the source ranged from $\simeq10^{39}$ down to several$\times10^{34}$ erg/s. Moreover, we discovered its dramatic drop by a factor of more than 100 below the limiting value of $L_{\rm lim}\simeq4\times10^{36}$ erg/s, which can be interpreted as a transition to the propeller regime. These measurements make SMC X-2 the sixth pulsating X-ray source where such a transition is observed and allow us to estimate the magnetic field of the neutron star in the system $B\simeq3\times10^{12}$ G, which is in agreement with independent results of the spectral analysis.

We report on the first focused hard X-ray view of the absorbed supergiant system IGRJ18027-2016 performed with the NuSTAR observatory. The pulsations are clearly detected with a period of $P_{spin}=139.866(1)$ s and a pulse fraction of about 50-60% at energies from 3 to 80 keV. The source demonstrates an approximately constant X-ray luminosity on a time scale of more than dozen years with an average spin-down rate of $dP/dt\simeq6\times10^{-10}$ s/s. This behaviour of the pulsar can be explained in terms of the wind accretion model in the settling regime. The detailed spectral analysis at energies above 10 keV was performed for the first time and revealed a possible cyclotron absorption feature at energy ~23 keV. This energy corresponds to the magnetic field $B\simeq3\times10^{12}$ G at the surface of the neutron star, which is typical for X-ray pulsars.

We present results of long NuSTAR (200 ks) and XMM-Newton (100 ks) observations of the Arches stellar cluster, a source of bright thermal (kT~2 keV) X-rays with prominent Fe XXV K_alpha 6.7 keV line emission and a nearby molecular cloud, characterized by an extended non-thermal hard X-ray continuum and fluorescent Fe K_alpha 6.4 keV line of a neutral or low ionization state material around the cluster. Our analysis demonstrates that the non-thermal emission of the Arches cloud underwent a dramatic change, with its homogeneous morphology, traced by fluorescent Fe K_alpha line emission, vanishing after 2012, revealing three bright clumps. The declining trend of the cloud emission, if linearly fitted, is consistent with half-life decay time of ~8 years. Such strong variations have been observed in several other molecular clouds in the Galactic Centre, including the giant molecular cloud Sgr B2, and point toward a similar propagation of illuminating fronts, presumably induced by the past flaring activity of Sgr A*.

We present a study of the evolution of the optical and X-ray fluxes during outbursts of two short-period cataclysmic variables, the confirmed intermediate polar CC Scl and the intermediate polar candidate FS Aur. We found that the X-ray and optical light curves are well correlated in both objects, although the amplitudes of outbursts in X-rays are smaller than those in the optical. The ratio of the outburst amplitudes in X-rays and the optical in both objects is close to ~0.6. This is significantly higher than was observed during the outburst of the non-magnetic dwarf nova U Gem, in which this ratio was only ~0.03. The obtained data also suggest that the dependence between the X-ray and optical fluxes must steepen significantly toward very low accretion rates and very low fluxes. Similarities in the behaviour of CC Scl and FS Aur indicate strongly the magnetic nature of the white dwarf in FS Aur.

We describe the evolution of optical and X-ray temporal characteristics during the outburst decline of the black hole X-ray binary SWIFT J1753.5-0127. The optical/X-ray cross-correlation function demonstrates a single positive correlation at the outburst peak, then it has multiple dips and peaks during the decline stage, which are then replaced by the precognition dip plus peak structure in the outburst tail. Power spectral densities and phase lags show a complex evolution, revealing the presence of intrinsically connected optical and X-ray quasi-periodic oscillations. For the first time, we quantitatively explain the evolution of these timing properties during the entire outburst within one model, the essence of which is the expansion of the hot accretion flow towards the tail of the outburst. The pivoting of the spectrum produced by synchrotron Comptonization in the hot flow is responsible for the appearance of the anti-correlation with the X-rays and for the optical quasi-periodic oscillations. Our model reproduces well the cross-correlation and phase lag spectrum during the decline stage, which could not be understood with any model proposed before.

In this paper we investigated the long-term evolution of the pulse-period in the high-mass X-ray binary LMC X-4 by taking advantage of more than 43~yrs of measurements in the X-ray domain. Our analysis revealed for the first time that the source is displaying near-periodical variations of its spin period on a time scale of roughly 6.8~yrs, making LMC X-4 one of the known binary systems showing remarkable long term spin torque reversals. We discuss different scenarios to interpret the origin of these torque reversals.

We present results of the spectral and timing analysis of the X-ray pulsar LMC X-4 with the NuSTAR observatory in the broad energy range 3-79 keV. Along with the detailed analysis of the averaged source spectrum, the high-precision pulse phased-resolved spectra were obtained for the first time. It has been shown that the comptonization model gives the best approximation of the obtained spectra. The evolution of its parameters was traced depending on the pulse phase as well. The search for the possible cyclotron absorption line was performed for all energy spectra in the 5-55 keV energy range. The obtained upper limit for the depth of the cyclotron absorption line $\tau\simeq0.15$ ($3\sigma$) indicates no cyclotron absorption line in this energy range, which provides an estimate of the magnitude of the magnetic field on the surface of the neutron star: $B < 3 \times 10^{11}$ G or $B > 6.5 \times 10^{12}$ G. The latter one is agree with the estimate of the magnetic field obtained from the analysis of the power spectrum of the pulsar: $B \sim 3 \times 10^{13}$ G. Based on results of the pulse phase-resolved spectroscopy we revealed a delay between maxima of the source emission and the equivalent width of the fluorescent iron line. This delay can be apparently associated with the travel time of photons between the emitting regions in the vicinity of the neutron star and the relatively cold regions where this emission is reflected (presumably, at the inflowing stream or at the place of an interaction of the stream and the outer edge of the accretion disk).

We report the results of a comprehensive analysis of X-ray (Chandra and Swift observatories), optical (Southern African Large Telescope, SALT) and near-infrared (the VVV survey) observations of the Be/X-ray binary pulsar 2S1553-542. Accurate coordinates for the X-ray source are determined and are used to identify the faint optical/infrared counterpart for the first time. Using VVV and SALTICAM photometry, we have constructed the spectral energy distribution (SED) for this star and found a moderate NIR excess that is typical for Be stars and arises due to the presence of circumstellar material (disk). A comparison of the SED with those of known Be/X-ray binaries has allowed us to estimate the spectral type of the companion star as B1-2V and the distance to the system as $>15$ kpc. This distance estimation is supported by the X-ray data and makes 2S1553-542 one of the most distant X-ray binaries within the Milky Way, residing on the far side in the Scutum-Centaurus arm or even further.

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.

Aims. The black hole binary SWIFT J1753.5-0127 is providing a unique data set to study accretion flows. Various investigations of this system and of other black holes have not, however, led to an agreement on the accretion flow geometry or on the seed photon source for Comptonization during different stages of X-ray outbursts. We place constraints on these accretion flow properties by studying long-term spectral variations of this source. Methods. We performed phenomenological and self-consistent broad band spectral modeling of Swift J1753.5-0127 using quasi-simultaneous archived data from INTEGRAL/ISGRI, Swift/UVOT/XRT/BAT, RXTE/PCA/HEXTE and MAXI/GSC instruments. Results. We identify a critical flux limit, F \sim 1.5 \times 10^{-8} erg/cm^2/s, and show that the spectral properties of SWIFT J1753.5-0127 are markedly different above and below this value. Above the limit, during the outburst peak, the hot medium seems to intercept roughly 50 percent of the disk emission. Below it, in the outburst tail, the contribution of the disk photons reduces significantly and the entire spectrum from the optical to X-rays can be produced by a synchrotron-self-Compton mechanism. The long-term variations in the hard X-ray spectra are caused by erratic changes of the electron temperatures in the hot medium. Thermal Comptonization models indicate unreasonably low hot medium optical depths during the short incursions into the soft state after 2010, suggesting that non-thermal electrons produce the Comptonized tail in this state. The soft X-ray excess, likely produced by the accretion disk, shows peculiarly stable temperatures for over an order of magnitude changes in flux. Conclusions. The long-term spectral trends of SWIFT J1753.5-0127 are likely set by variations of the truncation radius and a formation of a hot, quasi-spherical inner flow in the vicinity of the black hole. (abridged)

We investigate the structure of the O-type binary system HD 48099 by measuring linear polarization that arises due to light scattering process. Linear polarization measurements of HD 48099 in the B, V and R passbands with the high-precision Dipol-2 polarimeter have been carried out. The data have been obtained with the 60 cm KVA (Observatory Roque de los Muchachos, La Palma, Spain) and T60 (Haleakala, Hawaii, USA) remotely controlled telescopes during 31 observing nights. Polarimetry in the optical wavelengths has been complemented by observations in the X-rays with the SWIFT space observatory. Optical polarimetry revealed small intrinsic polarization in HD 48099 with 0.1% peak to peak variation over the orbital period of 3.08 days. The variability pattern is typical for binary systems, showing strong second harmonic of the orbital period. We apply our model code for the electron scattering in the circumstellar matter to put constraints on the system geometry. A good model fit is obtained for scattering of light on a cloud produced by the colliding stellar winds. The geometry of the cloud, with a broad distribution of scattering particles away from the orbital plane, helps in constraining the (low) orbital inclination. We derive from the polarization data the inclination i = 17 deg +/- 2deg and the longitude of the ascending node Omega = 82 deg +/- 1 deg of the binary orbit. The available X-ray data provide additional evidence for the existence of the colliding stellar winds in the system. Another possible source of the polarized light could be scattering from the stellar photospheres. The models with circumstellar envelopes, or matter confined to the orbital plane, do not provide good constraints on the low inclination, better than i < 27 deg, as is already suggested by the absence of eclipses.

We report the results of the deepest Galactic Plane (|b| < 17.5 deg) survey in the 67.9 and 78.4 keV nuclear de-excitation lines of titanium-44 (44Ti) performed using the data acquired with the IBIS/ISGRI instrument onboard the INTEGRAL satellite during 12 years of operation. The peak sensitivity of our survey reached an unprecedented level of $4.8\times10^{-6}$ erg s$^{-1}$ cm$^{-2}$ (3 sigma) that improves the sensitivity of the survey done by CGRO/COMPTEL by a factor of ~5. As a result, constraining upper limits for all sources from the catalog of Galactic supernova remnants (SNRs; Green 2014) are derived. These upper limits can be used to estimate the exposure needed to detect 44Ti emission from any known SNR using existing and prospective X- and gamma-ray telescopes. Among the youngest Galactic SNRs, only Cas A shows significant 44Ti emission flux in good agreement with the NuSTAR measurements. We did not detect any other sources of titanium emission in the Galactic Plane at significance level higher than 5 sigma confirming previous claims of the rarity of such 44Ti-producing SNRs.

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.

The progenitors of hydrogen-poor core-collapse supernovae (SNe) of types Ib, Ic and IIb are believed to have shed their outer hydrogen envelopes either by extremely strong stellar winds, characteristic of classical Wolf-Rayet stars, or by binary interaction with a close companion star. The exact nature of the progenitors and the relative importance of these processes are still open questions. One relatively unexplored method to constrain the progenitors is to search for high-mass X-ray binaries (HMXB) at SN locations in pre-explosion X-ray observations. In a HMXB, one star has already exploded as a core-collapse SN, producing a neutron star or a stellar-mass black hole. It is likely that the second star in the system will also explode as a supernova, which should cause a detectable long-term change in the system's X-ray luminosity. In particular, a pre-explosion detection of a HMXB coincident with a SN could be informative about the progenitor's nature. In this paper we analyze pre-explosion ACIS observations of 18 nearby type Ib, Ic and IIb supernovae from the Chandra X-ray observatory public archive. Two sources that could potentially be associated with the supernova are identified in the sample. Additionally we make similar post-explosion measurements for 46 SNe. Although our modelling indicates that progenitor systems with compact binary companions are probably quite rare, studies of this type can in the future provide more stringent constraints as the number of discovered nearby SNe and suitable pre-explosion X-ray data are both increasing.

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 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.

We report results of a spectral and timing analysis of the poorly studied transient X-ray pulsar 2S 1553-542 using data collected with the NuSTAR and Chandra observatories and the Fermi/GBM instrument during an outburst in 2015. Properties of the source at high energies (>30 keV) are studied for the first time and the sky position had been essentially improved. The source broadband spectrum has a quite complicated shape and can be reasonably described by a composite model with two continuum components - a black body emission with the temperature about 1 keV at low energies and a power law with an exponential cutoff at high energies. Additionally an absorption feature at $\sim23.5$ keV is discovered both in phase-averaged and phase-resolved spectra and interpreted as the cyclotron resonance scattering feature corresponding to the magnetic field strength of the neutron star $B\sim3\times10^{12}$ G. Based on the Fermi/GBM data the orbital parameters of the system were substantially improved, that allowed us to determine the spin period of the neutron star P = 9.27880(3) s and a local spin-up $\dot P \simeq -7.5\times10^{-10}$ s s$^{-1}$ due to the mass accretion during the NuSTAR observations. Assuming accretion from the disk and using standard torque models we have estimated the distance to the system $d=20\pm4$ kpc.

We present the updated orbital solution for the transient Be X-ray binary V 0332+53 comple- menting historical measurements with the data from the gamma-ray burst monitor onboard Fermi obtained during the outburst in June-October 2015. We model the observed changes in the spin- frequency of the pulsar and deduce the orbital parameters of the system. We significantly improve existing constrains and show that contrary to the previous findings no change in orbital parameters is required to explain the spin evolution of the source during the outbursts in 1983, 2005 and 2015. The reconstructed intrinsic spin-up of the neutron star during the latest outburst is found to be comparable with previosly observed values and predictions of the accretion torque theory.

We report on the results of the extensive multi-wavelength campaign from optical to GeV gamma-rays of the 2014 periastron passage of PSR B1259-63, which is a unique high-mass gamma-ray emitting binary system with a young pulsar companion. Observations demonstrate the stable nature of the post-periastron GeV flare and prove the coincidence of the flare with the start of rapid decay of the H$\alpha$ equivalent width, usually interpreted as a disruption of the Be stellar disk. Intensive X-ray observations reveal changes in the X-ray spectral behaviour happening at the moment of the GeV flare. We demonstrate that these changes can be naturally explained as a result of synchrotron cooling of monoenergetic relativistic electrons injected into the system during the GeV flare.

This paper is devoted to determining the nature of the hard X-ray source IGR J17463-2854 located toward the Galactic bulge. Using data from the INTEGRAL and Chandra X-ray observatories, we show that five point X-ray sources with approximately identical fluxes in the 2-10 keV energy band are detected in the error circle of the object under study. In addition, significant absorption at low energies has been detected in the spectra of all these sources. Based on data from the VVV (VISTA/ESO) infrared Galactic Bulge Survey, we have unambiguously identified three of the five sources, determined the J, H and K magnitudes of the corresponding stars, and obtained upper limits on the fluxes for the remaining two sources. Analysis of the color-magnitude diagrams has shown that one of these objects most likely belongs to a class of rarely encountered objects, symbiotic binary systems (several tens are known with certainty), i.e., low-mass binary systems consisting of a white dwarf and a red giant. Note that all our results were obtained using improved absorption values and an extinction law differing in this direction from the standard one.

We present spectral and timing analysis of NuSTAR observations of the accreting X-ray pulsar 2RXP J130159.6-635806. The source was serendipitously observed during a campaign focused on the gamma-ray binary PSR B1259-63 and was later targeted for a dedicated observation. The spectrum has a typical shape for accreting X-ray pulsars, consisting of a simple power law with an exponential cutoff starting at ~7 keV with a folding energy of E_fold=~18 keV. There is also an indication of the presence of a 6.4 keV iron line in the spectrum at the ~3 sigma significance level. NuSTAR measurements of the pulsation period reveal that the pulsar has undergone a strong and steady spin-up for the last 20 years. The pulsed fraction is estimated to be ~80%, and is constant with energy up to 40 keV. The power density spectrum shows a break towards higher frequencies relative to the current spin period. This, together with steady persistent luminosity, points to a long-term mass accretion rate high enough to bring the pulsar out of spin equilibrium.

High-mass X-ray binaries are fundamental in the study of stellar evolution, nucleosynthesis, structure and evolution of galaxies and accretion processes. Hard X-rays observations by INTEGRAL and Swift have broad- ened significantly our understanding in particular for the super-giant systems in the Milky Way, which number has increased by almost a factor of three. INTEGRAL played a crucial role in the discovery, study and understanding of heavily obscured systems and of fast X-ray transients. Most super-giant systems can now be classified in three categories: classical/obscured, eccentric and fast transient. The classical systems feature low eccentricity and variability factor of about 1000, mostly driven by hydrodynamic phenomena occurring on scales larger than the accretion radius. Among them, systems with short orbital periods and close to Roche-Lobe overflow or with slow winds, appear highly obscured. In eccentric systems, the variability amplitude can reach even higher factors, because of the contrast of the wind density along the orbit. Four super-giant systems, featuring fast outbursts, very short orbital periods and anomalously low accretion rates, are not yet understood. Simulations of the accretion processes on relatively large scales have pro- gressed and reproduce parts of the observations. The combined effects of wind clumps, magnetic fields, neutron star rotation and eccentricity ought to be included in future modelling work. Observations with INTEGRAL in combination with other observatories were also important for detecting cyclotron resonant scattering features in spectra of X-ray pulsars, probing their variations and the geometry of the accretion column and emission regions. Finally, the unique characteristics of INTEGRAL and its long life time played a fundamental role for building a complete catalogue of HXMBs, to study the different populations of these systems ...

We present the results of the identification of six objects from the ASCA Galactic center and Galactic plane surveys: AXJ173548-3207, AXJ173628-3141, AXJ1739.5-2910, AXJ1740.4-2856, AXJ1740.5-2937, AXJ1743.9-2846. Chandra, XMM-Newton, and XRT/Swift X-ray data have been used to improve the positions of the optical counterparts to these sources. Thereafter, we have carried out a series of spectroscopic observations of the established optical counterparts at the RTT-150 telescope. Analysis of X-ray and optical spectra as well as photometric measurements in a wide wavelength range based on optical and infrared catalogs has allowed the nature of the program sources to be determined. Two X-ray objects have been detected in the error circle of AXJ173628-3141: one is a coronally active G star and the other may be a symbiotic star, a red giant with an accreting white dwarf. Three sources (AXJ1739.5-2910, AXJ1740.5-2937, AXJ1743.9-2846) have turned out to be active G-K stars, presumably RS CVn objects, one (AXJ1740.4-2856) is an M dwarf, and another one (AXJ173548-3207) may be a low-mass X-ray binary in its low state. The distances and corresponding luminosities of the sources in the soft X-ray band have been estimated; analysis of deep INTEGRAL Galactic Center observations has not revealed a statistically significant flux at energies higher 20 keV from any of them.

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.

We present results of all sky survey, performed with data acquired by the IBIS telescope onboard the INTEGRAL observatory over eleven years of operation, at energies above 100 keV. A catalogue of detected sources includes 132 objects. The statistical sample detected on the time-averaged 100-150 keV map at a significance above 5 sigma contains 88 sources: 28 AGNs, 38 LMXBs, 10 HMXBs and 12 rotation-powered young X-ray pulsars. The catalogue includes also 15 persistent sources, which were registered with the significance 4 sigma < S/N < 5 in hard X-rays, but at the same time were firmly detected (>12 sigma) in the 17-60 keV energy band. All sources from these two groups are known X-ray emitters, that means that the catalogue has 100% purity in respect to them. Additionally, 29 sources were found in different time intervals. In the context of the survey we present a hardness ratio of galactic and extragalactic sources, a LMXBs longitudinal asymmetry and a number-flux relation for non-blazar AGNs. At higher energies, in the 150-300 keV energy band, 25 sources have been detected with a signal-to-noise ratio S/N > 5 sigma, including 7 AGNs, 13 LMXBs, 3 HMXBs, and 2 rotation-powered pulsars. Among LMXBs and HMXBs we identified 12 black hole candidates (BHC) and 4 neutron star (NS) binaries.

The black hole (BH) candidate SWIFT J1753.5-0127 has remained active since the onset of its 2005 outburst. Emission lines in the optical spectrum were observed at the very beginning of the outburst, but since then the spectrum has been featureless making a precise BH mass estimation impossible. Here we present results from our optical and UV observations of SWIFT J1753.5-0127 taken in 2012-2013. Our new observations show extremely broad, double-peaked emission lines in the optical and UV spectra. The optical data also show narrow absorption and emission features with nearly synchronous and significant Doppler motions. A radial velocity study of these lines which we associate with the secondary star, yields a semi-amplitude of K_2=382 km/s. A time-series analysis of the spectral and photometric data revealed a possible orbital periodicity of 2.85 h, significantly shorter than the reported 3.2 h periodic signal by Zurita et al. The observed variability properties argue against a low orbital inclination angle and we present several observational arguments in favour of the BH interpretation. However, the measured radial velocity semi-amplitude of the donor star and the short orbital period imply that SWIFT J1753.5-0127 has one of the lowest measured mass function for a BH in a low-mass X-ray binary. We show that the compact object mass in excess of 5 Msun is highly improbable. Thus, SWIFT J1753.5-0127 is a BH binary that has one of the shortest orbital period and hosts probably one of the smallest stellar-mass BH found to date.

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$.

We study hard X-ray emission of the brightest accreting neutron star Sco X-1 with INTEGRAL observatory. Up to now INTEGRAL have collected ~4 Msec of deadtime corrected exposure on this source. We show that hard X-ray tail in time average spectrum of Sco X-1 has a power law shape without cutoff up to energies ~200-300 keV. An absence of the high energy cutoff does not agree with the predictions of a model, in which the tail is formed as a result of Comptonization of soft seed photons on bulk motion of matter near the compact object. The amplitude of the tail varies with time with factor more than ten with the faintest tail at the top of the so-called flaring branch of its color-color diagram. We show that the minimal amplitude of the power law tail is recorded when the component, corresponding to the innermost part of optically thick accretion disk, disappears from the emission spectrum. Therefore we show that the presence of the hard X-ray tail may be related with the existence of the inner part of the optically thick disk. We estimate cooling time for these energetic electrons and show that they can not be thermal. We propose that the hard X-ray tail emission originates as a Compton upscattering of soft seed photons on electrons, which might have initial non-thermal distribution.

Accretion onto magnetized neutron stars is expected to be centrifugally inhibited at low accretion rates. Several sources, however, are known to pulsate in quiescence at luminosities below the theoretical limit predicted for the onset of the centrifugal barrier. The source 1A 0535+262 is one of them. Here we present the results of an analysis of a ~50 ks long XMM-Newton observation of 1A 0535+262 in quiescence. At the time of the observation, the neutron star was close to apastron, and the source had remained quiet for two orbital cycles. In spite of this, we detected a pulsed X-ray flux of ~3e-11 erg/cm2/s . Several observed properties, including the power spectrum, remained similar to those observed in the outbursts. Particularly, we have found that the frequency of the break detected in the quiescent noise power spectrum follows the same correlation with flux observed when the source is in outburst. This correlation has been associated with the truncation of the accretion disk at the magnetosphere boundary. We argue that our result, along with other arguments previously reported in the literature, suggests that the accretion in quiescence also proceeds from an accretion disk around the neutron star. The proposed scenario consistently explains the energy of the cyclotron line observed in 1A 0535+262, and the timing properties of the source including the spin frequency evolution within and between the outbursts, and the frequency of the break in power spectrum.

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.

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.

Four hard X-ray sources from the INTEGRAL and Swift catalogs have been identified. X-ray and optical spectra have been obtained for each of the objects being studied by using data from the INTEGRAL, Swift, ROSAT, and Chandra X-ray observatories as well as observations with the RTT-150 and AZT-33IK optical telescopes. Two sources (SWIFT J1553.6+2606 and SWIFT J1852.2+8424) are shown to be extragalactic in nature: the first is a quasar, while the registered X-ray flux from the second is the total emission from two Seyfert 1 galaxies at redshifts 0.1828 and 0.2249. The source IGR J22534+6243 resides in our Galaxy and is an X-ray pulsar with a period of ~46.674 s that is a member of a high-mass X-ray binary, probably with a Be star. The nature of yet another Galactic source, SWIFT J1852.8+3002, is not completely clear and infrared spectroscopy is needed to establish it.

An overview of the results of observations for the transient X-ray pulsar 4U 0115+63, a member of a binary system with a Be star, since its discovery to the present day (~40 years) based on data from more than dozen observatories and instruments is presented. A overall light curve and the history of change in the spin frequency of the neutron star over the entire history of its observations, which also includes the results of recent measurements made by the INTEGRAL observatory during the 2004, 2008, and 2011 outbursts, are provided. The source's energy spectra have also been constructed from the INTEGRAL data obtained during the 2011 outburst for a dynamic range of its luminosities (~10^{37} - 7 x 10^{37} erg/s). We show that apart from the fundamental harmonic of the cyclotron absorption line at energy ~11 keV, its four higher harmonics at energies ~24, 35.6, 48.8, and 60.7 keV are detected in the spectrum. We have performed a detailed analysis of the source's spectra in the 4-28 keV energy band based on all of the available RXTE archival data obtained during bright outbursts in 1995-2011. We have confirmed that modifying the source's continuum model can lead to the disappearance of the observed anticorrelation between the energy of the fundamental harmonic of the cyclotron absorption line and the source's luminosity. Thus, the question about the evolution of the cyclotron absorption line energy with the luminosity of the X-ray pulsar 4U 0115+63 remains open and a physically justified radiation model for X-ray pulsars is needed to answer it.

We present results of the study of persistent high mass X-ray binaries (HMXBs) in the Milky Way, obtained from the deep INTEGRAL Galactic plane survey. This survey provides us a new insight into the population of high mass X-ray binaries because almost half of the whole sample consists of sources discovered with INTEGRAL. It is demonstrated for the first time that the majority of persistent HMXBs have supergiant companions and their luminosity function steepens somewhere around ~2x10^{36} erg/s. We show that the spatial density distribution of HMXBs correlates well with the star formation rate distribution in the Galaxy. The vertical distribution of HMXBs has a scale-height h~85 pc, that is somewhat larger than the distribution of young stars in the Galaxy. We propose a simple toy model, which adequately describes general properties of HMXBs in which neutron stars accrete a matter from the wind of the its companion (wind-fed NS-HMXBs population). Using the elaborated model we argue that a flaring activity of so-called supergiant fast X-ray transients, the recently recognized sub-sample of HMXBs, is likely related with the magnetic arrest of their accretion.

It is assumed that the radioactive decay of Ti-44 powers the infrared, optical and UV emission of supernova remnants after the complete decay of Co-56 and Co-57 (the isotopes that dominated the energy balance during the first three to four years after the explosion) until the beginning of active interaction of the ejecta with the surrounding matter. Simulations show that the initial mass of Ti-44 synthesized in core-collapse supernovae is (0.02-2.5) x 10^{-4} solar masses (M_sun). Hard X-rays and gamma-rays from the decay of this Ti-44 have been unambiguously observed from Cassiopeia A only, leading to the suggestion that the values of the initial mass of Ti-44 near the upper bound of the predictions occur only in exceptional cases. For the remnant of supernova 1987A, an upper limit to the initial mass of Ti-44 of < 10^{-3} M_sun has been obtained from direct X-ray observations, and an estimate of (1-2) x 10^{-4} M_sun has been made from infrared light curves and ultraviolet spectra by complex model-dependent computations. Here we report observations of hard X-rays from the remnant of supernova 1987A in the narrow band containing two direct-escape lines of Ti-44 at 67.9 and 78.4 keV. The measured line fluxes imply that this decay provided sufficient energy to power the remnant at late times. We estimate that the initial mass of Ti-44 was (3.1+/-0.8) x 10^{-4} M_sun, which is near the upper bound of theoretical predictions.

Results of the deep survey of the Large Magellanic Cloud (LMC), performed with the INTEGRAL observatory, are presented. The long exposure (~7 Ms) allowed us to detect twenty one sources in this sky region: ten belonging to the LMC itself (7 HMXBs, 2 PSRs, 1 LMXB), six of extragalactic origin and three belonging to other galaxies from the Local Group - the Milky Way (2 sources) and Small Magellanic Cloud (1 source). Four new hard X-ray sources of these 21 ones were discovered during the survey in addition to IGR J05414-6858 reported earlier; two of them were identified with extragalactic objects. We report also for the first time the detection of a hard X-ray emission from the Crab-like pulsar PSR J0537-6910 and identification of the hard X-ray source IGR J05305-6559 with the high-mass X-ray binary EXO 053109-6609.

Context. The INTEGRAL observatory operating in a hard X-ray/gamma domain has gathered a large observational data set over nine years starting in 2003. Most of the observing time was dedicated to the Galactic source population study, making possible the deepest Galactic survey in hard X-rays ever compiled. Aims. We aim to perform a Galactic survey that can be used as the basis of Galactic source population studies, and perform mapping of the Milky Way in hard X-rays over the maximum exposure available at |b|<17.5 deg. Methods. We used sky reconstruction algorithms especially developed for the high quality imaging of INTEGRAL/IBIS data. Results. We present sky images, sensitivity maps, and catalogs of detected sources in the three energy bands 17-60, 17-35, and 35-80 keV in the Galactic plane at |b|<17.5 deg. The total number of sources in the reference 17-60 keV band includes 402 objects exceeding a 4.7 sigma detection threshold on the nine-year time-averaged map. Among the identified sources with known and tentatively identified natures, 253 are Galactic objects (108 low-mass X-ray binaries, 82 high-mass X-ray binaries, 36 cataclysmic variables, and 27 are of other types), and 115 are extragalactic objects, including 112 active galactic nuclei (AGNs) and 3 galaxy clusters. The sample of Galactic sources with S/N>4.7 sigma has an identification completeness of ~92%, which is valuable for population studies. Since the survey is based on the nine-year sky maps, it is optimized for persistent sources and may be biased against finding transients.

In 2003-2012, the INTEGRAL observatory has performed long-term observations of the Large Magellanic Cloud (LMC). At present, this is one of the deepest hard X-ray (20-60 keV) surveys of extragalactic fields in which more than 20 sources of different natures have been detected. We present the results of a statistical analysis of the population of high-mass X-ray binaries in the LMC and active galactic nuclei (AGNs) observed in its direction. The hard X-ray luminosity function of high-mass X-ray binaries is shown to be described by a power law with a slope alpha~1.8, that in agreement with the luminosity function measurements both in the LMC itself, but made in the soft X-ray energy band, and in other galaxies. At the same time, the number of detected AGNs toward the LMC turns out to be considerably smaller than the number of AGNs registered in other directions, in particular, toward the source 3C 273. The latter confirms the previously made assumption that the distribution of matter in the local Universe is nonuniform.

We present results of a comprehensive analysis of the X-ray pulsar X Persei over the period 1996 to 2011, encompassing the quite low state and subsequent strong outburst activity. Using data from the RXTE and Swift observatories we detected several consecutive outbursts, during which the source luminosity increased by a factor of ~5 up to L_x~1.2x10^{35} erg/s. Previously the source was observed in a high state only once. The source spectrum in a standard energy band (4-25 keV) is independent of the flux change and can be described by a model that includes both thermal and non-thermal components. With the help of the INTEGRAL observatory data we registered the highly significant cyclotron absorption line in the source spectrum and, for the first time, significantly detected a hard X-ray emission from the pulsar up to ~160 keV. We also report drastic changes of the pulse period during the outburst activity: a long episode of spin-down was changed to spin-up with a rate of Pdot/P = -(3-5)x10^{-4} yr^{-1}, that is several times higher than previous rates of spin-up and spin-down. To search for a correlation between the X-ray and optical lightcurves we took data from the AAVSO International Database. No significant correlation between optical and X-ray fluxes at any time lag from dozens of days to years was found.

We present results of a study of the Galactic ridge X-ray emission (GRXE) in hard X-rays with the IBIS telescope on board INTEGRAL in the region near the Galactic Anticenter (GA) at l=155 deg. We assumed a conservative 2 sigma upper limit on the flux from the GA in the 25-60 keV energy band of 1.25E-10 erg/s/cm^2 (12.8 mCrab) per IBIS field of view, or 6.6E-12 erg/s/cm^2 (0.7 mCrab) per degree longitude in the 135 deg. < l < 175 deg. region. This upper limit exceeds the expected GRXE intensity in the GA direction by an order of magnitude, given the near-infrared (NIR) surface brightness of the Milky Way in this region and the standard hard X-ray-to-NIR intensity ratio for the GRXE, assuming stellar origin. Based on the CGRO/EGRET surface brightness of the Galaxy above 100 MeV as a tracer of the cosmic-ray (CR) induced gamma-ray background, the expected GRXE flux in GA exceeds the measured 2 sigma upper limit by a factor of 4. Therefore, the non-detection of hard X-ray emission from the GA does not contradict the stellar nature of the GRXE, but is inconsistent with CR origin.

We report the results of an analysis of data obtained with the INTEGRAL, Swift and RXTE observatories during the 2010 April and September outbursts of the X-ray pulsar RX J0440.9+4431. The temporal and spectral properties of the pulsar in a wide energy band (0.6-120 keV) were studied for the first time. We discovered a ~32 keV cyclotron resonant scattering feature in the source spectrum, that allowed us to estimate the magnetic field strength of the neutron star as B~3.2 x 10^12 G. The estimate of the magnetic field strength was confirmed by a comprehensive analysis of the noise power spectrum of the source. Based on the recurrence time between Type I outbursts the orbital period of the binary system can be estimated as ~155 days. We have shown that the pulse profile has a sinusoidal-like single-peaked shape and has practically no dependence on the source luminosity or energy band.

We analyse SPI/INTEGRAL data on the 511 keV line from the Galactic Centre, accumulated over ~6 years of observations. We decompose the X-ray and soft gamma-ray emission of the central part of the Milky Way into a relatively compact "Bulge" and a more extended "Disk" components and report their spectral properties. The Bulge component shows a prominent 511 keV line and essentially no flux at 1.8 MeV, while the Disk component on the contrary contains a prominent 1.8 MeV line and a very weak annihilation line. We show that the spectral shape of the annihilation radiation (the narrow 511 keV line and the associated othro-positronium continuum) is surprisingly well described by a model of annihilation of hot positrons in a radiatively cooling interstellar medium (ISM). The model assumes that positrons are initially injected into a hot ($\sim 10^6$~K), volume filling ISM, which is allowed to freely cool via radiative losses. The annihilation time in such a medium is longer than the cooling time for temperatures higher than a few $10^{4}$~K. Thus, most of the positrons annihilate only after the gas has cooled down to $\sim 10^5$~K, giving rise to annihilation emission characteristic of a warm, ionized ISM.

This paper is the second in a series devoted to the hard X-ray (17-60 keV) whole sky survey performed by the INTEGRAL observatory over seven years. Here we present a catalog of detected sources which includes 521 objects, 449 of which exceed a 5 sigma detection threshold on the time-averaged map of the sky, and 53 were detected in various subsamples of exposures. Among the identified sources with known and suspected nature, 262 are Galactic (101 low-mass X-ray binaries, 95 high-mass X-ray binaries, 36 cataclysmic variables, and 30 of other types) and 219 are extragalactic, including 214 active galactic nuclei (AGNs), 4 galaxy clusters, and galaxy ESO 389-G 002. The extragalactic (|b|>5 deg) and Galactic (|b|<5 deg) persistently detected source samples are of high identification completeness (respectively ~96% and ~94%) and valuable for population studies.

This paper is the first in a series devoted to the hard X-ray whole sky survey performed by the INTEGRAL observatory over seven years. Here we present an improved method for image reconstruction with the IBIS coded mask telescope. The main improvements are related to the suppression of systematic effects which strongly limit sensitivity in the region of the Galactic Plane (GP), especially in the crowded field of the Galactic Center (GC). We extended the IBIS/ISGRI background model to take into account the Galactic Ridge X-ray Emission (GRXE). To suppress residual systematic artifacts on a reconstructed sky image we applied nonparametric sky image filtering based on wavelet decomposition. The implemented modifications of the sky reconstruction method decrease the systematic noise in the ~20 Ms deep field of GC by ~44%, and practically remove it from the high-latitude sky images. New observational data sets, along with an improved reconstruction algorithm, allow us to conduct the hard X-ray survey with the best currently available minimal sensitivity 3.7E-12 erg/s/cm2 ~0.26 mCrab in the 17-60 keV band at a 5 sigma detection level. The survey covers 90% of the sky down to the flux limit of 6.2E-11 erg/s/cm2 (~4.32 mCrab) and 10% of the sky area down to the flux limit of 8.6E-12 erg/s/cm2 (~0.60 mCrab).

We present the results of a systematic search for outbursts in the narrow positron annihilation line on various time scales (5x10^4 - 10^6 s) based on the SPI/INTEGRAL data obtained from 2003 to 2008. We show that no outbursts were detected with a statistical significance higher than ~6 sigma for any of the time scales considered over the entire period of observations. We also show that, given the large number of independent trials, all of the observed spikes could be associated with purely statistical flux fluctuations and, in part, with a small systematic prediction error of the telescope's instrumental background. Based on the exposure achieved in ~6 yr of INTEGRAL operation, we provide conservative upper limits on the rate of outbursts with a given duration and flux in different parts of the sky.

We present here the preliminary results of the study of the fluorescent iron line emission from X-ray pulsars with Be companions. We propose to use properties of this emission to investigate the spatial distribution and physical conditions of the matter around the compact object as well as in the binary system as a whole. Using data of the RXTE observatory the iron line behavior in the transient X-ray pulsar V 0332+53 spectrum was studied during the powerful type II outburst in 2004-2005. Particularly, we investigated a variability of the iron line equivalent width on different time scales (pulse period, orbital period, outburst phase) and searched for its correlation with the continuum flux, spectral parameters, etc.

Analysis of the data obtained with the RXTE observatory during a powerful outburst of the X-ray pulsar V0332+53 in 2004-2005 is presented. Observational data covering the outburst brightening phase are analysed in detail for the first time. A comparison of source parameters and their evolution during the brightening and fading phases shows no evidence for any hysteresis behaviour. It is found that the dependences of the energy of the cyclotron absorption line on the luminosity during the brightening and fading phases are almost identical. The complete data sequence including the outburst brightening and fading phases makes it possible to impose the more stringent constraints on the magnetic field in the source. The pulse profile and pulsed fraction are studied as functions of the luminosity and photon energy.

We study power density spectra (PDS) of X-ray flux variability in binary systems where the accretion flow is truncated by the magnetosphere. PDS of accreting X-ray pulsars where the neutron star is close to the corotation with the accretion disk at the magnetospheric boundary, have a distinct break/cutoff at the neutron star spin frequency. This break can naturally be explained in the "perturbation propagation" model, which assumes that at any given radius in the accretion disk stochastic perturbations are introduced to the flow with frequencies characteristic for this radius. These perturbations are then advected to the region of main energy release leading to a self-similar variability of X-ray flux P~f^{-1...-1.5}. The break in the PDS is then a natural manifestation of the transition from the disk to magnetospheric flow at the frequency characteristic for the accretion disk truncation radius (magnetospheric radius). The proximity of the PDS break frequency to the spin frequency in corotating pulsars strongly suggests that the typical variability time scale in accretion disks is close to the Keplerian one. In transient accreting X-ray pulsars characterized by large variations of the mass accretion rate during outbursts, the PDS break frequency follows the variations of the X-ray flux, reflecting the change of the magnetosphere size with the accretion rate. Above the break frequency the PDS steepens to ~f^{-2} law which holds over a broad frequency range. These results suggest that strong f^{-1...-1.5} aperiodic variability which is ubiquitous in accretion disks is not characteristic for magnetospheric flows.

We review the results of a timing analysis of the observations for ten bright X-ray pulsars (with fluxes >100 mCrab in the 20-100 keV energy band) that fell within the INTEGRAL field of view from 2003 to 2007. The dependence of the pulse profile on the energy and intrinsic source luminosity has been investigated; particular attention has been paid to searching for changes in the pulse profile near the cyclotron frequency. The dependence of the pulsed fraction for X-ray pulsars on their luminosity and energy band has been studied in detail for the first time.

Recent results of the spectral and timing analysis of X-ray pulsars in hard X-rays with the INTEGRAL observatory are reviewed. The evolution of the cyclotron line energy with the source luminosity was studied in detail for the first time for several sources. It was shown that for V0332+53 this dependence is linear, but for 4U0115+63 and A0535+262 it is more complicated. There are some evidences of the "reverse" evolution for GX301-2 and Her X-1, and no evolution was found for Vela X-1, Cen X-3, etc. A strong dependence of the pulse fraction on the energy and source luminosity was revealed and studied in detail. A prominent feature in the pulse fraction dependence on the energy was discovered near the cyclotron frequency for several bright sources. The obtained results are compared with results of observations in standard X-rays and briefly discussed in terms of current models; some preliminary explanations are proposed.

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.

We are reporting a discovery of X-ray pulsations from the source AX J1749.1-2733 with the period of ~132 s based on the XMM-Newton data obtained in March 2007. The observed pulse profile has a double-peaked structure with the pulse fraction of about 25-30 % in the 3-10 keV energy band. We have also found that a periodicity with practically the same period has been detected from the source by the IBIS telescope onboard the INTEGRAL observatory during an outburst on Sept. 9, 2003 in the 20-50 keV energy band. Due to the double-peaked pulse profile, there is an additional peak on both periodograms of nearly ~66 s, therefore we have also investigated the possibility that the last value is the true pulse period. The source spectrum obtained by the XMM-Newton observatory in the soft energy band is being heavily absorbed ($N_H\simeq2\times10^{23}$ cm$^{-2}$) due to a strong intrinsic absorption in the binary system that leads to the conclusion that AX J1749.1-2733 is a new transient X-ray pulsar in the high mass X-ray binary system.

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 analyze the observations of the transient X-ray pulsar 4U 0115+63 with the RXTE and INTEGRAL observatories in a wide X-ray (3-100 keV) energy band during its intense outbursts in 1999 and 2004. The energy of the fundamental harmonic of the cyclotron resonance absorption line near the maximum of the X-ray flux from the source (luminosity range 5x10^{37} - 2x10^{38} erg/s) is ~11 keV. When the pulsar luminosity falls below ~5x10^{37} erg/s, the energy of the fundamental harmonic is displaced sharply toward the high energies, up to ~16 keV. Under the assumption of a dipole magnetic field configuration, this change in cyclotron harmonic energy corresponds to a decrease in the height of the emitting region by ~2 km, while other spectral parameters, in particular, the cutoff energy, remain essentially constant. At a luminosity ~7x10^{37} erg/s, four almost equidistant cyclotron line harmonics are clearly seen in the spectrum. This suggests that either the region where the emission originates is compact or the emergent spectrum from different (in height) segments of the accretion column is uniform. We have found significant pulse profile variations with energy, luminosity, and time. In particular, we show that the profile variations from pulse to pulse are not reduced to a simple modulation of the accretion rate specified by external conditions.

We studied spectra for 34 accretion-powered X-ray and one millisecond pulsars that were within the field of view of the INTEGRAL observatory over two years (December 2002 - January 2005) of its in-orbit operation and that were detected by its instruments at a statistically significant level (>8 sigma in the energy range 18-60 keV). There are seven recently discovered objects of this class among the pulsars studied: 2RXP J130159.6-635806, IGR/AX J16320-4751, IGR J16358-4726, AX J163904-4642, IGR J16465-4507, SAX/IGR J18027-2017 and AX J1841.0-0535. We analyze the evolution of spectral parameters as a function of the intensity of the sources and compare these with the results of previous studies.

We present the results of broad band (3-100 keV) observations of several X-ray pulsars with the INTEGRAL and RXTE observatories. We concentrate on the luminosity and energy dependence of the pulse profile and the variations of the cyclotron line energy. In V0332+53 the line energy changes nearly linearly with the source luminosity, while in 4U0115+63 its behavior is more complicated. Strong variations of the pulse profile with the energy and source intensity were found for both of pulsars; in V0332+53 the changes of the pulse profile near the cyclotron line are especially drastic. The preliminary results obtained for Her X-1 and GX 301-2 in a high intensity state show the absence of significant pulse profile changes with the energy. Results and possible emission mechanisms are briefly discussed in terms of theoretical models of accreting pulsars.

We present results of observations of the transient X-ray pulsar V0332+53 performed during a very powerful outburst in Dec, 2004 -- Feb, 2005 with the INTEGRAL and RXTE observatories in a wide (3-100 keV) energy band. A cyclotron resonance scattering line at an energy of ~26 keV has been detected in the source spectrum together with its two higher harmonics at ~50 and ~73 keV, respectively. We show that the energy of the line is not constant but linearly changes with the source luminosity. Strong pulse profile variations, especially near the cyclotron line, are revealed for different levels of the source intensity. We discuss the obtained results in terms of the theoretical models of X-ray pulsars.

We present spectra for 34 accretion-powered X-ray and one millisecond pulsars that were within the field of view of the INTEGRAL observatory over two years (December 2002 - January 2005) of its in-orbit operation and that were detected by its instruments at a statistically significant level (> 8 sigma in the energy range 18--60 keV). There are seven recently discovered objects of this class among the pulsars studied: 2RXP J130159.6-635806, IGR/AX J16320-4751, IGR J16358-4726, AX J163904-4642, IGR J16465-4507, SAX/IGR J18027-2017 and AX J1841.0-0535. We have also obtained hard X-ray (> 20 keV) spectra for the accretion-powered pulsars A 0114+650, RX J0146.9+6121, AX J1820.5-1434, AX J1841.0-0535 and the millisecond pulsar XTE J1807-294 for the first time. We analyze the evolution of spectral parameters as a function of the intensity of the sources and compare these with the results of previous studies.

We analyze the observations of the X-ray pulsar LMC X-4 performed by the INTEGRAL observatory and the All-Sky Monitor (ASM) of the RXTE observatory over a wide energy range. The observed hard X-ray flux from the source is shown to change by more than a factor of 50 (from ~70 mCrab in the high state to ~1.3 mCrab in the low state) on the time scale of the accretion-disk precession period, whose mean value for 1996-2004 was determined with a high accuracy, Pprec=30.275+/-0.004 days. In the low state, a flare about 10 h in duration was detected from the source; the flux from the source increased by more than a factor of 4 during this flare. The shape of the pulsar's broadband spectrum is essentially invariable with its intensity; no statistically significant features associated with the possible resonance cyclotron absorption line were found in the spectrum of the source.

We analyze the observations of the X-ray pulsar KS 1947+300 performed by the INTEGRAL and RXTE observatories over a wide (3-100 keV) X-ray energy range. The shape of the pulse profile was found to depend on the luminosity of the source. Based on the model of a magnetized neutron star, we study the characteristics of the pulsar using the change in its spin-up rate. We estimated the magnetic field strength of the pulsar and the distance to the binary.

We present the first results of the observations of the X-ray pulsars LMC X-4, 4U0352+309 and EXO1722-363 performed with the INTEGRAL observatory. The LMC X-4 was investigated during the whole superorbital cycle ($\sim$30 days) and it was found that its period was not stable at this time scale. We detected a variable X-ray flux (18-60 keV) from the pulsar EXO1722-363, which could be connected with the orbital motion in the binary system. A more accurate position and the estimate of the orbital period for this source are reported. We also investigated a hard X-ray spectrum of 4U0352+309 (X Persei) measured with INTEGRAL and report the detection of the cyclotron absorption line at about 29 keV.

The variability of the X-ray flux from the pulsar GX 301-2 is analyzed by using data from the ART-P telescope of the Granat Observatory. The intensity variations with time scales of several thousand seconds are studied at various orbital phases. The high-state flux from the source exceeds its low-state flux by as much as a factor of 10. The hardness and spectrum of the source are shown to change greatly with its intensity. These intensity variations are most likely caused by substantial inhomogeneities in the stellar wind from the companion star.

We present the observations of the pulsar SMC X-1 with the ART-P telescope onboard the Granat observatory. We investigate the variability of the flux from the source on time scales of several tens of days. The intensity variation of the pulsar are shown to be consistent with the presence of a periodicity in the system with a characteristic time scale of $\sim61$ days. The precession of an inclined accretion disk, as indirectly confirmed by the absence of low-state pulsations, may be responsible for the observed variability. The spectrum of the source is well described by a power-law energy dependence of the photon flux density with a slope of $\sim1.5$ and an exponential cutoff at energies above $\sim14-18$ keV. We estimated the inclinations between the planes of the orbit and the accretion disk and the magnetic field of the neutron star.