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.

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.

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.

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.

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

The Be X-ray binary pulsar EXO 2030+375, first detected in 1985, has shown a significant detected X-ray outburst at nearly every periastron passage of its 46-day orbit for the past ~25 years, with one low state accompanied by a torque reversal in the 1990s. In early 2015 the outbursts progressively became fainter and less regular while the monotonic spin-up flattened. At the same time a decrease in the H$\alpha$ line equivalent width was reported, indicating a change in the disk surrounding the mass donor. In order to explore the source behaviour in the poorly explored low-flux state with a possible transition to a state of centrifugal inhibition of accretion we have undertaken an observing campaign with Swift/XRT, NuSTAR and the Nordic Optical Telescope (NOT). This conference contribution reports the preliminary results obtained from our campaign.

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.

In July 2015, the high-mass X-ray binary V 0332+53 underwent a giant outburst, a decade after the previous one. This system hosts a strongly magnetized neutron star. During the 2004--2005 outburst, an anti-correlation between the centroid energy of its fundamental cyclotron resonance scattering features (CRSFs) and the X-ray luminosity was observed. The long ($\approx 100$ d) and bright ($L_{\rm x} \approx 10^{38}$ erg s$^{-1}$) 2015 outburst offered the opportunity to study during another outburst the unique properties of the fundamental CRSF and its dependence on the X-ray luminosity. The source was observed by the INTEGRAL satellite for $\sim 330$ ks. We exploit the spectral resolution at high energies of the SPectrometer on INTEGRAL (SPI) and the Joint European X-ray Monitors to characterize its spectral properties, focusing in particular on the CRSF-luminosity dependence. We complement the data of the 2015 outburst with those collected by SPI in 2004--2005 and left unpublished so far. We find a highly significant anti-correlation of the centroid energy of the fundamental CRSF and the $3-100$ keV luminosity of $E_1 \propto -0.095(8)L_{37}$ keV. This trend is observed for both outbursts. We confirm the correlation between the width of the fundamental CRSF and the X-ray luminosity previously found in the JEM-X and IBIS dataset of the 2004--2005 outburst. By exploiting the RXTE/ASM and Swift/BAT monitoring data we also report on the detection of a $\sim 34$ d modulation superimposed on the mean profiles and roughly consistent with the orbital period of the pulsar. We discuss possible interpretations of such variability.

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.

Parallel computing and steadily increasing computation speed have led to a new tool for analyzing multiple datasets and datatypes: fitting several datasets simultaneously. With this technique, physically connected parameters of individual data can be treated as a single parameter by implementing this connection into the fit directly. We discuss the terminology, implementation, and possible issues of simultaneous fits based on the X-ray data analysis tool Interactive Spectral Interpretation System (ISIS). While all data modeling tools in X-ray astronomy allow in principle fitting data from multiple data sets individually, the syntax used in these tools is not often well suited for this task. Applying simultaneous fits to the transient X-ray binary GRO J1008-57, we find that the spectral shape is only dependent on X-ray flux. We determine time independent parameters such as, e.g., the folding energy E_fold, with unprecedented precision.

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

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.

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.

We present a scheme to determine the spectral state of the canonical black hole Cyg X-1 using data from previous and current X-ray all sky monitors (RXTE-ASM, Swift-BAT, MAXI, and Fermi-GBM). State determinations of the hard/intermediate and soft state agree to better than 10% between different monitors, facilitating the determination of the state and its context for any observation of the source, potentially over the lifetimes of different individual monitors. A separation of the hard and the intermediate state, which strongly differ in their spectral shape and short-term timing behavior, is only possible when monitor data in the soft X-rays (<5 keV) are available. A statistical analysis of the states confirms the different activity patterns of the source (e.g., months to years long hard state periods or phases during which numerous transitions occurs). It also shows the hard and soft states to be stable, with the probability of Cyg X-1 remaining in a given state for at least one week to be larger than 85% for the hard state and larger than 75%, for the soft state. Intermediate states are short lived, with a probability of 50% that the source leaves the intermediate state within three days. A reliable detection of these potentially short-lived events is only possible with monitor data with a time resolution of better than 1 d.

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.