A summary is given of the present state of our knowledge of High-Mass X-ray Binaries (HMXBs), their formation and expected future evolution. Among the HMXB-systems that contain neutron stars, only those that have orbital periods upwards of one year will survive the Common-Envelope (CE) evolution that follows the HMXB phase. These systems may produce close double neutron stars with eccentric orbits. The HMXBs that contain black holes do not necessarily evolve into a CE phase. Systems with relatively short orbital periods will evolve by stable Roche-lobe overflow to short-period Wolf-Rayet (WR) X-ray binaries containing a black hole. Two other ways for the formation of WR X-ray binaries with black holes are identified: CE-evolution of wide HMXBs and homogeneous evolution of very close systems. In all three cases, the final product of the WR X-ray binary will be a double black hole or a black hole neutron star binary.

Understanding the complex behavior of High Mass X-ray binaries (HMXBs) is not possible without detailed information about their donor stars. While crucial, this turns out to be a challenge on multiple fronts. First, multi-wavelength spectroscopy is vital. As such systems can be highly absorbed, this is often already hard to accomplish. Secondly, even if the spectroscopic data is available, the determination of reliable stellar parameters requires sophisticated model atmospheres that accurately describe the outermost layers and the wind of the donor star.

For early-type donors, the stellar wind is radiatively driven and there is a smooth transition between the outermost layers of the star and the wind. The intricate non-LTE conditions in the winds of hot stars complicate the situation even further, as proper model atmospheres need to account for a multitude of physics to accurately provide stellar and wind parameters. The latter are especially crucial for the so-called “wind-fed” HXMBs, where the captured wind of the supergiant donor is the only source for the material accreted by the compact object.

In this review I will briefly address the different approaches for treating stellar winds in the analysis of HMXBs. The fundamentals of stellar atmosphere modeling will be discussed, also addressing the limitations of modern models. Examples from recent analysis results for particular HMXBs will be outlined. Furthermore, the path for the next generation of stellar atmosphere models will be outlined, where models can be used not only for measurement purposes, but also to make predictions and provide a laboratory for theoretical conclusions. Stellar atmospheres are a key tool in understanding HMXBs, e.g. by providing insights about the accretion of stellar winds onto the compact object, or by placing the studied systems in the correct evolutionary context in order to identify potential gravitational wave (GW) progenitors.

Winds of hot stars are driven by the radiative force due to absorption of light in the lines of heavier elements. Consequently, the mass-loss rate and the wind velocity depend on the ionization state of the wind. As a result of this, there is a feedback between the ionizing X-ray source and the stellar wind in HMXBs powered by wind accretion. We study the influence of the small-scale wind structure (clumping) on this feedback using our NLTE hydrodynamical wind models. We find that clumping weakens the effect of X-ray irradiation. Moreover, we show that the observed X-ray luminosities of HMXBs can not be explained by wind accretion scenario without introducing the X-ray feedback. Taking into account the feedback, the observed and estimated X-ray luminosities nicely agree. We identify two cases of X-ray feedback with low and high X-ray luminosities that can explain the dichotomy between SFXTs and sgXBs.

Supergiant X-ray Binaries host a compact object, generally a neutron star, orbiting an evolved O/B star. Mass transfer proceeds through the intense radiatively-driven wind of the stellar donor, a fraction of which is captured by the gravitational field of the neutron star. The subsequent accretion process onto the neutron star is responsible for the abundant X-ray emission from those systems. They also display variations in time of the X-ray flux by a factor of a few 10, along with changes in the hardness ratios believed to be due to varying absorption along the line-of-sight. We used the most recent results on the inhomogeneities (aka clumps) in the non-stationary wind of massive hot stars to evaluate their impact on the time-variable accretion process. We ran three-dimensional simulations of the wind in the vicinity of the accretor to witness the formation of the bow shock and follow the inhomogeneous flow over several spatial orders of magnitude, down to the neutron star magnetosphere. In particular, we show that the impact of the clumps on the time-variability of the intrinsic mass accretion rate is severely damped by the crossing of the shock, compared to the purely ballistic Bondi-Hoyle-Lyttleton estimation. We also account for the variable absorption due to clumps passing by the line-of-sight and estimate the final effective variability of the mass accretion rate for different orbital separations. These results are confronted to recent analysis of Vela X-1 observations with Chandra by Grinberg et al. (2017). It shows that clumps account well for time-variability at low luminosity but can not generate, per se, the high luminosity activity observed.

Previous circular polarization observations obtained with the ESO FOcal Reducer low dispersion spectrograpgh at the VLT in 2007–2008 revealed the presence of a weak longitudinal magnetic field on the surface of the optical component of the X-ray binary Cyg X-1, which contains a black hole and an O9.7Iab supergiant on a 5.6 d orbit. In this contribution we report on recently acquired FORS 2 spectropolarimetric observations of Cyg X-1 along with measurements of a few additional high-mass X-ray binaries.

A small subset of Galactic O-stars possess surface magnetic fields that alter the outflowing stellar wind by magnetically confining it. Key to the magnetic confinement is that it induces rotational modulation of spectral lines over the full EM domain; this allows us to infer basic quantities, e.g., mass-loss rate and magnetic geometry. Here, we present an empirical study of the Hα line in Galactic magnetic O-stars to constrain the mass fed from the stellar base into the magnetosphere, using realistic multi-dimensional magnetized wind models, and compare with theoretical predictions. Our results suggest that it may be reasonable to use mass-feeding rates from non-magnetic wind theory if the absolute mass-loss rate is scaled down according to the amount of wind material falling back upon the stellar surface. This provides then some empirical support to the proposal that such magnetic O-stars might evolve into heavy stellar-mass black holes (Petit et al.2017).

INTEGRAL has been observing the γ-ray sky for 15 years and has discovered many high-energy sources of various nature. Among them, active galactic nuclei (AGN), low or high-mass X-ray binaries (LMXB and HMXB) and cataclysmic variables (CV) are rather difficult to differentiate from one another at high energies and require further optical or near-infrared observations to constrain their exact nature. Using near-infrared photometric and spectroscopic data from ESO VLT/ISAAC, we aim to reveal the nature of 14 high-energy INTEGRAL sources and improve the census of X-ray binaries. By comparing their spectral features to stellar spectra atlases, we identified 5 new CVs, 2 low or intermediate mass X-ray binaries, 2 HMXBs and 5 AGNs.

The progenitors of many core-collapse supernovae (CCSNe) are expected to be in binary systems. By performing a series of three-dimensional hydrodynamical simulations, we investigate how CCSN explosions affect their binary companion. We find that the amount of removed stellar mass, the resulting impact velocity, and the chemical contamination of the companion that results from the impact of the SN ejecta, strongly increases with decreasing binary separation and increasing explosion energy. Also, it is foud that the impact effects of CCSN ejecta on the structure of main-sequence (MS) companions, and thus their long term post-explosion evolution, are in general not dramatic.

GX 301–2, a bright high-mass X-ray binary with an orbital period of 41.5 days, exhibits stable periodic orbital intensity modulations with a strong pre-periastron X-ray flare. Several models have been proposed to explain the accretion at different orbital phases. In Islam & Paul (2014), we presented results from an orbital resolved spectroscopic study of GX 301–2 using data from MAXI Gas Slit Camera. We have found a strong orbital dependence of the absorption column density and equivalent width of the iron emission line. A very large equivalent width of the iron line along with a small value of the column density in the orbital phase range 0.1–0.3 after the periastron passage indicates the presence of high density accretion stream. We aim to further investigate the characteristics of the accretion stream with an AstroSat observation of the system.

Eta Carinae is the most massive active binary within 10,000 light-years. While famous for the largest non-terminal stellar explosion ever recorded, observations reveal a supermassive (∼120 M⊙) binary consisting of an LBV and either a WR or extreme O star in a very eccentric orbit (e=0.9) with a 5.54-year period. Dramatic changes across multiple wavelengths are routinely observed as the stars move about in their highly elliptical orbits, especially around periastron when the hot (∼40 kK) companion star delves deep into the denser and much cooler (∼15 kK) extended wind photosphere of the LBV primary. Many of these changes are due to a dynamic wind-wind collision region (WWCR) that forms between the stars, plus expanding radiation-illuminated fossil WWCRs formed one, two, and three 5.54-year orbital cycles ago. These fossil WWCRs have been spatially and spectrally resolved by the Hubble Space Telescope/Space Telescope Imaging Spectrograph (HST/STIS) at multiple epochs, resulting in data cubes that spatially map Eta Carinae’s innermost WWCRs and follow temporal changes in several forbidden emission lines (e.g. [Fe iii] 4659 Å, [Fe ii] 4815 Å) across the 5.54-year cycle. We present initial results of 3D time-dependent hydrodynamical and radiative-transfer simulations of the Eta Carinae binary and its WWCRs with the goal of producing synthetic data cubes of forbidden emission lines for comparison to the available HST/STIS observations. Comparison of the theoretical models to the observations reveals important details about the binary’s orbital motion, photoionization properties, and recent (5–15year) mass loss history. Such an analysis also provides a baseline for following future changes in Eta Carinae, essential for understanding the late-stage evolution of a nearby supernova progenitor. Our modeling methods can also be adapted to a number of other colliding wind binary systems (e.g. WR 140) that are scheduled to be studied with future observatories (e.g. the James Webb Space Telescope).

Many high-mass X-ray binaries (HMXBs) are runaways. Stellar wind and radiation of donor stars in HMXBs along with outflows and jets from accretors interact with the local interstellar medium and produce curious circumstellar structures. Several such structures are presented and discussed in this contribution.

We propose that a significant fraction of the wide massive binaries in the field are formed as a result of the disintegration of multiple systems of trapezium type. As examples we discuss here the binaries formed from the evolution of the mini-cluster associated with the B component of the Orion Trapezium, from that of the Orion Trapezium itself, and from 10 additional massive trapezia for which we found reliable data in the literature.

Wolf-Rayet (WR) stars are stripped stellar cores that form through strong stellar wind or binary mass transfer. It is proposed that binary evolution plays a vital role in the formation of WR stars in low metallicity environments due to the metallicity dependance of stellar winds. However observations indicate a similar binary fraction of WR stars in the Small Magellanic Cloud (SMC) compared to the Milky Way. There are twelve WR stars in the SMC and five of them are members of binary systems. One of them (SMC AB8) harbors a WO type star. In this work we explore possible formation channels of this binary. We use the MESA code to compute large grids of binary evolution models, and then use least square fitting to compare our models with the observations. In order to reproduce the key properties of SMC AB8, we require efficient semiconvection to produce a sufficiently large convective core, as well as a longer He-burning lifetime. We also need a high mass loss rate during the WN stage to assist the removal of the outer envelope. In this way, we can reproduce the observed properties of AB8, except for the surface carbon to oxygen ratio, which requires further investigation.

Massive stars lose a considerable amount of mass during their lifetime. When the star explodes as a supernova (SN), the resulting shock wave expands in the medium created by the stellar mass-loss. Thermal X-ray emission from the SN depends on the square of the density of the ambient medium, which in turn depends on the mass-loss rate (and velocity) of the progenitor wind. The emission can therefore be used to probe the stellar mass-loss in the decades or centuries before the star’s death.

We have aggregated together data available in the literature, or analysed by us, to compute the X-ray lightcurves of almost all young supernovae detectable in X-rays. We use this database to explore the mass-loss rates of massive stars that collapse to form supernovae. Mass-loss rates are lowest for the common Type IIP supernovae, but increase by several orders of magnitude for the highest luminosity X-ray SNe.

The stellar winds of hot stars have an important impact on both stellar and galactic evolution, yet their structure and internal processes are not fully understood in detail. One of the best nearby laboratories for studying such massive stellar winds is the O4I(n)fp star ζ Pup. After briefly discussing existing X-ray observations from Chandra and XMM, we present a simulation of X-ray emission line profile measurements for the upcoming 840 kilosecond Chandra HETGS observation. This simulation indicates that the increased S/N of this new observation will allow several major steps forward in the understanding of massive stellar winds. By measuring X-ray emission line strengths and profiles, we should be able to differentiate between various stellar wind models and map the entire wind structure in temperature and density. This legacy X-ray spectrum of ζ Pup will be a useful benchmark for future X-ray missions.

Using high resolution 3D hydrodynamical simulations we quantify the amount of mass accreted onto the secondary star of the binary system η Carinae during periastron passage on its highly eccentric orbit. The accreted mass is responsible for the spectroscopic event occurring every orbit close to periastron passage, during which many lines vary and the x-ray emission associated with the destruction wind collision structure declines. The system is mainly known for its giant eruptions that occurred in the nineteenth century. The high mass model of the system, M1=170M⊙ and M2=80M⊙, gives Macc≍ 3×10−6M⊙ compatible with the amount required for explaining the reduction in secondary ionization photons during the spectroscopic event, and also matches its observed duration. As accretion occurs now, it surely occurred during the giant eruptions. This implies that mass transfer can have a huge influence on the evolution of massive stars.

In order to test massive star evolution above 25 M⊙, we perform spectral analysis on a sample of massive stars in the Small Magellanic Cloud that includes both O stars as well as more evolved Wolf-Rayet stars. We present a grid of non-LTE stellar atmospheres that has been calculated using the cmfgen code, in order to have a systematic and homogeneous approach. We obtain stellar and wind parameters for O stars, spectral types ranging from O2 to O6, and the complete sample of known Wolf-Rayet stars. We discuss the evolutionary status of both the O and WR stars and the links between them, as well as the most likely evolutionary path towards black hole formation in a low metallicity environment, including testing theoretical predictions for mass-loss rates at low metallicities.

Rapidly rotating B-type stars with gaseous mass-loss disks in Keplerian rotation are common central objects in X-Ray binaries. These disks are physically well understood in the framework of the viscous decretion disk, and their typical parameters have been established for a large number of single Be stars in the recent years. According to the current observational evidence, the Be stars and disks found in BeXRBs are well within the boundaries known from single Be stars, i.e., they are normal Be stars. New results have also been obtained on the orbital disk truncation and other tidal effects of the companion objects on the disk.

Because of their small angular size < 1 mas, spatial information on High-mass X-ray binaries (HMXB) has typically been inferred from photometry or spectroscopy. Optical interferometry offers the possibility to spatially resolve such systems, but has been traditionally limited to bright targets or low spectral resolution. The VLTI instrument GRAVITY, working in the near-infrared K band, achieves unprecedented precision in differential interferometric quantities at high spectral resolution, allowing to study HMXBs through the lens of optical interferometry for the first time. We present GRAVITY observations on two X-ray binaries: the microquasar SS 433 and the supergiant HMXB BP Cru. The former is the only known steady super-Eddington accretor in the Galaxy and is in a unique stage of binary evolution, with probable ties to at least part of the ULX population. With GRAVITY, we resolve its massive winds and optical baryonic jets for the first time, finding evidence for powerful equatorial outflows and photoionization as the main heating process along the jets. BP Cru harbors an X-ray pulsar accreting from the wind of its early-blue hypergiant companion Wray 977. The GRAVITY observations resolve the inner parts of the stellar wind and allow probing the influence of the orbiting pulsar on the circumstellar environment.

We present a case study of GG Carinae (GG Car), a Galactic B[ e ] supergiant binary having significant eccentricity (0.28), based on Global Jet Watch spectroscopy data which has been collecting high-time-sampled optical spectra since early 2015. GG Car has so far not been observed in the X-ray band, however it is of similar phenomenology to known X-ray binaries and may therefore be an obscured X-ray source. We have discovered that the absorption component of the H-alpha line displays a ∼62-478-day period in both equivalent width and wavelength centroid indicating cycles in the dynamics of the circumstellar environment, such as precession of the circumbinary or circumprimary disk. Circumbinary disk precession is an as-of-yet underexplored origin of super-orbital variations in the X-ray flux of X-ray binaries, since the rate of precession is generally much longer than the orbital period of the inner binary.