Search

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., Comment: 48 pages, 4 figures, accepted for publication in MNRAS

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., Comment: 7 pages, 2 figures, Proceedings of the 8th INTEGRAL Workshop, Dublin; http://pos.sissa.it/archive/conferences/115/015/INTEGRAL%202010_015.pdf

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., Comment: 7 pages, 3 figures; Proceedings of the 8th INTEGRAL Workshop "The Restless Gamma-ray Universe" (INTEGRAL 2010). September 27-30 2010. Dublin Castle, Dublin, Ireland. Published online at http://pos.sissa.it/cgi-bin/reader/conf.cgi?confid=115., id.16

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., Comment: 6 pages, 1 figure. Accepted for publication in A&A

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., Comment: 27 pages, 16 colour figures, this version is translated by authors, so it differs from that, which is published in Astronomy Letters

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., Comment: 6 pages, 2 figures, Accepted for publication in Astrophysics & Space Science. Largely extended version of ArXiv:0912.5216

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., Comment: 6 pages, 3 figures. Accepted for publication on PoS (contribution PoS(extremesky2009)054), proceedings of "The Extreme sky: Sampling the Universe above 10 keV", held in Otranto (Italy) in October 2009

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., Comment: This paper has been withdrawn by the authors due to the appearance of the largely extended version in Astrophysics & Space Science and corresponding publication policy. The new arXiv identifier is 1004.5115.

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., Comment: Proceedings of 7th INTEGRAL Workshop in Copenhagen

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., Comment: 2 pages, 1 figure

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., Comment: Accepted for publication in Astronomy and Astrophysics

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., Comment: v.2, LATEX, 25 pages, inc. 7 ps figures, Astron. Lett., in press. Typos corrected, reference added

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., Comment: 9 pages, 1 figure, generally matches the MNRAS version

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., Comment: 7 pages, 5 figures, accepted for publication in Astronomy & Astrophysics

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., Comment: 23 pages, 15 figures, MNRAS accepted

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., Comment: MNRAS accepted. 13 pages, 15 figures, 1 table

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., Comment: 20 pages, 8 figures; Submitted to MNRAS

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., Comment: 6 pages, 5 figures, accepted for publication by A&A

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., Comment: 12 pages, 10 figures, Proc. 7th INTEGRAL Workshop, Copenhagen, 8-11 September 2008

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)., Comment: 8 pages, Astron. Lett. in press

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., Comment: 9 pages, 10 figures, accepted for publication in A&A on 29/01/2008

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., Comment: 4 pages, 4 figures, accepted for publication in A&A Letters. To be published in parallel to Postnov et al. 2008

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., Comment: 4 pages, 1 figure, accepted for publication in A&A Letters. To be published in parallel to Caballero et al. 2008

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., Comment: 7 pages, 5 figures, accepted for publiction in MNRAS

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., Comment: 3 pages, no figures, Revised version with the new estimate of mERB rate by Lorimer et al. strengthening the discrepancy with the binary neutron star merging rate. References added, typos corrected

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., Comment: To be published in the proceedings of the 6th INTEGRAL Workshop "The Obscured Universe" (3-7 July 2006, Moscow)

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, Comment: 32 pages, 3 figures. Corrected typos

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., Comment: LaTeX, 4 pages, 3 figures; talk given at XLIId Rencontres de Moriond "Gravitaional Waves and Experimental Gravitation", March 11-18, La Thuile; to be published in the Proceedings; references added

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., Comment: 21 pages, 5 figures. Accepted to MNRAS

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., Comment: 6 pages, 4 figures, In: The transient Milky Way: a perspective for MIRAX, AIP Conf. Proc. 840 (2006) 65-70

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., Comment: 4 pages, 2 figures, accepted by A&A Letters

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., Comment: Proceedings of the 363. WE-Heraeus Seminar on: Neutron Stars and Pulsars (Posters and contributed talks) Physikzentrum Bad Honnef, Germany, May.14-19, 2006, eds. W.Becker, H.H.Huang, MPE Report 291, pp.212-215

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., Comment: 5 pages, 10 figures

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}., Comment: 9 pages, 9 figures, accepted for publication in Astronomy and Astrophysics

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., Comment: 9 pages, 8 figures, accepted for publication in PAZH

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., Comment: Published in Proceedings of the 13th Young Scientists' Conference on Astronomy and Space Physics, held in Kyiv, Ukraine, April 25-29, 2006, Eds.: Golovin, A.; Ivashchenko, G.; Simon, A

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., Comment: LATEX, 21 pages, 8 gif figures, Astronomy Letters, in press

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., Comment: 16 pages, 13 figures, accepted for publication in A&A

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., Comment: 15 pages, AA style, 17 figures, accepted for publication in AA; replaced to match the accepted version

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., Comment: 8 pages, 2 figures. Submitted to A&A

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., Comment: 15 pages, latex, 5 figures, to appear in Astronomy Letters

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., Comment: rvised version: typos corrected, important references added. 11 pages, no figures, LaTeX, to be published in Proc. Int. Workshop QUARKS-2004, Ed. V.A.Rubakov et al

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., Comment: 7 pages, 5 figures, accepted for publication in A&A

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., Comment: 4 pages, 3 figures, Submitted to Astronomy and Astrophysics Letters

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$, Comment: 8 pages, 11 figures, Proceedings of the V INTEGRAL Workshop, Munich 16-20 February 2004. ESA SP-552

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, Comment: 5 pages, 2 figures. Accepted for publication in Astronomy Letters

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., Comment: 11 pages, 2 figures