Variability in high-mass X-ray binaries

Strongly magnetized, accreting neutron stars show periodic and aperiodic variability over a wide range of time scales. By obtaining spectral and timing information on these different time scales, we can have a closer look into the physics of accretion close to the neutron star and the properties of the accreted material. One of the most prominent time scales is the strong pulsation, i.e., the rotation period of the neutron star itself. Over one rotation, our view of the accretion column and the X-ray producing region changes significantly. This allows us to sample different physical conditions within the column but at the same time requires that we have viewing-angle-resolved models to properly describe them. In wind-fed high-mass X-ray binaries, the main source of aperiodic variability is the clumpy stellar wind, which leads to changes in the accretion rate (i.e., luminosity) as well as absorption column. This variability allows us to study the behavior of the accretion column as a function of luminosity, as well as to investigate the structure and physical properties of the wind, which we can compare to winds in isolated stars.Comment: 6 pages, 4 figures, accepted for publication in Astronomische Nachrichten (proceedings of the XMM-Newton Workshop 2019

Probing large-scale wind structures in Vela X-1 using off-states with INTEGRAL

Vela X-1 is the prototype of the class of wind-fed accreting pulsars in high mass X-ray binaries hosting a supergiant donor. We have analyzed in a systematic way ten years of INTEGRAL data of Vela X-1 (22-50 keV) and we found that when outside the X-ray eclipse, the source undergoes several luminosity drops where the hard X-rays luminosity goes below 3x10^35 erg/s, becoming undetected by INTEGRAL. These drops in the X-ray flux are usually referred to as "off-states" in the literature. We have investigated the distribution of these off-states along the Vela X-1 ~8.9 d orbit, finding that their orbital occurrence displays an asymmetric distribution, with a higher probability to observe an off-state near the pre-eclipse than during the post-eclipse. This asymmetry can be explained by scattering of hard X-rays in a region of ionized wind, able to reduce the source hard X-ray brightness preferentially near eclipse ingress. We associate this ionized large-scale wind structure with the photoionization wake produced by the interaction of the supergiant wind with the X-ray emission from the neutron star. We emphasize that this observational result could be obtained thanks to the accumulation of a decade of INTEGRAL data, with observations covering the whole orbit several times, allowing us to detect an asymmetric pattern in the orbital distribution of off-states in Vela X-1.Comment: Accepted for publication in Monthly Notices of the Royal Astronomical Society (5 pages, 3 figures). A few typos fixed to match the published versio

Coupling hydrodynamics with comoving frame radiative transfer: II. Stellar wind stratification in the high-mass X-ray binary Vela X-1

CONTEXT: Vela X-1, a prototypical high mass X-ray binary (HMXB), hosts a neutron star (NS) in a close orbit around an early-B supergiant donor star. Accretion of the donor star's wind onto the NS powers its strong X-ray luminosity. To understand the physics of HMXBs, detailed knowledge about the donor star winds is required. AIMS: To gain a realistic picture of the donor star in Vela X-1, we constructed a hydrodynamically consistent atmosphere model describing the wind stratification while properly reproducing the observed donor spectrum. To investigate how X-ray illumination affects the stellar wind, we calculated additional models for different X-ray luminosity regimes. METHODS: We use the recently updated version of the PoWR code to consistently solve the hydrodynamic equation together with the statistical equations and the radiative transfer. RESULTS: The wind flow in Vela X-1 is driven by ions from various elements with Fe III and S III leading in the outer wind. The model-predicted mass-loss rate is in line with earlier empirical studies. The mass-loss rate is almost unaffected by the presence of the accreting NS in the wind. The terminal wind velocity is confirmed at $v_\infty \approx 600$ km/s. On the other hand, the wind velocity in the inner region where the NS is located is only $\approx 100$ km/s, which is not expected on the basis of a standard $\beta$-velocity law. In models with an enhanced level of X-rays, the velocity field in the outer wind can be altered. If the X-ray flux is too high, the acceleration breaks down because the ionization increases. CONCLUSIONS: Accounting for radiation hydrodynamics, our Vela X-1 donor atmosphere model reveals a low wind speed at the NS location, and it provides quantitative information on wind driving in this important HMXB.Comment: 19 pages, 10 figures, accepted for publication in Astronomy & Astrophysic

A large spin-up rate measured with INTEGRAL in the High Mass X-ray Binary Pulsar SAXJ2103.5+4545

The High Mass X-ray Binary Pulsar SAXJ2103.5+4545 has been observed with INTEGRAL several times during the last outburst in 2002-2004. We report a comprehensive study of all INTEGRAL observations, allowing a study of the pulse period evolution during the recent outburst. We measured a very rapid spin-up episode, lasting 130days, which decreased the pulse period by 1.8s. The spin-up rate, pdot=-1.5e-7 s/s, is the largest ever measured for SAXJ2103.5+4545, and it is among the fastest for an accreting pulsar. The pulse profile shows evidence for temporal variability, apparently not related to the source flux or to the orbital phase. The X-ray spectrum is hard and there is significant emission up to 150keV. A new derivation of the orbital period, based on RXTE data, is also reported.Comment: 8 pages, 7 figures, accepted for publication in A&

Study of the cyclotron feature in MXB 0656-072

We have monitored a type II outburst of the Be/X-ray binary MXB 0656−072 in a series of pointed RXTE observations during October through December 2003. The source spectrum shows a cyclotron resonance scattering feature at 32.8 +0.5 −0.4 keV, corresponding to a magnetic field strength of 3.67 +0.06 −0.04 × 10 12 G and is stable through the outburst and over the pulsar spin phase. The pulsar, with an average pulse period of 160.4 ± 0.4s,shows a spin-up of 0.45 s over the duration of the outburst. From optical data, the source distance is estimated to be 3.9 ± 0.1 kpc and this is used to estimate the X-ray luminosity and a theoretical prediction of the pulsar spin-up during the outburst

Formation of phase lags at the cyclotron energies in the pulse profiles of magnetized, accreting neutron stars

Context: Accretion-powered X-ray pulsars show highly energy-dependent and complex pulse-profile morphologies. Significant deviations from the average pulse profile can appear, in particular close to the cyclotron line energies. These deviations can be described as energy-dependent phase lags, that is, as energy-dependent shifts of main features in the pulse profile. Aims: Using a numerical study we explore the effect of cyclotron resonant scattering on observable, energy-resolved pulse profiles. Methods: We generated the observable emission as a function of spin phase, using Monte Carlo simulations for cyclotron resonant scattering and a numerical ray-tracing routine accounting for general relativistic light-bending effects on the intrinsic emission from the accretion columns. Results: We find strong changes in the pulse profile coincident with the cyclotron line energies. Features in the pulse profile vary strongly with respect to the average pulse profile with the observing geometry and shift and smear out in energy additionally when assuming a non-static plasma. Conclusions: We demonstrate how phase lags at the cyclotron energies arise as a consequence of the effects of angular redistribution of X-rays by cyclotron resonance scattering in a strong magnetic field combined with relativistic effects. We also show that phase lags are strongly dependent on the accretion geometry. These intrinsic effects will in principle allow us to constrain a system's accretion geometry.Comment: 4 pages, 4 figures; updated reference lis