We report the discovery of a nearby massive white dwarf with He-H atmosphere. The white dwarf is located at a distance of $74.5\pm0.9$ pc. Its radius, mass, effective temperature, H/He ratio and age are $R=2500\pm100$ km, $M=1.33\pm0.01$ $\rm M_{\odot}$, $T_{\mathrm{eff}}=31200\pm 1200$ K, ${\mathrm{H/He}}\sim0.1$ and $330\pm40$ Myr, respectively. The observed spectrum is redshifted by $V_{\mathrm{r}}=+240\pm15$ km s$^{-1},$ which is mostly attributed to the gravitational redshift. The white dwarf shows a regular stable photometric variability with amplitude $\Delta g\approx0.06^{\mathrm{m}}$ and period $P=353.456$ s suggesting rapid rotation. This massive, hot and rapidly rotating white dwarf is likely to originate from the merging of close binary white dwarf system that avoided explosion in a thermonuclear type Ia supernova at the Carboniferous Period of the Earth history.

Seventeen years of hard X-ray observations with the instruments of the INTEGRAL observatory, with a focus on the Milky Way and in particular on the Galactic Centre region, have provided a unique database for exploration of the Galactic population of low-mass X-ray binaries (LMXBs). Our understanding of the diverse energetic phenomena associated with accretion of matter onto neutron stars and black holes has greatly improved. We review the large variety of INTEGRAL based results related to LMXBs. In particular, we discuss the spatial distribution of LMXBs over the Galaxy and their X-ray luminosity function as well as various physical phenomena associated with Atoll and Z sources, bursters, symbiotic X-ray binaries, ultracompact X-ray binaries and persistent black hole LMXBs. We also present an up-to-date catalogue of confirmed LMXBs detected by INTEGRAL, which comprises 166 objects. Last but not least, the long-term monitoring of the Galactic Centre with INTEGRAL has shed light on the activity of Sgr A* in the recent past, confirming previous indications that our supermassive black hole experienced a major accretion episode just ~100 years ago. This exciting topic is covered in this review too.

It is shown that a mechanism of PBH formation from high-baryon bubbles with log-normal mass spectrum naturally leads to the central mass of the PBH distribution close to ten solar masses independently of the model details. This result is in good agreement with observations.

Possible models for the generation of electromagnetic (EM) radiation during the coalescence of neutron star-black hole binaries are considered. The mass of the remnant disk around the black hole during the coalescence of neutron stars and black holes is calculated by taking into account the equation of state for neutron stars and the rotation of the binary components before the coalescence. The parameters of binary systems before the coalescence (the mass ratio, the component rotation, the neutron star magnetic field) are calculated by the population synthesis method. The derived mass of the remnant disk around the black hole after the coalescence is used to estimate the kinetic energy of the relativistic jet launched by the Blandford-Znajek mechanism. A disk mass of more than $\sim 0.05 M_\odot$ required for the formation of short gamma-ray bursts is shown to be obtained in no more than 1-10\% of the coalescences (depending on the equation of state). Less efficient common envelopes (a large parameter $\alpha_{CE}$) lead to a noticeably larger percentage of events with astrophysically interesting EM energy release. For binaries with a large mass ratio, in which a magnetized neutron star is not subjected to tidal disruption before the coalescence, the possibility of the formation of an electrically charged rotating black hole (Wald charge) is considered and estimates of the maximum EM power released by such a black hole after the coalescence are made. The conversion of the emitted gravitational waves into electromagnetic ones in the relativistic lepton plasma generated in coalescing pulsar-black hole binaries at the pre-coalescence stage is discussed.

The available data on the chirp mass distribution of the black holes in the coalescing binaries in O1-O3 LIGO/Virgo runs are analyzed and compared with theoretical expectations based on the hypothesis that these black holes are primordial with log-normal mass spectrum. The inferred best-fit mass spectrum parameters, $M_0=17 M_\odot$ and $\gamma=0.9$, fall within the theoretically expected range and shows excellent agreement with observations. On the opposite, binary black hole models based on massive binary star evolution require additional adjustments to reproduce the observed chirp mass distribution.

We model the Galactic population of detached binaries that harbor black holes with (0.5-1.7) solar mass companions -- the remnants of case B mass-exchange that rapidly cross Hertzsprung gap after the termination of the Roche-lobe overflow or as He-shell burning stars. Several such binaries can be currently present in the Galaxy. The range of black hole masses in them is about 4 to 10 solar ones, the orbital periods are tens to hundreds day. The unique BH-binary LB-1 fits well into this extremely rare class of double stars.

We report on the results of a NuSTAR observation of the Supergiant Fast X-ray Transient pulsar IGRJ11215-5952 during the peak of its outburst in June 2017. IGRJ11215-5952 is the only SFXT undergoing strictly periodic outbursts, every 165 days. NuSTAR caught several X-ray flares, spanning a dynamic range of 100, and detected X-ray pulsations at 187.0 s, consistent with previous measurements. The spectrum from the whole observation is well described by an absorbed power-law (with a photon index of 1.4) modified, above 7 keV, by a cutoff with an e-folding energy of 24 keV. A weak emission line is present at 6.4 keV, consistent with Kalpha emission from cold iron in the supergiant wind. The time-averaged flux is 1.5E-10 erg/cm2/s (3-78 keV, corrected for the absorption), translating into an average luminosity of about 9E35 erg/s (1-100 keV, assuming a distance of 6.5 kpc). The NuSTAR observation allowed us to perform the most sensitive search for cyclotron resonant scattering features in the hard X-ray spectrum, resulting in no significant detection in any of the different spectral extractions adopted (time-averaged, temporally-selected, spin-phase-resolved and intensity-selected spectra). The pulse profile showed an evolution with both the energy (3-12 keV energy range compared with 12-78 keV band) and the X-ray flux: a double peaked profile was evident at higher fluxes (and in both energy bands), while a single peaked, sinusoidal profile was present at the lowest intensity state achieved within the NuSTAR observations (in both energy bands). The intensity-selected analysis allowed us to observe an anti-correlation of the pulsed fraction with the X-ray luminosity. The pulse profile evolution can be explained by X-ray photon scattering in the accreting matter above magnetic poles of a neutron star at the quasi-spherical settling accretion stage.

It is widely accepted that dark matter contributes about a quarter of the critical mass-energy density in our Universe. The nature of dark matter is currently unknown, with the mass of possible constituents spanning nearly one hundred orders of magnitude. The ultralight scalar field dark matter, consisting of extremely light bosons with $m \sim 10^{-22}$ eV and often called "fuzzy" dark matter, provides intriguing solutions to some challenges at sub-Galactic scales for the standard cold dark matter model. As shown by Khmelnitsky and Rubakov, such a scalar field in the Galaxy would produce an oscillating gravitational potential with nanohertz frequencies, resulting in periodic variations in the times of arrival of radio pulses from pulsars. The Parkes Pulsar Timing Array (PPTA) has been monitoring 20 millisecond pulsars at two to three weeks intervals for more than a decade. In addition to the detection of nanohertz gravitational waves, PPTA offers the opportunity for direct searches for fuzzy dark matter in an astrophysically feasible range of masses. We analyze the latest PPTA data set which includes timing observations for 26 pulsars made between 2004 and 2016. We perform a search in this data set for evidence of ultralight dark matter in the Galaxy using Bayesian and Frequentist methods. No statistically significant detection has been made. We therefore place upper limits on the local dark matter density. Our limits, improving on previous searches by a factor of two to five, constrain the dark matter density of ultralight bosons with $m \leq 10^{-23}$ eV to be below $6\,\text{GeV}\,\text{cm}^{-3}$ with 95\% confidence in the Earth neighborhood. Finally, we discuss the prospect of probing the astrophysically favored mass range $m \gtrsim 10^{-22}$ eV with next-generation pulsar timing facilities.

Modern astrophysical methods of determination of spins of rotating stellar-mass black hole in close binaries and of supermassive black holes in active galactic nuclei are briefly discussed. Effective spins of coalescing binary black holes derived from LIGO/Virgo gravitational wave observations are specially addressed. We consider three types of coalescing binaries: double black holes, black hole-neutron star binaries and primordial double black holes. The effective spins of coalescing astrophysical binary black holes and black holes with neutron stars are calculated for two plausible models of black hole formations from stellar core collapses (without or with additional fallback from the stellar envelope) taking into account the stellar metallicity and star formation rate evolution in the Universe. The calculated distributions do not contradict the reported LIGO/Virgo observations. The effective spins of primordial coalescing stellar-mass black holes can reach a few per cent due to accretion spin-up in the cold external medium.

Primordial stellar-mass black holes, which may contribute to dark matter and to the observed LIGO binary black hole coalescences, are expected to be born with very low spins. Here we show that accretion mass gain by the components of a primordial black hole binary from the surrounding matter could lead to noticeable spins of the components prior to the coalescence provided high initial orbital eccentricities.

INTEGRAL IBIS/ISGRI 18-60 keV observations of SS433 performed in 2003-2011 enabled the hard X-ray phase-resolved orbital and precessional light curves and spectra to be constructed. The spectra can be fitted by a power-law with photon index $\simeq 3.8$ and remain almost constant while the X-ray flux varies by a factor of a few. This suggests that the hard X-ray emission is produced in an extended quasi-isothermal hot 'corona' surrounding central parts of a supercritical accretion disc. A joint analysis of the broadband 18-60 keV orbital and precessional light curves was performed in the model assuming a significant Roche lobe overfilling by the optical star, up to its filling the outer Lagrangian surface enabling mass loss through the outer Lagrangian L$_2$ point. From this modeling, the relativistic-to-optical component mass ratio $q=M_x/M_v\gtrsim0.4\div 0.8$ is estimated. An analysis of the observed long-term stability of the orbital period of SS433 with an account of the recent observations of SS433 by the VLTI GRAVITY interferometer enabled an independent mass ratio estimate $q>0.6$. This estimate in combination with the radial velocity semi-amplitude for stationary He II emission, $K_x=168\pm 18$ km/s (Hillwig et al 2004) suggests the optical component mass in SS433 $M_v>12 M_\odot$. Thus, the mass of the relativistic component in SS433 is $M_x>7 M_\odot$, which is close to the mean mass of black holes in X-ray binaries ($\sim 8 M_\odot$). The large binary mass ratio in SS433 allows us to understand why there is no common envelope in this binary at the secondary mass transfer evolutionary stage and the system remains semi-detached (van den Heuvel et al. 2017). We also discuss unsolved issues and outline prospects for further study of SS433.

The low luminosity, X-ray flaring activity, of the sub-class of high mass X-ray binaries called Supergiant Fast X-ray Transients, has been investigated using XMM-Newton public observations, taking advantage of the products made publicly available by the EXTraS project. One of the goals of EXTraS was to extract from the XMM-Newton public archive information on the aperiodic variability of all sources observed in the soft X-ray range with EPIC (0.2-12 keV). Adopting a Bayesian block decomposition of the X-ray light curves of a sample of SFXTs, we picked out 144 X-ray flares, covering a large range of soft X-ray luminosities (1e32-1e36 erg/s). We measured temporal quantities, like the rise time to and the decay time from the peak of the flares, their duration and the time interval between adjacent flares. We also estimated the peak luminosity, average accretion rate and energy release in the flares. The observed soft X-ray properties of low-luminosity flaring activity from SFXTs is in qualitative agreement with what is expected by the application of the Rayleigh-Taylor instability model in accreting plasma near the neutron star magnetosphere. In the case of rapidly rotating neutron stars, sporadic accretion from temporary discs cannot be excluded.

We revisit the determination of binary mass ratio in the Galactic microquasar SS433 based on recent GRAVITY VLTI measurements of mass and angular momentum outflow through a circumbinary disc. The new observations combined with the constancy of the binary orbital period over $\sim 30$ yrs confirm that the mass ratio in SS433 is $q=M_\mathrm{x}/M_\mathrm{v}\gtrsim 0.6$. For the assumed optical star mass $M_\mathrm{v}$ ranging from $\sim 8$ to 15 $M_\odot$ such a mass ratio suggests a low limit of the compact object mass of $M_\mathrm{x}\sim 5-9 M_\odot$, placing the compact object in SS433 as a stellar-mass black hole.

We present a new model of the population of symbiotic X-ray binaries (SyXBs) that takes into account non-stationary character of quasi-spherical sub-sonic accretion of the red giant's stellar wind onto slowly rotating neutron stars. Updates of the earlier models are given, which include more strict criteria of slow NS rotation for plasma entry into the NS magnetosphere via Rayleigh-Taylor instability, as well as more strict conditions for settling accretion for slow stellar winds, with an account of variations in the specific angular momentum of captured stellar wind in eccentric binaries. These modifications enabled a more adequate description of the distributions of observed systems over binary orbital periods, NS spin periods and their X-ray luminosity in the $\sim 10^{32}-10^{36}$~erg s$^{-1}$ range and brought their model Galactic number into reasonable agreement with the observed one. Reconciliation of the model and observed orbital periods of SyXBs requires a low efficiency of matter expulsion from common envelopes during the evolution that results in the formation of NS-components of symbiotic X-ray systems.

Cyclotron lines, also called cyclotron resonant scattering features (CRSF) are spectral features, generally appearing in absorption, in the X-ray spectra of objects containing highly magnetized neutron stars, allowing the direct measurement of the magnetic field strength in these objects. Cyclotron features are thought to be due to resonant scattering of photons by electrons in the strong magnetic fields. The main content of this contribution focusses on electron cyclotron lines as found in accreting X-ray binary pulsars (XRBP) with magnetic fields on the order of several 1012 Gauss. Also, possible proton cyclotron lines from single neutron stars with even stronger magnetic fields are briefly discussed.With regard to electron cyclotron lines, we present an updated list of XRBPs that show evidence of such absorption lines. The first such line was discovered in a 1976 balloon observation of the accreting binary pulsar Hercules X-1, it is considered to be the first direct measurement of the magnetic field of a neutron star. As of today (mid 2018), we list 36 XRBPs showing evidence of one ore more electron cyclotron absorption line(s). A few have been measured only once and must be confirmed (several more objects are listed as candidates). In addition to the Tables of objects, we summarize the evidence of variability of the cyclotron line as a function of various parameters (especially pulse phase, luminosity and time), and add a discussion of the different observed phenomena and associated attempts of theoretical modeling. We also discuss our understanding of the underlying physics of accretion onto highly magnetized neutron stars. For proton cyclotron lines, we present tables with seven neutron stars and discuss their nature and the physics in these objects.

A fraction of high-mass X-ray binaries are supergiant fast X-ray transients. These systems have on average low X-ray luminosities, but display short flares during which their X-ray luminosity rises by a few orders of magnitude. The leading model for the physics governing this X-ray behaviour suggests that the winds of the donor OB supergiants are magnetized. In agreement with this model, the first spectropolarimetric observations of the SFXT IGR J11215-5952 using the FORS2 instrument at the Very Large Telescope indicate the presence of a kG longitudinal magnetic field. Based on these results, it seems possible that the key difference between supergiant fast X-ray transients and other high-mass X-ray binaries are the properties of the supergiant's stellar wind and the physics of the wind's interaction with the neutron star magnetosphere.

The possible formation mechanisms of massive close binary black holes (BHs) that can merge in the Hubble time to produce powerful gravitational wave bursts detected during advanced LIGO O1 and O2 science runs include the evolution from field low-metallicity massive binaries, the dynamical formation in dense stellar clusters and primordial BHs. Different formation channels produce different source distributions of total masses ${M_\mathrm{tot}}$ and effective spins $\chi_\mathrm{eff}$ of coalescing binary BHs. Using a modified \textsc{bse} code, we carry out extensive population synthesis calculations of the expected effective spin and total mass distributions from the standard field massive binary formation channel for different metallicities of BH progenitors (from zero-metal Population III stars up to solar metal abundance), different initial rotations of the binary components, stellar wind mass loss prescription, different BH formation models and a range of common envelope efficiencies. The stellar rotation is treated in two-zone (core-envelope) approximation using the effective core-envelope coupling time and with an account of the tidal synchronization of stellar envelope rotation during the binary system evolution. The results of our simulations, convolved with the metallicity-dependent star-formation history, show that the total masses and effective spins of the merging binary black holes detected during LIGO O1-O2 runs but the heaviest one (GW170729) can be simultaneously reproduced by the adopted BH formation models. Noticeable effective spin of GW170729 requires additional fallback from the rotating stellar envelope.

Different accretion regimes onto magnetized NSs in HMXBs are considered: wind-fed supersonic (Bondi) regime at high accretion rates $\dot M\gtrsim 4\times 10^{16}$~[g s$^{-1}$] , subsonic settling regime at lower $\dot M$ and supercritical disc accretion during Roche lobe overflow. In wind-fed stage, NSs in HMXBs reach equilibrium spin periods $P^*$ proportional to binary orbital period $P_b$. At supercritical accretion stage, the system may appear as a pulsating ULX. Population synthesis of Galactic HMXBs using standard assumptions on the binary evolution and NS formation is presented. Comparison of the model $P^*$ -- $P_b$ (the Corbet diagram), $P^*$ -- $L_x$ and $P_b$ -- $L_x$ distributions with those for the observed HMXBs (including Be X-ray binaries) and pulsating ULXs suggests the importance of the reduction of $P^*$ in non-circular orbits, explaining the location of Be X-ray binaries in the model Corbet diagram, and the universal parameters of pulsating ULXs depending only on the NS magnetic fields.

A huge optical luminosity of the supercritical accretion disc and powerful stellar wind in the high-mass X-ray binary SS433 make it difficult to reliably estimate the mass ratio of the binary components from spectroscopic observations. We analyze different indirect methods of the mass ratio estimate. We show that with an account of the possible Roche lobe overflow by the optical star, the analysis of X-ray eclipses in the standard and hard X-ray bands suggests the estimate $q=M_\mathrm{x}/M_\mathrm{v}\gtrsim 0.3$. We argue that the double-peak hydrogen Brackett lines in SS433 should form not in the accretion disc but in a circumbinary envelope, suggesting a total mass of $M_\mathrm{v}+M_\mathrm{x}\gtrsim 40 M_\odot$. The observed long-term stability of the orbital period in SS433 $|\dot P_b/P_b|\le 1.793\times 10^{-14}$~s$^{-1}$ over $\sim 28$ year period is used to place an independent constraint of $q\gtrsim 0.6$ in SS433, confirming its being a Galactic microquasar hosting a superaccreting black hole.

Based on our new NuSTAR X-ray telescope data, we rule out any cyclotron line up to 60 keV in the spectra of the high mass X-ray binary 4U2206+54. In particular, we do not find any evidence of the previously claimed line around 30 keV, independently of the source flux, along the spin pulse. The spin period has increased significantly, since the last observation, up to $5750\pm10$ s, confirming the rapid spin down rate $\dot \nu=-1.8\times 10^{-14}$~Hz s$^{-1}$. This behaviour might be explained by the presence of a strongly magnetized neutron star ($B_{\rm s}>$ several times $10^{13}$ G) accreting from the slow wind of its main sequence O9.5 companion.

The phenomenon of fast radio bursts (FRBs) was discovered in 2007. These are powerful (0.1-100 Jy) single radio pulses with durations of several milliseconds, large dispersion measures, and record high brightness temperatures suggesting coherent emission mechanism. As of time of writing, 32 FRBs were recorded. There is also one repeating source from which already hundreds of bursts have been detected. The FRB rate is estimated to amount up to several thousand per day over the sky, and their isotropic sky distribution likely suggests a cosmological origin. Since the discovery, different hypotheses on the possible FRB nature have been proposed, however up to now the origin of these transient events remains obscure. The most prospective models treat them as being related to a bursting emission from magnetars -- neutron stars whose activity is due to their magnetic field dissipation, -- or as being analogs of giant pulses observed from several radio pulsars . Future increase in the statistics of the observed FRBs and an improvement upon the FRB population characteristics will make possible to use them as a new tool to probe the intergalactic medium and to test fundamental physical theories.

We briefly review main observational properties of fast radio bursts (FRBs) and discuss two most popular hypothesis for the explanation of these enigmatic intense millisecond radio flashes. FRBs most probably originate on extragalactic distances, and their rate on the sky is about a few thousand per day with fluences above $\sim$~1~Jy~ms (or with fluxes larger than few tenths of Jy). Two leading scenarios describing these events include strong flares of magnetars and supergiant pulses of young radio pulsars with large rotational energy losses, correspondingly. At the moment, it is impossible to choose between these models. However, new telescopes can help to solve the puzzle of FRBs in near future.

A significant fraction of high-mass X-ray binaries are supergiant fast X-ray transients (SFXTs). The prime model for the physics governing their X-ray behaviour suggests that the winds of donor OB supergiants are magnetized. To investigate if magnetic fields are indeed present in the optical counterparts of such systems, we acquired low-resolution spectropolarimetric observations of the two optically brightest SFXTs, IGR J08408-4503 and IGR J11215-5952 with the ESO FORS2 instrument during two different observing runs. No field detection at a significance level of 3sigma was achieved for IGR J08408-4503. For IGR J11215-5952, we obtain 3.2sigma and 3.8sigma detections (<B_z>_hydr = -978+-308G and <B_z>_hydr = 416+-110G) on two different nights in 2016. These results indicate that the model involving the interaction of a magnetized stellar wind with the neutron star magnetosphere can indeed be considered to characterize the behaviour of SFXTs. We detected long-term spectral variability in IGR J11215-5952, while for IGR J08408-4503 we find an indication of the presence of short-term variability on a time scale of minutes.

Based on new Chandra X-ray telescope data, we present empirical evidence of plasma Compton cooling during a flare in the non pulsating massive X-ray binary 4U1700-37. This behaviour might be explained by quasispherical accretion onto a slowly rotating magnetised neutron star. In quiescence, the neutron star in 4U1700-37 is surrounded by a hot radiatively cooling shell. Its presence is supported by the detection of mHz quasi periodic oscillations likely produced by its convection cells. The high plasma temperature and the relatively low X-ray luminosity observed during the quiescence, point to a small emitting area about 1 km, compatible with a hot spot on a NS surface. The sudden transition from a radiative to a significantly more efficient Compton cooling regime triggers an episode of enhanced accretion resulting in a flare. During the flare, the plasma temperature drops quickly. The predicted luminosity for such transitions, Lx = 3 x 10^35 erg s-1, is very close to the luminosity of 4U1700-37 during quiescence. The transition may be caused by the accretion of a clump in the stellar wind of the donor star. Thus, a magnetised NS nature of the compact object is strongly favoured.

The transformation of powerful gravitational waves, created by the coalescence of massive black hole binaries, into electromagnetic radiation in external magnetic fields is revisited. In contrast to the previous calculations of the similar effect, we study the realistic case of the gravitational radiation frequency below the plasma frequency of the surrounding medium. The gravitational waves propagating in the plasma constantly create electromagnetic radiation dragging it with them, despite the low frequency. The plasma heating by the unattenuated electromagnetic wave may be significant in a hot rarefied plasma with strong magnetic field and can lead to a noticeable burst of electromagnetic radiation with higher frequency. The graviton-to-photon conversion effect in plasma is discussed in the context of possible electromagnetic counterparts of GW150914 and GW170104.

4U0114$+$65 is one of the slowest known X-ray pulsars. We present an analysis of a pointed observation by the XMM-Newton X-ray telescope. The energy-resolved light curve, the time-resolved X-ray spectra provided by the EPIC cameras on board and the first high-resolution spectrum of this source, provided by the Reflection Grating Spectrometer, were analysed. An X-ray pulse of $9350\pm 160$ s was measured, confirming the secular spin up of this source. We also fit the pulse-phase-resolved spectra with Comptonisation models. These models imply a very small ($r\sim 3$ km) and hot ($kT\sim 2-3$ keV) emitting region and therefore point to a hot spot over the neutron star (NS) surface as the most reliable explanation for the X-ray pulse. The long NS spin period, the spin-up rate, and persistent X-ray emission can be explained within the theory of quasi-spherical settling accretion, which may indicate that the magnetic field is in the magnetar range. Thus, 4U 0114$+$65 could be a wind-accreting magnetar. We also observed two episodes of low luminosity. The first can be explained as an absorption by an over-dense structure in the wind of the donor. The second episode may be due to temporal cessation of accretion onto one magnetic pole caused by non-spherical matter capture from the structured stellar wind. The light curve displays two types of dips seen during the high-flux intervals. The short dips, with durations of tens of seconds, are produced through absorption by wind clumps. The long dips might be associated with the rarefied interclump medium. We also found evidence of emission lines in the X-ray photoionised wind of the B1Ia donor. The Fe K$\alpha$ line is highly variable and much weaker than in other X-ray binaries with supergiant donors. The degree of wind clumping, measured through the covering fraction, is much lower than in supergiant donor stars with earlier spectral types.

Pulsations from the high mass X-ray binary AXJ1910.7+0917 were discovered during Chandra observations performed in 2011 (Israel et al. 2016). We report here more details on this discovery and discuss the source nature. The period of the X-ray signal is P=36200+/-110s, with a pulsed fraction, PF, of 63+/-4%. Given the association with a massive B-type companion star, we ascribe this long periodicity to the rotation of the neutron star, making AXJ1910.7+0917 the slowest known X-ray pulsar. We report also on the spectroscopy of XMM-Newton observations that serendipitously covered the source field, resulting in an highly absorbed (column density almost reaching 1e23cm-2), power law X-ray spectrum. The X-ray flux is variable on a timescale of years, spanning a dynamic range >60. The very long neutron star spin period can be explained within a quasi-spherical settling accretion model, that applies to low luminosity, wind-fed, X-ray pulsars.

Primordial black holes (PBHs) that form in the early Universe in the modified Affleck-Dine (AD) mechanism of baryogenesis should have intrinsic log-normal mass distribution of PBHs. We show that the parameters of this distribution adjusted to provide the required spatial density of massive seeds ($\ge 10^4 M_\odot$) for early galaxy formation and not violating the dark matter density constraints, predict the existence of the population of intermediate-mass PBHs with a number density of $\sim 1000$~Mpc$^{-3}$. We argue that the population of intermediate-mass AD PBHs can also seed the formation of globular clusters in galaxies. In this scenario, each globular cluster should host an intermediate-mass black hole with a mass of a few thousand solar masses, and should not obligatorily be immersed in a massive dark matter halo.

Possible formation mechanisms of massive close binary black holes that can merge in the Hubble time to produce powerful gravitational wave bursts detected during advanced LIGO O1 science run are briefly discussed. The pathways include the evolution from field low-metallicity massive binaries, the dynamical formation in globular clusters and primordial black holes. Low effective black hole spins inferred for LIGO GW150914 and LTV151012 events are discussed. Population synthesis calculations of the expected spin and chirp mass distributions from the standard field massive binary formation channel are presented for different metallicities (from zero-metal Population III stars up to solar metal abundance). We conclude that that merging binary black holes can contain systems from different formation channels, discrimination between which can be made with increasing statistics of mass and spin measurements from ongoing and future gravitational wave observations.

X-ray spectra of accreting pulsars are generally observed to vary with their X-ray luminosity. In particular, the hardness of the X-ray continuum is found to depend on luminosity. In a few sources, the correlation between the energy of the cyclotron resonance scattering feature (CRSF) and the luminosity is clear. Different types (signs) of the correlation are believed to reflect different accretion modes. We analyse two NuSTAR observations of the transient accreting pulsar Cep X-4 during its 2014 outburst. Our analysis is focused on a detailed investigation of the dependence of the CRSF energy and of the spectral hardness on X-ray luminosity, especially on short timescales. To investigate the spectral changes as a function of luminosity within each of the two observations, we used the intrinsic variability of the source on the timescale of individual pulse cycles (tens of seconds), the so-called pulse-to-pulse variability. We find that the NuSTAR spectrum of Cep X-4 contains two CRSFs: the fundamental line at ~30 keV and its harmonic at ~55 keV. We find for the first time that the energy of the fundamental CRSF increases and the continuum becomes harder with increasing X-ray luminosity not only between the two observations, that is, on the long timescale, but also within an individual observation, on the timescale of a few tens of seconds. We investigate these dependencies in detail including their non-linearity. We discuss a possible physical interpretation of the observed behaviour in the frame of a simple one-dimensional model of the polar emitting region with a collisionless shock formed in the infalling plasma near the neutron star surface. With this model, we are able to reproduce the observed variations of the continuum hardness ratio and of the CRSF energy with luminosity.

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

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

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

Features and applications of quasi-spherical settling accretion onto rotating magnetized neutron stars in high-mass X-ray binaries are discussed. The settling accretion occurs in wind-fed HMXBs when the plasma cooling time is longer than the free-fall time from the gravitational capture radius, which can take place in low-luminosity HMXBs with $L_x\lesssim 4\times 10^{36}$ erg/s. We briefly review the implications of the settling accretion, focusing on the SFXT phenomenon, which can be related to instability of the quasi-spherical convective shell above the neutron star magnetosphere due to magnetic reconnection from fast temporarily magnetized winds from OB-supergiant. If a young neutron star in a wind-fed HMXB is rapidly rotating, the propeller regime in a quasi-spherical hot shell occurs. We show that X-ray spectral and temporal properties of enigmatic $\gamma$ Cas Be-stars are consistent with failed settling accretion regime onto a propelling neutron star. The subsequent evolutionary stage of $\gamma$ Cas and its analogs should be the X Per-type binaries comprising low-luminosity slowly rotating X-ray pulsars.

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.

The black hole binary properties inferred from the LIGO gravitational wave signal GW150914 posed several serious problems. The high masses and low effective spin of black hole binary can be explained if they are primordial (PBH) rather than the products of the stellar binary evolution. Such PBH properties are postulated ad hoc but not derived from fundamental theory. We show that the necessary features of PBHs naturally follow from the slightly modified Affleck-Dine (AD) mechanism of baryogenesis. The log-normal distribution of PBHs, predicted within the AD paradigm, is adjusted to provide an abundant population of low-spin stellar mass black holes. The same distribution gives a sufficient number of quickly growing seeds of supermassive black holes observed at high redshifts and may comprise an appreciable fraction of Dark Matter which does not contradict any existing observational limits. Testable predictions of this scenario are discussed.

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

The enigmatic X-ray emission from the bright optical star, $\gamma$ Cassiopeia, is a long-standing problem. $\gamma$ Cas is known to be a binary system consisting of a Be-type star and a low-mass ($M\sim 1\,M_\odot$) companion of unknown nature orbiting in the Be-disk plane. Here we apply the quasi-spherical accretion theory onto a compact magnetized star and show that if the low-mass companion of $\gamma$ Cas is a fast spinning neutron star, the key observational signatures of $\gamma$ Cas are remarkably well reproduced. Direct accretion onto this fast rotating neutron star is impeded by the propeller mechanism. In this case, around the neutron star magnetosphere a hot shell is formed that emits thermal X-rays in qualitative and quantitative agreement with observed properties of the X-ray emission from $\gamma$ Cas. We suggest that $\gamma$ Cas and its analogs constitute a new subclass of Be-type X-ray binaries hosting rapidly rotating neutron stars formed in supernova explosions with small kicks. The subsequent evolutionary stage of $\gamma$ Cas and its analogs should be the X Per-type binaries comprising low-luminosity slowly rotating X-ray pulsars. The model explains the enigmatic X-ray emission from $\gamma$ Cas, and also establishes evolutionary connections between various types of rotating magnetized neutron stars in Be-binaries.

Rotating proto-neutron stars can be important sources of gravitational waves to be searched for by present-day and future interferometric detectors. It was demonstrated by Imshennik that in extreme cases the rapid rotation of a collapsing stellar core may lead to fission and formation of a binary proto-neutron star which subsequently merges due to gravitational wave emission. In the present paper, we show that such dynamically unstable collapsing stellar cores may be the product of a former merger process of two stellar cores in a common envelope. We applied population synthesis calculations to assess the expected fraction of such rapidly rotating stellar cores which may lead to fission and formation of a pair of proto-neutron stars. We have used the BSE population synthesis code supplemented with a new treatment of stellar core rotation during the evolution via effective core-envelope coupling, characterized by the coupling time, $\tau_c$. The validity of this approach is checked by direct MESA calculations of the evolution of a rotating 15 $M_\odot$ star. From comparison of the calculated spin distribution of young neutron stars with the observed one, reported by Popov and Turolla, we infer the value $\tau_c \simeq 5 \times 10^5$ years. We show that merging of stellar cores in common envelopes can lead to collapses with dynamically unstable proto-neutron stars, with their formation rate being $\sim 0.1-1\%$ of the total core collapses, depending on the common envelope efficiency.

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

Gamma-ray haloes can exist around galaxies due to the interaction of escaping galactic cosmic rays with the surrounding gas. We have searched for such a halo around the nearby giant spiral Andromeda galaxy M31 using almost 7 years of Fermi LAT data at energies above 300 MeV. The presence of a diffuse gamma-ray halo with total photon flux $2.6\pm 0.6\times 10^{-9}$ cm$^{-2}$~s$^{-1}$, corresponding to a luminosity (0.3-100 GeV) of $(3.2\pm 0.6)\times 10^{38}$ erg s$^{-1}$ (for a distance of 780~kpc) was found at a 5.3$\sigma$ confidence level. The halo form does not correspond to the extended baryonic HI disc of M31, as would be expected in hadronic production of gamma photons from cosmic ray interaction, nor it is spherically symmetric, as could be in the case of dark matter annihilation. The best-fit halo template corresponds to two 6-7.5 kpc bubbles symmetrically located perpendicular to the M31 galactic disc, similar to the 'Fermi bubbles' found around the Milky Way centre, which suggests the past activity of the central supermassive black hole or a star-formation burst in M31.

We have characterized the typical temporal behaviour of the bright X-ray flares detected from the three Supergiant Fast X-ray Transients showing the most extreme transient behaviour (XTEJ1739-302, IGRJ17544-2619, SAXJ1818.6-1703). We focus here on the cumulative distributions of the waiting-time (time interval between two consecutive X-ray flares), and the duration of the hard X-ray activity (duration of the brightest phase of an SFXT outburst), as observed by INTEGRAL/IBIS in the energy band 17-50 keV. Adopting the cumulative distribution of waiting-times, it is possible to identify the typical timescale that clearly separates different outbursts, each composed by several single flares at ks timescale. This allowed us to measure the duration of the brightest phase of the outbursts from these three targets, finding that they show heavy-tailed cumulative distributions. We observe a correlation between the total energy emitted during SFXT outbursts and the time interval covered by the outbursts (defined as the elapsed time between the first and the last flare belonging to the same outburst as observed by INTEGRAL). We show that temporal properties of flares and outbursts of the sources, which share common properties regardless different orbital parameters, can be interpreted in the model of magnetized stellar winds with fractal structure from the OB-supergiant stars.

Theories of galaxy formation predict the existence of extended gas halo around spiral galaxies. If there are 10-100 nG magnetic fields at several ten kpc distances from the galaxies, extended galactic cosmic ray (CR) haloes could also exist. Galactic CRs can interact with the tenuous hot halo gas to produce observable gamma-rays. We have performed search for a gamma-ray halo around the M31 galaxy -- the closest large spiral galaxy. Our analysis of almost 7 years of the Fermi LAT data revealed the presence of a spatially extended diffuse emission excess around M31. The data can be fitted using the simplest morphology of a uniformly bright circle. The best fit gave a 4.7$\sigma$ significance for a $0.9^{\circ}$ (12 kpc) halo with a photon flux of $\sim (3.2\pm1.0)\times 10^{-9} ~\mathrm{cm^{-2}s^{-1}}$ and a luminosity of $(4.0\pm1.5)\times 10^{38} ~\mathrm{erg~s^{-1}}$ in the energy range 0.3--100 GeV. Our results also imply a low level of the flux from the disc of the M31 galaxy $(3.3 \pm 1.0) \times 10^{-10}~\mathrm{cm^{-2}s^{-1}}$. The corresponding gamma-ray luminosity, $5\times10^{37} ~\mathrm{erg~s^{-1}}$ is several times smaller than the luminosity of the Milky Way. This difference could be explained by a lower star formation rate in M31: there are less CRs and the level of the ISM turbulence is lower, which in turn leads to a shorter time of CR containment.

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

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

Context: The magnetic field is a crucial ingredient of neutron stars. It governs the physics of accretion and of the resulting high-energy emission in accreting pulsars. Studies of the cyclotron resonant scattering features (CRSFs) seen as absorption lines in the X-ray spectra of the pulsars permit direct measuremets of the field strength. Aims: From an analysis of a number of pointed observations with different instruments, the energy of CRSF, Ecyc, has recently been found to decay in Her X-1, which is one of the best-studied accreting pulsars. We present our analysis of a homogeneous and almost uninterrupted monitoring of the line energy with Swift/BAT. Methods: We analyzed the archival Swift/BAT observations of Her X-1 from 2005 to 2014. The data were used to measure the CRSF energy averaged over several months. Results: The analysis confirms the long-term decay of the line energy. The downward trend is highly significant and consistent with the trend measured with the pointed observations: dEcyc/dt ~-0.3 keV per year. Conclusions: The decay of Ecyc either indicates a local evolution of the magnetic field structure in the polar regions of the neutron star or a geometrical displacement of the line-forming region due to long-term changes in the structure of the X-ray emitting region. The shortness of the observed timescale of the decay, -Ecyc/(dEcyc/dt) ~ 100 yr, suggests that trend reversals and/or jumps of the line energy might be observed in the future.

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

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

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

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

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

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

We consider the sensitivity of the pulsar timing array (PTA) technique to specific kind of narrow-band stochastic signals in nano-Hz frequency range. Specifically, we examine the narrow-band signal produced by oscillating gravitational scalar potentials in the Galaxy (Gravitational Potential Background), which arise if an ultralight massive scalar field is the galactic dark matter. We have performed a Bayesian analysis of publicly available data on 12 pulsars obtained by the NANOGrav project. In the monochromatic approximation, the upper limit on the variable gravitational potential amplitude is $\Psi_c<1.14 \times 10^{-15}$, corresponding to the dimensionless strain amplitude $h_c=2\sqrt{3}\Psi_c < 4\times 10^{-15}$ at frequency $f=1.75\times 10^{-8} Hz$. In the narrow-band approximation, the upper limit on the energy density of GPB is found to be $\Omega_{GPB}<1.27 \times 10^{-9}$ at $f=6.2\times 10^{-9} Hz$. These limits are an order of magnitude higher than the theoretically expected values, if the ultralight scalar field with a mass of $\sim 10^{-23} eV$ is assumed to be the galactic dark matter with a local density of $\sim 0.3 GeV cm^{-3}$.

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 consider consequences of hypothetical existence of baryo-dense stars created in the very early universe within an extension of Affleck-Dine scenario of baryogenesis. New constraints on the possible number of compact antimatter objects are derived. The contemporary observational data do not exclude significant amount of antimatter in the Galaxy (and in other galaxies) in the form of the baryo-dense low-massive stars.

We performed a Bayesian analysis of pulsar timing residuals from the NANOGrav pulsar timing array to search for a specific form of stochastic narrow-band signal produced by oscillating gravitational potential in the Galactic halo. Such oscillations arise in models of warm dark matter composed of an ultralight massive scalar field (m\sim 10^{-23} eV), recently considered by Khmelnitsky and Rubakov (2014). In the monochromatic approximation, the stringent upper limit (95% c.l.) on the variable gravitational potential amplitude is found to be Psi_c<1.14 10^{-15}, corresponding to the characteristic strain h_c=2\sqrt{3}\Psi_c < 4 10^{-15} at frequency f=1.75 10^{-8} Hz. In the narrow-band approximation, the upper limit of this background energy density is \Omega_{GPB}<1.27 10^{-9} at frequency f=6.2 10^{-9} Hz. These limits are an order of magnitude higher than the expected signal amplitude if the galactic dark matter consists of the ultralight scalar field. The applied analysis of pulsar timing residuals can be used to search for any narrow-band stochastic signals with different correlation properties. As a by-product, parameters of the red noise present in four NANOGrav pulsars (J1713+0747, J2145-0750, B1855+09, J1744-1134) have been evaluated.

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

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

We review the formation and evolution of compact binary stars consisting of white dwarfs (WDs), neutron stars (NSs), and black holes (BHs). Mergings of compact binary stars are expected to be the most important sources for the forthcoming gravitational-wave (GW) astronomy. In the first part of the review, we discuss observational manifestations of close binary stars with NS and/or black components and their merger rate, crucial points in the formation and evolution of compact stars in binary systems, including the treatment of the natal kicks which NSs and BHs acquire during the core collapse of massive stars and the common envelope phase of binary evolution, which are most relevant to the merging rates of NS-NS, NS-BH and BH-BH binaries. The second part of the review is devoted mainly to formation and evolution of binary WDs and their observational manifestations, including their role as progenitors of cosmologically important thermonuclear SN Ia. We also consider AM CVn-stars which are thought to be the best verification binary GW sources for future low-frequency GW space interferometers.

Based on multiyear INTEGRAL observations of SS433 in 2003-2011, a composite IBIS/ISGRI 18-60 keV orbital light curve is constructed around zero precessional phases psi_{pr}= 0 at the maximim accretion disk opening angle. It shows a peculiar shape with significant excess near the orbital phase phi_orb= 0.25, which is not seen in the softer 2-10 keV energy band. The 40-60 keV orbital light curve demonstrates two almost equal humps at phases \sim 0.25 and \sim 0.75, most likely due to nutation effects of the accretion disk. The nutational variability of SS433 in 15-50 keV with a period of 6.290 days is independently found from timing analysis of Swift/BAT data. The change of the off-eclipse 18-60 keV X-ray flux with the precessional phase shows a double-wave form with strong primary maximum at psi_{pr}= 0 and weak but significant secondary maximum at psi_{pr}= 0.6. A weak variability of the 18-60 keV flux in the middle of the orbital eclipse correlated with the disk precessional phase is also observed. The joint analysis of the broadband 18-60 keV orbital and precessional light curves confirms the presence of a hot extended corona in the central parts of the supercritical accretion disk and constrains the binary mass ratio in SS433 in the range 0.5>q>0.3, suggesting the black hole nature of the compact object.

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

Properties of the population of millisecond extragalactic radio bursts discovered by Thornton et al. (2013) are in good correspondence with the hypothesis that such events are related to hyperflares of magnetars, as was proposed by us after the first observation of an extragalactic millisecond radio burst by Lorimer et al. (2007). We also point that some of multiple millisecond radio bursts from M31 discovered by Rubio-Herrera et al. (2013) also can be related to weaker magnetar bursts.

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

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

We study thermal emission from circumstellar structures heated by gamma-ray burst (GRB) radiation and ejecta and calculate its contribution to GRB optical and X-ray afterglows using the modified radiation hydro-code small STELLA. It is shown that thermal emission originating in heated dense shells around the GRB progenitor star can reproduce X-ray plateaus (like observed in GRB 050904, 070110) as well as deviations from a power law fading observed in optical afterglows of some GRBs (e.g. 020124, 030328, 030429X, 050904). Thermal radiation pressure in the heated circumburst shell dominates the gas pressure, producing rapid expansion of matter similar to supenova-like explosions close to opacity or radiation flux density jumps in the circumburst medium. This phenomenon can be responsible for so-called supernova bumps in optical afterglows of several GRBs. Such a `quasi-supernova' suggests interpretation of the GRB-SN connection which does not directly involve the explosion of the GRB progenitor star.

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

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

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

Based on multiyear INTEGRAL observations of SS433, a composite IBIS/ISGRI 18-60 keV orbital light curve is constructed around zero precessional phase $\psi_{pr}= 0$. It shows a peculiar shape characterized by a significant excess near the orbital phase $\phi_{orb}= 0.25$, which is not seen in the softer 2-10 keV energy band. Such a shape is likely to be due to a complex asymmetric structure of the funnel in a supercritical accretion disk in SS433. The orbital light curve at 40-60 keV demonstrates two almost equal bumps at phases $\sim 0.25$ and $\sim 0.75$, most likely due to nutation effects of the accretion disk. The change of the off-eclipse 18-60 keV X-ray flux with the precessional phase shows a double-wave form with strong primary maximum at $\psi_{pr}= 0$ and weak but significant secondary maximum at $\psi_{pr}= 0.6$. A weak variability of the 18-60 keV flux in the middle of the orbital eclipse correlated with the disk precessional phase is also observed. The joint analysis of the broadband (18-60 keV) orbital and precessional light curves obtained by INTEGRAL confirms the presence of a hot extended corona in the central parts of the supercritical accretion disk and constrain the binary mass ratio in SS433 in the range $0.5\gtrsim q\gtrsim 0.3$, confirming the black hole nature of the compact object. Orbital and precessional light curves in the hardest X-ray band 40-60 keV, which is free from emission from thermal X-ray jets, are also best fitted by the same geometrical model with hot extended corona at $q\sim 0.3$, stressing the conclusions of the modeling of the broad-band X-ray orbital and precessional light curves.

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

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

Symbiotic X-ray binaries (SyXBs) comprise a rare class of low-mass X-ray binaries. We study the Galactic SyXBs, which we consider as detached binaries composed of low-mass giants and wind-fed neutron star companions, by simulation of the interaction of a magnetized neutron star (NS) with its environment and utilizing a population synthesis code. We focus mainly on the parameters that influence observational appearance of the SyXB: the donor wind velocity (vw) and the angular momentum distribution in the shell of matter settling onto NS. We estimate the birthrate of SyXB as $\sim 4.1\times 10^{-5}$ yr$^{-1}$ to $ \sim 6.6\times 10^{-6}$ yr$^{-1}$ and their number in the Galaxy as $\sim$(100 -- 1000). Assumed stellar wind velocity from cool giants is the input parameter that influences the model SyXBs population most. Among known SyXBs or candidate systems, 4U 1954+31 and IGR J16358-4724 in which NS have very long spin periods may host quasi-spherically accreting NSs. GX 1+4 has a peculiar long-term spin behaviour and it may also be a quasi-spherical wind-accreting source. We cannot identify whether there are wind-fed accretion disks in 4U 1700+24, Sct X-1, IRXS J180431.1-273932 and 2XMM J174016.0-290337.

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

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

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

Context: Episodic flaring activity is a common feature of X-ray pulsars in HMXBs. In some Be/X-ray binaries flares were observed in quiescence or prior to outbursts. EXO 2030+375 is a Be/X-ray binary showing "normal" outbursts almost every ~46 days, near periastron passage of the orbital revolution. Some of these outbursts were occasionally monitored with the INTEGRAL observatory. Aims: The INTEGRAL data revealed strong quasi-periodic flaring activity during the rising part of one of the system's outburst. Such activity has previously been observed in EXO 2030+375 only once, in 1985 with EXOSAT. (Some indications of single flares have also been observed with other satellites.) Methods: We present the analysis of the flaring behavior of the source based on INTEGRAL data and compare it with the flares observed in EXO 2030+375 in 1985. Results: Based on the observational properties of the flares, we argue that the instability at the inner edge of the accretion disk is the most probable cause of the flaring activity.

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

We present an analysis of the highly variable accreting X-ray pulsar 3A 1954+319 using 2005-2009 monitoring data obtained with INTEGRAL and Swift. This considerably extends the pulse period history and covers flaring episodes in 2005 and 2008. In 2006 the source was identified as one of only a few known symbiotic X-ray binaries (SyXBs), i.e., systems composed of a neutron star accreting from the inhomogeneous medium around an M-giant star. The extremely long pulse period of 5.3 hr is directly visible in the 2008 INTEGRAL-ISGRI outburst light curve. The pulse profile is double peaked and generally not significantly energy dependent although there is an indication of possible softening during the main pulse. During the outburst a strong spin-up of -1.8 10^(-4) hr hr^(-1) occurred. Between 2005 and 2008 a long-term spin-down trend of 2.1 10^-5 hr hr^(-1) was observed for the first time for this source. The 3-80 keV pulse peak spectrum of 3A 1954+319 during the 2008 flare could be well described by a thermal Comptonization model. We interpret the results within the framework of a recently developed quasi-spherical accretion model for SyXBs.

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

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

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

We analyze a flux-limited sample of persistent and bright (with 2-10 keV fluxes exceeding 1.4e-10 erg/s/cm2) low-mass X-ray binaries (LMXBs) in our Galaxy. It is demonstrated that the majority of binary systems with X-ray luminosities below logL(erg/sec)~37.3 have unevolved secondary companions (except for those with white dwarf donors), while systems with higher X-ray luminosity predominantly harbor giant donors. Mass transfer in binary systems with giants significantly shortens their life time thus steepening the X-ray luminosity function of LMXBs at high luminosity. We argue that this is the reason why the LMXB luminosity function constructed in the last years from observations of sources in our and distant galaxies demonstrates a break at logL(erg/sec)~37.3.

We present methodology and results of numerical modeling of the interaction of GRB prompt emission with the circumburst medium using a modified version of the multi-group radiation hydrocode STELLA. The modification includes the nonstationary photoionization, the photoionization heating and the Compton heating along with the hydrodynamics and radiation transfer. The lightcurves and spectra of the outcoming gamma-ray, X-ray and optical emission are calculated for a set of models (shells) of the circumburst environment, which differ in dimensions, density, density profile, composition, temperature. In some cases total bolometric and optical luminosities can reach 10^47 and 10^43 erg/s respectively. These effects can be responsible for irregularities which are seen on lightcurves of some GRB's X-ray and optical afterglows.

We discuss possible radio bursts which can be generated during binary neutron stars mergers associated with short gamma-ray bursts. Low-frequency radio band appear to be advantageous due to the time delay of a radio signal propagating in the intergalactic medium, which makes it possible to use short gamma-ray burst alerts to search for specific regions on the sky by low-frequency radio instruments, especially LOFAR. The LOFAR sensitivity will allow significant detection of such bursts up to redshifts z~0.3 with a rate of ~25 events per year.

We discuss a possible generation of radio bursts preceding final stages of binary neutron star mergings which can be accompanied by short gamma-ray bursts. Detection of such bursts appear to be advantageous in the low-frequency radio band due to a time delay of ten to several hundred seconds required for radio signal to propagate in the ionized intergalactic medium. This delay makes it possible to use short gamma-ray burst alerts to promptly monitor specific regions on the sky by low-frequency radio facilities, especially by LOFAR. To estimate the strength of the radio signal, we assume a power-law dependence of the radio luminosity on the total energy release in a magnetically dominated outflow, as found in millisecond pulsars. Based on the planned LOFAR sensitivity at 120 MHz, we estimate that the LOFAR detection rate of such radio transients could be about several events per month from redshifts up to $z\sim1.3$ in the most optimistic scenario. The LOFAR ability to detect such events would crucially depend on exact efficiency of low-frequency radio emission mechanism.

We propose a novel method for observing the gravitational wave signature of super-massive black hole (SMBH) mergers. This method is based on detection of a specific type of gravitational waves, namely gravitational wave burst with memory (BWM), using pulsar timing. We study the unique signature produced by BWM in anomalous pulsar timing residuals. We show that the present day pulsar timing precision allows one to detect BWM due to SMBH mergers from distances up to 1 Gpc (for case of equal mass 10^8 Msun SMBH). Improvements in precision of pulsar timing together with the increase in number of observed pulsars should eventually lead to detection of a BWM signal due to SMBH merger, thereby making the proposed technique complementary to the capabilities of the planned LISA mission.

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

The Be/X-ray binary A 0535+26 showed a normal (type I) outburst in August/September 2005, which reached a maximum X-ray flux of 400mCrab in the 5-100keV range. The outburst was observed by INTEGRAL and RXTE. The energy of the fundamental cyclotron line has been measured with INTEGRAL and RXTE at ~45keV. Flaring activity was observed during the rise to the peak of the outburst. RXTE observations during one of these flares found the energy of the fundamental cyclotron line shifted to a significantly higher position than during the rest of the outburst, where it remains constant. Also, the energy-dependent pulse profiles during the flare differ significantly from the rest of the outburst. These differences have been interpreted with the presence of magnetospheric instabilities at the onset of the accretion. A decomposition method is applied to A 0535+26 pulse profiles. Basic assumptions of the method are that the asymmetry observed in the pulse profiles is caused by a distorted magnetic dipole field, and that the emission regions have axisymmetric beam patterns. Using pulse profiles obtained from RXTE observations, the contribution of the two emission regions has been disentangled. Constraints on geometry of the pulsar and a possible solution of the beam pattern are given. First results of the comparison of the reconstructed beam pattern with a geometrical model that includes relativistic light deflection are presented.

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

Coalescing neutron star binaries are believed to be the most reliable sources for ground-based detectors of gravitational waves and likely progenitors of short gamma-ray bursts. In the process of coalescence, magnetic fields of neutron stars can induce interesting observational manifestations and affect the form of gravitational wave signal. In this papaer we use the population synthesis method to model the expected distribution of neutron star magnetic fields during the coalescence under different assumptions on the initial parameters of neutron stars and their magnetic field evolution. We discuss possible elecotrmagnetic phenomena preceding the coalescence of magnetized neutron star binaries and the effect of magnetic field on the gravitational wave signal. We find that a log-normal (Gaussian in logarithms) distribution of the initial magnetic fields of neutron stars, which agrees with observed properties of radio pulsars, produces the distribution of the magnetic field energy during the coalescence that adequately describes the observed luminosity function of short gamma-ray bursts under different assumptions on the field evolution and initial parameters of neutron stars. This agreement lends further support to the model of coalescing neutron star binaries as progenitors of gamma-ray bursts.

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.

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

The analysis of hard X-ray INTEGRAL observations (2003-2008) of superaccreting galactic microquasar SS433 at precessional phases of the source with the maximum disk opening angle is carried out. It is found that the shape and width of the primary X-ray eclipse is strongly variable suggesting additional absorption in dense stellar wind and gas outflows from the optical A7I-component and the wind-wind collision region. The independence of the observed hard X-ray spectrum on the accretion disk precessional phase suggests that hard X-ray emission (20-100 keV) is formed in an extended, hot, quasi-isothermal corona, probably heated by interaction of relativistic jet with inhomogeneous wind outflow from the precessing supercritical accretion disk. A joint modeling of X-ray eclipsing and precessional hard X-ray variability of SS433 revealed by INTEGRAL by a geometrical model suggests the binary mass ratio $q=m_x/m_v\simeq 0.25\div 0.5$. The absolute minimum of joint orbital and precessional $\chi^2$ residuals is reached at $q\simeq 0.3$. The found binary mass ratio range allows us to explain the substantial precessional variability of the minimum brightness at the middle of the primary optical eclipse. For the mass function of the optical star $f_v=0.268 M_\odot$ as derived from Hillwig & Gies data, the obtained value of $q\simeq 0.3$ yields the masses of the components $m_x\simeq 5.3 M_\odot$, $m_v\simeq 17.7 M_\odot$, confirming the black hole nature of the compact object in SS433.

We develop a Monte-Carlo technique based on L.B. Lucy's indivisible photon packets method to calculate X-ray continuum spectra of comptonized thermal plasma in arbitrary geometry and apply it to describe the broadband X-ray continuum of the galactic superaccreting microquasar SS433 observed by INTEGRAL. A physical model of the X-ray emitting region is proposed that includes thermal emission from the accretion disk, jets and hot corona where the photons of different origin are comptonized. From comparison with INTEGRAL observations, we estimate physical parameters of the complex X-ray emitting region in SS433 and present model spectra for different viewing angles of the object.

The role of binary progenitors of neutron stars in the apparent distribution of space velocities and spin-velocity alignment observed in young pulsars is studied. A Monte-Carlo synthesis of pulsar population from single and binary stars with different assumptions about the NS natal kick model (direction distribution, amplitude, and kick reduction in binary progenitors which experienced mass exchange due to Roche lobe overflow with initial masses on the main sequence from the range 8-11 $M_\odot$) is performed. The calculated spin-velocity alignment distributions are compared with observational data obtained from radio polarization measurements. The observed space velocity of pulsars is found to be mostly shaped by the natal kick velocity form and its amplitude; the fraction of binaries is not important here for reasonably large kicks. The distribution of kick direction relative to the spin axis during the formation of a NS is found to affect strongly the spin-velocity correlation of pulsars. Comparison with observed pulsar spin-velocity angles favours a sizeable fraction of binary progenitors and the kick-spin angle $\sim 5-20^\circ$. The form of the initial binary mass ratio distribution does not affect our results.

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

The effect of a narrow-band isotropic stochastic GW background on pulsar timing and astrometric measurements is studied. Such a background appears in some theories of gravity. We show that the existing millisecond pulsar timing accuracy ($\sim 0.2 \rm{\mu s}$) strongly constrains possible observational consequences of theory of massive gravity with spontaneous Lorentz braking \cite{dtt:2005}, essentially ruling out significant contribution of massive gravitons to the local dark halo density. The present-day accuracy of astrometrical measurements ($\sim 100 \rm{\mu as}$) sets less stringent constraints on this theory.

The analysis of INTEGRAL observations (2003-2008) of superaccreting galactic microquasar SS433 at precessional phases with the maximum disk opening angle is carried out. The shape and width of the primary X-ray eclipse is found to be strongly variable suggesting additional absorption in dense stellar wind and gas outflows from the optical A7I-component. The joint modeling of X-ray eclipse and precessional X-ray variability by a geometrical model suggests the binary mass ratio q=m_x/m_v=0.3, allowing an explnation of peculiarities of the optical variability of SS433, in particular, the substantial precessional variability at the primary optical eclipse minimum. For the mass function of the optical star f_v=0.268 M_\odot as derived from Hillwig and Gies (2008) data, the obtained q yields the masses of the components m_x=5 M_\odot, m_v=15 M_\odot, confirming the black hole nature of the relativistic object in SS433. The independence of the observed hard X-ray spectrum on the precession phase suggests that hard X-ray emission is formed in an extended hot corona. The Monte-Carlo simulations of the broadband X-ray spectrum of SS433 at the maximum disk opening precessional phases allowed us to determine physical parameters of the corona (temperature T_{cor}=20 keV, Thomson optical depth \tau=0.2), and to estimate the jet mass outflow rate \dot M_j=3\times 10^{19} g/s yielding the kinetic power of the jets \sim 10^{39} erg/s.

Impact of supernova explosion on the neutron star magnetosphere in a massive binary system is considered. The supernova shock striking the NS magnetosphere filled with plasma can lead to the formation of a magnetospheric tail with significant magnetic energy. The magnetic field reconnection in the current sheet formed can convert the magnetic energy stored in the tail into kinetic energy of accelerated charged particles. Plasma instabilities excited by beams of relativistic particles can lead to the formation of a short pulse of coherent radio emission with parameters similar to those of the observed bright extragalactic millisecond radio burst (Lorimer et al. 2007).

In this work, analyzing the propagation of electromagnetic waves in the field of gravitational waves, we show the presence and significance of the so called surfing effect for pulsar timing measurements. It is shown that, due to the transverse nature of gravitational waves, the surfing effect leads to enormous pulsar timing residuals if the speed of gravitational waves is smaller than speed of light. This fact allows to place significant constraints on parameter $\epsilon$, which characterizes the relative deviation of the speed of gravitational waves from the speed of light. We show that the existing constraints from pulsar timing measurements already place stringent limits on $\epsilon$ and consequently on the mass of graviton $m_g$. These limits on $m_g$ are three orders of magnitude stronger than the current constraints from Solar System tests. The current constraints also allow to rule out massive gravitons as possible candidates for cold dark matter in galactic halo. In the near future, the gravitational wave background from extragalactic super massive black hole binaries, along with the expected sub-microsecond pulsar timing accuracy, will allow to achieve constrains of $\epsilon\lesssim 0.4\%$ and possibly stronger.

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 study the spectral and temporal behavior of the High Mass X-ray Binary A 0535+26 during a `pre-outburst flare' which took place ~5 d before the peak of a normal (type I) outburst in August/September 2005. We compare the studied behavior with that observed during the outburst. We analyse RXTE observations that monitored A 0535+26 during the outburst. We complete spectral and timing analyses of the data. We study the evolution of the pulse period, present energy-dependent pulse profiles both at the initial pre-outburst flare and close to outburst maximum, and measure how the cyclotron resonance-scattering feature (hereafter CRSF) evolves. We present three main results: a constant period P=103.3960(5)s is measured until periastron passage, followed by a spin-up with a decreasing period derivative of Pdot=(-1.69+/-0.04)x10^(-8)s/s at MJD 53618, and P remains constant again at the end of the main outburst. The spin-up provides evidence for the existence of an accretion disk during the normal outburst. We measure a CRSF energy of Ecyc~50kev during the pre-outburst flare, and Ecyc~46kev during the main outburst. The pulse shape, which varies significantly during both pre-outburst flare and main outburst, evolves strongly with photon energy.

We argue that X-ray flaring variability observed in the transient X-ray pulsar A 0535+26 is due to low-mode magnetospheric instability. This instability develops at the onset of accretion, in the thin boundary layer between the accretion disk and neutron star magnetosphere. As a result, the matter collected in the boundary layer can rapidly fall onto the NS surface close to the magnetic poles, but not exactly along the field lines by which the stationary accretion proceeds. This explains the shift in cyclotron line energy measured using RXTE data in a pre-outburst spike, with respect to the line energy observed during the main outburst. Furthermore, the instability can account for the difference in pulse profiles, and their energy evolution that is different in the pre-outburst flare and main outburst.

We propose that the strong millisecond extragalactic radio burst (mERB) discovered by Lorimer et al. (2007) may be related to a hyperflare from an extragalactic soft gamma-ray repeater. The expected rate of such hyperflares, $\sim$ 20 - 100 d$^{-1}$ Gpc$^{-3}$, is in good correspondence with the value estimated by Lorimer et al. The possible mechanism of radio emission can be related to the tearing mode instability in the magnetar magnetosphere as discussed by Lyutikov (2002), and can produce the radio flux corresponding to the observed $\sim$ 30 Jy from the mERB using a simple scaling of the burst energy.

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

We study the effect of the neutron star spin -- kick velocity alignment observed in young radio pulsars on the coalescence rate of binary neutron stars. Two scenarios of the neutron star formation are considered: when the kick is always present and when it is small or absent if a neutron star is formed in a binary system due to electron-capture degenerate core collapse. The effect is shown to be especially strong for large kick amplitudes and tight alignments, reducing the expected galactic rate of binary neutron star coalescences compared to calculations with randomly directed kicks. The spin-kick correlation also leads to a much narrower NS spin-orbit misalignment.

We present here an updated description of the "Scenario Machine" code. This tool is used to carry out a population synthesis of binary stars. Previous version of the description can be found at http://xray.sai.msu.ru/~mystery//articles/review/contents.html

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 study the effect of the kick velocity -- neutron star spin alignment observed in young radio pulsars on the coalescence rate of binary neutron stars. The effect is shown to be especially strong for large kick amplitudes and tight alignments, reducing the galactic rate of binary neutron star coalescences up to an order of magnitude with respect to the rates calculated for random kicks. The spin-kick correlation also leads to much narrower NS spin-orbit misalignments compared to random kicks.

We compute the probability that intermediate mass black holes (IMBHs) capture companions due to dynamical interactions and become accreting sources, and explore the possibility that the accreting IMBHs would appear as ultra-luminous X-ray sources (ULXs). We focus on IMBHs originating from low-metallicity Population III stars. Two channels of IMBH formation are considered: from primordial halos in the framework of hierarchical clustering, and from non-mixed, zero-metallicity primeval gas in galactic discs. IMBHs can form binary systems due to tidal captures of single stars and exchange interactions with existing binary systems in galactic discs. We find that neither formation mechanism of the accreting IMBH binary is able to provide enough sources to explain the observed population of ULXs. Even at sub-ULX luminosity, the total number of accreting IMBHs with $L> 10^{36}$ erg s$^{-1}$ with dynamically captured companions is found to be $< 0.01$ per galaxy.

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

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

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

We review the formation and evolution of compact binary stars consisting of white dwarfs (WDs), neutron stars (NSs), and black holes (BHs). Binary NSs and BHs are thought to be the primary astrophysical sources of gravitational waves (GWs) within the frequency band of ground-based detectors, while compact binaries of WDs are important sources of GWs at lower frequencies to be covered by space interferometers (LISA). Major uncertainties in the current understanding of properties of NSs and BHs most relevant to the GW studies are discussed, including the treatment of the natal kicks which compact stellar remnants acquire during the core collapse of massive stars and the common envelope phase of binary evolution. We discuss the coalescence rates of binary NSs and BHs and prospects for their detections, the formation and evolution of binary WDs and their observational manifestations. Special attention is given to AM CVn-stars -- compact binaries in which the Roche lobe is filled by another WD or a low-mass partially degenerate helium-star, as these stars are thought to be the best LISA verification binary GW sources.

Hard X-ray INT observations of SS 433 carried out during 2003-2005 years with an analysis of precessional and orbital variability is presented. The width of X-ray eclipse in the 25-50 keV range at the precessional phase $\psi=0.1$ (accretion disk is open to observer) is higher than that in the Ginga 18.4-27.6 keV range. This fact suggests existance the presence of hot extended corona around the supercritical accretion disk. Spectrum of hard X-rays in the range 10-200 keV does not change with the precessional phase which also suggests that hard X-ray flux is generated in the hot extended corona around the accretion disk. The parameters of this hot corona are: kT=23-25 keV, \tau = 1.8-2.8. Mass ratio estimated from the analysis of the ingress part of the eclipse light curve is in the range q=m_x/m_v=0.3-0.5.

We present our analysis of the extensive monitoring of SS433 by the RXTE observatory collected over the period 1996-2005. The difference between energy spectra taken at different precessional and orbital phases shows the presence of strong photoabsorption (N_H>10^{23}cm^{-2}) near the optical star, probably due to its powerful, dense wind. Therefore the size of the secondary deduced from analysis of X-ray orbital eclipses might be significantly larger than its Roche lobe size, which must be taken into account when evaluating the mass ratio from analysis of X-ray eclipses. Assuming that a precessing accretion disk is geometrically thick, we recover the temperature profile in the X-ray emitting jet that best fits the observed precessional variations in the X-ray emission temperature. The hottest visible part of the X-ray jet is located at a distance of l_0/a~0.06-0.09, or ~2-3*10^{11}cm from the central compact object, and has a temperature of about T_{max}~30 keV. We discovered appreciable orbital X-ray eclipses at the ``crossover'' precessional phases (jets are in the plane of the sky, disk is edge-on), which under model assumptions put a lower limit on the size of the optical component R/a>0.5 and an upper limit on a mass ratio of binary companions q=M_x/M_{opt}<0.3-0.35, if the X-ray opaque size of the star is not larger than 1.2R_{Roche, secondary}.

We present analysis of extensive monitoring of SS433 by the RXTE observatory collected over the period 1996-2005. The difference between energy spectra taken at different precessional and orbital phases shows the presence of a strong photoabsorption near the optical star, probably due to its powerful dense wind. Assuming that a precessing accretion disk is thick, we recover the temperature profile in the X-ray emitting jet that best fits the observed precessional variations of the X-ray emission temperature. The hottest visible part of the X-ray jet is located at a distance of $l_0/a\sim0.06-0.09$, or $\sim2-3\times10^{11}$cm from the central compact object and has a temperature of about $T_{\rm max}\sim30$ keV. We discovered appreciable orbital X-ray eclipses at the ``crossover'' precessional phases (jets are in the plane of the sky, disk is edge-on) which put a lower limit on the size of the optical component $R/a\ga0.5$ and an upper limit on a mass ratio of binary companions $q=M_{\rm x}/M_{\rm opt}\la0.3-0.35$. The size of the eclipsing region can be larger than secondary's Roche lobe because of substantial photoabsorption by dense stellar wind. This must be taken into account when evaluating the mass ratio from analysis of X-ray eclipses.

Population synthesis is used to model the number of neutron stars in globular clusters that are observed as LMXBs and millisecond PSRs. The dynamical interaction between binary and single stars in a GC are assumed to take place with a permanently replenished "background" of single stars whose density distribution keeps track with the cluster evolution as a whole and evolution of single stars. We use the hypothesis (Podsiadlowski et al) that NS forming in binary systems from components with initial masses \sim 8-12 M_\odot during the electron-capture collapse of the degenerate O-Ne-Mg core do not acquire a high space velocities (kicks). The remaining NSs (i.e. from single stars with M>8 M_\odot or binary comonents with M>12 M_\odot) are assumed to be born with high kicks, as found from obsrevations of single pulsars (Hobbs et al. 2005). Under this assumption, a sizeable fraction of NSs remain in GCs (about 1000 NSs in a GC with a mass of 5\times 10^5 M_\odot). The number of ms PSRs formed in the cluster via accretion spin-up in binaries is then about 10, which is consistent with observations. Our modelling reproduces the observed shape of the X-ray luminosity function for accreting NSs in binaries with normal and degenerate components and the distribution of spin periods of ms PSRs in GCs under the assumption of accretion-driven magnetic field decay of NSs up to a bottom value of 10^8 G. The number of LMXBs and ms PSRs dynamically expelling from GCs is also calculated.

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

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

We present broad band power spectra of variations of SS433 in radio, optical and X-ray spectral bands. We show that at frequencies lower than 1e-5 Hz the source demonstrates the same variability pattern in all these bands. The broad band power spectrum can be fitted by one power law down to frequencies ~1e-7 Hz with flattening afterwards. Such a flattening means that on time scales longer than ~1e7 sec the source variability becomes uncorrelated. This naturally leads to the appearance of quasi-poissonian flares in the source light curve, which have been regularly observed in radio and optical spectral bands. The radio flux power spectrum appears to have a second break at Fourier frequencies ~1e-5 Hz which can be caused by the smearing of the intrinsic radio variability on timescale of the light-crossing time of the radio emitting region. We find a correlation of the radio and optical fluxes of SS433 and the radio flux is delayed by about ~2 days with respect to the optical one. Power spectra of optical and X-ray variabilities continue with the same power law from 1e-7 Hz up to ~0.01-0.05 Hz. The broad band power spectrum of SS433 can be interpreted in terms of self-similar accretion rate modulations in the accretion disk proposed by Lyubarskii (1997) and elaborated by Churazov et al. (2001). We discuss a viscous time-scale in the accretion disk of SS433 in implication to the observed broad band power spectrum.

X-ray luminosity function derived from observations of X-ray sources in galactic bulges can be explained by prinicipal evolutionary relations for mass accretion rate onto the compact object. The observed mean distribution of individual X-ray luminosities of galactic LMXB is satisfactorily described by a symmetric quasi-Lorentzian curve. The flux variance for bright sources is found to be proportional to the mean luminosity. Such a distribution does not change the slope of the power law luminosity function of the source population, which is expected from the dependence of the mass transfer rate on the mass of the Roche-lobe filling non-degenerate optical component.

Current understanding of core collapse and thermonuclear supernovae is reviewed. Recent progress in unveiling the nature of cosmic gamma-ray bursts (GRB) is discussed, with the focus on the apparent link of several GRBs with an energetic subclass of stellar explosions, type Ib/c core-collapse supernovae. This relation provides the strong case that the GRB phenomenon is connected with the final stages of massive star evolution and possibly with the formation of neutron stars and black holes.

In these proceedings, I discuss recent progress in understanding the nature of cosmic gamma-ray bursts (GRB), with the focus on the apparent relation of several GRBs with an energetic subclass of stellar explosions, type Ib/c core-collapse supernovae. This relation provides the strong case that the GRB phenomenon is connected with the final stages of massive star evolution and possibly with the formation of neutron stars and black holes. I speculate that intrinsically faint, apparently spherically symmetric nearby GRB 980425 and 031203 associated with bright hypernovae SN 1998bw and SN 2003lw, respectively, can signal the formation of a neutron star in the end of gravitational collapse, while the bulk of cosmological GRBs with a universal energy release of $\sim 10^{51}$ ergs in narrow-collimated jets are produced when a black hole is formed. In the former case, the energy source of GRB is the neutron star rotational energy; in the latter case the GRB energy is due to non-stationary accretion onto the black hole.

We present optical and NIR (near infrared) follow up observations of the GRB 001109 from 1 to 300 days after the burst. No transient emission was found at these wavelengths within this GRB's (Gamma Ray Burst) 50" radius BeppoSAX error box. Strong limits (3 sigma) are set with: R >~ 21, 10.2 hr after the GRB; I >~ 23, 11.4 hr after the GRB; H >~ 20.7, 9.9 hr after the GRB; and Ks >~ 20, 9.6 hours after the GRB. We discuss whether the radio source found in the GRB's error box (Taylor et al. 2000) might be related to the afterglow. We also present a multiwavelength study of a reddened starburst galaxy, found coincident with the potential radio and the X ray afterglow. We show that our strong I band upper limit makes of the GRB 001109 the darkest one localised by the BeppoSAX's NFI (Narrow Field Instrument), and it is one of the most constraining upper limits on GRB afterglows to date. Further to it, the implications of these observations in the context of dark GRBs are considered.

We present results of first simultaneous optical and X-ray observations of peculiar binary system SS433. For the first time, chaotic variability of SS433 in the optical spectral band (R band) on time scales as small as tens of seconds was detected. We find that the X-ray flux of SS433 is delayed with respect to the optical emission by approximately 80 sec. Such a delay can be interpreted as the travel time of mass accretion rate perturbations from the jet base to the observed X-ray emitting region. In this model, the length of the supercritical accretion disk funnel in SS433 is ~1e12 cm.

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

Analysis of 18-120 keV images of the Sagittarius Arm Tangent region (SATR) obtained by IBIS telescope onboard INTEGTRAL observatory during the spring of 2003 is performed. In the 18-60 keV energy range, 28 sources have been detected with a flux level above 1.4 mCrab. Of these sources, 16 were identified earlier as galactic X-ray binary systems, 3 as extragalactic objects, 2 as pulsars inside supernova remnants, and 7 has unknown nature. The analysis revealed the presence of three previously unknown sources. Fourteen sources show significant flux in the 60-120 keV energy range

We present the earliest BTA (SAO RAS 6-m telescope) spectroscopic observations of the GRB030329 optical transient (OT), which almost coincide in time with the first "break" (t ~ 0.5 day after the GRB) of the OT light curve. The BTA/MPFS (Multi Pupil Fiber Spectrograph) spectra are clearly not smooth: the OT spectra showed a continuum with several broad spectral features at about 4000, 4450, 5900AA. The OT spectrum showed some evolution starting from the first night after the burst, and the beginning of spectral changes are seen as early as ~10-12 hours after the GRB. The onset of the spectral changes for t < 1 day indicates that the contribution from Type Ic supernova (SN) into the OT optical flux can be detected earlier. So, summarizing our and other spectroscopic data on the OT we confirm the evolution of the OT spectrum and color pointed also by other authors (Hjorth et al. 2003; Matheson et al. 2003b). The early spectra of GRB030329/SN2003dh properties can be consistent with a shock moving into a stellar wind formed from the pre-SN. Such a behavior (similar to one near the UV shock breakout in SNe) can be explained by the existence of a dense matter in the immediate surroundings of GRB/SN progenitor.

We analyze the mass distribution of stellar black holes derived from the light and radial velocity curves of optical stars in close binary systems using dynamical methods. The systematic errors inherent in this approach are discussed. These are associated primarily with uncertainties in models for the contribution from gaseous structures to the optical brightness of the systems under consideration. The mass distribution is nearly flat in the range 4-15M_sun. This is compared with the mass distribution for black holes in massive close binaries, which can be manifest as ultraluminous X-ray sources (L_x > 10^39 erg/s) observed in other galaxies. If the X-ray luminosities of these objects correspond to the Eddington limit, the black-hole mass distribution should be described by a power law, which is incompatible with the flat shape derived dynamically from observations of close binaries in our Galaxy. One possible explanation of this discrepancy is the rapid evaporation of stellar-mass black holes predicted in recent multi-dimensional models of gravity. This hypothesis can be verifed by measuring the stellar black-hole mass spectrum or finding isolated or binary black holes with masses below 3M_sun.

The population synthesis method is used to study the possibility of explaining the appreciable fraction (20^+12_15%) of the intergalactic (no-host) type Ia supernovae observed in galaxy clusters (Gal-Yam ete al. 2003) by binary whote dwarf merginngs in the cores of globular clusters. In a typical globular cluster, the number of merging double white dwarfs is fount to be smaller than 10^{-13} per year per average cluster star during the entire evolution of the cluster, which is a factor of 3 higher than in a Milky-Way-type galaxy. From 5 to 30% of the merging white dwarfs are dynamically expelled from the cluster with barycenter velocities up to 150 km/s. SN Ia explosions during the mergers of binary white dwarfs in dense star clusters may account for \sim 1% of the total rate of SN Ia in the central parts of galaxy clusters if the baryon mass fraction in such star clusters is \sim 0.3%.

Time-dependent thermal effects should accompany standard non-thermal afterglows of GRB when gamma-rays pass through inhomogeneous surroundings of the GRB site. Thermal relaxation of an optically thin plasma is calculated using time-dependent collisional ionization of the plasma ion species. X-ray emission lines are similar to those found in the fading X-ray afterglow of GRB 011211. Thermal relaxation of clouds or shells around the GRB site could also contribute to the varying late optical GRB afterglows, such as in GRB021004 and GRB030329.

Observations of SS433 by INTEGRAL carried out in March -- May 2003 are presented. SS433 is evidently detected on the INTEGRAL images of the corresponding sky region in the energy bands 25-50 and 50-100 keV. The precessional variability of the hard X-ray flux is clearly seen. The X-ray eclipse caused by the binary orbital motion is also detected. A possible origin of the hard continuum is briefly discussed.

We discuss masses of stellar black holes found in binary systems and errors in their determination. The observed mass distribution has a broad shape within the range $4-16 M_\odot$ without visible concentration to some preferred value. On the other hand, the black hole mass funciton as inferred from observations of luminous X-ray sources in other galaxies exhibits a power-law form. We show that both dynamically obtained black hole mass function and derived from X-ray observations can be made consistent in the frame of the hypothesis of enhanced black hole evaporation on a RSII brane for reasonable values of the warp-factor of the Anti-de-Sitter bulk space.

Time-dependent thermal X-ray spectra are calculated from physically plausible conditions around GRB. It is shown that account for time-dependent ionization processes strongly affects the observed spectra of hot rarefied plasma. These calculations may provide an alternative explanation to the observed X-ray lines of early GRBs afterglows (such as GRB 011211). Our technique will allow one to obtain independent constraints on the GRB collimation angle and on the clumpiness of circumstellar matter.

The empirically determined universal power-law shape of X-ray luminosity function of high mass X-ray binaries in galaxies is explained by fundamental mass-luminosity and mass-radius relations for massive stars.

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

Limitations on the multi-dimensional TeV-scale quantum gravity model by Arkani-Hamed, Dimopoulos and Dvali (1998) and on model by Dvali, Gabadadze, Porrati (2000) are obtained from an analysis of gamma-ray bursts at cosmological distances (relativistic fireball model).

Maximum gravitational wave emission from binary stars in eccentric orbits occurs near the periastron passage. We show that for a stationary distribution of binary neutron stars in the Galaxy, several high-eccentricity systems with orbital periods in the range from tens of minutes to several days should exist that emit broad gravitational-wave pulses in the frequency range 1-100 mHz. The space interferometer LISA could register the pulsed signal from these system at a signal-to-noise ratio level $S/N>5\sqrt{5}$ in the frequency range $\sim 10^{-3}-10^{-1}$ Hz during one-year observational time. Some detection algorithms for such a signal are discussed.

The observed mass distribution for the compact remnants of massive stars (neutron stars and black holes) and its relationship to possible mechanisms for the ejection of the envelopes of type II and Ib/c supernovae is analyzed. The conclusion is drawn that this distribution can be obtained only by a magneto-rotational mechanism for the supernovae with sufficiently long time of the field amplification, and a soft equation of state for neutron stars with limiting masses \sim1.5-1.6M_\odot. Some consequences of this hypothesis are discussed.

Unresolved coalescing binary neutron stars in the Galaxy and beyond produce a stochastic gravitational wave background within the LISA frequency band, which can be potentially dangerous for possible detection of fundamentally more interesting relic cosmological backgrounds. Here we address the question what effects the unavoidable eccentricity of orbits of these systems should have on the properties of this backgrounds. In particularly, we examine starting from which frequency one-year observations with LISA can be secure from contamination by the noise produced by galactic binary neutron stars in eccentric orbits. We come to the conclusion that harmonics from such binaries do not contribute significantly above $\sim 10^{-3}$ Hz.

The recent progress in studies of gamma-ray bursts, their afterglows, and host galaxies is discussed. The emphasis is given to high-energy phenomena associated with gamma-ray burst explosions: high-energy cosmic rays, neutrinos, gravitational waves. We also show how the relativistic fireball model for GRBs can be used to constrain modern theories of large and infinite extra-dimensions. In particular, in the frame of 5D gravity with the Standard Model localized on 3D brane (Dvali et al. 2000), the very existence of relativistic fireballs of $\sim 10^{53}$ ergs puts the lower bound on the quantum gravity scale $\sim 0.1$ eV.

Stochastic gravitational wave background produced by a stationary coalescing population of binary neutron stars in the Galaxy is calculated. This background is found to constitute a confusion limit within the LISA frequency band up to a limiting frequency $\NUlim{}\sim 10^{-3}$ Hz, leaving the frequency window $\sim 10^{-3}$--$10^{-2}$ Hz open for the potential detection of cosmological stochastic gravitational waves and of signals involving massive black holes out to cosmological distances.

The first generation of long-baseline laser interferometric detectors of gravitational waves will start collecting data in 2001-2003. We carefully analyse their planned performance and compare it with the expected strengths of astrophysical sources. The scientific importance of the anticipated discovery of various gravitatinal wave signals and the reliability of theoretical predictions are taken into account in our analysis. We try to be conservative both in evaluating the theoretical uncertainties about a source and the prospects of its detection. After having considered many possible sources, we place our emphasis on (1) inspiraling binaries consisting of stellar mass black holes and (2) relic gravitational waves. We draw the conclusion that inspiraling binary black holes are likely to be detected first by the initial ground-based interferometers. We estimate that the initial interferometers will see 2-3 events per year from black hole binaries with component masses 10-15M_\odot, with a signal-to-noise ratio of around 2-3, in each of a network of detectors consisting of GEO, VIRGO and the two LIGOs. It appears that other possible sources, including coalescing neutron stars, are unlikely to be detected by the initial instruments. We also argue that relic gravitational waves may be discovered by the space-based interferometers in the frequency interval 2x10^{-3}-10^{-2} Hz, at the signal-to-noise ratio level around 3.

Observed properties of GRBs, WR stars and their CO-cores in the end of evolution are analyzed. A possible bimodality of the observed GRB energy distribution ($10^{48}$ erg for GRB9809425; $3\times 10^{51}\div 2\times 10^{54}$ erg for others) is in accord with the bimodal mass distribution of relativistic objects ($M_{NS}=(1.35\pm 0.15) M_\odot$; $M_{BH}=(4\div 15) M_\odot$). The peculiarity of SN1998bw can be related to the rotation of the collapsing CO-core. The expected galactic collapse rate of CO-cores of most compact WR stars of type WO is $\sim 10^{-5}$ per year, only by one and a half order of magnitude higher than the GRB rate. The allowance for a gamma-ray beaming or random outcome of the CO-core collapse due to some instabilities brings this rate in accordance GRB rate. We argue that WR stars (most probably, of type WO) can be considered as progenitors of cosmic gamma-ray bursts. Two types of GRBs are predicted in correspondence with the bimodal mass distribution of the relativistic objects. Three types of optical afterglows should appear depending on which CO-core is collapsing: of a single WR star, of a WR star in a WR+O or a hypothetic WR+(A-M) binary system. We briefly discuss a model of GRB as a transient phenomenon occurring at early stages of galactic evolution ($z>1$), when very massive ($M>100 M_\odot$) low-metallicity stars could form. WR-galaxies can be most probable candidates for GRB host galaxies.

After presenting a short history of gamma-ray burst (GRB) studies, we discuss the current constraints on GRB models which follows from astronomical observations. We concentrate on the energetics of the GRBs with known redshifts and the association of the well-localized GRBs with star forming regions in remote galaxies. We also discuss implications of the observed GRB rate. Arguments are given favoring possible connection of GRBs with core collapse of massive Wolf-Rayet stars. The possibility of GRB to be a transient phenomenon in the early history of galactic star formation related to evolution of very massive metal-free stars is briefly outlined.

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

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

Angular fluctuations of stochastic gravitational wave backgrounds (GWB) produced by extragalactic astrophysical sources are calculated. The angular properties of such backgrounds are determined by the large scale structure of Universe (galaxy clustering). The evolution of star formation rate with redshift is taken into account. Fluctuations of the metric strain amplitude associated with such noises at angular scales of about one degree are found to be of order 5-20% slowly growing toward smaller angular scales. This feature can be potentially used to separate astrophysical GWB from cosmological ones in future experiments.

We show that the distribution of observed energies of GRB with known redshifts can be explained by the hypothesis of the standard energy release E_0=5 10^51 ergs. Two situations are possible, either the beaming angle differs from burst to burst, or there is a universal emssion diagram in each burst, the observed difference being due to different viewing angles.

The formation and evolution of binary black holes (BH) is studied using the modern evolutionary scenario for very massive stars with high mass loss (Vanbeveren et al. 1998). Main sequence stars with masses $M>35 M_\odot$ are assumed to form a BH in the end of their nuclear evolution. The mass of BH formed is parametrized as $M_{bh}=k_{bh}\times M_{SN}$, where $M_{SN}$ is the mass of the pre-supernova star taken from evolutionary calculations, $k_{bh}\le 1$. The possibility is explored that a newly formed BH acquires a kick velocity 0-250 km/s. Binary BH are found to merge within the Hubble time at an appreciable rate only for non-zero kick velocities. We calculate the galactic merging rates of binary BH systems, their detection rate by the initial laser interferometers, and the distributions of merging binary BH over orbital eccentricities at different frequencies. The distribution of angles between BH spins and the orbital angular momentum is also presented.

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

In this article we use "Scenario Machine" - the population synthesis simulator- to calculate the evolution of the populations of the selected types of X-ray sources after the starformation burst with the total mass in binaries $1.5\cdot10^6 M_{\odot}$ during the first 20 Myr after the burst. Sources of the four types were calculated: transient sources- accreting neutron stars with Be- stars; accreting neutron stars in pair with supergiants; Cyg X-1-like sources- black holes with supergiants; superaccreting black holes. We used two values of the $\alpha - coefficient$ in the mass-function: 2.35 (Salpeter's function) and 1.01 ("flat spectrum"). The calculations were made for two values of the upper limit of the mass-function: 120 and 30 $M_{\odot}$. For the flat spectrum the number of sources of all types significantly increased. Decreasing of the upper mass limit below the critical mass of a black hole formation increase (for the "flat spectrum") the number of transient sources with neutron stars up to $\approx 300$. We give approximating formulae for the time dependence of source numbers.

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

We calculate stochastic gravitational wave background produced by extragalactic merging binary white dwarfs at the LISA frequencies $10^{-3}-10^{-2}$ Hz with account of a strong evolution of global star formation rate in the Universe recently established observationally. We show that for the observed global star formation history and modern cosmological models the extragalactic background is an order of magnitude smaller than the mean Galactic value. An early star formation burst at high redshifts can bring it at a higher level but still a few times lower than the mean Galactic one.

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

Gravitational wave noise associated with unresolved binary stars in the Galaxy is studied with the special aim of determining the upper frequency at which it stops to contribute at the rms noise level of the proposed space-born interferometer (LISA). The upper limit to this background is derived from the statistics of SN Ia explosions, part of which can be triggered by binary white dwarf coalescences. The upper limiting frequency at which binary stochastic noise crosses LISA rms sensitivity is found to lie within the range 0.03-0.07 Hz, depending on the galactic binary white dwarf coalescence rate. To be reliably detectable by LISA, the energy density of relic cosmological background per logarithmic frequency interval should be Omega_{GW}h_{100}^2>10^{-8} at f>0.03 Hz.

We use "Scenario Machine" -- the population synthesis simulator -- to calculate the evolution of populations of selected types of X-ray sources after a starformation burst with the total mass in binaries (1--1.5) \cdot 10^6 M_{\odot}$ during the first 20 Myr after a burst. We present here the results of two sets of runs of the program. In the first set we examined the following types of close binaries: transient sources-- neutron stars with Be- stars; ``X-ray pulsars''-- neutron stars in pairs with supergiants; Cyg X-1-like sources-- black holes with supergiants; ``SS443-like sources''-- superaccreting black holes. We used two values of the exponent $\alpha$ in the initial mass function: 2.35 (Salpeter's function) and 1 (``flat spectrum''). The calculations were made for the folowing values of the upper limit of the mass function: 120 and 30 $M_{\odot}$. For the ``flat spectrum'', suggested in (Contini et al, 1995), the number of sources of all types significantly increased. With the ``flat spectrum'' and with the upper mass limit 120 $M_{\odot}$ we obtained hundreds of sources of all calculated types. Decreasing of the upper mass limit below the critical mass of a black hole formation increase the number of transient sources with neutron stars up to $\approx 300$. In the second set we examined the evolution of of 12 other types of X-ray sources for $\alpha=1$, $\alpha=1.35$ and $\alpha=2.35$ and for three upper mass limits: 120 $M_{\odot}$, 60 $M_{\odot}$, 40 $M_{\odot}$ (see Perez-Olia & Colina 1995 for the reasons for such upper limits) on the same time scale 20 Myr after a star formation burst.

An energy deposition of \sim 10^{50} ergs into the exterior 10^{-3} M_\sun layers of a red giant is calculated to produce an optical phenomenon similar to afterglows of gamma-ray bursts (GRB) recently observed. This model can be realized if a GRB is generated by some mechanism in a close binary system. In contrast to a ``hypernova'' scenario for GRB recently proposed by Paczy\'nski (1997), this model does not require a huge kinetic energy of the expanding shell to explain optical afterglows of GRB.

Stochastic gravitational waves (GW) associated with unresolved astrophysical sources at frequency bands of the ongoing GW interferometers LIGO/VIRGO and LISA are studied. We show that GW noise from rotating galactic neutron stars with low magnetic fields may reach the advanced LIGO sensitivity level at frequency $f\sim 100$ Hz. Within LISA frequency band (10^{-4}-10^{-1} Hz), the GW background from galactic binary stars is shown to mainly contribute up to a frequency of 3.10^{-2} Hz, depending on the galactic rate of binary white dwarf mergers. To be detectable by LISA, relic GW backgrounds should be as high as \Omega_{GW}h_{100}^2>10^{-8} at 10^{-2} Hz.

Assuming binary neutron star mergings as a standard-candle model for GRBs, an independent estimate is obtained for the redshift of GRB970228 and GRB970508 with optical counterparts, using mean statistical properties of GRBs observed by BATSE. We derive z=0.7\pm 0.1 and z=1.9\pm 0.1 for GRB970228 and GRB970508 respectively, depending on the power-law index of the GRB spectrum s=-1.1\pm 0.3 and the value of the redshift of the initial star formation z_*=3-10 in a flat \Omega=1 Universe with a cosmological term \Lambda=0.7.

Merging rates of compact binaries (double neutron stars or black holes) are calculated based on the modern concept of binary stellar evolution. It is found that the initial laser interferometers with an rms-sensitivity of $10^{-21}$ at the frequency 100 Hz can detect 10-700 black holes and only \sim 1 neutron star coalescences in a 1-year integration time. Implications of the evolutionary effects to the cosmological origin of GRB are also discussed.

Galactic binary pulsars with different companion types (OB-star, white dwarf (WD), neutron star (NS), black hole (BH), or planet) assuming various phenomenological distributions for kick velocities of newborn NS are studied. The merging rates of relativistic binary stars (NS+NS, NS+BH, BH+BH) dependence on the kick velocity is specially adresed. The BH+BH merging is shown to occur, depending on parameters of BH formation, at a rate of one per 200,000 -- 500,000 years in a Milky Way type galaxy. The NS+NS merging rate $R_{ns}$ is found to be 1 per \sim 3,000 years for zero recoil, and decrease to one per 10,000 years even for the highest kick velocities of 400 km/s. That the merging rates derived from evolutionary calculations are by two order of magnitude higher than those based on binary pulsar statistics only, is suggested to be due to the fact that the observable binary pulsars in pairs with NS form only a fraction of the total number of binary NS systems. The expected detection rate of binary BH by a LIGO-type gravitational wave detector comparable with and even higher than the binary NS merging rate for a wide range of parameters.

Study of gravitational-radiation induced merging rates of relativistic binary stars (double neutron stars; neutron star + black hole; double black holes) shows that the first-generation gravitational wave interferometers with an rms-sensitivity of $10^{-21}$ at frequency 100 Hz can detect 10-700 black hole and only $\sim 1$ neutron star coalescences in a 1-year integration time in a wide range of stellar evolution parameters. It is notable that modern concepts of stellar evolution predict that the first detection of gravitational wave will independently discover black holes.

Assuming the observed pulsar velocities to originate during asymmetric collapse of stellar cores, we compute the amplitude of gravitational waves emitted during type II and Ib supernova explosions and their detection rate from within a distance of 30 Mpc. At the rms-level of advanced laser interferometers h ~ 1.0e-22 at frequencies 300-1000 Hz the expected rate is about 1 per year.

Based on evolutionary scenarios for binary stellar evolution we study the merging rates of relativistic binary stars (NS+NS, NS+BH, BH+BH) under different assumptions of BH formation. We find the BH+BH merging rate in the range one per 200,000 -- 500,000 year in a Milky-Way type galaxy, while the NS+NS merging rate $R_{ns}$ is approximately 10 times as high, which means that the expected event rate even for high mean kick velocities of NS up to 400 km/s is at least 30-50 binary NS mergings per year from within a distance of 200 Mpc. As typical BH is formed with a mass 3-10 times the NS mass (assumed 1.4 M$_\odot$), the rates obtained imply that the expected detection rate of binary BH by a LIGO-type gravitational wave detector is 10-100 times higher than the binary NS merging rate for a wide range of parameters.

We present a WWW-version of the {\it Scenario Machine} - a computer code designed to calculate the evolution of close binary stellar systems. The Internet users can directly access to the code and calculate binary evolutionary tracks with parameters at the user's will. The program is running on the {\it Pentium} server of the Division of the Relativistic Astrophysics of the Sternberg Astronimical Institute (http://xray.sai.msu.su/ ). The results are presented both in the form of tables and graphic diagrams. The work is always in progress. More possibilities for Internet users are intended to become available in the near future.

We model the evolution of the selected types of X-ray sources (both with BHs and NSs) for the star formation burst at the Galactic center. For the currently assumed starburst age of 4-7 millions years, our results are consistent with the recent {\it GRANAT} X-ray observations of the Galactic center and predict a substantial number of X-ray binaries with black holes. In addition we made a very simple model of the spatial distribution of the X-ray sources, which is also consistent with the {\it GRANAT} observations. The results can be applied to other starformation regions.

X-ray binary population at the Galactic center is proposed here to form as a result of a starburst, which is synthesized by using the ``Scenario Machine'' code of binary star evolution. For the currently assumed starburst age of 4-7 millions years, our results are consistent with the recent {\it GRANAT} X-ray observations of the Galactic center and predict a substantial number of X-ray binaries with black holes.

Calculation of gravitational radiation during binary inspiral leading to possible formation of a Thorne-Zytkow (TZ) object (a neutron star inside a supergiant core) is performed. The calculations were done for polytropic density distributions with different indeces $n$, as well as for realistic models of supergiants. A maximum frequency of the emitted gravitational waves during this process was found to range from a few to 300 initial keplerian orbital frequensies, that is from $10^{-5}$ to $\sim 0.1$ Hz, depending on the model. A dimensionless strain metric amplitude can reach $h\sim 10^{-23.5}$ for a source lying 10 kpc away from the Sun. We conclude that TZ objects forming at a rate of 1 per 500 yrs in the Galaxy could be potential astrophysical targets for space laser intereferometers. Analysis of gravitational waveforms emitted during TZ-object formation could bring a unique information about stellar structure.

Evolution of the coalescence rate of double neutron stars (NS) and neutron star -- black hole (BH) binaries are computed for model galaxies with different star formation rates. Assuming gamma-ray bursts (GRB) to originate from NS+NS or NS+BH merging in distant galaxies, theoretical logN--logS distributions and <V/V_max> tests of gamma-ray bursts (GRB) are calculated for the first time taking the computed merging rates into account. We use a flat cosmological model (Omega=1) with different values of the cosmological constant Lambda and under various assumptions about the star formation history in galaxies. The calculated source evolution predicts a 5-10 times increase of the source statistics at count rates 3-10 times lower than the exising BATSE sensitivity limit. The most important parameter in fitting the 2nd BATSE catalogue is the initial redshift of star formation, which is found to be z_*=2-5 depending on a poorly determined average spectral index of GRB.