We discovered 2.8 s pulsations in the X-ray emission of the ultraluminous X-ray source (ULX) M51 ULX-7 within the UNSEeN project, which was designed to hunt for new pulsating ULXs (PULXs) with XMM-Newton. The pulse shape is sinusoidal and large variations of its amplitude were observed even within single exposures (pulsed fraction from less than 5% to 20%). M51 ULX-7 is a variable source, generally observed at an X-ray luminosity between $10^{39}$ and $10^{40}$ erg s$^{-1}$, located in the outskirts of the spiral galaxy M51a at a distance of 8.6 Mpc. According to our analysis, the X-ray pulsar orbits in a 2-d binary with a projected semi-major axis $a_\mathrm{X} \sin i \simeq$ 28 lt-s. For a neutron star (NS) of 1.4 $M_{\odot}$, this implies a lower limit on the companion mass of 8 $M_{\odot}$, placing the system hosting M51 ULX-7 in the high-mass X-ray binary class. The barycentric pulse period decreased by $\simeq$0.4 ms in the 31 d spanned by our May -- June 2018 observations, corresponding to a spin-up rate $\dot{P} \simeq -1.5\times10^{-10}\text{s s}^{-1}$. In an archival 2005 XMM-Newton exposure, we measured a spin period of $\sim$3.3 s, indicating a secular spin-up of $\dot{P}_{\mathrm{sec}}\simeq -10^{-9}\text{ s s}^{-1}$, a value in the range of other known PULXs. Our findings suggest that the system consists of an OB giant and a moderately magnetic (dipole field component in the range $10^{12}$ G $\lesssim B_{\mathrm{dip}}\lesssim 10^{13}$G) accreting NS with weakly beamed emission ($1/12\lesssim b\lesssim1/4$).

We report the discovery of a third ULX in NGC 925 (ULX-3), detected in November 2017 by Chandra at a luminosity of $L_{\rm X} = (7.8\pm0.8)\times10^{39}$ erg s$^{-1}$. Examination of archival data for NGC 925 reveals that ULX-3 was detected by Swift at a similarly high luminosity in 2011, as well as by XMM-Newton in January 2017 at a much lower luminosity of $L_{\rm X} = (3.8\pm0.5)\times10^{38}$ erg s$^{-1}$. With an additional Chandra non-detection in 2005, this object demonstrates a high dynamic range of flux of factor >26. In its high-luminosity detections, ULX-3 exhibits a hard power-law spectrum with $\Gamma=1.6\pm0.1$, whereas the XMM-Newton detection is slightly softer, with $\Gamma=1.8^{+0.2}_{-0.1}$ and also well-fitted with a broadened disc model. The long-term light curve is sparsely covered and could be consistent either with the propeller effect or with a large-amplitude superorbital period, both of which are seen in ULXs, in particular those with neutron star accretors. Further systematic monitoring of ULX-3 will allow us to determine the mechanism by which ULX-3 undergoes its extreme variability and to better understand the accretion processes of ULXs.

We present a qualitative search for ultra-fast outflows (UFOs) in excess variance spectra of radio-quiet active galactic nuclei (AGN). We analyse 42 sources from the Tombesi et al. (2010) spectroscopic UFO detection sample, and an additional 22 different sources from the Kara et al. (2016) variability sample. A total of 58 sources have sufficient observational data from XMM-Newton EPIC-pn and variability for an excess variance spectrum to be calculated. We examine these spectra for peaks corresponding to variable blue-shifted H- and He-like ion absorption lines from UFOs. We find good evidence for such outflows in 28% of the AGN sample and weak evidence in a further 31%, meaning that $\sim$ 30-60% of the AGN sample hosts such UFOs. The mean and median blue-shifted velocity is found to be $\sim$ 0.14c and 0.12c, respectively. Current variability methods allow for a fast, model-independent determination of UFOs, however, further work needs to be undertaken to better characterize the statistical significance of the peaks in these spectra by more rigorous modelling. Detecting good evidence for variable UFO lines in a large number of sources also lays the groundwork for detailed analysis of the variability timescales of the absorbers. This will allow us to probe their densities and hence distances from the central super-massive black hole.

We present a detailed analysis of the spectral properties of the black hole transient GRS 1716-249, based on the archival Swift and NuSTAR observations taken during the outburst of this source in 2016-2017. The first six NuSTAR observations show that the source is in a canonical hard state, where the spectrum is dominated by a power-law continuum. The seventh NuSTAR observation is taken during the intermediate state where both a disc thermal component and a power-law continuum are shown. All of our observations show a broad emission line feature in the iron band and a Compton hump above 10 keV. We model the broad band spectra using a high density disc reflection model, where the soft X-ray emission in the hard state is interpreted as part of the disc reflection component. This model enables us to constrain the disc density parameter of GRS 1716-249 in the range of $10^{19}$-$10^{20}$ cm$^{-3}$. We only obtain an upper limit of the inner disc radius using high density disc reflection spectroscopy and the results indicate either a non-truncated disc or a slightly truncated disc with $R_{\rm in}<20r_{\rm g}$.

We present an analysis of X-ray spectra of the high mass X-ray binary 4U 0115+634 as observed with Suzaku and RXTE in 2011 July, during the fading phase of a giant X-ray outburst. We used a continuum model consisting of an absorbed cutoff power-law and an ad-hoc Gaussian emission feature centered around 8.5 keV, which we discuss to be due to cyclotron emission. Our results are consistent with a fundamental cyclotron absorption line centered at ${\sim}10.2$ keV for all observed flux ranges. At the same time we rule out significant influence of the 8.5 kev Gaussian on the CRSF parameters, which are not consistent with the cyclotron line energies and depths of previously reported flux-dependent descriptions. We also show that some continuum models can lead to artificial line-like residuals in the analyzed spectra, which are then misinterpreted as unphysically strong cyclotron lines. Specifically, our results do not support the existence of a previously claimed additional cyclotron feature at ${\sim}15$ keV. Apart from these features, we find for the first time evidence for a He-like Fe XXV emission line at ${\sim}6.7$ keV and weak H-like Fe XXVI emission close to ${\sim}7.0$ keV.

We present a broad-band X-ray study of the effect of superorbital periods on X-ray spectra and pulse profiles in the neutron star X-ray binaries LMC X-4 and SMC X-1. These two sources display periodic or quasi-periodic variations in luminosity on the order of tens of days which are known to be superorbital, and are attributed to warped, precessing accretion disks. Using joint NuSTAR and XMM-Newton observations that span a complete superorbital cycle, we examine the broad-band spectra of these sources and find the shape to be well described by an absorbed power law with a soft blackbody component. Changes in spectral shape and pulse profile shape are periodic with superorbital period, as expected from a precessing disk. We perform X-ray tomography using the changes in pulse profiles to model the geometry and kinematics of the inner accretion disk. Our simple beam and inner disk geometric model indicates that the long term changes in soft pulse shape and phase are consistent with reprocessed emission from a precessing inner disk.

We present the results from a monitoring campaign made with the Neil Gehrels Swift Observatory of the M51 galaxies, which contain several variable ultraluminous X-ray sources (ULXs). The ongoing campaign started in May 2018, and we report here on $\sim1.5$ years of observations. The campaign, which consists of 105 observations, has a typical cadence of 3--6 days, and has the goal of determining the long-term X-ray variability of the ULXs. Two of the most variable sources were ULX7 and ULX8, both of which are known to be powered by neutron stars that are exceeding their isotropic Eddington luminosities by factors of up to 100. This is further evidence that neutron star powered ULXs are the most variable. Our two main results are, first, that ULX7 exhibits a periodic flux modulation with a period of 38 days varying over a magnitude and a half in flux from peak to trough. Since the orbital period of the system is known to be 2 days, the modulation is super-orbital, which is a near-ubiquitous property of ULX pulsars. Secondly we identify a new transient ULX, M51 XT-1, the onset of which occurred during our campaign, reaching a peak luminosity of $\sim10^{40}$ erg s$^{-1}$, before gradually fading over the next $\sim200$ days until it slipped below the detection limit of our observations. Combined with the high-quality Swift/XRT lightcurve of the transient, serendipitous observations made with Chandra and XMM-Newton provide insights into the onset and evolution of a likely super-Eddington event.

M82 X-2 is the first pulsating ultraluminous X-ray source (PULX)to be identified. Since the discovery in 2014, NuSTAR has observed the M82 field 15 times throughout 2015 and 2016. In this paper, we report the results of pulsation searches in all these datasets, and find only one new detection. This new detection allows us to refine the orbital period of the source and measure an average spin down rate between 2014 and 2016 of 5x10^-11 Hz/s, which is in contrast to the strong spin up seen during the 2014 observations and represent the first detection of spin down in a PULX system. Thanks to the improved orbital solution allowed by this new detection, we are also able to detect pulsations in additional segments of the original 2014 dataset. We find a glitch superimposed on the very strong and variable spin-up already reported, the first positive glitch identified in a PULX system. We discuss the new findings in the context of current leading models for PULXs.

We present a spectral study of the ultraluminous Be/X-ray transient pulsar Swift J0243.6+6124 using Neutron Star Interior Composition Explorer (NICER) observations during the system's 2017--2018 giant outburst. The 1.2--10~keV energy spectrum of the source can be approximated with an absorbed cut-off power law model. We detect strong, luminosity-dependent emission lines in the 6--7 keV energy range. A narrow 6.42 keV line, observed in the sub-Eddington regime, is seen to evolve into a broad Fe-line profile in the super-Eddington regime. Other features are found at 6.67 and 6.97 keV in the Fe-line complex. An asymmetric broad line profile, peaking at 6.67 keV, is possibly due to Doppler effects and gravitational redshift. The 1.2--79 keV broadband spectrum from NuSTAR and NICER observations at the outburst peak is well described by an absorbed cut-off power law plus multiple Gaussian lines and a blackbody component. Physical reflection models are also tested to probe the broad iron line feature. Depending on the mass accretion rate, we found emission sites that are evolving from ~5000 km to a range closer to the surface of the neutron star. Our findings are discussed in the framework of the accretion disk and its implication on the magnetic field, the presence of optically thick accretion curtain in the magnetosphere, jet emission, and the massive, ultra-fast outflow expected at super-Eddington accretion rates. We do not detect any signatures of a cyclotron absorption line in the NICER or NuSTAR data.

MAXI J1820+070 (optical counterpart ASASSN-18ey) is a black hole candidate discovered through its recent very bright outburst. The low extinction column and long duration at high flux allow detailed measurements of the accretion process to be made. In this work, we compare the evolution of X-ray spectral and timing properties through the initial hard state of the outburst. We show that the inner accretion disc, as measured by relativistic reflection, remains steady throughout this period of the outburst. Nevertheless, subtle spectral variability is observed, which is well explained by a change in coronal geometry. However, characteristic features of the temporal variability - low-frequency roll-over and QPO frequency - increase drastically in frequency, as the outburst proceeds. This suggests that the variability timescales are governed by coronal conditions rather than solely by the inner disc radius. We also find a strong correlation between X-ray luminosity and coronal temperature. This can be explained by electron pair production with a changing effective radius and a non-thermal electron fraction of ~20%.

SN2010da/NGC 300 ULX-1 was first detected as a supernova impostor in May 2010 and was recently discovered to be a pulsating ultraluminous X-ray source. In this letter, we present VLT/X-shooter spectra of this source obtained in October 2018, covering the wavelength range 350-2300 nm. The $J$- and $H$-bands clearly show the presence of a red supergiant donor star that is best matched by a MARCS stellar atmosphere with $T_{\rm eff} = 3650 - 3900$ K and $\log(L_{\rm bol}/L_{\odot}) = 4.25\pm0.10$, which yields a stellar radius $R = 310 \pm 70 R_{\odot}$. To fit the full spectrum, two additional components are required: a blue excess that can be fitted either by a hot blackbody (T $\gtrsim 20,000$ K) or a power law (spectral index $\alpha \approx 4$) and is likely due to X-ray emission reprocessed in the outer accretion disk or the donor star; and a red excess that is well fitted by a blackbody with a temperature of $\sim 1100$ K, and is likely due to warm dust in the vicinity of SN2010da. The presence of a red supergiant in this system implies an orbital period of at least 0.8-2.1 years, assuming Roche lobe overflow. Given the large donor-to-compact object mass ratio, orbital modulations of the radial velocity of the red supergiant are likely undetectable. However, the radial velocity amplitude of the neutron star is large enough (up to 40-60 km s$^{-1}$) to potentially be measured in the future, unless the system is viewed at a very unfavorable inclination.

Hypervelocity stars are rare objects, mostly main-sequence (MS) B stars, traveling so fast that they will eventually escape from the Milky Way. Recently, it has been shown that the popular Hills mechanism, in which a binary system is disrupted via a close encounter with the supermassive black hole at the Galactic center, may not be their only ejection mechanism. The analyses of Gaia data ruled out a Galactic center origin for some of them, and instead indicated that they are extreme disk runaway stars ejected at velocities exceeding the predicted limits of classical scenarios (dynamical ejection from star clusters or binary supernova ejection). We present the discovery of a new extreme disk runaway star, PG 1610$+$062, which is a slowly pulsating B star bright enough to be studied in detail. A quantitative analysis of spectra taken with ESI at the Keck Observatory revealed that PG 1610$+$062 is a late B-type MS star of 4-5 $M_\odot$ with low projected rotational velocity. Abundances (C, N, O, Ne, Mg, Al, Si, S, Ar, and Fe) were derived differentially with respect to the normal B star HD 137366 and indicate that PG 1610$+$062 is somewhat metal rich. A kinematic analysis, based on our spectrophotometric distance (17.3 kpc) and on proper motions from Gaia's second data release, shows that PG 1610$+$062 was probably ejected from the Carina-Sagittarius spiral arm at a velocity of $550\pm40$ km s$^{-1}$, which is beyond the classical limits. Accordingly, the star is in the top five of the most extreme MS disk runaway stars and is only the second among the five for which the chemical composition is known.

Vela X-1 is an eclipsing high mass X-ray binary (HMXB) consisting of a 283s accreting X-ray pulsar in a close orbit of 8.964 days around the B0.5Ib supergiant HD77581 at a distance of just 2.4 kpc. The system is considered a prototype of wind-accreting HMXB and it has been used as a baseline in different theoretical or modelling studies. We discuss the observational properties of the system and the use of the observational data as laboratory to test recent developments in modelling the accretion process in High-Mass X-ray Binaries (e.g., Sander et al. 2018; El Mellah et al. 2018), which range from detailed descriptions of the wind acceleration to modelling of the structure of the flow of matter close to the neutron star and its variations.

Two recent observations of the nearby galaxy NGC 6946 with NuSTAR, one simultaneous with an XMM-Newton observation, provide an opportunity to examine its population of bright accreting sources from a broadband perspective. We study the three known ultraluminous X-ray sources (ULXs) in the galaxy, and find that ULX-1 and ULX-2 have very steep power-law spectra with $\Gamma=3.6^{+0.4}_{-0.3}$ in both cases. Their properties are consistent with being super-Eddington accreting sources with the majority of their hard emission obscured and down-scattered. ULX-3 (NGC 6946 X-1) is significantly detected by both XMM-Newton and NuSTAR at $L_{\rm X}=(6.5\pm0.1)\times10^{39}$ erg s$^{-1}$, and has a power-law spectrum with $\Gamma=2.51\pm0.05$. We are unable to identify a high-energy break in its spectrum like that found in other ULXs, but the soft spectrum likely hinders our ability to detect one. We also characterise the new source, ULX-4, which is only detected in the joint XMM-Newton and NuSTAR observation, at $L_{\rm X}=(2.27\pm0.07)\times10^{39}$ erg s$^{-1}$, and is absent in a Chandra observation ten days later. It has a very hard cut-off power-law spectrum with $\Gamma=0.7\pm0.1$ and $E_{\rm cut}=11^{+9}_{-4}$ keV. We do not detect pulsations from ULX-4, but its transient nature can be explained either as a neutron star ULX briefly leaving the propeller regime or as a micro-tidal disruption event induced by a stellar-mass compact object.

Studying the physical processes occurring in the region just above the magnetic poles of strongly magnetized, accreting binary neutron stars is essential to our understanding of stellar and binary system evolution. Perhaps more importantly, it provides us with a natural laboratory for studying the physics of high temperature and high density plasmas exposed to extreme radiation, gravitational, and magnetic fields. Observations over the past decade have shed new light on the manner in which plasma falling at velocities near the speed of light onto a neutron star surface is halted. Recent advances in modeling these processes have resulted in direct measurement of the magnetic fields and plasma properties. On the other hand, numerous physical processes have been identified that challenge our current picture of how the accretion process onto neutron stars works. Observation and theory are our essential tools in this regime because the extreme conditions cannot be duplicated on Earth. This white paper gives an overview of the current theory, the outstanding theoretical and observational challenges, and the importance of addressing them in contemporary astrophysics research.

With the recent discoveries of massive and highly luminous quasars at high redshifts ($z\sim7$; e.g. Mortlock et al. 2011), the question of how black holes (BHs) grow in the early Universe has been cast in a new light. In order to grow BHs with $M_{\rm BH} > 10^9$ M$_{\odot}$ by less than a billion years after the Big Bang, mass accretion onto the low-mass seed BHs needs to have been very rapid (Volonteri & Rees, 2005). Indeed, for any stellar remnant seed, the rate required would need to exceed the Eddington limit. This is the point at which the outward force produced by radiation pressure is equal to the gravitational attraction experienced by the in-falling matter. In principle, this implies that there is a maximum luminosity an object of mass $M$ can emit; assuming spherical accretion and that the opacity is dominated by Thompson scattering, this Eddington luminosity is $L_{\rm{E}} = 1.38 \times 10^{38} (M/M_{\odot})$ erg s$^{-1}$. In reality, it is known that this limit can be violated, due to non-spherical geometry or various kinds of instabilities. Nevertheless, the Eddington limit remains an important reference point, and many of the details of how accretion proceeds above this limit remain unclear. Understanding how this so-called super-Eddington accretion occurs is of clear cosmological importance, since it potentially governs the growth of the first supermassive black holes (SMBHs) and the impact this growth would have had on their host galaxies (`feedback') and the epoch of reionization, as well as improving our understanding of accretion physics more generally.

We report on NuSTAR observations of transient pulsations in the neutron star X-ray binary SMC X-1. The transition from non-pulsing to pulsing states was not accompanied by a large change in flux. Instead, both pulsing and non-pulsing states were observed in a single observation during the low-flux super-orbital state. During the high-state, we measure a pulse period of $P = 0.70117(9)\,\mathrm{s}$ at $T_{ref} = 56145\,\mathrm{MJD}$. Spectral analysis during non-pulsing and pulsing states reveals that the observations can be consistently modeled by an absorbed power law with a phenomenological cutoff resembling a Fermi-Dirac distribution, or by a partially obscured cutoff power law. The shapes of the underlying continua show little variability between epochs, while the covering fraction and column density vary between super-orbital states. The strength of pulsations also varies, leading us to infer that the absence and reemergence of pulsations are related to changing obscuration, such as by a warped accretion disk. SMC X-1 is accreting near or above its Eddington limit, reaching an unabsorbed X-ray luminosity of $L_{\rm X}({\rm 2-10~keV}) \approx 5 \times 10^{38}\, {\rm erg}\, {\rm s}^{-1}$. This suggests that SMC X-1 may be a useful local analog to ultraluminous X-ray pulsars (ULXPs), which likewise exhibit strong variability in their pulsed fractions, as well as flux variability on similar timescales. In particular, the gradual pulse turn-on which has been observed in M82 X-2 is similar to the behavior we observe in SMC X-1. Thus we propose that pulse fraction variability of ULXPs may also be due to variable obscuration.

We have performed a comprehensive spectral and timing analysis of the first NuSTAR observation of the high-mass X-ray binary 4U 1538-522. The observation covers the X-ray eclipse of the source, plus the eclipse ingress and egress. We use the new measurement of the mid-eclipse time to update the orbital parameters of the system and find marginally-significant evolution in the orbital period, with $\dot{P}_{\rm orb}/P_{\rm orb} = \left(-0.95 \pm 0.37\right) \times 10^{-6}$ yr$^{-1}$. The cyclotron line energy is found approximately 1.2 keV higher than RXTE measurements from 1997--2003, in line with the increased energy observed by Suzaku in 2012 and strengthening the case for secular evolution of 4U 1538-522's CRSF. We additionally characterize the behavior of the iron fluorescence and emission lines and line-of-sight absorption as the source moves into and out of eclipse.

Ultraluminous X-ray (ULX) pulsars are a new class of object powered by apparent super-critical accretion onto magnetized neutron stars. Three sources in this class identified so far; M82 X-2, NGC 5907 ULX-1 and NGC 7793 P13, have been found to have two properties in common; $\sim1$-s spin periods, and for NGC 5907 ULX-1 and NGC 7793 P13 periodic X-ray flux modulations on timescales of $\sim60-80$ days. M82 X-2 resides in a crowded field that includes the ULX M82 X-1 separated from X-2 by 5", as well as other bright point sources. A 60-day modulation has been observed from the region but the origin has been difficult to identify; both M82 X-1 and X-2 have been suggested as the source. In this paper we present the analysis of a systematic monitoring campaign by Chandra, the only X-ray telescope capable of resolving the crowded field. From a simple Lomb-Scargle periodogram analysis and a more sophisticated Gaussian Process analysis we find that only X-2 exhibits a periodic signal around 60 days supporting previous claims that it is the origin. We also construct a phase-averaged flux profile of the modulations from higher cadence Swift/XRT data and find that the flux variations in the Chandra data are fully consistent with the flux profile. Since the orbit of the neutron star and its companion is known to be 2.5 days, the $\sim60$-day period must be super-orbital in origin. The flux of the modulations varies by a factor of $\sim$100 from minimum to maximum, with no evidence for spectral variations, making the origin difficult to explain.

We report on our search for the optical counterparts of two ultraluminous X-ray pulsars with known orbital periods, M82 X-2 and NGC 5907 X-1, in new and archival HST observations, in an effort to characterize the donor stars in these systems. We detect five near-infrared sources consistent with the position of M82 X-2 that are too bright to be single stars. We also detect seven sources in the WFC3/UVIS F336W image whose photometry matches that of 10-15 M$_\odot$ stars turning off the main sequence. Such stars have densities consistent with the properties of the donor star of M82 X-2 as inferred from X-ray timing analysis, although it is also possible that the donor is a lower mass star below our detection limit or that there is a significant contribution from the accretion disc to the optical emission. We detect three candidate counterparts to NGC 5907 X-1 in the near-infrared. All of these are too bright to be the donor star of the ULX, which based on its orbital period is a red giant. The high background at the location of NGC 5907 X-1 precludes us from detecting this expected donor star. The recently discovered NGC 5907 ULX-2 also falls within the field of view of the near-infrared imaging; we detect four sources in the error circle, with photometry that matches AGB stars. The star suggested to be the counterpart of NGC 5907 ULX-2 by Pintore et al. (2018) falls outside our 2-$\sigma$ error circle.

In this paper we present the science potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies of strongly magnetized objects. We will focus on the physics and astrophysics of strongly magnetized objects, namely magnetars, accreting X-ray pulsars, and rotation powered pulsars. We also discuss the science potential of eXTP for QED studies. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Sciences, the eXTP mission is expected to be launched in the mid 2020s.

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.

The high mass X-ray binary Cep X-4, during its 2014 outburst, showed evidence for an asymmetric cyclotron line in its hard X-ray spectrum. The 2014 spectrum provides one of the clearest cases of an asymmetric line profile among all studied sources with Cyclotron Resonance Scattering Features (CRSF). We present a phase-resolved analysis of NuSTAR and Suzaku data taken at the peak and during the decline phases of this outburst. We find that the pulse-phased resolved spectra are well-fit by a single, symmetric cyclotron feature. The fit parameters vary strongly with pulse phase: most notably the central energy and depth of the cyclotron feature, the slope of the power-law component, and the absorbing column density. We synthesise a phase averaged spectrum using the best-fit parameters for these individual pulse phases, and find that this combined model spectrum has a similar asymmetry in the cyclotron features as discovered in phase-averaged data. We conclude that the pulse phase resolved analysis with simple symmetric line profiles when combined can explain the asymmetry detected in the phase-averaged data.

Two \textit{NuSTAR} observations of the luminous X-ray pulsar LMC X-4 in October and November 2015 captured several bright accretion flares from this source, which has a long history of stable pulse and superorbital behavior. We present a timing analysis of these data in which we detect a rapid pulse turn-on in association with the accretion flares, during which the source reaches super-Eddington luminosities. Pulsations, which are normally seen from this source, are found to only occur for approximately one hour before and during the bright flares. Beyond one hour before and after the flares, we find pulsations to be weak or nonexistent, with fractional RMS amplitudes of less than 0.05. At the onset of the flare, the pulse profiles exhibit a phase shift of 0.25 cycles that could be associated with a change in the emission geometry. This increase in pulse strength occurring well before the flare cannot be explained by the propeller effect, and potentially offers a connection between the magnetic properties of pulsars that accrete close to their Eddington limits and ULX pulsars.

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.

Ultraluminous pulsars are a definite proof that persistent super-Eddington accretion occurs in nature. They support the scenario according to which most Ultraluminous X-ray Sources (ULXs) are super-Eddington accretors of stellar mass rather than sub-Eddington intermediate mass black holes. An important prediction of theories of supercritical accretion is the existence of powerful outflows of moderately ionized gas at mildly relativistic speeds. In practice, the spectral resolution of X-ray gratings such as RGS onboard XMM-Newton is required to resolve their observational signatures in ULXs. Using RGS, outflows have been discovered in the spectra of 3 ULXs (none of which are currently known to be pulsars). Most recently, the fourth ultraluminous pulsar was discovered in NGC 300. Here we report detection of an ultrafast outflow (UFO) in the X-ray spectrum of the object, with a significance of more than 3{\sigma}, during one of the two simultaneous observations of the source by XMM-Newton and NuSTAR in December 2016. The outflow has a projected velocity of 65000 km/s (0.22c) and a high ionisation factor with a log value of 3.9. This is the first direct evidence for a UFO in a neutron star ULX and also the first time that this its evidence in a ULX spectrum is seen in both soft and hard X-ray data simultaneously. We find no evidence of the UFO during the other observation of the object, which could be explained by either clumpy nature of the absorber or a slight change in our viewing angle of the accretion flow.

We present an X-ray spectral and timing analysis of two $NuSTAR$ observations of the transient Be X-ray binary SAX J2103.5+4545 during its April 2016 outburst, which was characterized by the highest flux since $NuSTAR$'s launch. These observations provide detailed hard X-ray spectra of this source during its bright precursor flare and subsequent fainter regular outburst for the first time. In this work, we model the phase-averaged spectra for these observations with a negative and positive power law with an exponential cut-off (NPEX) model and compare the pulse profiles at different flux states. We found that the broad-band pulse profile changes from a three peaked pulse in the first observation to a two peaked pulse in the second observation, and that each of the pulse peaks has some energy dependence. We also perform pulse-phase spectroscopy and fit phase-resolved spectra with NPEX to evaluate how spectral parameters change with pulse phase. We find that while the continuum parameters are mostly constant with pulse phase, a weak absorption feature at ~12 keV that might, with further study, be classified as a cyclotron line, does show strong pulse phase dependence.

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.

Bright and eclipsing, the high-mass X-ray binary Vela X-1 offers a unique opportunity to study accretion onto a neutron star from clumpy winds of O/B stars and to disentangle the complex accretion geometry of these systems. In Chandra-HETGS spectroscopy at orbital phase ~0.25, when our line of sight towards the source does not pass through the large-scale accretion structure such as the accretion wake, we observe changes in overall spectral shape on timescales of a few kiloseconds. This spectral variability is, at least in part, caused by changes in overall absorption and we show that such strongly variable absorption cannot be caused by unperturbed clumpy winds of O/B stars. We detect line features from high and low ionization species of silicon, magnesium and neon whose strengths and presence depend on the overall level of absorption. They imply a co-existence of cool and hot gas phases in the system that we interpret as a highly variable, structured accretion flow close to the compact object such as has been recently seen in simulations of wind accretion in high-mass X-ray binaries.

Recent observations of Her X-1 with NuSTAR and INTEGRAL in 2016 have provided evidence that the 20-year decay of the cyclotron line energy found between 1996 and 2015 has ended and that an inversion with a new increase, possibly similar to the one observed around 1990-1993, has started. We consider this a strong motivation for further observations and for enhanced efforts to significantly improve our theoretical understanding of the accretion process in binary X-ray pulsars. We speculate about the physics behind the long-term decay and its inversion, a possible cyclic behavior, and correlations with other variable observables.

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.

We present a comprehensive spectral analysis of the BeXRB GRO J1008-57 over a luminosity range of three orders of magnitude using NuSTAR, Suzaku and RXTE data. We find significant evolution of the spectral parameters with luminosity. In particular the photon index hardens with increasing luminosity at intermediate luminosities between $10^{36}$ $-$ $10^{37}$ erg s$^{-1}$. This evolution is stable and repeatedly observed over different outbursts. However, at the extreme ends of the observed luminosity range, we find that the correlation breaks down, with a significance level of at least $3.7\sigma$. We conclude that these changes indicate transitions to different accretion regimes, which are characterized by different deceleration processes, such as Coulomb or radiation breaking. We compare our observed luminosity levels of these transitions to theoretical predications and discuss the variation of those theoretical luminosity values with fundamental neutron star parameters. Finally, we present detailed spectroscopy of the unique "triple peaked" outburst in 2014/15 which does not fit in the general parameter evolution with luminosity. The pulse profile on the other hand is consistent with what is expected at this luminosity level, arguing against a change in accretion geometry. In summary, GRO J1008-57 is an ideal target to study different accretion regimes due to the well constrained evolution of its broad-band spectral continuum over several orders of magnitude in luminosity.

We report on two NuSTAR observations of the HMXB A 0535+26 taken toward the end of its normal 2015 outburst at very low $3-50$ keV luminosities of ${\sim}1.4\times10^{36}$ erg/s and ${\sim}5\times10^{35}$ erg/s which are complemented by 9 Swift observations. The data clearly confirm indications seen in earlier data that the source's spectral shape softens as it becomes fainter. The smooth, exponential rollover at high energies present in the first observation evolves to a much more abrupt steepening of the spectrum at $20-30$ keV. The continuum evolution can be well described with emission from a magnetized accretion column, modeled using the compmag model modified by an additional Gaussian emission component for the fainter observation. Between the two observations, the optical depth changes from $0.75\pm0.04$ to $0.56^{+0.01}_{-0.04}$, the electron temperature remains constant, and there is an indication that the column decreases in radius. Since the energy resolved pulse profiles remain virtually unchanged in shape between the two observations, the emission properties of the accretion column, however, reflect the same accretion regime. This conclusion is also confirmed by our result that the energy of the cyclotron resonant scattering feature (CRSF) at ${\sim}45$ keV is independent of the luminosity, implying that the magnetic field in the region in which the observed radiation is produced is the same in both observations. Finally, we also constrain the evolution of the continuum parameters with rotational phase of the neutron star. The width of the CRSF could only be constrained for the brighter observation. Based on Monte-Carlo simulations of CRSF formation in single accretion columns, its pulse phase dependence supports a simplified fan beam emission pattern. The evolution of the CRSF width is very similar to that of the CRSF depth, which is in disagreement with expectations.

We report on the RXTE detection of a sudden increase in the absorption column density, $N_\mathrm{H}$, during the 2011 May outburst of GX 304-1. The $N_\mathrm{H}$ increased up to ${\sim}16\times 10^{22}$ atoms cm$^{-2}$, which is a factor of 3-4 larger than what is usually measured during the outbursts of GX 304-1 as covered by RXTE. Additionally, an increase in the variability of the hardness ratio as calculated from the energy resolved RXTE-PCA light curves is measured during this time range. We interpret these facts as an occultation event of the neutron star by material in the line of sight. Using a simple 3D model of an inclined and precessing Be disk around the Be type companion, we are able to qualitatively explain the $N_\mathrm{H}$ evolution over time. We are able to constrain the Be-disk density to be on the order of $10^{-11}$ g cm$^{-3}$. Our model strengthens the idea of inclined Be disks as origin of double-peaked outbursts as the derived geometry allows accretion twice per orbit under certain conditions.

IGR J18214-1318, a Galactic source discovered by the International Gamma-Ray Astrophysics Laboratory, is a high-mass X-ray binary (HMXB) with a supergiant O-type stellar donor. We report on the XMM-Newton and NuSTAR observations that were undertaken to determine the nature of the compact object in this system. This source exhibits high levels of aperiodic variability, but no periodic pulsations are detected with a 90% confidence upper limit of 2% fractional rms between 0.00003-88 Hz, a frequency range that includes the typical pulse periods of neutron stars (NSs) in HMXBs (0.1-10$^3$ s). Although the lack of pulsations prevents us from definitively identifying the compact object in IGR J18214-1318, the presence of an exponential cutoff with e-folding energy $\lesssim30$ keV in its 0.3-79 keV spectrum strongly suggests that the compact object is an NS. The X-ray spectrum also shows a Fe K$\alpha$ emission line and a soft excess, which can be accounted for by either a partial-covering absorber with $N_{\mathrm{H}}\approx10^{23}$ cm$^{-2}$ which could be due to the inhomogeneous supergiant wind, or a blackbody component with $kT=1.74^{+0.04}_{-0.05}$ keV and $R_{BB}\approx0.3$ km, which may originate from NS hot spots. Although neither explanation for the soft excess can be excluded, the former is more consistent with the properties observed in other supergiant HMXBs. We compare IGR J18214-1318 to other HMXBs that lack pulsations or have long pulsation periods beyond the range covered by our observations.

Important tracers for the dark matter halo of the Galaxy are hypervelocity stars (HVSs), which are faster than the local escape velocity of the Galaxy and their slower counterparts, the high-velocity stars in the Galactic halo. Such HVSs are believed to be ejected from the Galactic Centre (GC) through tidal disruption of a binary by the super-massive black hole (Hills mechanism). The Hyper-MUCHFUSS survey aims at finding high-velocity potentially unbound hot subdwarf stars. We present the spectroscopic and kinematical analyses of a He-sdO as well as three candidates among the sdB stars using optical Keck/ESI and VLT (Xshooter, FORS) spectroscopy. Proper motions are determined by combining positions from early-epoch photographic plates with those derived from modern digital sky surveys. The Galactic rest frame velocities range from 203 km\,s$^{-1}$ to 660 km\,s$^{-1}$, indicating that most likely all four stars are gravitationally bound to the Galaxy. With $T_\text{eff}=47000$ K and a surface gravity of $\log g = 5.7$, SDSS J205030.39$-$061957.8 (J2050) is a spectroscopic twin of the hypervelocity He-sdO US 708. As for the latter, the GC is excluded as a place of origin based on the kinematic analysis. Hence, the Hills mechanism can be excluded for J2050. The ejection velocity is much more moderate ($385\pm79$ km\,s$^{-1}$) than that of US 708 ($998\pm68$ km\,s$^{-1}$). The binary thermonuclear supernova scenario suggested for US 708 would explain the observed properties of J2050 very well without pushing the model parameters to their extreme limits, as required for US 708. Accordingly, the star would be the surviving donor of a type Ia supernova. Three sdB stars also showed extreme kinematics; one could be a HVS ejected from the GC, whereas the other two could be ejected from the Galactic disk through the binary supernova mechanism. Alternatively, they might be extreme halo stars.

Massive stars, at least $\sim$ 10 times more massive than the Sun, have two key properties that make them the main drivers of evolution of star clusters, galaxies, and the Universe as a whole. On the one hand, the outer layers of massive stars are so hot that they produce most of the ionizing ultraviolet radiation of galaxies; in fact, the first massive stars helped to re-ionize the Universe after its Dark Ages. Another important property of massive stars are the strong stellar winds and outflows they produce. This mass loss, and finally the explosion of a massive star as a supernova or a gamma-ray burst, provide a significant input of mechanical and radiative energy into the interstellar space. These two properties together make massive stars one of the most important cosmic engines: they trigger the star formation and enrich the interstellar medium with heavy elements, that ultimately leads to formation of Earth-like rocky planets and the development of complex life. The study of massive star winds is thus a truly multidisciplinary field and has a wide impact on different areas of astronomy. [...] This detailed review summarises the current knowledge on the theory and observations of winds from massive stars, as well as on observations and accretion processes in wind-fed high mass X-ray binaries. The aim is to combine in the near future all available theoretical diagnostics and observational measurements to achieve a unified picture of massive star winds in isolated objects and in binary systems.

Strong winds from massive stars are a topic of interest to a wide range of astrophysical fields. In High-Mass X-ray Binaries the presence of an accreting compact object on the one side allows to infer wind parameters from studies of the varying properties of the emitted X-rays; but on the other side the accretor's gravity and ionizing radiation can strongly influence the wind flow. Based on a collaborative effort of astronomers both from the stellar wind and the X-ray community, this presentation attempts to review our current state of knowledge and indicate avenues for future progress.

Cyclotron resonant scattering features (CRSFs) are formed by scattering of X-ray photons off quantized plasma electrons in the strong magnetic field (of the order 10^12 G) close to the surface of an accreting X-ray pulsar. The line profiles of CRSFs cannot be described by an analytic expression. Numerical methods such as Monte Carlo (MC) simulations of the scattering processes are required in order to predict precise line shapes for a given physical setup, which can be compared to observations to gain information about the underlying physics in these systems. A versatile simulation code is needed for the generation of synthetic cyclotron lines. Sophisticated geometries should be investigatable by making their simulation possible for the first time. The simulation utilizes the mean free path tables described in the first paper of this series for the fast interpolation of propagation lengths. The code is parallelized to make the very time consuming simulations possible on convenient time scales. Furthermore, it can generate responses to mono-energetic photon injections, producing Green's functions, which can be used later to generate spectra for arbitrary continua. We develop a new simulation code to generate synthetic cyclotron lines for complex scenarios, allowing for unprecedented physical interpretation of the observed data. An associated XSPEC model implementation is used to fit synthetic line profiles to NuSTAR data of Cep X-4. The code has been developed with the main goal of overcoming previous geometrical constraints in MC simulations of CRSFs. By applying this code also to more simple, classic geometries used in previous works, we furthermore address issues of code verification and cross-comparison of various models. The XSPEC model and the Green's function tables are available online at http://www.sternwarte.uni-erlangen.de/research/cyclo .

We report on FORS2 optical spectroscopy of the black hole X-ray binary V404 Cygni, performed at the very beginning of its 2015 outburst decay, complemented by quasi-simultaneous $Swift$ X-ray and ultra-violet as well as REM near-infrared observations. Its peculiar spectrum is dominated by a wealth of emission signatures of HI, HeI, and higher ionisation species, in particular FeII. The spectral features are divided between broad red-shifted and narrow stationary varieties, the latter being emitted in the outer regions. Continuum and line variability at short time scale is high and we find Baldwin effect-like anti-correlations between the full-widths at half-maximum and equivalent widths of the broad lines with their local continua. The Balmer decrement H{\alpha}/H{\beta} is also abnormally large at $4.61\pm0.62$. We argue that these properties hint at the broad lines being optically thick and arising within a circumbinary component in which shocks between faster optically thick and slower optically thin regions may occur. We associate it to a nova-like nebula formed by the cooling remnant of strong accretion disc winds that turned off when the mass-accretion rate dropped following the last major flare. The FeII lines likely arise from the overlap region between this nebula and the companion star winds, whereas we favour the shocks within the nebula as responsible for the optical continuum via self-absorbed optically thin bremsstrahlung. The presence of a near-infrared excess also points towards the contribution of a strongly variable compact jet or a dusty component.

We report new spectral modeling of the accreting X-ray pulsar Hercules X- 1. Our radiation-dominated radiative shock model is an implementation of the analytic work of Becker & Wolff on Comptonized accretion flows onto magnetic neutron stars. We obtain a good fit to the spin-phase averaged 4 to 78 keV X-ray spectrum observed by the Nuclear Spectroscopic Telescope Array during a main- on phase of the Her X-1 35-day accretion disk precession period. This model allows us to estimate the accretion rate, the Comptonizing temperature of the radiating plasma, the radius of the magnetic polar cap, and the average scattering opacity parameters in the accretion column. This is in contrast to previous phenomenological models that characterized the shape of the X-ray spectrum but could not determine the physical parameters of the accretion flow. We describe the spectral fitting details and discuss the interpretation of the accretion flow physical parameters.

M82 X-1 is one of the brightest ultraluminous X-ray sources (ULXs) known, which, assuming Eddington-limited accretion and other considerations, makes it one of the best intermediate-mass black hole (IMBH) candidates. However, the ULX may still be explained by super-Eddington accretion onto a stellar-remnant black hole. We present simultaneous NuSTAR, Chandra and Swift/XRT observations during the peak of a flaring episode with the aim of modeling the emission of M82 X-1 and yielding insights into its nature. We find that thin-accretion disk models all require accretion rates at or above the Eddington limit in order to reproduce the spectral shape, given a range of black hole masses and spins. Since at these high Eddington ratios the thin-disk model breaks down due to radial advection in the disk, we discard the results of the thin-disk models as unphysical. We find that the temperature profile as a function of disk radius ($T(r)\propto r^{-p}$) is significantly flatter ($p=0.55^{+ 0.07}_{- 0.04}$) than expected for a standard thin disk ($p=0.75$). A flatter profile is instead characteristic of a slim disk which is highly suggestive of super-Eddington accretion. Furthermore, radiation hydrodynamical simulations of super-Eddington accretion have shown that the predicted spectra of these systems are very similar to what we observe for M82 X-1. We therefore conclude that M82 X-1 is a super-Eddington accretor. Our mass estimates inferred from the inner disk radius imply a stellar-remnant black hole ($M_{\rm BH}=26^{+9}_{-6} M_{\odot}$) when assuming zero spin, or an IMBH ($M_{\rm BH}=125^{+45}_{-30} M_{\odot}$) when assuming maximal spin.

In a previous work, we introduced a tool for analyzing multiple datasets simultaneously, which has been implemented into ISIS. This tool was used to fit many spectra of X-ray binaries. However, the large number of degrees of freedom and individual datasets raise an issue about a good measure for a simultaneous fit quality. We present three ways to check the goodness of these fits: we investigate the goodness of each fit in all datasets, we define a combined goodness exploiting the logical structure of a simultaneous fit, and we stack the fit residuals of all datasets to detect weak features. These tools are applied to all RXTE-spectra from GRO 1008-57, revealing calibration features that are not detected significantly in any single spectrum. Stacking the residuals from the best-fit model for the Vela X-1 and XTE J1859+083 data evidences fluorescent emission lines that would have gone undetected otherwise.

New spectroscopic observations of the halo hyper-velocity star candidate SDSS J121150.27+143716.2 ($V=17.92$ mag) revealed a cool companion to the hot subdwarf primary. The components have a very similar radial velocity and their absolute luminosities are consistent with the same distance, confirming the physical nature of the binary, which is the first double-lined hyper-velocity candidate. Our spectral decomposition of the Keck/ESI spectrum provided an sdB+K3V pair, analogous to many long-period subdwarf binaries observed in the Galactic disk. We found the subdwarf atmospheric parameters: $T_{\rm eff}=30\,600\pm500$ K, $\log{g}=5.57\pm0.06$ cm s$^{-2}$ and He abundance $\log(n{\rm He}/n{\rm H})=-3.0\pm0.2$. Oxygen is the most abundant metal in the hot subdwarf atmosphere, and Mg and Na lines are the most prominent spectral features of the cool companion, consistent with a metallicity of $[{\rm Fe}/{\rm H}]=-1.3$. The non-detection of radial velocity variations suggest the orbital period to be a few hundred days, in agreement with similar binaries observed in the disk. Using the SDSS-III flux calibrated spectrum we measured the distance to the system $d=5.5\pm0.5$ kpc, which is consistent with ultraviolet, optical, and infrared photometric constraints derived from binary spectral energy distributions. Our kinematic study shows that the Galactic rest-frame velocity of the system is so high that an unbound orbit cannot be ruled out. On the other hand, a bound orbit requires a massive dark matter halo. We conclude that the binary either formed in the halo or it was accreted from the tidal debris of a dwarf galaxy by the Milky Way.

We have performed a full time- and luminosity-resolved spectral analysis of the high-mass X-ray binary 4U 1538-522 using the available RXTE, INTEGRAL, and Suzaku data, examining both phase-averaged and pulse-phase-constrained datasets and focusing on the behavior of the cyclotron resonance scattering feature (CRSF). No statistically significant trend between the energy of the CRSF and luminosity is observed in the combined dataset. However, the CRSF energy appears to have increased by ~1.5 keV in the ~8.5 years between the RXTE and Suzaku measurements, with Monte Carlo simulations finding the Suzaku measurement 4.6$\sigma$ above the RXTE points. Interestingly, the increased Suzaku CRSF energy is much more significant and robust in the pulse-phase-constrained spectra from the peak of the main pulse, suggesting a change that is limited to a single magnetic pole. The 7 years of RXTE measurements do not show any strongly-significant evolution with time on their own. We discuss the significance of the CRSF's behavior with respect to luminosity and time in the context of historical observations of this source as well as recent observational and theoretical work concerning the neutron star accretion column, and suggest some mechanisms by which the observed change over time could occur.

We report on the timing and spectral analysis of two Suzaku observations with different flux levels of the high-mass X-ray binary KS 1947+300 during its 2013 outburst. In agreement with simultaneous NuSTAR observations, the continuum is well described by an absorbed power law with a cut-off and an additional black body component. In addition we find fluorescent emission from neutral, He-like, and even H-like iron. We determine a pulse period of ~18.8 s with the source showing a spin-up between the two observations. Both Suzaku observations show a very similar behavior of the pulse profile, which is strongly energy dependent, with an evolution from a profile with one peak at low energies to a profile with two peaks of different widths towards higher energies seen in both, the Suzaku and NuSTAR data. Such an evolution to a more complex profile at higher energies is rarely seen in X-ray pulsars, most cases show the opposite behavior. Pulse phase-resolved spectral analysis shows a variation in the absorbing column density, NH , over pulse phase. Spectra taken during the pulse profile minima are intrinsically softer compared to the pulse phase-averaged spectrum.

We present results from the spectral fitting of the candidate black hole X-ray binary Swift J1753.5-0127 in an accretion state previously unseen in this source. We fit the 0.7-78 keV spectrum with a number of models, however the preferred model is one of a multi-temperature disk with an inner disk temperature $\mathrm{k}T_\mathrm{in}=0.252\pm0.003$ keV scattered into a steep power-law with photon index $\Gamma=6.39^{+0.08}_{-0.02}$ and an additional hard power law tail ($\Gamma=1.79\pm0.02$). We report on the emergence of a strong disk-dominated component in the X-ray spectrum and we conclude that the source has entered the soft state for the first time in its ~10 year prolonged outburst. Using reasonable estimates for the distance to the source ($3$ kpc) and black hole mass ($5M_{\odot}$), we find the unabsorbed luminosity (0.1-100 keV) to be $\approx0.60$% of the Eddington luminosity, making this one of the lowest luminosity soft states recorded in X-ray binaries. We also find that the accretion disk extended towards the compact object during its transition from hard to soft, with the inner radius estimated to be $R_{\mathrm{in}}=28.0^{+0.7}_{-0.4} R_g$ or ~$12R_g$, dependent on the boundary condition chosen, assuming the above distance and mass, a spectral hardening factor $f=1.7$ and a binary inclination $i=55^{\circ}$.

We present a timing and spectral analysis of the X-ray pulsar XTE J1946+274 observed with Suzaku during an outburst decline in 2010 October and compare with previous results. XTE J1946+274 is a transient X-ray binary consisting of a Be-type star and a neutron star with a 15.75 s pulse period in a 172 d orbit with 2-3 outbursts per orbit during phases of activity. We improve the orbital solution using data from multiple instruments. The X-ray spectrum can be described by an absorbed Fermi-Dirac cutoff power law model along with a narrow Fe K line at 6.4 keV and a weak Cyclotron Resonance Scattering Feature (CRSF) at ~35 keV. The Suzaku data are consistent with the previously observed continuum flux versus iron line flux correlation expected from fluorescence emission along the line of sight. However, the observed iron line flux is slightly higher, indicating the possibility of a higher iron abundance or the presence of non-uniform material. We argue that the source most likely has only been observed in the subcritical (non-radiation dominated) state since its pulse profile is stable over all observed luminosities and the energy of the CRSF is approximately the same at the highest (~$5 \times 10^{37} $erg s$^{-1}$) and lowest (~$5 \times 10^{36} $erg s$^{-1}$) observed 3-60 keV luminosities.

We present results from a 244\,ks \textit{NuSTAR} observation of 3C\,273 obtained during a cross-calibration campaign with the \textit{Chandra}, \textit{INTEGRAL}, \textit{Suzaku}, \textit{Swift}, and \textit{XMM-Newton} observatories. We show that the spectrum, when fit with a power-law model using data from all observatories except \textit{INTEGRAL} over the 1--78\,keV band, leaves significant residuals in the \textit{NuSTAR} data between 30--78\,keV. The \nustar\ 3--78\,keV spectrum is well-described by an exponentially cutoff power-law ($\Gamma = 1.646 \pm 0.006$, E$_\mathrm{cutoff} = 202_{-34}^{+51}$\,keV) with a weak reflection component from cold, dense material. There is also evidence for a weak ($EW = 23 \pm 11$ eV) neutral iron line. We interpret these features as arising from coronal emission plus reflection off an accretion disk or distant material. Beyond 80\,keV \textit{INTEGRAL} data show clear excess flux relative to an extrapolation of the AGN model fit to \nustar. This high-energy power-law is consistent with the presence of a beamed jet, which begins to dominate over emission from the inner accretion flow at 30-40 keV. Modeling the jet locally (in the \textit{NuSTAR} + \textit{INTEGRAL} band) as a power-law, we find the coronal component is fit by $\Gamma_\mathrm{AGN} = 1.638 \pm 0.045$, $E_\mathrm{cutoff} = 47 \pm 15$\,keV, and jet photon index by $\Gamma_\mathrm{jet} = 1.05 \pm 0.4$. We also consider \textit{Fermi}/LAT observations of 3C\,273 and here the broad-band spectrum of the jet can be described by a log-parabolic model, peaking at $\sim 2$\,MeV. Finally, we investigate the spectral variability in the \textit{NuSTAR} band and find an inverse correlation between flux and $\Gamma$.

We report on NuSTAR and Swift observations of a soft state of the neutron star low-mass X-ray binary GS 1826-24, commonly known as the "clocked" burster. The transition to the soft state was recorded in 2014 June through an increase of the 2-20 keV source intensity measured by MAXI, simultaneous with a decrease of the 15-50 keV intensity measured by Swift/BAT. The episode lasted approximately two months, after which the source returned to its usual hard state. We analyze the broad-band spectrum measured by Swift/XRT and NuSTAR, and estimate the accretion rate during the soft episode to be about 13% of Eddington, within the range of previous observations. However, the best fit spectral model, adopting the double Comptonization used previously, exhibits significantly softer components. We detect seven type-I X-ray bursts, all significantly weaker (and with shorter rise and decay times) than observed previously. The burst profiles and recurrence times vary significantly, ruling out the regular bursts that are typical for this source. One burst exhibited photospheric radius expansion, and we estimate the source distance at about (5.7 / xi_b^1/2) kpc, where xi_b parameterizes the possible anisotropy of the burst emission. Interpreting the soft state as a transition from an optically thin inner flow to an optically thick flow passing through a boundary layer, as is commonly observed in similar systems, is contradicted by the lower optical depth measured for the double Comptonization model we find for this soft state. The effect of a change in disk geometry on the burst behavior remains unclear.

Between July 2012 and February 2013, NuSTAR and XMM-Newton performed four long-look joint obser- vations of the type 1.8 Seyfert, NGC 1365. We have analyzed the variable absorption seen in these observations in order to characterize the geometry of the absorbing material. Two of the observations caught NGC 1365 in an unusually low absorption state, revealing complexity in the multi-layer absorber which had previously been hidden. We find the need for three distinct zones of neutral absorption in addition to the two zones of ionized absorption and the Compton-thick torus previously seen in this source. The most prominent absorber is likely associated with broad line region clouds with column densities of around $\sim\,$10$^{23}$ cm$^{-2}$ and a highly clumpy nature as evidenced by an occultation event in February 2013. We also find evidence of a patchy absorber with a variable column around $\sim\,10^{22}$ cm$^{-2}$ and a line of sight covering fraction of 0.3$-$0.9 which responds directly to the intrinsic source flux, possibly due to a wind geometry. A full-covering, constant absorber with a low column density of $\sim\,1 \times$ 10$^{22}$ cm$^{-2}$ is also present, though the location of this low density haze is unknown.

Hypervelocity stars (HVS) travel with velocities so high, that they exceed the escape velocity of the Galaxy. Several acceleration mechanisms have been discussed. Only one HVS (US 708, HVS 2) is a compact helium star. Here we present a spectroscopic and kinematic analysis of US\,708. Travelling with a velocity of $\sim1200\,{\rm km\,s^{-1}}$, it is the fastest unbound star in our Galaxy. In reconstructing its trajectory, the Galactic center becomes very unlikely as an origin, which is hardly consistent with the most favored ejection mechanism for the other HVS. Furthermore, we discovered US\,708 to be a fast rotator. According to our binary evolution model it was spun-up by tidal interaction in a close binary and is likely to be the ejected donor remnant of a thermonuclear supernova.

Measurements of the high-energy cut-off in the coronal continuum of active galactic nuclei have long been elusive for all but a small number of the brightest examples. We present a direct measurement of the cut-off energy in the nuclear continuum of the nearby Seyfert 1.9 galaxy MCG -05-23-016 with unprecedented precision. The high sensitivity of NuSTAR up to 79 keV allows us to clearly disentangle the spectral curvature of the primary continuum from that of its reflection component. Using a simple phenomenological model for the hard X-ray spectrum, we constrain the cut-off energy to $116_{-5}^{+6}$ keV with 90% confidence. Testing for more complex models and nuisance parameters that could potentially influence the measurement, we find that the cut-off is detected robustly. We further use simple Comptonized plasma models to provide independent constraints for both the kinetic temperature of the electrons in the corona and its optical depth. At the 90% confidence level, we find $kT_e=29\pm2$ keV and $\tau_e=1.23\pm0.08$ assuming a slab (disk-like) geometry, and $kT_e=25\pm2$ keV and $\tau_e=3.5\pm0.2$ assuming a spherical geometry. Both geometries are found to fit the data equally well and their two principal physical parameters are correlated in both cases. With the optical depth in the optically thick regime, the data are pushing the currently available theoretical models of the Comptonized plasma to the limits of their validity. Since the spectral features and variability arising from the inner accretion disk have been observed previously in MCG -05-23-016, the inferred high optical depth implies that a spherical or disk-like corona cannot be homogeneous.

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.

Ultraluminous X-ray sources (ULX) are off-nuclear point sources in nearby galaxies whose X-ray luminosity exceeds the theoretical maximum for spherical infall (the Eddington limit) onto stellar-mass black holes. Their luminosity ranges from $10^{40}$ erg s$^{-1} < L_X$(0.5 - 10 keV) $<10^{40}$ erg s$^{-1}$. Since higher masses imply less extreme ratios of the luminosity to the isotropic Eddington limit theoretical models have focused on black hole rather than neutron star systems. The most challenging sources to explain are those at the luminous end ($L_X$ > $10^{40}$ erg s$^{-1}$), which require black hole masses MBH >50 solar masses and/or significant departures from the standard thin disk accretion that powers bright Galactic X-ray binaries. Here we report broadband X-ray observations of the nuclear region of the galaxy M82, which contains two bright ULXs. The observations reveal pulsations of average period 1.37 s with a 2.5-day sinusoidal modulation. The pulsations result from the rotation of a magnetized neutron star, and the modulation arises from its binary orbit. The pulsed flux alone corresponds to $L_X$(3 - 30 keV) = $4.9 \times 10^{39}$ erg s$^{-1}$. The pulsating source is spatially coincident with a variable ULX which can reach $L_X$ (0.3 - 10 keV) = $1.8 \times 10^{40}$ erg s$^{-1}$. This association implies a luminosity ~100 times the Eddington limit for a 1.4 solar mass object, or more than ten times brighter than any known accreting pulsar. This finding implies that neutron stars may not be rare in the ULX population, and it challenges physical models for the accretion of matter onto magnetized compact objects.

We present spectral and timing analysis of NuSTAR observations of RX J0520.5$-$6932 in the 3-79 keV band collected during its outburst in January 2014. The target was observed on two epochs and we report the detection of a cyclotron resonant scattering feature with central energies of $E_\mathrm{CRSF} = 31.3_{-0.7}^{+0.8}$ keV and $31.5_{-0.6}^{+0.7}$ keV during the two observations, respectively, corresponding to a magnetic field of $B \approx 2 \times10^{12}$ G. The 3-79 keV luminosity of the system during the two epochs assuming a nominal distance of 50 kpc was $3.667\pm0.007\times 10^{38}\,\mathrm{erg\,s^{-1}}$ and $3.983\pm0.007\times10^{38}\,\mathrm{erg\,s^{-1}}$. Both values are much higher than the critical luminosity of $\approx1.5\times10^{37}\,\mathrm{erg\,s^{-1}}$ above which a radiation dominated shock front may be expected. This adds a new object to the sparse set of three systems that have a cyclotron line observed at luminosities in excess of $10^{38}\,\mathrm{erg\,s^{-1}}$. A broad ($\sigma\approx0.45$ keV) Fe emission line is observed in the spectrum at a central energy of $6.58_{-0.05}^{+0.05}$ keV in both epochs. The pulse profile of the pulsar was observed to be highly asymmetric with a sharply rising and slowly falling profile of the primary peak. We also observed minor variations in the cyclotron line energy and width as a function of the rotation phase.% As in observations of other cyclotron absorption line sources, there is a small ($\Delta\phi\lesssim0.1$) phase difference between the peak of the cyclotron energy variation and the peak of the flux variation.

Timing of high-count rate sources with the NuSTAR Small Explorer Mission requires specialized analysis techniques. NuSTAR was primarily designed for spectroscopic observations of sources with relatively low count-rates rather than for timing analysis of bright objects. The instrumental dead time per event is relatively long (~2.5 msec), and varies by a few percent event-to-event. The most obvious effect is a distortion of the white noise level in the power density spectrum (PDS) that cannot be modeled easily with the standard techniques due to the variable nature of the dead time. In this paper, we show that it is possible to exploit the presence of two completely independent focal planes and use the cross power density spectrum to obtain a good proxy of the white noise-subtracted PDS. Thereafter, one can use a Monte Carlo approach to estimate the remaining effects of dead time, namely a frequency-dependent modulation of the variance and a frequency-independent drop of the sensitivity to variability. In this way, most of the standard timing analysis can be performed, albeit with a sacrifice in signal to noise relative to what would be achieved using more standard techniques. We apply this technique to NuSTAR observations of the black hole binaries GX 339-4, Cyg X-1 and GRS 1915+105.

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.

We present analysis of 100 ks contiguous XMM-Newton data of the prototypical wind accretor Vela X-1. The observation covered eclipse egress between orbital phases 0.134 and 0.265, during which a giant flare took place, enabling us to study the spectral properties both outside and during the flare. This giant flare with a peak luminosity of $3.92^{+0.42}_{-0.09} \times 10^{37}$ erg s$^{-1}$ allows estimates of the physical parameters of the accreted structure with a mass of $\sim$ $10^{21}$ g. We have been able to model several contributions to the observed spectrum with a phenomenological model formed by three absorbed power laws plus three emission lines. After analysing the variations with orbital phase of the column density of each component, as well as those in the Fe and Ni fluorescence lines, we provide a physical interpretation for each spectral component. Meanwhile, the first two components are two aspects of the principal accretion component from the surface of the neutron star, and the third component seems to be the \textit{X-ray light echo} formed in the stellar wind of the companion.

We present the results of NuSTAR and XMM-Newton observations of the two ultraluminous X-ray sources (ULX) NGC 1313 X-1 and X-2. The combined spectral bandpass of the two satellites enables us to produce the first spectrum of X-1 between 0.3 and 30 keV, while X-2 is not significantly detected by NuSTAR above 10 keV. The NuSTAR data demonstrate that X-1 has a clear cutoff above 10 keV, whose presence was only marginally detectable with previous X-ray observations. This cutoff rules out the interpretation of X-1 as a black hole in a standard low/hard state, and it is deeper than predicted for the downturn of a broadened iron line in a reflection-dominated regime. The cutoff differs from the prediction of a single-temperature Comptonization model. Further, a cold disk-like black body component at ~0.3 keV is required by the data, confirming previous measurements by XMM-Newton only. We observe a spectral transition in X-2, from a state with high luminosity and strong variability to a lower-luminosity state with no detectable variability, and we link this behavior to a transition from a super-Eddington to a sub-Eddington regime.

The NuSTAR hard X-ray telescope observed the transient Be/X-ray binary GS 0834-430 during its 2012 outburst - the first active state of this system observed in the past 19 years. We performed timing and spectral analysis, and measured the X-ray spectrum between 3-79keV with high statistical significance. We find the phase-averaged spectrum to be consistent with that observed in many other magnetized accreting pulsars. We fail to detect cyclotron resonance scattering features that would allow us to constrain the pulsar's magnetic field in either phase-averaged or phase-resolved spectra. Timing analysis shows a clearly detected pulse period of ~12.29s in all energy bands. The pulse profiles show a strong, energy-dependent hard phase lag of up to 0.3 cycles in phase, or about 4s. Such dramatic energy-dependent lags in the pulse profile have never before been reported in high-mass X-ray binary (HMXB) pulsars. Previously reported lags have been significantly smaller in phase and restricted to low-energies (E<10keV). We investigate the possible mechanisms that might produce this energy-dependent pulse phase shift. We find the most likely explanation for this effect to be a complex beam geometry.

A study of archival RXTE, Swift, and Suzaku pointed observations of the transient high mass X-ray binary GRO J1008-57 is presented. A new orbital ephemeris based on pulse arrival timing shows times of maximum luminosities during outbursts of GRO J1008-57 to be close to periastron at orbital phase -0.03. This makes the source one of a few for which outburst dates can be predicted with very high precision. Spectra of the source in 2005, 2007, and 2011 can be well described by a simple power law with high energy cutoff and an additional black body at lower energies. The photon index of the power law and the black body flux only depend on the 15-50 keV source flux. No apparent hysteresis effects are seen. These correlations allow to predict the evolution of the pulsar's X-ray spectral shape over all outbursts as a function of just one parameter, the source's flux. If modified by an additional soft component, this prediction even holds during GRO J1008-57's 2012 type II outburst.

The knowledge of the spectral state of a black hole is essential for the interpretation of data from black holes in terms of their emission models. Based on pointed observations of Cyg X-1 with the Rossi X-ray timing Explorer (RXTE) that are used to classify simultaneous RXTE-ASM observations, we develop a scheme based on RXTE -ASM colors and count rates that can be used to classify all observations of this canonical black hole that were performed between 1996 and 2011. We show that a simple count rate criterion, as used previously, leads to a significantly higher fraction of misclassified observations. This scheme enables us to classify single INTEGRAL-IBIS science windows and to obtain summed spectra for the soft, intermediate and hard state with low contamination by other states.

The Be/X-ray binary A 0535+26 showed a normal (type I) outburst in August 2009. It is the fourth in a series of normal outbursts associated with the periastron, but is unusual by presenting a double-peaked light curve. The two peaks reached a flux of ~450 mCrab in the 15-50 keV range. We present results of the timing and spectral analysis of INTEGRAL, RXTE, and Suzaku observations of the outburst. The energy dependent pulse profiles and their evolution during the outburst are studied. No significant differences with respect to other normal outbursts are observed. The centroid energy of the fundamental cyclotron line shows no significant variation during the outburst. A spectral hardening with increasing luminosity is observed. We conclude that the source is accreting in the sub-critical regime. We discuss possible explanations for the double-peaked outburst.

We report on an outburst of the high mass X-ray binary 4U 0115+63 with a pulse period of 3.6s in 2008 March/April as observed with RXTE and INTEGRAL. During the outburst the neutron star's luminosity varied by a factor of 10 in the 3--50\,keV band. In agreement with earlier work we find evidence for five cyclotron resonance scattering features at ~10.7, 21.8, 35.5, 46.7, and 59.7keV. Previous work had found an anti-correlation between the fundamental cyclotron line energy and the X-ray flux. We show that this apparent anti-correlation is probably due to the unphysical interplay of parameters of the cyclotron line with the continuum models used previously, e.g., the negative and positive exponent power law (NPEX). For this model, we show that cyclotron line modeling erroneously leads to describing part of the exponential cutoff and the continuum variability, and not the cyclotron lines. When the X-ray continuum is modeled with a simple exponentially cutoff power law modified by a Gaussian emission feature around 10keV, the correlation between the line energy and the flux vanishes and the line parameters remain virtually constant over the outburst. We therefore conclude that the previously reported anti-correlation is an artifact of the assumptions adopted in the modeling of the continuum.

We report on a series of outbursts of the high mass X-ray binary XTE J1946+274 in 2010/2011 as observed with INTEGRAL, RXTE, and Swift. We discuss possible mechanisms resulting in the extraordinary outburst behavior of this source. The X-ray spectra can be described by standard phenomenological models, enhanced by an absorption feature of unknown origin at about 10 keV and a narrow iron K alpha fluorescence line at 6.4 keV, which are variable in flux and pulse phase. We find possible evidence for the presence of a cyclotron resonance scattering feature at about 25 keV at the 93% level. The presence of a strong cyclotron line at 35 keV seen in data from the source's 1998 outburst and confirmed by a reanalysis of these data can be excluded. This result indicates that the cyclotron line feature in XTE J1946+274 is variable between individual outbursts.

Seventeen accreting neutron star pulsars, mostly high mass X-ray binaries with half of them Be-type transients, are known to exhibit Cyclotron Resonance Scattering Features (CRSFs) in their X-ray spectra, with characteristic line energies from 10 to 60 keV. To date about two thirds of them, plus a few similar systems without known CRSFs, have been observed with Suzaku. We present an overview of results from these observations, including the discovery of a CRSF in the transient 1A 1118-61 and pulse phase resolved spectroscopy of GX 301-2. These observations allow for the determination of cyclotron line parameters to an unprecedented degree of accuracy within a moderate amount of observing time. This is important since these parameters vary - e.g., with orbital phase, pulse phase, or luminosity - depending on the geometry of the magnetic field of the pulsar and the properties of the accretion column at the magnetic poles. We briefly introduce a spectral model for CRSFs that is currently being developed and that for the first time is based on these physical properties. In addition to cyclotron line measurements, selected highlights from the Suzaku analyses include dip and flare studies, e.g., of 4U 1907+09 and Vela X-1, which show clumpy wind effects (like partial absorption and/or a decrease in the mass accretion rate supplied by the wind) and may also display magnetospheric gating effects.

We present a timing and spectral analysis of the variable X-ray source 3A 1954+319. Our analysis is mainly based on an outburst serendipitously observed during INTEGRAL Key Program observations of the Cygnus region in 2008 fall and on the Swift/BAT longterm light curve. Previous observations, though sparse, have identified the source to be one of only nine known symbiotic X-ray binaries, i.e., systems composed of an accreting neutron star orbiting in a highly inhomogeneous medium around an M-giant companion. The spectrum of 3A 1954+319 above 20 keV can be best described by a broken power law model. The extremely long pulse period of ~5.3 hours is clearly visible in the INTEGRAL/ISGRI light curve and confirmed through an epoch folding period search. Furthermore, the light curve allows us to determine a very strong spin up of -2x10^-4 h/h during the outburst. This spin up is confirmed by the pulse period evolution calculated from Swift/BAT data. The Swift/BAT data also show a long spin-down trend prior to the 2008 outburst, which is confirmed in archival INTEGRAL/ISGRI data. We discuss possible accretion models and geometries allowing for the transfer of such large amounts of angular momentum and investigate the harder spectrum of this outburst compared to previously published results.

We investigate the structure of the wind in the neutron star X-ray binary system Vela X-1 by analyzing its flaring behavior. Vela X-1 shows constant flaring, with some flares reaching fluxes of more than 3.0 Crab between 20-60 keV for several 100 seconds, while the average flux is around 250 mCrab. We analyzed all archival INTEGRAL data, calculating the brightness distribution in the 20-60 keV band, which, as we show, closely follows a log-normal distribution. Orbital resolved analysis shows that the structure is strongly variable, explainable by shocks and a fluctuating accretion wake. Analysis of RXTE ASM data suggests a strong orbital change of N_H. Accreted clump masses derived from the INTEGRAL data are on the order of 5 x 10^19 -10^21 g. We show that the lightcurve can be described with a model of multiplicative random numbers. In the course of the simulation we calculate the power spectral density of the system in the 20-100 keV energy band and show that it follows a red-noise power law. We suggest that a mixture of a clumpy wind, shocks, and turbulence can explain the measured mass distribution. As the recently discovered class of supergiant fast X-ray transients (SFXT) seems to show the same parameters for the wind, the link between persistent HMXB like Vela X-1 and SFXT is further strengthened.

We present results from two observations of the wind-accreting X-ray pulsar 4U 1907+09 using the Suzaku observatory. The broadband time-averaged spectrum allows us to examine the continuum emission of the source and the cyclotron resonance scattering feature at ~19 keV. Additionally, using the narrow CCD response of Suzaku near 6 keV allows us to study in detail the Fe K bandpass and to quantify the Fe K beta line for this source for the first time. The source is absorbed by fully-covering material along the line of sight with a column density of NH ~2e22 /cm^2, consistent with a wind accreting geometry, and a high Fe abundance (~3-4 x solar). Time and phase-resolved analyses allow us to study variations in the source spectrum. In particular, dips found in the 2006 observation which are consistent with earlier observations occur in the hard X-ray bandpass, implying a variation of the whole continuum rather than occultation by intervening material, while a dip near the end of the 2007 observation occurs mainly in the lower energies implying an increase in NH along the line of sight, perhaps indicating clumpiness in the stellar wind.

We present a detailed analysis of all archival INTEGRAL data of the accreting X-ray pulsar Vela X-1. We extracted lightcurves in several energy bands from 20 keV up to 60 keV. The lightcurves show that the source was found in very active as well as quiet states. During the active states several giant flares were detected. For these states spectra between 5 keV and 120 keV were obtained. The spectra of the active states were found to be significantly softer than those from the quiet states. We performed a statistical analysis of the flaring behavior. The resulting log-normal distribution of the intensity of Vela X-1 shows that the source spends most of the time at an average flux level of 300 mCrab but also that the distribution extends well up to more than 2.0 Crab.

The evolution of the particle background at an altitude of ~540 km during the time interval between 1996 and 2007 is studied using the particle monitor of the High Energy X-ray Timing Experiment on board NASA's Rossi X-ray Timing Explorer. A special emphasis of this study is the location and strength of the South Atlantic Anomaly (SAA). The size and strength of the SAA are anti-correlated with the the 10.7 cm radio flux of the Sun, which leads the SAA strength by ~1 year reflecting variations in solar heating of the upper atmosphere. The location of the SAA is also found to drift westwards with an average drift rate of about 0.3 deg/yr following the drift of the geomagnetic field configuration. Superimposed to this drift rate are irregularities, where the SAA suddenly moves eastwards and where furthermore the speed of the drift changes. The most prominent of these irregularities is found in the second quarter of 2003 and another event took place in 1999. We suggest that these events are previously unrecognized manifestations of the geomagnetic jerks of the Earth's magnetic field.