Wolf-Rayet stars are advanced evolutionary stages of massive stars. Despite their large mass-loss rates and high wind velocities, none of them display a bow shock, although a fraction of them are classified as runaway. Our 2.5-D numerical simulations of circumstellar matter around a 60Mo runaway star show that the fast Wolf-Rayet stellar wind is released into a wind-blown cavity filled with various shocks and discontinuities generated throughout the precedent evolutionary phases. The resulting fast-wind slow-wind interaction leads to the formation of spherical shells of swept-up dusty material similar to those observed in near infrared 24 micron with Spitzer, and which appear to be co-moving with the runaway massive stars, regardless of their proper motion and/or the properties of the local ambient medium. We interpret bright infrared rings around runaway Wolf-Rayet stars in the Galactic plane, like WR138a, as indication of their very high initial masses and a complex evolutionary history. Stellar-wind bow shocks become faint as stars run in diluted media, therefore, our results explain the absence of detected bow shocks around Galactic Wolf-Rayet stars such as the high-latitude, very fast-moving objects WR71, WR124 and WR148. Our results show that the absence of a bow shock is consistent with a runaway nature of some Wolf-Rayet stars. This questions the in-situ star formation scenario of high-latitude Wolf-Rayet stars in favor of dynamical ejection from birth sites in the Galactic plane.

The mass, origin and evolutionary stage of the binary system LB-1 has been the subject of intense debate, following the claim that it hosts an $\sim$70$M_{\odot}$ black hole, in stark contrast with the expectations for stellar remnants in the Milky Way. We conducted a high-resolution, phase-resolved spectroscopic study of the near-infrared Paschen lines in this system, using the 3.5-m telescope at Calar Alto Observatory. We find that Pa$\beta$ and Pa$\gamma$ (after proper subtraction of the stellar absorption component) are well fitted with a standard double-peaked model, typical of disk emission. We measured the velocity shifts of the red and blue peaks at 28 orbital phases: the line center has an orbital motion in perfect antiphase with the stellar motion, and the radial velocity amplitude ranges from 8 to 13 km/s for different choices of lines and profile modelling. We interpret this curve as proof that the disk is tracing the orbital motion of the primary, ruling out the circumbinary disk and the hierarchical triple scenarios. The phase-averaged peak-to-peak half-separation (proxy for the projected rotational velocity of the outer disk) is $\sim$70 km s$^{-1}$, larger than the stellar orbital velocity and also inconsistent with a circumbinary disk. From those results, we infer a primary mass 4--8 times higher than the secondary mass. Moreover, we show that the ratio of the blue and red peaks (V/R intensity ratio) has a sinusoidal behaviour in phase with the secondary star, which can be interpreted as the effect of external irradiation by the secondary star on the outer disk. Finally, we briefly discuss our findings in the context of alternative scenarios recently proposed for LB-1. Definitive tests between alternative solutions will require further astrometric data from $Gaia$.

We have obtained a very deep exposure (813 ks) of $\zeta\,$Puppis (O4 supergiant) with the Chandra/HETG Spectrometer. Here we report on analysis of the 1-9 \r{A} region, especially well suited for Chandra, which has a significant contribution from continuum emission between well separated emission lines from high-ionization species. These data allow us to study the hottest plasma present through the continuum shape and emission line strengths. Assuming a powerlaw emission measure distribution which has a high-temperature cut-off, we find that the emission is consistent with a thermal spectrum having a maximum temperature of 12 MK. This implies an effective wind shock velocity of $900\,\mathrm{km\,s^{-1}}$, well below the wind terminal speed of $2250\,\mathrm{km\,s^{-1}}$. For X-ray emission which forms close to the star, the speed and X-ray flux are larger than can be easily reconciled with strictly self-excited line-deshadowing-instability models, suggesting a need for a fraction of the wind to be accelerated extremely rapidly right from the base. This is not so much a dynamical instability as a nonlinear response to changing boundary conditions.

We present a multiwavelength study of the iconic Bubble Nebula (NGC 7635) and its ionising star BD$+$60$^{\circ}$2522. We obtained XMM-Newton EPIC X-ray observations to search for extended X-ray emission as in other similar wind-blown bubbles around massive stars. We also obtained San Pedro M\'{a}rtir spectroscopic observations with the Manchester Echelle Spectrometer to study the dynamics of the Bubble Nebula. Although our EPIC observations are deep, we do not detect extended X-ray emission from this wind-blown bubble. On the other hand, BD$+$60$^{\circ}$2522 is a bright X-ray source similar to other O stars. We used the stellar atmosphere code PoWR to characterise BD$+$60$^{\circ}$2522 and found that this star is a young O-type star with stellar wind capable of producing a wind-blown bubble that in principle could be filled with hot gas. We discussed our findings in line with recent numerical simulations proposing that the Bubble Nebula has been formed as the result of the fast motion of BD$+$60$^{\circ}$2522 through the medium. Our kinematic study shows that the Bubble Nebula is composed by a series of nested shells, some showing blister-like structures, but with little signatures of hydrodynamical instabilities that would mix the material producing diffuse X-ray emission as seen in other wind-blown bubbles. Its morphology seems to be merely the result of projection effects of these different shells.

We present new ATCA multi-wavelength radio measurements (range 2.1-21.2 GHz) of the early-type magnetic star rho Oph A, performed in March 2019 during 3 different observing sessions. These new ATCA observations evidence a clear rotational modulation of the stellar radio emission and the detection of coherent auroral radio emission from rho Oph A at 2.1 GHz. We collected high-resolution optical spectra of rho Oph A acquired by several instruments over a time span of about ten years. We also report new magnetic field measurements of rho Oph A that, together with the radio light curves and the temporal variation of the equivalent width of the HeI line (lambda=5015 Angstrom), were used to constrain the rotation period and the stellar magnetic field geometry. The above results have been used to model the stellar radio emission, modelling that allowed us to constrain the physical condition of rho Oph A's magnetosphere. Past XMM measurements showed periodic X-ray pulses from rho Oph A. We correlate the X-ray light curve with the magnetic field geometry of rho Oph A. The already published XMM data have been re-analyzed showing that the X-ray spectra of rho Oph A are compatible with the presence of a non-thermal X-ray component. We discuss a scenario where the emission phenomena occurring at the extremes of the electromagnetic spectrum, radio and X-ray, are directly induced by the same plasma process. We interpret the observed X-ray and radio features of rho Oph A as having an auroral origin.

High-mass X-ray binaries (HMXBs) are exceptional astrophysical laboratories that offer a rare glimpse into the physical processes that govern accretion on compact objects, massive-star winds, and stellar evolution. In a subset of the HMXBs, the compact objects accrete matter solely from winds of massive donor stars. These so-called wind-fed HMXBs are divided in persistent HMXBs and supergiant fast X-ray transients (SFXTs) according to their X-ray properties. While it has been suggested that this dichotomy depends on the characteristics of stellar winds, they have been poorly studied. With this investigation, we aim to remedy this situation by systematically analyzing donor stars of wind-fed HMXBs that are observable in the UV, concentrating on those with neutron star (NS) companions. We obtained Swift X-ray data, HST UV spectra, and additional optical spectra for all our targets. Our multi-wavelength approach allows us to provide stellar and wind parameters for six donor stars (four wind-fed systems and two OBe X-ray binaries). The wind properties are in line with the predictions of the line-driven wind theory. Three of the donor stars are in an advanced evolutionary stage, while for some of the stars, the abundance pattern indicates that processed material might have been accreted. When passing by the NS in its tight orbit, the donor star wind has not yet reached its terminal velocity but it is still significantly slower; its speed is comparable with the orbital velocity of the NS companion. There are no systematic differences between the two types of wind-fed HMXBs (persistent versus transients) with respect to the donor stars. For the SFXTs in our sample, the orbital eccentricity is decisive for their transient X-ray nature. Based on the orbital parameters and the further evolution of the donor stars, the investigated HMXBs will presumably form Thorne-\.Zytkow objects in the future.

The supergiant ionized shell SMC-SGS 1 (DEM 167), located in the outer Wing of the Small Magellanic Cloud (SMC), resembles structures that originate from an energetic star-formation event and later stimulate star formation as they expand into the ambient medium. However, stellar populations within and surrounding SMC-SGS 1 tell a different story. We present a photometric study of the stellar population encompassed by SMC-SGS 1 in order to trace the history of this structure and its potential influence on star formation within the low-density, low-metallicity SMC Wing. For a stellar population that is physically associated with SMC-SGS 1, we combined near-ultraviolet (NUV) photometry from the Galaxy Evolution Explorer (GALEX) with archival optical (V-band) photometry from the ESO Danish 1.54m Telescope. Given their colors and luminosities, we estimated stellar ages and masses by matching observed photometry to theoretical stellar isochrone models. We find that the investigated region supports an active, extended star-formation event spanning $\sim$ 25 - 40 Myr ago, as well as continued star formation into the present. Using a standard initial mass function (IMF), we infer a lower bound on the stellar mass from this period of $\sim 3 \times 10^4 M_{\odot}$, corresponding to a star-formation intensity of $\sim$ 6 $\times$ 10$^{-3}$ M$_{\odot}$ kpc$^{-2}$ yr$^{-1}$. The spatial and temporal distributions of young stars encompassed by SMC-SGS 1 imply a slow, consistent progression of star formation over millions of years. Ongoing star formation along the edge of and interior to SMC-SGS 1 suggests a combined stimulated and stochastic mode of star formation within the SMC Wing. A slow expansion of the shell within this low-density environment may preserve molecular clouds within the volume of the shell, leaving them to form stars even after nearby stellar feedback expels local gas and dust.

Massive Wolf-Rayet (WR) stars dominate the radiative and mechanical energy budget of galaxies and probe a critical phase in the evolution of massive stars prior to core-collapse. It is not known whether core He-burning WR stars (classical WR, cWR) form predominantly through wind-stripping (w-WR) or binary stripping (b-WR). With spectroscopy of WR binaries so-far largely avoided due to its complexity, our study focuses on the 44 WR binaries / binary candidates of the Large Magellanic Cloud (LMC, metallicity Z~0.5 Zsun), identified on the basis of radial velocity variations, composite spectra, or high X-ray luminosities. Relying on a diverse spectroscopic database, we aim to derive the physical and orbital parameters of our targets, confronting evolution models of evolved massive stars at sub-solar metallicity, and constraining the impact of binary interaction in forming them. Spectroscopy is performed using the Potsdam Wolf-Rayet (PoWR) code and cross-correlation techniques. Disentanglement is performed using the code Spectangular or the shift-and-add algorithm. Evolutionary status is interpreted using the Binary Population and Spectral Synthesis (BPASS) code, exploring binary interaction and chemically-homogeneous evolution. No obvious dichotomy in the locations of apparently-single and binary WN stars on the Hertzsprung-Russell diagram is apparent. According to commonly used stellar evolution models (BPASS, Geneva), most apparently-single WN stars could not have formed as single stars, implying that they were stripped by an undetected companion. Otherwise, it must follow that pre-WR mass-loss/mixing (e.g., during the red supergiant phase) are strongly underestimated in standard stellar evolution models.

X-ray radiation from accreting compact objects is an important part of stellar feedback. The metal-poor galaxy ESO 338-4 has experienced vigorous starburst during the last 40 Myr and contains some of the most massive super star clusters in the nearby Universe. Given its starburst age and its star-formation rate, ESO 338-4 is one of the most efficient nearby manufactures of neutron stars and black holes, hence providing an excellent laboratory for feedback studies. We compared X-ray images and spectra obtained by XMM-Newton and Chandra telescopes with integral field spectroscopic VLT MUSE observations in the optical to constrain the nature of strong X-ray emitters. X-ray observations uncover three ultraluminous X-ray sources (ULXs) in ESO 338-4. The brightest among them, ESO 338~X-1, has X-ray luminosity in excess of 10^{40} erg/s. We speculate that ESO 338-4 is powered by accretion on an intermediate-mass (~300Msun) black hole. We show that X-ray radiation from ULXs and hot superbubbles strongly contributes to HeII ionization and general stellar feedback in this template starburst galaxy.

Stars which start their lives with spectral types O and early-B are the progenitors of core-collapse supernovae, long gamma-ray bursts, neutron stars, and black holes. These massive stars are the primary sources of stellar feedback in star-forming galaxies. At low metallicities, the properties of massive stars and their evolution are not yet fully explored. Here we report a spectroscopic study of 320 OB stars in the Small Magellanic Cloud. The data, which we obtained with the ESO Very Large Telescope, were analyzed using state-of-the-art stellar atmosphere models. We find that stellar winds of our sample stars are much weaker than theoretically expected. The stellar rotation rates show a bi-modal distribution. The well-populated upper Hertzsprung-Russell diagram including our sample OB stars from SMC Wing as well as additional evolved stars all over SMC from the literature shows a strict luminosity limit. The comparison with single-star evolutionary tracks suggests a dichotomy in the fate of massive stars in the SMC. Only stars with Minit<30M$_{\odot}$ seem to evolve from the main sequence to the cool side of the HRD to become a red supergiant and to explode as type II-P supernova. In contrast, stars with Minit>30M$_{\odot}$ appear to stay always hot and might evolve quasi chemically homogeneously, finally collapsing to relatively massive black holes. However, we find no indication that chemical mixing is correlated with rapid rotation. We report extended star-formation episodes in a quiescent low-density region of the Wing, which is progressing stochastically. We measure the key parameters of stellar feedback and establish the links between the rates of star formation and supernovae. Our study reveals that in metal-poor environments the stellar feedback is dominated by core-collapse supernovae in combination with winds and ionizing radiation supplied by a few of the most massive stars.

Comprehensive spectral analyses of the Galactic Wolf-Rayet stars of the nitrogen sequence (i.e.\ the WN subclass) have been performed in a previous paper. However, the distances of these objects were poorly known. Distances have a direct impact on the "absolute" parameters, such as luminosities and mass-loss rates. The recent Gaia Data Release (DR2) of trigonometric parallaxes includes nearly all WN stars of our Galactic sample. In the present paper, we apply the new distances to the previously analyzed Galactic WN stars and rescale the results accordingly. On this basis, we present a revised catalog of 55 Galactic WN stars with their stellar and wind parameters. The correlations between mass-loss rate and luminosity show a large scatter, for the hydrogen-free WN stars as well as for those with detectable hydrogen. The slopes of the $\log L - \log \dot{M}$ correlations are shallower than found previously. The empirical Hertzsprung-Russell diagram (HRD) still shows the previously established dichotomy between the hydrogen-free early WN subtypes that are located on the hot side of the zero-age main sequence (ZAMS), and the late WN subtypes, which show hydrogen and reside mostly at cooler temperatures than the ZAMS (with few exceptions). However, with the new distances, the distribution of stellar luminosities became more continuous than obtained previously. The hydrogen-showing stars of late WN subtype are still found to be typically more luminous than the hydrogen-free early subtypes, but there is a range of luminosities where both subclasses overlap. The empirical HRD of the Galactic single WN stars is compared with recent evolutionary tracks. Neither these single-star evolutionary models nor binary scenarios can provide a fully satisfactory explanation for the parameters of these objects and their location in the HRD.

Over the last two decades, the discovery of exoplanets has fundamentally changed our perception of the universe and humanity's place within it. Recent work indicates that a solar system's X-ray and high energy particle environment is of fundamental importance to the formation and development of the atmospheres of close-in planets such as hot Jupiters, and Earth-like planets around M stars. X-ray imaging and spectroscopy provide powerful and unique windows into the high energy flux that an exoplanet experiences, and X-ray photons also serve as proxies for potentially transfigurative coronal mass ejections. Finally, if the host star is a bright enough X-ray source, transit measurements akin to those in the optical and infrared are possible and allow for direct characterization of the upper atmospheres of exoplanets. In this brief white paper, we discuss contributions to the study of exoplanets and their environs which can be made by X-ray data of increasingly high quality that are achievable in the next 10--15 years.

It is now recognized that energetic stellar photon and particle radiation evaporates and erodes planetary atmospheres and controls upper atmospheric chemistry. Key exoplanet host stars will be too faint at X-ray wavelengths for accurate characterization using existing generation and future slated X-ray telescopes. Observation of stellar coronal mass ejections and winds are also beyond current instrumentation. In line with theCommittee on an Exoplanet Science Strategy recognition that holistic observational approaches are needed, we point out here that a full understanding of exoplanet atmospheres, their evolution and determination of habitability requires a powerful high-resolution X-ray imaging and spectroscopic observatory. This is the only capability that can: (1) characterize by proxy the crucial, difficult to observe, EUV stellar flux, its history and its variations for planet hosting stars; (2) observe the stellar wind; (3) detect the subtle Doppler signatures of coronal mass ejections.

We report the discovery of a circular shell centred on the Be X-ray binary (BeXB) SXP 1323 in the Small Magellanic Cloud (SMC). The shell was detected in an Halpha image obtained with the Very Large Telescope (VLT). Follow-up spectroscopy with the Southern African Large Telescope (SALT) showed that the shell expands with a velocity of $\approx$ 100 km/s and that its emission is due to shock excitation. We suggest that this shell is the remnant of the supernova explosion that led to the formation of the SXP 1323's neutron star $\approx$ 40 000 yr ago. SXP 1323 represents the second known case of a BeXB associated with a supernova remnant (the first one is SXP 1062). Interestingly, both these BeXBs harbour long period pulsars and are located in a low-metallicity galaxy.

We present new, contemporaneous HST STIS and XMM observations of the O7 III(n)((f)) star xi Per. We supplement the new data with archival IUE spectra, to analyze the variability of the wind lines and X-ray flux of xi Pper. The variable wind of this star is known to have a 2.086 day periodicity. We use a simple, heuristic spot model which fits the low velocity (near surface) IUE wind line variability very well, to demonstrate that the low velocity absorption in the new STIS spectra of N IV 1718 and Si IV 1402 vary with the same 2.086 day period. It is remarkable that the period and amplitude of the STIS data agree with those of the IUE spectra obtained 22 years earlier. We also show that the time variability of the new XMM fluxes are also consistent with the 2.086 day period. Thus, our new, multi-wavelength coordinated observations demonstrate that the mechanism which causes the UV wind line variability is also responsible for a significant fraction of the X-rays in single O stars. The sequence of events for the multi-wavelength light curve minima is: Si IV 1402, N IV 1718, and X-ray flux, each separated by a phase of about 0.06 relative to the 2.086 day period. Analysis of the X-ray fluxes shows that they become softer as they weaken. This is contrary to expectations if the variability is caused by periodic excess absorption. Further, the high resolution X-ray spectra suggest that the individual emission lines at maximum are more strongly blue shifted. If we interpret the low velocity wind line light curves in terms of our model, it implies that there are two bright regions, i.e., regions with less absorption, separated by 180 deg, on the surface of the star. We note that the presence and persistent of two spots separated by 180 deg suggests that a weak dipole magnetic field is responsible for the variability of the UV wind line absorption and X-ray flux in xi Per.

HD 93129A was classified as the earliest O-type star in the Galaxy (O2~If*) and is considered as the prototype of its spectral class. However, interferometry shows that this object is a binary system, while recent observations even suggest a triple configuration. None of the previous spectral analyses of this object accounted for its multiplicity. With new high-resolution UV and optical spectra, we have the possibility to reanalyze this key object, taking its binary nature into account for the first time. We aim to derive the fundamental parameters and the evolutionary status of HD 93129A, identifying the contributions of both components to the composite spectrum. We analyzed UV and optical observations acquired with the Hubble Space Telescope and ESO's Very Large Telescope. A multiwavelength analysis of the system was performed using the latest version of the Potsdam Wolf-Rayet model atmosphere code. Despite the similar spectral types of the two components, we are able to find signatures from each of the components in the combined spectrum, which allows us to estimate the parameters of both stars. We derive $\log (L/L_\odot) = 6.15$, $T_{\textrm{eff}}=52$ kK, and $\log \dot{M}=-4.7 [M_\odot\text{yr}^{-1}]$ for the primary Aa, and $\log (L/L_\odot)=5.58$, $T_{\textrm{eff}}=45$ kK, and $\log\dot{M}=-5.8 [M_\odot\text{yr}^{-1}]$ for the secondary Ab. Even when accounting for the binary nature, the primary of HD 93129A is found to be one of the hottest and most luminous O stars in our Galaxy. Based on the theoretical decomposition of the spectra, we assign spectral types O2~If* and O3~III(f*) to components Aa and Ab, respectively. While we achieve a good fit for a wide spectral range, specific spectral features are not fully reproduced. The data are not sufficient to identify contributions from a hypothetical third component in the system.

Wolf-Rayet stars of the carbon sequence (WC stars) are an important cornerstone in the late evolution of massive stars before their core collapse. As core-helium burning, hydrogen-free objects with huge mass-loss, they are likely the last observable stage before collapse and thus promising progenitor candidates for type Ib/c supernovae. Their strong mass-loss furthermore provides challenges and constraints to the theory of radiatively driven winds. Thus, the determination of the WC star parameters is of major importance for several astrophysical fields. With Gaia DR2, for the first time parallaxes for a large sample of Galactic WC stars are available, removing major uncertainties inherent to earlier studies. In this work, we re-examine the sample from Sander et al. (2012) to derive key properties of the Galactic WC population. All quantities depending on the distance are updated, while the underlying spectral analyses remain untouched. Contrasting earlier assumptions, our study yields that WC stars of the same subtype can significantly vary in absolute magnitude. With Gaia DR2, the picture of the Galactic WC population becomes more complex: We obtain luminosities ranging from log L = 4.9 to 6.0 with one outlier having log L = 4.7. This indicates that the WC stars are likely formed from a broader initial mass range than previously assumed. We obtain mass-loss rates ranging between log Mdot = -5.1 and -4.1, with Mdot propto L^0.68 and a linear scaling of the modified wind momentum with luminosity. We discuss the implications for stellar evolution, including unsolved issues regarding the need of envelope inflation to address the WR radius problem, and the open questions in regard to the connection of WR stars with Gamma-ray bursts. WC and WO stars are progenitors of massive black holes, collapsing either silently or in a supernova that most-likely has to be preceded by a WO stage.

Among different types of massive stars in advanced evolutionary stages is the enigmatic WN8h type. There are only a few Wolf-Rayet (WR) stars with this spectral type in our Galaxy. It has long been suggested that WN8h-type stars are the products of binary evolution that may harbor neutron stars (NS). One of the most intriguing WN8h stars is the runaway WR124 surrounded by its magnificent nebula M1-67. We test the presence of an accreting NS companion in WR124 using 100 ks long observations by the Chandra X-ray observatory. The hard X-ray emission from WR124 with a luminosity of $L_\mathrm{X}$= 10$^{31}$ erg s$^{-1}$ is marginally detected. We use the non-LTE stellar atmosphere code PoWR to estimate the WR wind opacity to the X-rays. The wind of a WN8-type star is effectively opaque for X-rays, hence the low X-ray luminosity of WR124 does not rule out the presence of an embedded compact object. We suggest that, in general, high opacity WR winds could prevent X-ray detections of embedded NS, and be an explanation for the apparent lack of WR+NS systems.

The study of massive stars in different metallicity environments is a central topic of current stellar research. The spectral analysis of massive stars requires adequate model atmospheres. The computation of such models is difficult and time-consuming. Therefore, spectral analyses are greatly facilitated if they can refer to existing grids of models. Here we provide grids of model atmospheres for OB-type stars at metallicities corresponding to the Small and Large Magellanic Clouds, as well as to solar metallicity. In total, the grids comprise 785 individual models. The models were calculated using the state-of-the-art Potsdam Wolf-Rayet (PoWR) model atmosphere code. The parameter domain of the grids was set up using stellar evolution tracks. For all these models, we provide normalized and flux-calibrated spectra, spectral energy distributions, feedback parameters such as ionizing photons, Zanstra temperatures, and photometric magnitudes. The atmospheric structures (the density and temperature stratification) are available as well. All these data are publicly accessible through the PoWR website.

We report the detection of the auroral radio emission from the early-type magnetic star HD142301. New VLA observations of HD142301 detected highly polarized amplified emission occurring at fixed stellar orientations. The coherent emission mechanism responsible for the stellar auroral radio emission amplifies the radiation within a narrow beam, making the star where this phenomenon occurs similar to a radio lighthouse. The elementary emission process responsible for the auroral radiation mainly amplifies one of the two magneto-ionic modes of the electromagnetic wave. This explains why the auroral pulses are highly circularly polarized. The auroral radio emission of HD142301 is characterized by a reversal of the sense of polarization as the star rotates. The effective magnetic field curve of HD142301 is also available making it possible to correlate the transition from the left to the right-hand circular polarization sense (and vice-versa) of the auroral pulses with the known orientation of the stellar magnetic field. The results presented in this letter have implications for the estimation of the dominant magneto-ionic mode amplified within the HD142301 magnetosphere.

The X-ray emission of O-type stars was first discovered in the early days of the Einstein satellite. Since then many different surveys have confirmed that the ratio of X-ray to bolometric luminosity in O-type stars is roughly constant, but there is a paucity of studies that account for detailed information on spectral and wind properties of O-stars. Recently a significant sample of O stars within our Galaxy was spectroscopically identified and presented in the Galactic O-Star Spectroscopic Survey (GOSS). At the same time, a large high-fidelity catalog of X-ray sources detected by the XMM-Newton X-ray telescope was released. Here we present the X-ray catalog of O stars with known spectral types and investigate the dependence of their X-ray properties on spectral type as well as stellar and wind parameters. We find that, among the GOSS sample, 127 O-stars have a unique XMM-Newton source counterpart and a Gaia data release 2 (DR2) association. Terminal velocities are known for a subsample of 35 of these stars. We confirm that the X-ray luminosities of dwarf and giant O stars correlate with their bolometric luminosity. For the subsample of O stars with measure terminal velocities we find that the X-ray luminosities of dwarf and giant O stars also correlate with wind parameters. However, we find that these correlations break down for supergiant stars. Moreover, we show that supergiant stars are systematically harder in X-rays compared to giant and dwarf O-type stars. We find that the X-ray luminosity depends on spectral type, but seems to be independent of whether the stars are single or in a binary system. Finally, we show that the distribution of log(Lx/Lbol) in our sample stars is non-Gaussian, with the peak of the distribution at log(Lx/Lbol) around -6.6.

We present X-ray observations of CU Vir performed with XMM-Newton and Chandra. With Lx = 3 x 10^28 erg/s the source is moderately X-ray bright, but its spectrum is extremely hard compared to other Ap stars. Spectral modelling requires multi-component models with predominant hot plasma at temperatures of about Tx = 25 MK or, alternatively, a nonthermal spectral component. The Chandra observations was performed six years later than the one by XMM-Newton, yet the source has similar X-ray flux and spectrum, suggesting a steady and persistent X-ray emission. To explain its full X-ray properties, a generating mechanism beyond standard explanations like the presence of a low-mass companion or magnetically confined wind-shocks is required. Magnetospheric activity might be present or, as proposed for fast rotating strongly magnetic Bp stars, the X-ray emission of CU Vir is predominantly auroral in nature.

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.

We present the results from Chandra X-ray observations, and near- and mid-infrared analysis, using VISTA/VVV and Spitzer/GLIMPSE catalogs, of the high-mass star-forming region IRAS 16562-3959, which contains a candidate for a high mass protostar. We detected 249 X-ray sources within the ACIS-I field-of-view. The majority of the X-ray sources have low count rates (< 0.638 cts/ks), and hard X-ray spectra. The search for YSOs in the region using VISTA/VVV and Spitzer/GLIMPSE catalogs resulted in a total of 799 YSOs, with 66 Class I and 733 Class II YSOs. The search for near- and mid-infrared counterparts of the X-ray sources led to a total of 165 VISTA/VVV counterparts, and a total of 151 Spitzer/GLIMPSE counterparts. The infrared analysis of the X-ray counterparts allowed us to add an extra 123 YSOs associated with the region. We conclude that a total of 922 YSOs are associated with the region, with 67 Class I, 764 Class II, and 91 Class III YSOs. We also found that the region is composed of 11 subclusters. In the vicinity of the high-mass protostar, the stellar distribution has a core-halo structure. The subcluster containing the high-mass protostar is the densest and the youngest in the region, and the high-mass protostar is located at its center. The YSOs in this cluster appear to be substantially older than the high mass protostar.

Classic massive binary evolutionary scenarios predict that a transitional common-envelope (CE) phase could be preceded as well as succeeded by the evolutionary stage when a binary consists of a compact object and a massive star, that is, a high-mass X-ray binary (HMXB). The observational manifestations of common envelope are poorly constrained. We speculate that its ejection might be observed in some cases as a transient event at mid-infrared (IR) wavelengths. We estimate the expected numbers of CE ejection events and HMXBs per star formation unit rate, and compare these theoretical estimates with observations. We compiled a list of 85 mid-IR transients of uncertain nature detected by the Spitzer Infrared Intensive Transients Survey and searched for their associations with X-ray, optical, and UV sources. Confirming our theoretical estimates, we find that only one potential HMXB may be plausibly associated with an IR-transient and tentatively propose that X-ray source NGC 4490-X40 could be a precursor to the SPIRITS 16az event. Among other interesting sources, we suggest that the supernova remnant candidate [BWL2012] 063 might be associated with SPIRITS 16ajc. We also find that two SPIRITS events are likely associated with novae, and seven have potential optical counterparts. We conclude that the massive binary evolutionary scenarios that involve CE events do not contradict currently available observations of IR transients and HMXBs in star-forming galaxies.

SMC AB 6 is the shortest-period (6.5d) Wolf-Rayet (WR) binary in the Small Magellanic Cloud, and is therefore crucial for the study of binary interaction and formation of WR stars at low metallicity. The WR component in AB 6 was previously found to be very luminous (logL=6.3[Lsun]) compared to its reported orbital mass (8Msun), placing it significantly above the Eddington limit. Through spectroscopy and orbital analysis of newly acquired optical data taken with UVES, we aim to understand the peculiar results reported for this system and explore its evolutionary history. Results: We find that AB 6 contains at least four stars. The 6.5d period WR binary comprises the WR primary (WN3:h, star A) and a rather rapidly rotating early O-type companion (O5.5 V, star B). Static N and He lines suggest the presence of an emission line star (O5.5 I(f), star C). Finally, narrow absorption lines portraying a long-term radial velocity variation show the existence of a fourth star (O7.5 V, star D). Star D appears to form a second 140d period binary together with a fifth stellar member, which is a B-type dwarf or a black hole. It is not clear that these additional components are bound to the WR binary. The WR star is found to be less luminous than previously thought (logL = 5.9[Lsun]) and, adopting 41Msun for star B, more massive (18Msun). Correspondingly, the WR star does not exceed the Eddington limit. We derive the initial masses of 60 and 40Msun for stars A and B and an age of 3.9 Myr for the system. The WR binary likely experienced nonconservative mass transfer in the past supported by the relatively rapid rotation of star B. Conclusion: Our study shows that AB 6 is a multiple -- probably quintuple -- system. This finding resolves the previously reported puzzle of the WR primary exceeding the Eddington limit and suggests that the WR star exchanged mass with its companion in the past.

Clusters or associations of early-type stars are often associated with a 'superbubble' of hot gas. The formation of such superbubbles is caused by the feedback from massive stars. The complex N206 in the Large Magellanic Cloud exhibits a superbubble and a rich massive star population. We observed these massive stars using the FLAMES multi-object spectrograph at ESO-VLT. Available UV spectra from HST, IUE, and FUSE are also used. The spectral analysis is performed with Potsdam Wolf-Rayet (PoWR) model atmospheres. We present the stellar and wind parameters of the OB stars and the two WR binaries in the N206 complex. Twelve percent of the sample show Oe/Be type emission lines, although most of them appear to rotate far below critical. We found eight runaway stars based on their radial velocity. The wind-momentum luminosity relation of our OB sample is consistent with the expectations. The HRD of the OB stars reveals a large age spread (1-30 Myr), suggesting different episodes of star formation in the complex. The youngest stars are concentrated in the inner part of the complex, while the older OB stars are scattered over outer regions. We derived the present day mass function for the entire N206 complex as well as for the cluster NGC2018. Three very massive Of stars are found to dominate the feedback among 164 OB stars in the sample. The two WR winds alone release about as much mechanical luminosity as the whole OB star sample. The cumulative mechanical feedback from all massive stellar winds is comparable to the combined mechanical energy of the supernova explosions that likely occurred in the complex. Accounting also for the WR wind and supernovae, the mechanical input over the last five Myr is ~$2.3\times10^{52}$ erg, which exceeds the current energy content of the complex by more than a factor of five. The morphology of the complex suggests a leakage of hot gas from the superbubble.

We present new radio/millimeter measurements of the hot magnetic star HR5907 obtained with the VLA and ALMA interferometers. We find that HR5907 is the most radio luminous early type star in the cm-mm band among those presently known. Its multi-wavelength radio light curves are strongly variable with an amplitude that increases with radio frequency. The radio emission can be explained by the populations of the non-thermal electrons accelerated in the current sheets on the outer border of the magnetosphere of this fast rotating magnetic star. We classify HR5907 as another member of the growing class of strongly magnetic fast rotating hot stars where the gyro-synchrotron emission mechanism efficiently operates in their magnetospheres. The new radio observations of HR5907 are combined with archival X-ray data to study the physical condition of its magnetosphere. The X-ray spectra of HR5907 show tentative evidence for the presence of non-thermal spectral component. We suggest that non-thermal X-rays originate a stellar X-ray aurora due to streams of non-thermal electrons impacting on the stellar surface. Taking advantage of the relation between the spectral indices of the X-ray power-law spectrum and the non-thermal electron energy distributions, we perform 3-D modeling of the radio emission for HR5907. The wavelength-dependent radio light-curves probe magnetospheric layers at different heights above the stellar surface. A detailed comparison between simulated and observed radio light-curves leads us to conclude that the stellar magnetic field of HR5907 is likely non-dipolar, providing further indirect evidence of the complex magnetic field topology of HR5907.

Circumstantial evidence suggests that magnetism and enhanced X-ray emission are likely correlated in early B-type stars: similar fractions of them ($\sim$ 10 %) are strong and hard X-ray sources and possess strong magnetic fields. It is also known that some B-type stars have spots on their surface. Yet up to now no X-ray activity associated with spots on early-type stars was detected. In this Letter we report the detection of a magnetic field on the B2V star $\rho$ Oph A. Previously, we assessed that the X-ray activity of this star is associated with a surface spot, herewith we establish its magnetic origin. We analyzed FORS2 ESO VLT spectra of $\rho$ Oph A taken at two epochs and detected a longitudinal component of the magnetic field of order of $\sim500$ G in one of the datasets. The detection of the magnetic field only at one epoch can be explained by stellar rotation which is also invoked to explain observed periodic X-ray activity. From archival HARPS ESO VLT high resolution spectra we derived the fundamental stellar parameters of $\rho$ Oph A and further constrained its age. We conclude that $\rho$ Oph A provides strong evidence for the presence of active X-ray emitting regions on young magnetized early type stars.

SXP 1062 is a Be X-ray binary located in the Small Magellanic Cloud. It hosts a long-period X-ray pulsar and is likely associated with the supernova remnant MCSNR J0127-7332. In this work we present a multi-wavelength view on SXP 1062 in different luminosity regimes. We consider monitoring campaigns in optical (OGLE survey) and X-ray (SWIFT telescope). During these campaigns a tight coincidence of X-ray and optical outbursts is observed. We interpret this as typical Type I outbursts as often detected in Be X-ray binaries at periastron passage of the neutron star. To study different X-ray luminosity regimes in depth, during the source quiescence we observed it with XMM-Newton while Chandra observations followed an X-ray outburst. Nearly simultaneously with Chandra observations in X-rays, in optical the RSS/SALT telescope obtained spectra of SXP 1062. On the basis of our multi-wavelength campaign we propose a simple scenario where the disc of the Be star is observed face-on, while the orbit of the neutron star is inclined with respect to the disc. According to the model of quasi-spherical settling accretion our estimation of the magnetic field of the pulsar in SXP 1062 does not require an extremely strong magnetic field at the present time.

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.

The first directly detected gravitational waves (GW 150914) were emitted by two coalescing black holes (BHs) with masses of ~36Msun and ~29Msun. Several scenarios have been proposed to put this detection into an astrophysical context. The evolution of an isolated massive binary system is among commonly considered models. Various groups have performed detailed binary-evolution calculations that lead to BH merger events. However, the question remains open as to whether binary systems with the predicted properties really exist. The aim of this paper is to help observers to close this gap by providing spectral characteristics of massive binary BH progenitors during a phase where at least one of the companions is still non-degenerate. Stellar evolution models predict fundamental stellar parameters. Using these as input for our stellar atmosphere code (PoWR), we compute a set of models for selected evolutionary stages of massive merging BH progenitors at different metallicities. The synthetic spectra obtained from our atmosphere calculations reveal that progenitors of massive BH merger events start their lives as O2-3V stars that evolve to early-type blue supergiants before they undergo core-collapse during the Wolf-Rayet phase. When the primary has collapsed, the remaining system will appear as a wind-fed high-mass X-ray binary. We provide feedback parameters, broad band magnitudes, and spectral templates that should help to identify such binaries in the future. Comparisons of empirically determined mass-loss rates with those assumed by evolution calculations reveal significant differences. The consideration of the empirical mass-loss rates in evolution calculations will possibly entail a shift of the maximum in the predicted binary-BH merger rate to higher metallicities, that is, more candidates should be expected in our cosmic neighborhood than previously assumed.

AH Cep and CW Cep are both early B-type binaries with short orbital periods of 1.8~d and 2.7~d, respectively. All four components are B0.5V types. The binaries are also double-lined spectroscopic and eclipsing. Consequently, solutions for orbital and stellar parameters make the pair of binaries ideal targets for a study of the colliding winds between two B~stars. {\em Chandra} ACIS-I observations were obtained to determine X-ray luminosities. AH~Cep was detected with an unabsorbed X-ray luminosity at a 90\% confidence interval of $(9-33)\times 10^{30}$ erg s$^{-1}$, or $(0.5-1.7)\times 10^{-7} L_{\rm Bol}$, relative to the combined Bolometric luminosities of the two components. While formally consistent with expectations for embedded wind shocks, or binary wind collision, the near-twin system of CW~Cep was a surprising non-detection. For CW~Cep, an upper limit was determined with $L_X/L_{\rm Bol} < 10^{-8}$, again for the combined components. One difference between these two systems is that AH~Cep is part of a multiple system. The X-rays from AH~Cep may not arise from standard wind shocks nor wind collision, but perhaps instead from magnetism in any one of the four components of the system. The possibility could be tested by searching for cyclic X-ray variability in AH~Cep on the short orbital period of the inner B~stars.

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.

The long-term X-ray lightcurves of classical supergiant X-ray binaries and supergiant fast X-ray transients show relatively similar super-orbital modulations, which are still lacking a sound interpretation. We propose that these modulations are related to the presence of corotating interaction regions (CIRs) known to thread the winds of OB supergiants. To test this hypothesis, we couple the outcomes of 3-D hydrodynamic models for the formation of CIRs in stellar winds with a simplified recipe for the accretion onto a neutron star. The results show that the synthetic X-ray light curves are indeed modulated by the presence of the CIRs. The exact period and amplitude of these modulations depend on a number of parameters governing the hydrodynamic wind models and on the binary orbital configuration. To compare our model predictions with the observations, we apply the 3-D wind structure previously shown to well explain the appearance of discrete absorption components in the UV time series of a prototypical B0.5I-type supergiant. Using the orbital parameters of IGRJ16493-4348 which has the same B0.5I donor spectral type, the period and modulations in the simulated X-ray light-curve are similar to the observed ones, thus providing support to our scenario. We propose, that the presence of CIRs in donor star winds should be considered in future theoretical and simulation efforts of wind-fed X-ray binaries.

Massive stars are the key agents of feedback. Consequently, quantitative analysis of massive stars are required to understand how the feedback of these objects shapes/ creates the large scale structures of the ISM. The giant HII region N206 in the Large Magellanic Cloud contains an OB association that powers a X-ray superbubble, serving as an ideal laboratory in this context. We obtained optical spectra with the muti-object spectrograph FLAMES at the ESO-VLT. When possible, the optical spectroscopy was complemented by UV spectra from the HST, IUE, and FUSE archives. Detailed spectral classifications are presented for our sample Of-type stars. For the quantitative spectroscopic analysis we use the Potsdam Wolf-Rayet (PoWR) model atmosphere code. The physical parameters and nitrogen abundances of our sample stars are determined by fitting synthetic spectra to the observations. The stellar and wind parameters of nine Of-type stars are used to construct wind momentum,luminosity relationship. We find that our sample follows a relation close to the theoretical prediction, assuming clumped winds. The most massive star in the N206 association is an Of supergiant which has a very high mass-loss rate. Two objects in our sample reveal composite spectra, showing that the Of primaries have companions of late O subtype. All stars in our sample have an evolutionary age less than 4 million years, with the O2-type star being the youngest. All these stars show a systematic discrepancy between evolutionary and spectroscopic masses. All stars in our sample are nitrogen enriched. Nitrogen enrichment shows a clear correlation with increasing projected rotational velocities. The mechanical energy input from the Of stars alone is comparable to the energy stored in the N206 superbubble as measured from the observed X-ray and H alpha emission.

Supergiant fast X-ray transients (SFXTs) are characterized by a remarkable variability in the X-ray domain, widely ascribed to the accretion from a clumpy stellar wind. In this paper we performed a systematic and homogeneous analysis of sufficiently bright X-ray flares from the SFXTs observed with XMM-Newton to probe spectral variations on timescales as short as a few hundred of seconds. Our ultimate goal is to investigate if SFXT flares and outbursts are triggered by the presence of clumps and eventually reveal whether strongly or mildly dense clumps are required. For all sources, we employ a technique developed by our group, making use of an adaptive rebinned hardness ratio to optimally select the time intervals for the spectral extraction. A total of twelve observations performed in the direction of five SFXTs are reported. We show that both strongly and mildly dense clumps can trigger these events. In the former case, the local absorption column density may increase by a factor of >>3, while in the latter case, the increase is only by a factor of 2-3 (or lower). Overall, there seems to be no obvious correlation between the dynamic ranges in the X-ray fluxes and absorption column densities in SFXTs, with an indication that lower densities are recorded at the highest fluxes. This can be explained by the presence of accretion inhibition mechanism(s). We propose a classification of the flares/outbursts from these sources to drive future observational investigations. We suggest that the difference between the classes of flares/outbursts is related to the fact that the mechanism(s) inhibiting accretion can be overcome more easily in some sources compared to others. We also investigate the possibility that different stellar wind structures, rather than clumps, could provide the means to temporarily overcome the inhibition of accretion in SFXTs.

Context: HD 54879 (O9.7 V) is one of a dozen O-stars for which an organized atmospheric magnetic field has been detected. To gain insights into the interplay between atmospheres, winds, and magnetic fields of massive stars, we acquired UV and X-ray data of HD 54879 using the Hubble Space Telescope and the XMM-Newton satellite. In addition, 35 optical amateur spectra were secured to study the variability of HD 54879. A multiwavelength (X-ray to optical) spectral analysis is performed using the Potsdam Wolf-Rayet (PoWR) model atmosphere code and the xspec software. Results: The photospheric parameters are typical for an O9.7 V star. The microturbulent, macroturbulent, and projected rotational velocities are lower than previously suggested (<4 km/s). An initial mass of 16$\,M_\odot$ and an age of 5 Myr are inferred from evolutionary tracks. We derive a mean X-ray emitting temperature of $\log T_{\rm X} = 6.7\,$[K] and an X-ray luminosity of $\log L_\text{X} = 32\,$[erg/s]. Short- and long-scale variability is seen in the H-alpha line, but only a very long period of $P \approx 5\,$yr could be estimated. Assessing the circumstellar density of HD 54879 using UV spectra, we can roughly estimate the mass-loss rate HD 54879 would have in the absence of a magnetic field as $\log \dot{M}_{B=0}\approx -9.0\,[{M_\odot}/{\rm yr}]$. The magnetic field traps the stellar wind up to the Alfv\'en radius > $12\,R_\odot$, implying that its true mass-loss rate is $\log \dot{M}< -10.2\,[{M_\odot}/{\rm yr}]$. Hence, density enhancements around magnetic stars can be exploited to estimate mass-loss rates of non-magnetic stars of similar spectral types, essential for resolving the weak wind problem. Conclusions: Our study confirms that strongly magnetized stars lose little or no mass, and supplies important constraints on the weak-wind problem of massive main sequence stars.

The blue hypergiant Cyg OB2-12 (B3Ia+) is a representative member of the class of very massive stars in a poorly understood evolutionary stage. We obtained its high-resolution X-ray spectrum using Chandra observatory. PoWR model atmospheres were calculated to provide realistic wind opacities and to establish the wind density structure. We find that collisional de-excitation is the dominant mechanism de-populating the metastable upper levels of the forbidden lines of the He-like ions SiXIV and MgXII. Comparison between the model and observations reveals that X-ray emission is produced in a dense plasma, which could reside only at the photosphere or in a colliding wind zone between binary components. The observed X-ray spectra are well fitted by thermal plasma models, with average temperatures in excess of 10 MK. The wind speed in Cyg OB2-12 is not high enough to power such high temperatures, but the collision of two winds in a binary system can be sufficient. We used archival data to investigate the X-ray properties of other blue hypergiants. In general, stars of this class are not detected as X-rays sources. We suggest that our new Chandra observations of Cyg OB2-12 can be best explained if Cyg OB2-12 is a colliding wind binary possessing a late O-type companion. This makes Cyg OB2-12 only the second binary system among the 16 known Galactic hypergiants. This low binary fraction indicates that the blue hypergiants are likely products of massive binary evolution during which they either accreted a significant amount of mass or already merged with their companion.

4U0114$+$65 is one of the slowest known X-ray pulsars. We present an analysis of a pointed observation by the XMM-Newton X-ray telescope. The energy-resolved light curve, the time-resolved X-ray spectra provided by the EPIC cameras on board and the first high-resolution spectrum of this source, provided by the Reflection Grating Spectrometer, were analysed. An X-ray pulse of $9350\pm 160$ s was measured, confirming the secular spin up of this source. We also fit the pulse-phase-resolved spectra with Comptonisation models. These models imply a very small ($r\sim 3$ km) and hot ($kT\sim 2-3$ keV) emitting region and therefore point to a hot spot over the neutron star (NS) surface as the most reliable explanation for the X-ray pulse. The long NS spin period, the spin-up rate, and persistent X-ray emission can be explained within the theory of quasi-spherical settling accretion, which may indicate that the magnetic field is in the magnetar range. Thus, 4U 0114$+$65 could be a wind-accreting magnetar. We also observed two episodes of low luminosity. The first can be explained as an absorption by an over-dense structure in the wind of the donor. The second episode may be due to temporal cessation of accretion onto one magnetic pole caused by non-spherical matter capture from the structured stellar wind. The light curve displays two types of dips seen during the high-flux intervals. The short dips, with durations of tens of seconds, are produced through absorption by wind clumps. The long dips might be associated with the rarefied interclump medium. We also found evidence of emission lines in the X-ray photoionised wind of the B1Ia donor. The Fe K$\alpha$ line is highly variable and much weaker than in other X-ray binaries with supergiant donors. The degree of wind clumping, measured through the covering fraction, is much lower than in supergiant donor stars with earlier spectral types.

The Advanced Telescope for High ENergy Astrophysics (Athena) is the X-ray observatory mission selected by ESA within its Cosmic Vision 2015-2025 programme to address the Hot and Energetic Universe scientific theme. The ESO-Athena Synergy Team (EAST) has been tasked to single out the potential scientific synergies between Athena and optical/near-infrared (NIR) and sub/mm ground based facilities, in particular those of ESO (i.e., the VLT and ELT, ALMA and APEX), by producing a White Paper to identify and develop the: 1. needs to access ESO ground-based facilities to achieve the formulated Athena science objectives; 2. needs to access Athena to achieve the formulated science objectives of ESO facilities contemporary to Athena; 3. science areas where the synergetic use of Athena and ESO facilities in the late 2020s will result in scientific added value. Community input to the process happened primarily via a dedicated ESO - Athena Synergy Workshop that took place on Sept. 14 - 16, 2016 at ESO, Garching. This White Paper presents the results of the EAST's work, sorted by synergy area, and deals with the following topics: 1. the Hot Universe: Early groups and clusters and their evolution, Physics of the Intracluster medium, Missing baryons in cosmic filaments; 2. the Energetic Universe: Supermassive black hole (SMBH) history, SMBH accretion disks, Active Galactic Nuclei feedback - Molecular outflows, Ultra-fast outflows, Accretion Physics, Transient Science; 3. Observatory Science: Star Formation, Stars. It then discusses the optical-NIR-sub-mm perspective by providing details on VLT/MOONS, the E-ELT instruments, in particular the MOS, VISTA/4MOST, the ESO and ALMA archives, future ALMA and ESO developments, and finally the (likely) ESO - Athena astronomical scene in the 2020s. (abridged)

We present the results of a 140 ks XMM-Newton observation of the B2 star $\rho$ Ophiuchi A. The star has exhibited strong X-ray variability: a cusp-shaped increase of rate, similar to that which we partially observed in 2013, and a bright flare. These events are separated in time by about 104 ks, which likely corresponds to the rotational period of the star (1.2 days). Time resolved spectroscopy of the X-ray spectra shows that the first event is caused by an increase of the plasma emission measure, while the second increase of rate is a major flare with temperatures in excess of 60 MK ($kT\sim5$ keV). From the analysis of its rise, we infer a magnetic field of $\ge300$ G and a size of the flaring region of $\sim1.4-1.9\times10^{11}$ cm, which corresponds to $\sim25\%-30\%$ of the stellar radius. We speculate that either an intrinsic magnetism that produces a hot spot on its surface or an unknown low mass companion are the source of such X-rays and variability. A hot spot of magnetic origin should be a stable structure over a time span of $\ge$2.5 years, and suggests an overall large scale dipolar magnetic field that produces an extended feature on the stellar surface. In the second scenario, a low mass unknown companion is the emitter of X-rays and it should orbit extremely close to the surface of the primary in a locked spin-orbit configuration, almost on the verge of collapsing onto the primary. As such, the X-ray activity of the secondary star would be enhanced by its young age, and the tight orbit as in RS Cvn systems and $\rho$ Ophiuchi would constitute an extreme system that is worthy of further investigation.

Theoretical models predict that the compressed interstellar medium around runaway O stars can produce high-energy non-thermal diffuse emission, in particular, non-thermal X-ray and $\gamma$-ray emission. So far, detection of non-thermal X-ray emission was claimed for only one runaway star AE Aur. We present a search for non-thermal diffuse X-ray emission from bow shocks using archived XMM-Newton observations for a clean sample of 6 well-determined runaway O stars. We find that none of these objects present diffuse X-ray emission associated to their bow shocks, similarly to previous X-ray studies toward $\zeta$ Oph and BD$+$43$^{\circ}$3654. We carefully investigated multi-wavelength observations of AE Aur and could not confirm previous findings of non-thermal X-rays. We conclude that so far there is no clear evidence of non-thermal extended emission in bow shocks around runaway O 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.

Observations with powerful X-ray telescopes, such as XMM-Newton and Chandra, significantly advance our understanding of massive stars. Nearly all early-type stars are X-ray sources. Studies of their X-ray emission provide important diagnostics of stellar winds. High-resolution X-ray spectra of O-type stars are well explained when stellar wind clumping is taking into account, providing further support to a modern picture of stellar winds as non-stationary, inhomogeneous outflows. X-ray variability is detected from such winds, on time scales likely associated with stellar rotation. High-resolution X-ray spectroscopy indicates that the winds of late O-type stars are predominantly in a hot phase. Consequently, X-rays provide the best observational window to study these winds. X-ray spectroscopy of evolved, Wolf-Rayet type, stars allows to probe their powerful metal enhanced winds, while the mechanisms responsible for the X-ray emission of these stars are not yet understood.

The recent discovery of a gravitational wave from the merging of two black holes of about 30 solar masses each challenges our incomplete understanding of massive stars and their evolution. Critical ingredients comprise mass-loss, rotation, magnetic fields, internal mixing, and mass transfer in close binary systems. The imperfect knowledge of these factors implies large uncertainties for models of stellar populations and their feedback. In this contribution we summarize our empirical studies of Wolf-Rayet populations at different metallicities by means of modern non-LTE stellar atmosphere models, and confront these results with the predictions of stellar evolution models. At the metallicity of our Galaxy, stellar winds are probably too strong to leave remnant masses as high as 30 solar masses, but given the still poor agreement between evolutionary tracks and observation even this conclusion is debatable. At the low metallicity of the Small Magellanic Cloud, all WN stars which are (at least now) single are consistent with evolving quasi-homogeneously. O and B-type stars, in contrast, seem to comply with standard evolutionary models without strong internal mixing. Close binaries which avoided early merging could evolve quasi-homogeneously and lead to close compact remnants of relatively high masses that merge within a Hubble time.

In this paper we investigate the multiwavelengths properties of the magnetic early B-type star HR7355. We present its radio light curves at several frequencies, taken with the Jansky Very Large Array, and X-ray spectra, taken with the XMM X-ray telescope. Modeling of the radio light curves for the Stokes I and V provides a quantitative analysis of the HR7355 magnetosphere. A comparison between HR7355 and a similar analysis for the Ap star CUVir, allows us to study how the different physical parameters of the two stars affect the structure of the respective magnetospheres where the non-thermal electrons originate. Our model includes a cold thermal plasma component that accumulates at high magnetic latitudes that influences the radio regime, but does not give rise to X-ray emission. Instead, the thermal X-ray emission arises from shocks generated by wind stream collisions close to the magnetic equatorial plane. The analysis of the X-ray spectrum of HR7355 also suggests the presence of a non-thermal radiation. Comparison between the spectral index of the power-law X-ray energy distribution with the non-thermal electron energy distribution indicates that the non-thermal X-ray component could be the auroral signature of the non-thermal electrons that impact the stellar surface, the same non-thermal electrons that are responsible for the observed radio emission. On the basis of our analysis, we suggest a novel model that simultaneously explains the X-ray and the radio features of HR7355 and is likely relevant for magnetospheres of other magnetic early type stars.

Features and applications of quasi-spherical settling accretion onto rotating magnetized neutron stars in high-mass X-ray binaries are discussed. The settling accretion occurs in wind-fed HMXBs when the plasma cooling time is longer than the free-fall time from the gravitational capture radius, which can take place in low-luminosity HMXBs with $L_x\lesssim 4\times 10^{36}$ erg/s. We briefly review the implications of the settling accretion, focusing on the SFXT phenomenon, which can be related to instability of the quasi-spherical convective shell above the neutron star magnetosphere due to magnetic reconnection from fast temporarily magnetized winds from OB-supergiant. If a young neutron star in a wind-fed HMXB is rapidly rotating, the propeller regime in a quasi-spherical hot shell occurs. We show that X-ray spectral and temporal properties of enigmatic $\gamma$ Cas Be-stars are consistent with failed settling accretion regime onto a propelling neutron star. The subsequent evolutionary stage of $\gamma$ Cas and its analogs should be the X Per-type binaries comprising low-luminosity slowly rotating X-ray pulsars.

The "B fields in OB stars" (BOB) collaboration is based on an ESO Large Programme, to study the occurrence rate, properties, and ultimately the origin of magnetic fields in massive stars. In the framework of this programme, we carried out low-resolution spectropolarimetric observations of a large sample of massive stars using FORS2 installed at the ESO VLT 8-m telescope. We determined the magnetic field values with two completely independent reduction and analysis pipelines. Our in-depth study of the magnetic field measurements shows that differences between our two pipelines are usually well within 3sigma errors. From the 32 observations of 28 OB stars, we were able to monitor the magnetic fields in CPD-57 3509 and HD164492C, confirm the magnetic field in HD54879, and detect a magnetic field in CPD-62 2124. We obtain a magnetic field detection rate of 6+-3% for the full sample of 69 OB stars observed with FORS2 within the BOB programme. For the pre-selected objects with a v sin i below 60 km/s, we obtain a magnetic field detection rate of 5+-5%. We also discuss X-ray properties and multiplicity of the objects in our FORS2 sample with respect to the magnetic field detections.

The enigmatic X-ray emission from the bright optical star, $\gamma$ Cassiopeia, is a long-standing problem. $\gamma$ Cas is known to be a binary system consisting of a Be-type star and a low-mass ($M\sim 1\,M_\odot$) companion of unknown nature orbiting in the Be-disk plane. Here we apply the quasi-spherical accretion theory onto a compact magnetized star and show that if the low-mass companion of $\gamma$ Cas is a fast spinning neutron star, the key observational signatures of $\gamma$ Cas are remarkably well reproduced. Direct accretion onto this fast rotating neutron star is impeded by the propeller mechanism. In this case, around the neutron star magnetosphere a hot shell is formed that emits thermal X-rays in qualitative and quantitative agreement with observed properties of the X-ray emission from $\gamma$ Cas. We suggest that $\gamma$ Cas and its analogs constitute a new subclass of Be-type X-ray binaries hosting rapidly rotating neutron stars formed in supernova explosions with small kicks. The subsequent evolutionary stage of $\gamma$ Cas and its analogs should be the X Per-type binaries comprising low-luminosity slowly rotating X-ray pulsars. The model explains the enigmatic X-ray emission from $\gamma$ Cas, and also establishes evolutionary connections between various types of rotating magnetized neutron stars in Be-binaries.

We present the first SB2 orbital solution and disentanglement of the massive Wolf-Rayet binary R145 (P = 159d) located in the Large Magellanic Cloud. The primary was claimed to have a stellar mass greater than 300Msun, making it a candidate for the most massive star known. While the primary is a known late type, H-rich Wolf-Rayet star (WN6h), the secondary could not be so far unambiguously detected. Using moderate resolution spectra, we are able to derive accurate radial velocities for both components. By performing simultaneous orbital and polarimetric analyses, we derive the complete set of orbital parameters, including the inclination. The spectra are disentangled and spectroscopically analyzed, and an analysis of the wind-wind collision zone is conducted. The disentangled spectra and our models are consistent with a WN6h type for the primary, and suggest that the secondary is an O3.5 If*/WN7 type star. We derive a high eccentricity of e = 0.78 and minimum masses of M1 sin^3 i ~ M2 sin^3 i ~ 13 +- 2 Msun, with q = M2 / M1 = 1.01 +- 0.07. An analysis of emission excess stemming from a wind-wind collision yields a similar inclination to that obtained from polarimetry (i = 39 +- 6deg). Our analysis thus implies M1 = 53^{+40}_{-20} and M2 = 54^{+40}_{-20} Msun, excluding M1 > 300Msun. A detailed comparison with evolution tracks calculated for single and binary stars, as well as the high eccentricity, suggest that the components of the system underwent quasi-homogeneous evolution and avoided mass-transfer. This scenario would suggest current masses of ~ 80 Msun and initial masses of Mi,1 ~ 105 and Mi,2 ~ 90Msun, consistent with the upper limits of our derived orbital masses, and would imply an age of ~2.2 Myr.

In this paper we report on a long multi-wavelength observational campaign of the supergiant fast X-ray transient prototype IGR J17544-2619. A 150 ks-long observation was carried out simultaneously with XMM-Newton and NuSTAR, catching the source in an initial faint X-ray state and then undergoing a bright X-ray outburst lasting about 7 ks. We studied the spectral variability during outburst and quiescence by using a thermal and bulk Comptonization model that is typically adopted to describe the X-ray spectral energy distribution of young pulsars in high mass X-ray binaries. Although the statistics of the collected X-ray data were relatively high we could neither confirm the presence of a cyclotron line in the broad-band spectrum of the source (0.5-40 keV), nor detect any of the previously reported tentative detection of the source spin period. The monitoring carried out with Swift/XRT during the same orbit of the system observed by XMM-Newton and NuSTAR revealed that the source remained in a low emission state for most of the time, in agreement with the known property of all supergiant fast X-ray transients being significantly sub-luminous compared to other supergiant X-ray binaries. Optical and infrared observations were carried out for a total of a few thousands of seconds during the quiescence state of the source detected by XMM-Newton and NuSTAR. The measured optical and infrared magnitudes were slightly lower than previous values reported in the literature, but compatible with the known micro-variability of supergiant stars. UV observations obtained with the UVOT telescope on-board Swift did not reveal significant changes in the magnitude of the source in this energy domain compared to previously reported values.

Gamma Cassiopeiae is an enigmatic Be star with unusually strong hard X-ray emission. The Suzaku observatory detected six rapid X-ray spectral hardening events called "softness dips" in a ~100 ksec duration observation in 2011. All the softness dip events show symmetric softness ratio variations, and some of them have flat bottoms apparently due to saturation. The softness dip spectra are best described by either ~40% or ~70% partial covering absorption to kT ~12 keV plasma emission by matter with a neutral hydrogen column density of ~2-8e21 cm-2, while the spectrum outside of these dips is almost free of absorption. This result suggests the presence of two distinct X-ray emitting spots in the gamma Cas system, perhaps on a white dwarf companion with dipole mass accretion. The partial covering absorbers may be blobs in the Be stellar wind, the Be disk, or rotating around the white dwarf companion. Weak correlations of the softness ratios to the hard X-ray flux suggest the presence of stable plasmas at kT ~0.9 and 5 keV, which may originate from the Be or white dwarf winds. The formation of a Be star and white dwarf binary system requires mass transfer between two stars; gamma Cas may have experienced such activity in the past.

This review briefly summarizes our knowledge of the X-ray emission from single WN, WC, and WO stars. These stars have relatively modest X-ray luminosities, typically not exceeding L_sun. The analysis of X-ray spectra usually reveals thermal plasma with temperatures reaching a few 10 MK. X-ray variability is detected in some WN stars. At present we don't fully understand how X-ray radiation in produced in WR stars, albeit there are some promising research avenues, such as the presence of CIRs in the winds of some stars. To fully understand WR stars we need to unravel mechanisms of X-ray production in their winds.

Nearly all types of massive stars with radiatively driven stellar winds are X-ray sources that can be observed by the presently operating powerful X-ray telescopes. In this review I briefly address recent advances in our understanding of stellar winds obtained from X-ray observations. The winds of OB dwarfs with subtypes later than O9V may be predominantly in a hot phase, and X-ray observations offer the best window for their studies. The X-ray properties of OB supergiants are largely determined by the effects of radiative transfer in their clumped stellar winds. The recently suggested method to directly measure mass-loss rates of O stars by fitting the shapes of X-ray emission lines is considered but its validity cannot be confirmed. To obtain robust quantitative information on stellar wind parameters from X-ray spectroscopy, a multiwavelength analysis by means of stellar atmosphere models is required. Independent groups are now performing such analyses with encouraging results. Joint analyses of optical, UV, and X-ray spectra of OB supergiants yield consistent mass-loss rates. Depending on the adopted clumping parameters, the empirically derived mass-loss rates are comparable (within a factor of a few) to those predicted by standard recipes (Vink et al. 2001). All sufficiently studied O stars display variable X-ray emission that might be related to corotating interaction regions in their winds. In the latest stages of stellar evolution, single red supergiants and luminous blue variable stars do not emit observable amounts of X-rays. On the other hand, nearly all types of Wolf-Rayet (WR) stars are X-ray sources. X-ray spectroscopy allows a sensitive probe of WR wind abundances and opacities.

Massive stars are extremely luminous and drive strong winds, blowing a large part of their matter into the galactic environment before they finally explode as a supernova. Quantitative knowledge of massive star feedback is required to understand our Universe as we see it. Traditionally, massive stars have been studied under the assumption that their winds are homogeneous and stationary, largely relying on the Sobolev approximation. However, observations with the newest instruments, together with progress in model calculations, ultimately dictate a cardinal change of this paradigm: stellar winds are highly inhomogeneous. Hence, we are now advancing to a new stage in our understanding of stellar winds. Using the foundations laid by V.V. Sobolev and his school, we now update and further develop the stellar spectral analysis techniques. New sophisticated 3-D models of radiation transfer in inhomogeneous expanding media elucidate the physics of stellar winds and improve classical empiric mass-loss rate diagnostics. Applications of these new techniques to multiwavelength observations of massive stars yield consistent and robust stellar wind parameters.

Massive WR stars are evolved massive stars characterized by strong mass-loss. Hypothetically, they can form either as single stars or as mass donors in close binaries. About 40% of the known WR stars are confirmed binaries, raising the question as to the impact of binarity on the WR population. By performing a spectral analysis of all multiple WR systems in the SMC, we obtain the full set of stellar parameters for each individual component. Mass-luminosity relations are tested, and the importance of the binary evolution channel is assessed. The spectral analysis is performed with the PoWR model atmosphere code by superimposing model spectra that correspond to each component. Evolutionary channels are constrained using the BPASS evolution tool. Significant Hydrogen mass fractions (0.1 - 0.4) are detected in all WN components. A comparison with mass-luminosity relations and evolutionary tracks implies that the majority of the WR stars in our sample are not chemically homogeneous. The WR component in the binary AB 6 is found to be very luminous (Log L ~ 6.3 [Lsun]) given its orbital mass (~10 Msun), presumably because of observational contamination by a third component. Evolutionary paths derived for our objects suggest that Roche lobe overflow had occurred in most systems, affecting their evolution. However, the implied initial masses are large enough for the primaries to have entered the WR phase, regardless of binary interaction. Together with the results for the putatively single SMC WR stars, our study suggests that the binary evolution channel does not dominate the formation of WR stars at SMC metallicity.

The origin of the magnetic field in massive O-type stars is still under debate. To model the physical processes responsible for the generation of O star magnetic fields, it is important to understand whether correlations between the presence of a magnetic field and stellar evolutionary state, rotation velocity, kinematical status, and surface composition can be identified. The O4Ief supergiant zeta Pup is a fast rotator and a runaway star, which may be a product of a past binary interaction, possibly having had an encounter with the cluster Trumper 10 some 2Myr ago. The currently available observational material suggests that certain observed phenomena in this star may be related to the presence of a magnetic field. We acquired spectropolarimetric observations of zeta Pup with FORS2 mounted on the 8-m Antu telescope of the VLT to investigate if a magnetic field is indeed present in this star. We show that many spectral lines are highly variable and probably vary with the recently detected period of 1.78d. No magnetic field is detected in zeta Pup, as no magnetic field measurement has a significance level higher than 2.4sigma. Still, we studied the probability of a single sinusoidal explaining the variation of the longitudinal magnetic field measurements.

The accretion of the stellar wind material by a compact object represents the main mechanism powering the X-ray emission in classical supergiant high mass X-ray binaries and supergiant fast X-ray transients. In this work we present the first attempt to simulate the accretion process of a fast and dense massive star wind onto a neutron star, taking into account the effects of the centrifugal and magnetic inhibition of accretion ("gating") due to the spin and magnetic field of the compact object. We made use of a radiative hydrodynamical code to model the non-stationary radiatively driven wind of an O-B supergiant star and then place a neutron star characterized by a fixed magnetic field and spin period at a certain distance from the massive companion. Our calculations follow, as a function of time (on a total time scale of several hours), the transition of the system through all different accretion regimes that are triggered by the intrinsic variations in the density and velocity of the non-stationary wind. The X-ray luminosity released by the system is computed at each time step by taking into account the relevant physical processes occurring in the different accretion regimes. Synthetic lightcurves are derived and qualitatively compared with those observed from classical supergiant high mass X-ray binaries and supergiant fast X-ray transients. Although a number of simplifications are assumed in these calculations, we show that taking into account the effects of the centrifugal and magnetic inhibition of accretion significantly reduces the average X-ray luminosity expected for any neutron star wind-fed binary. The present model calculations suggest that long spin periods and stronger magnetic fields are favoured in order to reproduce the peculiar behavior of supergiant fast X-ray transients in the X-ray domain.

To investigate if magnetic fields are present in Wolf-Rayet stars, we selected a few stars in the Galaxy and one in the Large Magellanic Cloud (LMC). We acquired low-resolution spectropolarimetric observations with the ESO FORS2 instrument during two different observing runs. During the first run in visitor mode, we observed the LMC Wolf-Rayet star BAT99 7 and the stars WR6, WR7, WR18, and WR23 in our Galaxy. The second run in service mode was focused on monitoring the star WR6. Linear polarization was recorded immediately after the observations of circular polarization. During our visitor observing run, the magnetic field for the cyclically variable star WR6 was measured at a significance level of 3.3sigma (<B_z> = 258+-78G). Among the other targets, the highest value for the longitudinal magnetic field, <B_z> = 327+-141G, was measured in the LMC star BAT99 7. Spectropolarimetric monitoring of the star WR6 revealed a sinusoidal nature of the <B_z> variations with the known rotation period of 3.77d, significantly adding to the confidence in the detection. The presence of the rotation-modulated magnetic variability is also indicated in our frequency periodogram. The reported field magnitude suffers from significant systematic uncertainties at the factor 2 level, in addition to the quoted statistical uncertainties, owing to the theoretical approach used to characterize it. Linear polarization measurements showed no line effect in the stars, apart from WR6. BAT99 7, WR7, and WR,23 do not show variability of the linear polarization over two nights.

Classical Supergiant X-ray Binaries (SGXBs) and Supergiant Fast X-ray Transients (SFXTs) are two types of High-mass X-ray Binaries (HMXBs) that present similar donors but, at the same time, show very different behavior in the X-rays. The reason for this dichotomy of wind-fed HMXBs is still a matter of debate. Among the several explanations that have been proposed, some of them invoke specific stellar wind properties of the donor stars. Only dedicated empiric analysis of the donors' stellar wind can provide the required information to accomplish an adequate test of these theories. However, such analyses are scarce. To close this gap, we perform a comparative analysis of the optical companion in two important systems: IGR J17544-2619 (SFXT) and Vela X-1 (SGXB). We analyse the spectra of each star in detail and derive their stellar and wind properties. We compare the wind parameters, giving us an excellent chance of recognizing key differences between donor winds in SFXTs and SGXBs. We find that the stellar parameters derived from the analysis generally agree well with the spectral types of the two donors: O9I (IGR J17544-2619) and B0.5Iae (Vela X-1). An important difference between the stellar winds of the two stars is their terminal velocities v_inf=1500km/s in IGR J17544-2619 and v_inf=700km/s in Vela~X-1, which has important consequences on the X-ray luminosity of these sources. Their specific combination of wind speed and pulsar spin favours an accretion regime with a persistently high luminosity in Vela X-1, while it favours an inhibiting accretion mechanism in IGR~J17544-2619. Our study demonstrates that the wind relative velocity is critical in the determination of the class of HMXBs hosting a supergiant donor, given that it may shift the accretion mechanism from direct accretion to propeller regimes when combined with other parameters.

Substantial progress has been achieved over the last decade in studies of stellar magnetism due to the improvement of magnetic field measurement methods. We review recent results on the magnetic field characteristics of early B- and O-type stars obtained by various teams using different measurement techniques.

Non-thermal radiation has been predicted within bow shocks around runaway stars by recent theoretical works. We present X-ray observations towards the runaway stars $\zeta$ Oph (Chandra and Suzaku) and BD+43$^{\circ}$3654 (XMM-Newton) to search for the presence of non-thermal X-ray emission. We found no evidence of non-thermal emission spatially coincident with the bow shocks, nonetheless, diffuse emission is detected in the vicinity of $\zeta$ Oph. After a careful analysis of its spectral characteristics we conclude that this emission has a thermal nature with a plasma temperature of $T \approx 2 \times10^{6}$ K. The cometary shape of this emission seems to be in line with recent predictions of radiation-hydrodynamic models of runaway stars. The case of BD+43$^{\circ}$3654 is puzzling as non-thermal emission has been reported in a previous work for this source.

The massive Wolf-Rayet type star WR 102c is located near the Quintuplet Cluster - one of the three massive star clusters in the Galactic Center region. Previous studies indicated that WR 102c may have a dusty circumstellar nebula and is among the main ionizing sources of the Sickle Nebula associated with the Quintuplet cluster. We obtained observations with the ESO's VLT integral field spectrograph SINFONI in the K-band, extracted the stellar spectra, and analyzed them by means of stellar atmosphere models. Our new analysis supersedes the results reported for WR 102c previously. We significantly revise down its bolometric luminosity and hydrogen content. We detect four early OB type stars close to WR 102c. These stars have radial velocities similar to that of WR 102c. We suggest that together with WR 102c these stars belong to a distinct star cluster with a total mass of about 1000 M_sun. We identify a new WR nebula around WR 102c in the SINFONI map of the diffuse Br$\gamma$ emission and in the HST's Pa$\alpha$ images. The Br$\gamma$ line at different locations is not significantly broadened and similar to the width of nebular emission elsewhere in the H II region around WR 102c. The massive star WR 102c located in the Galactic Center region resides in a star cluster containing further early type stars. The stellar parameters of WR 102c are typical for hydrogen-free WN6 stars. We newly identify a nebula surrounding WR 102c that has a morphology similar to other nebulae around hydrogen-free WR stars, and propose that the formation of this nebula is linked to interaction of the fast stellar wind with the matter ejected at previous evolutionary stage of WR 102c.

We report the detection of a magnetic field in the helium-strong star CPD-57 3509 (B2 IV), a member of the Galactic open cluster NGC3293, and characterise the star's atmospheric and fundamental parameters. Spectropolarimetric observations with FORS2 and HARPSpol are analysed using two independent approaches to quantify the magnetic field strength. A high-S/N FLAMES/GIRAFFE spectrum is analysed using a hybrid non-LTE model atmosphere technique. Comparison with stellar evolution models constrains the fundamental parameters of the star. We obtain a firm detection of a surface averaged longitudinal magnetic field with a maximum amplitude of about 1 kG. Assuming a dipolar configuration of the magnetic field, this implies a dipolar field strength larger than 3.3 kG. Moreover, the large amplitude and fast variation (within about 1 day) of the longitudinal magnetic field implies that CPD-57 3509 is spinning very fast despite its apparently slow projected rotational velocity. The star should be able to support a centrifugal magnetosphere, yet the spectrum shows no sign of magnetically confined material; in particular, emission in H{\alpha} is not observed. Apparently, the wind is either not strong enough for enough material to accumulate in the magnetosphere to become observable or, alternatively, some leakage process leads to loss of material from the magnetosphere. The quantitative spectroscopic analysis of the star yields an effective temperature and a logarithmic surface gravity of 23750+-250 K and 4.05+-0.10, respectively, and a surface helium fraction of 0.28+-0.02 by number. The surface abundances of C, N, O, Ne, S, and Ar are compatible with the cosmic abundance standard, whereas Mg, Al, Si, and Fe are depleted by about a factor of 2. This abundance pattern can be understood as the consequence of a fractionated stellar wind. CPD-57 3509 is one of the most evolved He-strong stars known.

With a deep Chandra/HETGS exposure of WR 6, we have resolved emission lines whose profiles show that the X-rays originate from a uniformly expanding spherical wind of high X-ray-continuum optical depth. The presence of strong helium-like forbidden lines places the source of X-ray emission at tens to hundreds of stellar radii from the photosphere. Variability was present in X-rays and simultaneous optical photometry, but neither were correlated with the known period of the system or with each other. An enhanced abundance of sodium revealed nuclear processed material, a quantity related to the evolutionary state of the star. The characterization of the extent and nature of the hot plasma in WR 6 will help to pave the way to a more fundamental theoretical understanding of the winds and evolution of massive stars.

We present an overview of four phase-constrained Chandra HETGS X-ray observations of Delta Ori A. Delta Ori A is actually a triple system which includes the nearest massive eclipsing spectroscopic binary, Delta Ori Aa, the only such object which can be observed with little phase-smearing with the Chandra gratings. Since the fainter star, Delta Ori Aa2, has a much lower X-ray luminosity than the brighter primary, Delta Ori A provides a unique system with which to test the spatial distribution of the X-ray emitting gas around Delta Ori Aa1 via occultation by the photosphere of and wind cavity around the X-ray dark secondary. Here we discuss the X-ray spectrum and X-ray line profiles for the combined observation, having an exposure time of nearly 500 ksec and covering nearly the entire binary orbit. Companion papers discuss the X-ray variability seen in the Chandra spectra, present new space-based photometry and ground-based radial velocities simultaneous with the X-ray data to better constrain the system parameters, and model the effects of X-rays on the optical and UV spectrum. We find that the X-ray emission is dominated by embedded wind shock emission from star Aa1, with little contribution from the tertiary star Ab or the shocked gas produced by the collision of the wind of Aa1 against the surface of Aa2. We find a similar temperature distribution to previous X-ray spectrum analyses. We also show that the line half-widths are about $0.3-0.5\times$ the terminal velocity of the wind of star Aa1. We find a strong anti-correlation between line widths and the line excitation energy, which suggests that longer-wavelength, lower-temperature lines form farther out in the wind. Our analysis also indicates that the ratio of the intensities of the strong and weak lines of \ion{Fe}{17} and \ion{Ne}{10} are inconsistent with model predictions, which may be an effect of resonance scattering

We present time-resolved and phase-resolved variability studies of an extensive X-ray high-resolution spectral dataset of the $\delta$ Orionis Aa binary system. The four observations, obtained with Chandra ACIS HETGS, have a total exposure time of ~479 ks and provide nearly complete binary phase coverage. Variability of the total X-ray flux in the range 5-25 $\AA$ is confirmed, with maximum amplitude of about +/-15% within a single ~125 ks observation. Periods of 4.76d and 2.04d are found in the total X-ray flux, as well as an apparent overall increase in flux level throughout the 9-day observational campaign. Using 40 ks contiguous spectra derived from the original observations, we investigate variability of emission line parameters and ratios. Several emission lines are shown to be variable, including S XV, Si XIII, and Ne IX. For the first time, variations of the X-ray emission line widths as a function of the binary phase are found in a binary system, with the smallest widths at phase=0.0 when the secondary $\delta$ Orionis Aa2 is at inferior conjunction. Using 3D hydrodynamic modeling of the interacting winds, we relate the emission line width variability to the presence of a wind cavity created by a wind-wind collision, which is effectively void of embedded wind shocks and is carved out of the X-ray-producing primary wind, thus producing phase-locked X-ray variability.

The number of magnetic stars detected among massive stars is small; nevertheless, the role played by the magnetic field in stellar evolution cannot be disregarded. Links between line profile variability, enhancements/depletions of surface chemical abundances, and magnetic fields have been identified for low-mass B-stars, but for the O-type domain this is almost unexplored. Based on FORS2 and HARPS spectropolarimetric data, we present the first detection of a magnetic field in HD54879, a single slowly rotating O9.7 V star. Using two independent and different techniques we obtained the firm detection of a surface average longitudinal magnetic field with a maximum amplitude of about 600 G, in modulus. A quantitative spectroscopic analysis of the star with the stellar atmosphere code FASTWIND results in an effective temperature and a surface gravity of 33000$\pm1000$ K and 4.0$\pm0.1$ dex. The abundances of carbon, nitrogen, oxygen, silicon, and magnesium are found to be slightly lower than solar, but compatible within the errors. We investigate line-profile variability in HD54879 by complementing our spectra with spectroscopic data from other recent OB-star surveys. The photospheric lines remain constant in shape between 2009 and 2014, although H$\alpha$ shows a variable emission. The H$\alpha$ emission is too strong for a standard O9.7 V and is probably linked to the magnetic field and the presence of circumstellar material. Its normal chemical composition and the absence of photospheric line profile variations make HD54879 the most strongly magnetic, non-variable single O-star detected to date.

Massive B-type stars with strong magnetic fields and fast rotation are very rare and provide a mystery for theories of both star formation and magnetic field evolution. Only two such stars, called sigma Ori E analogs, were previously known. Recently, a team involved in APOGEE, one of the Sloan Digital Sky Survey III programs, announced the discovery of two additional rigidly rotating magnetosphere stars, HD23478 and HD345439. The presence of magnetic fields in these newly discovered sigma Ori E analogs was not investigated in the past. In the framework of our ESO Large Programme, and one normal ESO programme, we carried out low-resolution FORS2 spectropolarimetric observations of HD23478 and HD345439. From the measurements using hydrogen lines, we discover a rather strong longitudinal magnetic field of the order of up to 1.5kG in HD23478, and up to 1.3kG using the entire spectrum. The analysis of HD345439 using four subsequent spectropolarimetric subexposures does not reveal the presence of a magnetic field at a significance level of 3sigma. On the other hand, the inspection of individual subexposures indicates that HD345439 may host a strong magnetic field, rapidly varying over 88 minutes. A hint at the fast rotation of HD345439 is also given by the behaviour of several metallic and He I lines in the low-resolution FORS2 spectra, showing profile variations already on such a short time scale.

We report on both high-precision photometry from the MOST space telescope and ground-based spectroscopy of the triple system delta Ori A consisting of a binary O9.5II+early-B (Aa1 and Aa2) with P = 5.7d, and a more distant tertiary (O9 IV P > 400 yrs). This data was collected in concert with X-ray spectroscopy from the Chandra X-ray Observatory. Thanks to continuous coverage for 3 weeks, the MOST light curve reveals clear eclipses between Aa1 and Aa2 for the first time in non-phased data. From the spectroscopy we have a well constrained radial velocity curve of Aa1. While we are unable to recover radial velocity variations of the secondary star, we are able to constrain several fundamental parameters of this system and determine an approximate mass of the primary using apsidal motion. We also detected second order modulations at 12 separate frequencies with spacings indicative of tidally influenced oscillations. These spacings have never been seen in a massive binary, making this system one of only a handful of such binaries which show evidence for tidally induced pulsations.

Beta Cephei-type variables are early B-type stars that are characterized by oscillations observable in their optical light curves. At least one Beta Cep-variable also shows periodic variability in X-rays. Here we study the X-ray light curves in a sample of beta Cep-variables to investigate how common X-ray pulsations are for this type of stars. We searched the Chandra and XMM-Newton X-ray archives and selected stars that were observed by these telescopes for at least three optical pulsational periods. We retrieved and analyzed the X-ray data for kappa Sco, beta Cru, and alpha Vir. The X-ray light curves of these objects were studied to test for their variability and periodicity. While there is a weak indication for X-ray variability in beta Cru, we find no statistically significant evidence of X-ray pulsations in any of our sample stars. This might be due either to the insufficient data quality or to the physical lack of modulations. New, more sensitive observations should settle this question.

Eclipsing systems of massive stars allow one to explore the properties of their components in great detail. We perform a multi-wavelength, non-LTE analysis of the three components of the massive multiple system $\delta$ Ori A, focusing on the fundamental stellar properties, stellar winds, and X-ray characteristics of the system. The primary's distance-independent parameters turn out to be characteristic for its spectral type (O9.5 II), but usage of the ${\rm \it Hipparcos}$ parallax yields surprisingly low values for the mass, radius, and luminosity. Consistent values follow only if $\delta$ Ori lies at about twice the ${\rm \it Hipparcos}$ distance, in the vicinity of the $\sigma$-Orionis cluster. The primary and tertiary dominate the spectrum and leave the secondary only marginally detectable. We estimate the V-band magnitude difference between primary and secondary to be $\Delta V \approx 2.\!\!^{\rm m}8$. The inferred parameters suggest the secondary is an early B-type dwarf ($\approx$ B1 V), while the tertiary is an early B-type subgiant ($\approx$ B0 IV). We find evidence for rapid turbulent velocities ($\sim 200$ km ${\rm s}^{-1}$) and wind inhomogeneities, partially optically thick, in the primary's wind. The bulk of the X-ray emission likely emerges from the primary's stellar wind ($\log L_{\text{X}} / L_{\text{Bol}} \approx -6.85$), initiating close to the stellar surface at $R_0 \sim 1.1\,R_*$. Accounting for clumping, the mass-loss rate of the primary is found to be $\log \dot{M} \approx -6.4\,[M_\odot\,{\rm yr}^{-1}]$, which agrees with hydrodynamic predictions, and provides a consistent picture along the X-ray, UV, optical and radio spectral domains.

We present a comprehensive analysis of the whole sample of available XMM-Newton observations of High Mass X-ray Binaries (HMXBs) until August, 2013, focusing on the FeK{\alpha} emission line. This line is a key tool to better understand the physical properties of the material surrounding the X-ray source within a few stellar radii (the circumstellar medium). We have collected observations from 46 HMXBs, detecting FeK{\alpha} in 21 of them. We have used the standard classification of HMXBs to divide the sample in different groups. We find that: (1) FeK{\alpha} is centred at a mean value of 6.42 keV. Considering the instrumental and fits uncertainties, this value is compatible with ionization states lower than FeXVIII. (2) The flux of the continuum is well correlated with the flux of the line, as expected. Eclipse observations show that the Fe fluorescence emission comes from an extended region surrounding the X-ray source. (3) FeK{\alpha} is narrow (width lower than 0.15keV), reflecting that the reprocessing material does not move at high speeds. We attempt to explain the broadness of the line in terms of three possible broadening phenomena: line blending, Compton scattering and Doppler shifts (with velocities of the reprocessing material V=1000-2000 km/s). (4) The equivalent hydrogen column (NH) directly correlates with the EW of FeK{\alpha}, displaying clear similarities to numerical simulations. It highlights the strong link between the absorbing and the fluorescent matter. The obtained results clearly point to a very important contribution of the donors wind in the FeK{\alpha} emission and the absorption when the donor is a supergiant massive star.

A rather strong mean longitudinal magnetic field of the order of a few hundred Gauss was detected a few years ago in the Of?p star CPD -28 2561 using FORS2 low-resolution spectropolarimetric observations. In this work we present additional low-resolution spectropolarimetric observations obtained during several weeks in 2013 December using FORS2 (FOcal Reducer low dispersion Spectrograph) mounted at the 8-m Antu telescope of the VLT. These observations cover a little less than half of the stellar rotation period of 73.41d mentioned in the literature. The behaviour of the mean longitudinal magnetic field is consistent with the assumption of a single-wave variation during the stellar rotation cycle, indicating a dominant dipolar contribution to the magnetic field topology. The estimated polar strength of the surface dipole B_d is larger than 1.15kG. Further, we compared the behaviour of the line profiles of various elements at different rotation phases associated with different magnetic field strengths. The strongest contribution of the emission component is observed at the phases when the magnetic field shows a negative or positive extremum. The comparison of the spectral behaviour of CPD -28 2561 with that of another Of?p star, HD148937 of similar spectral type, reveals remarkable differences in the degree of variability between both stars. Finally, we present new X-ray observations obtained with the Suzaku X-ray Observatory. We report that the star is X-ray bright with log L_X/L_bol ~= -5.7. The low resolution X-ray spectra reveal the presence of a plasma heated up to 24MK. We associate the 24MK plasma in CPD -28 2561 with the presence of a kG strong magnetic field capable to confine stellar wind.

Within the context of the "B fields in OB stars (BOB)" collaboration, we used the HARPSpol spectropolarimeter to observe the early B-type stars beta CMa (HD44743; B1 II/III) and epsilon CMa (HD52089; B1.5 II). For both stars, we consistently detected the signature of a weak (<30 G in absolute value) longitudinal magnetic field. We determined the physical parameters of both stars and characterise their X-ray spectrum. For beta CMa, our mode identification analysis led to determining a rotation period of 13.6+/-1.2 days and of an inclination angle of the rotation axis of 57.6+/-1.7 degrees, with respect to the line of sight. On the basis of these measurements and assuming a dipolar field geometry, we derived a best fitting obliquity of ~22 degrees and a dipolar magnetic field strength (Bd) of ~100 G (60<Bd<230 G within 1 sigma), below what is typically found for other magnetic massive stars. For epsilon CMa we could only determine a lower limit on the dipolar magnetic field strength of 13 G. For this star, we determine that the rotation period ranges between 1.3 and 24 days. Both stars are expected to have a dynamical magnetosphere. We also conclude that both stars are most likely core hydrogen burning and that they have spent more than 2/3 of their main sequence lifetime. A histogram of the distribution of the dipolar magnetic field strength for the magnetic massive stars known to date does not show the magnetic field "desert" observed instead for intermediate-mass stars. The biases involved in the detection of (weak) magnetic fields in massive stars with the currently available instrumentation and techniques imply that weak fields might be more common than currently observed. Our results show that, if present, even relatively weak magnetic fields are detectable in massive stars and that more observational effort is probably still needed to properly access the magnetic field incidence.

We present the XMM-Newton discovery of X-ray emission from the planetary nebula (PN) A78, the second born-again PN detected in X-rays apart from A30. These two PNe share similar spectral and morphological characteristics: They harbor diffuse soft X-ray emission associated with the interaction between the H-poor ejecta and the current fast stellar wind, and a point-like source at the position of the central star (CSPN). We present the spectral analysis of the CSPN, using for the first time a NLTE code for expanding atmospheres which takes line blanketing into account for the UV and optical spectra. The wind abundances are used for the X-ray spectral analysis of the CSPN and the diffuse emission. The X-ray emission from the CSPN in A78 can be modeled by a single C VI emission line, while the X-ray emission from its diffuse component is better described by an optically thin plasma emission model with temperature $kT$=0.088 keV ($T\approx$1.0$\times$10${^6}$ K). We estimate X-ray luminosities in the 0.2--2.0 keV energy band of $L_{\mathrm{X,CSPN}}$=(1.2$\pm$0.3)$\times$10$^{31}$ erg~s$^{-1}$ and $L_{\mathrm{X,DIFF}}$=(9.2$\pm$2.3)$\times$10$^{30}$ erg~s$^{-1}$ for the CSPN and diffuse components, respectively.

We present the results of a broad-band X-ray study of the enigmatic Be star Gamma Cassiopeiae (herein gamma Cas) based on observations made with both the Suzaku and INTEGRAL observatories. gamma Cas has long been recognized as the prototypical example of a small subclass of Be stars with moderately strong X-ray emission dominated by a hot thermal component in the 0.5-12 keV energy range Lx ~ 10^32 - 10^33 erg s^-1. This places them at the high end of the known luminosity distribution for stellar emission, but several orders of magnitude below typical accretion powered Be X-ray binaries. The INTEGRAL observations spanned an 8 year baseline and represent the deepest measurement to date at energies above ~50 keV. We find that the INTEGRAL data are consistent within statistics to a constant intensity source above 20 keV, with emission extending up to ~100 keV and that searches for all of the previously reported periodicities of the system at lower energies led to null results. We further find that our combined Suzaku and INTEGRAL spectrum, which we suggest is the most accurate broad-band X-ray measurement of gamma Cas to date, is fitted extremely well with a thermal plasma emission model with a single absorption component. We found no compelling need for an additional non-thermal high-energy component. We discuss these results in the context of a currently favored models for gamma Cas and its analogs.

We analyze the expansion of hydrogen-poor knots and filaments in the born-again planetary nebulae A30 and A78 based on Hubble Space Telescope (HST) images obtained almost 20 yr apart. The proper motion of these features generally increases with distance to the central star, but the fractional expansion decreases, i.e., the expansion is not homologous. As a result, there is not a unique expansion age, which is estimated to be 610-950 yr for A30 and 600-1140 yr for A78. The knots and filaments have experienced complex dynamical processes: the current fast stellar wind is mass loaded by the material ablated from the inner knots; the ablated material is then swept up until it shocks the inner edges of the outer, hydrogen-rich nebula. The angular expansion of the outer filaments shows a clear dependence on position angle, indicating that the interaction of the stellar wind with the innermost knots channels the wind along preferred directions. The apparent angular expansion of the innermost knots seems to be dominated by the rocket effect of evaporating gas and by the propagation of the ionization front inside them. Radiation-hydrodynamical simulations show that a single ejection of material followed by a rapid onset of the stellar wind and ionizing flux can reproduce the variety of clumps and filaments at different distances from the central star found in A30 and A78.

The B fields in OB stars (BOB) survey is an ESO large programme collecting spectropolarimetric observations for a large number of early-type stars in order to study the occurrence rate, properties, and ultimately the origin of magnetic fields in massive stars. As of July 2014, a total of 98 objects were observed over 20 nights with FORS2 and HARPSpol. Our preliminary results indicate that the fraction of magnetic OB stars with an organised, detectable field is low. This conclusion, now independently reached by two different surveys, has profound implications for any theoretical model attempting to explain the field formation in these objects. We discuss in this contribution some important issues addressed by our observations (e.g., the lower bound of the field strength) and the discovery of some remarkable objects.

A very sensitive X-ray investigation of the giant HII region N11 in the LMC was performed using the Chandra X-ray Observatory. The 300ks observation reveals X-ray sources with luminosities down to 10^32 erg/s, increasing by more than a factor of 5 the number of known point sources in the field. Amongst these detections are 13 massive stars (3 compact groups of massive stars, 9 O-stars and one early B-star) with log(Lx/Lbol)~-6.5 to -7, which may suggest that they are highly magnetic or colliding wind systems. On the other hand, the stacked signal for regions corresponding to undetected O-stars yields log(Lx/Lbol)~-7.3, i.e., an emission level comparable to similar Galactic stars despite the lower metallicity. Other point sources coincide with 11 foreground stars, 6 late-B/A stars in N11, and many background objects. This observation also uncovers the extent and detailed spatial properties of the soft, diffuse emission regions but the presence of some hotter plasma in their spectra suggests contamination by the unresolved stellar population.

X-ray emission from stars much more massive than the Sun was discovered only 35 years ago. Such stars drive fast stellar winds where shocks can develop, and it is commonly assumed that the X-rays emerge from the shock-heated plasma. Many massive stars additionally pulsate. However, hitherto it was neither theoretically predicted nor observed that these pulsations would affect their X-ray emission. Here we report the discovery of pulsating X-rays from the massive B-type star Xi1 Canis Majoris. This star is a variable of beta Cephei type and has a strong magnetic field. Our observations with the XMM-Newton telescope reveal X-ray pulsations with the same period as the fundamental stellar pulsation. This discovery challenges our understanding of stellar winds from massive stars, their X-ray emission, and their magnetism.

Aims: Following our comprehensive studies of the WR stars in the Milky Way, we now present spectroscopic analyses of almost all known WN stars in the LMC. Methods: For the quantitative analysis of the wind-dominated emission-line spectra, we employ the Potsdam Wolf-Rayet (PoWR) model atmosphere code. By fitting synthetic spectra to the observed spectral energy distribution and the available spectra (ultraviolet and optical), we obtain the physical properties of 107 stars. Results: We present the fundamental stellar and wind parameters for an almost complete sample of WN stars in the LMC. Among those stars that are putatively single, two different groups can be clearly distinguished. While 12% of our sample are more luminous than 10^6 Lsun and contain a significant amount of hydrogen, 88% of the WN stars, with little or no hydrogen, populate the luminosity range between log (L/Lsun) = 5.3...5.8. Conclusions: While the few extremely luminous stars (log (L/Lsun) > 6), if indeed single stars, descended directly from the main sequence at very high initial masses, the bulk of WN stars have gone through the red-supergiant phase. According to their luminosities in the range of log (L/Lsun) = 5.3...5.8, these stars originate from initial masses between 20 and 40 Msun. This mass range is similar to the one found in the Galaxy, i.e. the expected metallicity dependence of the evolution is not seen. Current stellar evolution tracks, even when accounting for rotationally induced mixing, still partly fail to reproduce the observed ranges of luminosities and initial masses. Moreover, stellar radii are generally larger and effective temperatures correspondingly lower than predicted from stellar evolution models, probably due to subphotospheric inflation.

We analyze a 162 ks HETG Chandra observation of the O7.5 III(n)((f)) star xi Per, together with contemporaneous H alpha observations. The X-ray spectrum of this star is similar to other single O stars, and not pathological in any way. Its UV wind lines are known to display cyclical time variability, with a period of 2.086 days, which is thought to be associated with co-rotating interaction regions (CIRs). We examine the Chandra and H alpha data for variability on this time scale. We find that the X-rays vary by about 15% over the course of the observations and that this variability is out of phase with variable absorption on the blue wing of the H alpha profiles (assumed to be a surrogate for the UV absorption associated with CIRs). While not conclusive, both sets of data are consistent with models where the CIRs are either a source of X-rays or modulate them.

Recent magnetic field surveys in O- and B-type stars revealed that about 10% of the core-hydrogen-burning massive stars host large-scale magnetic fields. The physical origin of these fields is highly debated. To identify and model the physical processes responsible for the generation of magnetic fields in massive stars, it is important to establish whether magnetic massive stars are found in very young star-forming regions or whether they are formed in close interacting binary systems. In the framework of our ESO Large Program, we carried out low-resolution spectropolarimetric observations with FORS2 in 2013 April of the three most massive central stars in the Trifid nebula, HD164492A, HD164492C, and HD164492D. These observations indicated a strong longitudinal magnetic field of about 500-600G in the poorly studied component HD164492C. To confirm this detection, we used HARPS in spectropolarimetric mode on two consecutive nights in 2013 June. Our HARPS observations confirmed the longitudinal magnetic field in HD164492C. Furthermore, the HARPS observations revealed that HD164492C cannot be considered as a single star as it possesses one or two companions. The spectral appearance indicates that the primary is most likely of spectral type B1-B1.5V. Since in both observing nights most spectral lines appear blended, it is currently unclear which components are magnetic. Long-term monitoring using high-resolution spectropolarimetry is necessary to separate the contribution of each component to the magnetic signal. Given the location of the system HD164492C in one of the youngest star formation regions, this system can be considered as a Rosetta Stone for our understanding of the origin of magnetic fields in massive stars.

Very massive stars, 100 times heavier than the sun, are rare. It is not yet known whether such stars can form in isolation or only in star clusters. The answer to this question is of fundamental importance. The central region of our Galaxy is ideal for investigating very massive stars and clusters located in the same environment. We used archival infrared images to investigate the surroundings of apparently isolated massive stars presently known in the Galactic Center. We find that two such isolated massive stars display apparent bow shocks and hence may be "runaways" from their birthplace. Thus, some isolated massive stars in the Galactic Center region might have been born in star clusters known in this region. However, no bow shock is detected around the isolated star WR102ka (Peony nebula star), which is one of the most massive and luminous stars in the Galaxy. This star is located at the center of an associated dusty circumstellar nebula. To study whether a star cluster may be "hidden" in the surroundings of WR102ka, to obtain new and better spectra of this star, and to measure its radial velocity, we obtained observations with the integral-field spectrograph SINFONI at the ESO's Very Large Telescope (VLT). Our observations confirm that WR102ka is one of the most massive stars in the Galaxy and reveal that this star is not associated with a star cluster. We suggest that WR102ka has been born in relative isolation, outside of any massive star cluster.

Recent studies of O-type stars demonstrated that discrepant mass-loss rates are obtained when different diagnostic methods are employed - fitting the unsaturated UV resonance lines (e.g. P v) gives drastically lower values than obtained from the H{\alpha} emission. Wind clumping may be the main cause for this discrepancy. In a previous paper, we have presented 3-D Monte-Carlo calculations for the formation of scattering lines in a clumped stellar wind. In the present paper we select five O-type supergiants (from O4 to O7) and test whether the reported discrepancies can be resolved this way. In the first step, the analyses start with simulating the observed spectra with Potsdam Wolf-Rayet (PoWR) non-LTE model atmospheres. The mass-loss rates are adjusted to fit best to the observed H{\alpha} emission lines. For the unsaturated UV resonance lines (i.e. P v) we then apply our 3-D Monte-Carlo code, which can account for wind clumps of any optical depths, a non-void inter-clump medium, and a velocity dispersion inside the clumps. The ionization stratifications and underlying photospheric spectra are adopted from the PoWR models. From fitting the observed resonance line profiles, the properties of the wind clumps are constrained. Our results show that with the mass-loss rates that fit H{\alpha} (and other Balmer and He II lines), the UV resonance lines (especially the unsaturated doublet of P v) can also be reproduced without problem when macroclumping is taken into account. There is no need to artificially reduce the mass-loss rates, nor to assume a sub-solar phosphorus abundance or an extremely high clumping factor, contrary to what was claimed by other authors. These consistent mass-loss rates are lower by a factor of 1.3 to 2.6, compared to the mass-loss rate recipe from Vink et al. Macroclumping resolves the previously reported discrepancy between H{\alpha} and P v mass-loss diagnostics.

We obtained four pointings of over 100 ks each of the well-studied Wolf-Rayet star WR 6 with the XMM-Newton satellite. With a first paper emphasizing the results of spectral analysis, this follow-up highlights the X-ray variability clearly detected in all four pointings. However, phased light curves fail to confirm obvious cyclic behavior on the well-established 3.766 d period widely found at longer wavelengths. The data are of such quality that we were able to conduct a search for "event clustering" in the arrival times of X-ray photons. However, we fail to detect any such clustering. One possibility is that X-rays are generated in a stationary shock structure. In this context we favor a co-rotating interaction region (CIR) and present a phenomenological model for X-rays from a CIR structure. We show that a CIR has the potential to account simultaneously for the X-ray variability and constraints provided by the spectral analysis. Ultimately, the viability of the CIR model will require both intermittent long-term X-ray monitoring of WR 6 and better physical models of CIR X-ray production at large radii in stellar winds.

The Be/X-ray binary SXP1062 is of especial interest owing to the large spin period of the neutron star, its large spin-down rate, and the association with a supernova remnant constraining its age. This makes the source an important probe for accretion physics. To investigate the long-term evolution of the spin period and associated spectral variations, we performed an XMM-Newton target-of-opportunity observation of SXP1062 during X-ray outburst. Spectral and timing analysis of the XMM-Newton data was compared with previous studies, as well as complementary Swift/XRT monitoring and optical spectroscopy with the SALT telescope were obtained. The spin period was measured to be P_s =(1071.01+-0.16) s on 2012 Oct 14. The X-ray spectrum is similar to that of previous observations. No convincing cyclotron absorption features, which could be indicative for a high magnetic field strength, are found. The high-resolution RGS spectra indicate the presence of emission lines, which may not completely be accounted for by the SNR emission. The comparison of multi-epoch optical spectra suggest an increasing size or density of the decretion disc around the Be star. SXP1062 showed a net spin-down with an average of (2.27+-0.44) s/yr over a baseline of 915 days.

Stars over a wide range of masses and evolutionary stages are nowadays known to emit X-rays. This X-ray emission is a unique probe of the most energetic phenomena occurring in the circumstellar environment of these stars, and provides precious insight on magnetic phenomena or hydrodynamic shocks. Owing to its large collecting area, Athena+ will open up an entirely new window on these phenomena. Indeed, Athena+ will not only allow us to study many more objects with an unprecedented spectral resolution, but will also pioneer the study of the dynamics of these objects via time-resolved high-resolution spectroscopy. In this way, Athena+ will be a unique tool to study accretion processes in TTauri stars, flaring activity in young stars, dynamos in ultra-cool dwarfs, small and large-scale structures in the winds of single massive stars, wind interactions in massive binary systems, hot bubbles in planetary nebula... All these studies will lead to a deeper understanding of yet poorly understood processes which have profound impact in star and planetary system formation as well as in feedback processes on Galactic scale.

This White Paper, submitted to the recent ESA call for science themes to define its future large missions, advocates the need for a transformational leap in our understanding of two key questions in astrophysics: 1) How does ordinary matter assemble into the large scale structures that we see today? 2) How do black holes grow and shape the Universe? Hot gas in clusters, groups and the intergalactic medium dominates the baryonic content of the local Universe. To understand the astrophysical processes responsible for the formation and assembly of these large structures, it is necessary to measure their physical properties and evolution. This requires spatially resolved X-ray spectroscopy with a factor 10 increase in both telescope throughput and spatial resolving power compared to currently planned facilities. Feedback from supermassive black holes is an essential ingredient in this process and in most galaxy evolution models, but it is not well understood. X-ray observations can uniquely reveal the mechanisms launching winds close to black holes and determine the coupling of the energy and matter flows on larger scales. Due to the effects of feedback, a complete understanding of galaxy evolution requires knowledge of the obscured growth of supermassive black holes through cosmic time, out to the redshifts where the first galaxies form. X-ray emission is the most reliable way to reveal accreting black holes, but deep survey speed must improve by a factor ~100 over current facilities to perform a full census into the early Universe. The Advanced Telescope for High Energy Astrophysics (Athena+) mission provides the necessary performance (e.g. angular resolution, spectral resolution, survey grasp) to address these questions and revolutionize our understanding of the Hot and Energetic Universe. These capabilities will also provide a powerful observatory to be used in all areas of astrophysics.

Recent studies have claimed the existence of very massive stars (VMS) up to 300 solar masses in the local Universe. As this finding may represent a paradigm shift for the canonical stellar upper-mass limit of 150 Msun, it is timely to discuss the status of the data, as well as the far-reaching implications of such objects. We held a Joint Discussion at the General Assembly in Beijing to discuss (i) the determination of the current masses of the most massive stars, (ii) the formation of VMS, (iii) their mass loss, and (iv) their evolution and final fate. The prime aim was to reach broad consensus between observers and theorists on how to identify and quantify the dominant physical processes.

Theories on the origin of magnetic fields in massive stars remain poorly developed, because the properties of their magnetic field as function of stellar parameters could not yet be investigated. To investigate whether magnetic fields in massive stars are ubiquitous or appear only in stars with a specific spectral classification, certain ages, or in a special environment, we acquired 67 new spectropolarimetric observations for 30 massive stars. Among the observed sample, roughly one third of the stars are probable members of clusters at different ages, whereas the remaining stars are field stars not known to belong to any cluster or association. Spectropolarimetric observations were obtained during four different nights using the low-resolution spectropolarimetric mode of FORS2 (FOcal Reducer low dispersion Spectrograph) mounted on the 8-m Antu telescope of the VLT. Furthermore, we present a number of follow-up observations carried out with the high-resolution spectropolarimeters SOFIN mounted at the Nordic Optical Telescope (NOT) and HARPS mounted at the ESO 3.6m between 2008 and 2011. To assess the membership in open clusters and associations, we used astrometric catalogues with the highest quality kinematic and photometric data currently available. The presence of a magnetic field is confirmed in nine stars previously observed with FORS1/2: HD36879, HD47839, CPD-282561, CPD-472963, HD93843, HD148937, HD149757, HD328856, and HD164794. New magnetic field detections at a significance level of at least 3sigma were achieved in five stars: HD92206c, HD93521, HD93632, CPD-468221, and HD157857. Among the stars with a detected magnetic field, five stars belong to open clusters with high membership probability. According to previous kinematic studies, five magnetic O-type stars in our sample are candidate runaway stars.

We report the discovery of extended X-ray emission within the young star cluster NGC 602 in the Wing of the Small Magellanic Cloud (SMC) based on observations obtained with the Chandra X-ray Observatory. X-ray emission is detected from the cluster core area with the highest stellar density and from a dusty ridge surrounding the HII region. We use a census of massive stars in the cluster to demonstrate that a cluster wind or wind-blown bubble is unlikely to provide a significant contribution to the X-ray emission detected from the central area of the cluster. We therefore suggest that X-ray emission at the cluster core originates from an ensemble of low- and solar-mass pre-main-sequence (PMS) stars, each of which would be too weak in X-rays to be detected individually. We attribute the X-ray emission from the dusty ridge to the embedded tight cluster of the new-born stars known in this area from infrared studies. Assuming that the levels of X-ray activity in young stars in the low-metallicity environment of NGC 602a are comparable to their Galactic counterparts, then the detected spatial distribution, spectral properties, and level of X-ray emission are largely consistent with those expected from low- and solar-mass PMS stars and young stellar objects (YSOs). This is the first discovery of X-ray emission attributable to PMS stars and YSOs in the SMC, which suggests that the accretion and dynamo processes in young, low-mass objects in the SMC resemble those in the Galaxy.

Stellar winds are a crucial component of massive stars, but their exact properties still remain uncertain. To shed some light on this subject, we have analyzed an exceptional set of X-ray observations of zeya Pup, one of the closest and brightest massive stars. The sensitive lightcurves that were derived reveal two major results. On the one hand, a slow modulation of the X-ray flux (with a relative amplitude of up to 15% over 16h in the 0.3--4.0keV band) is detected. Its characteristic timescale cannot be determined with precision, but amounts from one to several days. It could be related to corotating interaction regions, known to exist in zeta Pup from UV observations. Hour-long changes, linked to flares or to the pulsation activity, are not observed in the last decade covered by the XMM observations; the 17h tentative period, previously reported in a ROSAT analysis, is not confirmed either and is thus transient, at best. On the other hand, short-term changes are surprisingly small (<1% relative amplitude for the total energy band). In fact, they are compatible solely with the presence of Poisson noise in the data. This surprisingly low level of short-term variability, in view of the embedded wind-shock origin, requires a very high fragmentation of the stellar wind, for both absorbing and emitting features (>10^5 parcels, comparing with a 2D wind model). This is the first time that constraints have been placed on the number of clumps in an O-type star wind and from X-ray observations.

An increasing number of OB stars have been shown to possess magnetic fields. Although the sample remains small, it is surprising that the magnetic and X-ray properties of these stars appear to be far less correlated than expected. This contradicts model predictions, which generally indicate that the X-rays from magnetic stars to be harder and more luminous than their non-magnetic counterparts. Instead, the X-ray properties of magnetic OB stars are quite diverse. $\tau$ Sco is one example where the expectations are better met. This bright main sequence, early B star has been studied extensively in a variety of wavebands. It has a surface magnetic field of around 500 G, and Zeeman Doppler tomography has revealed an unusual field configuration. Furthermore, tau Sco displays an unusually hard X-ray spectrum, much harder than similar, non-magnetic OB stars. In addition, the profiles of its UV P Cygni wind lines have long been known to possess a peculiar morphology. Recently, two stars, HD 66665 and HD 63425, whose spectral types and UV wind line profiles are similar to those of $\tau$ Sco, have also been determined to be magnetic. In the hope of establishing a magnetic field - X-ray connection for at least a sub-set of the magnetic stars, we obtained XMM-Newton EPIC spectra of these two objects. Our results for HD 66665 are somewhat inconclusive. No especially strong hard component is detected; however, the number of source counts is insufficient to rule out hard emission. longer exposure is needed to assess the nature of the X-rays from this star. On the other hand, we do find that HD 63425 has a substantial hard X-ray component, thereby bolstering its close similarity to tau Sco.

The clumping of massive star winds is an established paradigm, which is confirmed by multiple lines of evidence and is supported by stellar wind theory. We use the results from time-dependent hydrodynamical models of the instability in the line-driven wind of a massive supergiant star to derive the time-dependent accretion rate on to a compact object in the Bondi-Hoyle-Lyttleton approximation. The strong density and velocity fluctuations in the wind result in strong variability of the synthetic X-ray light curves. Photoionization of inhomogeneous winds is different from the photoinization of smooth winds. The degree of ionization is affected by the wind clumping. The wind clumping must also be taken into account when comparing the observed and model spectra of the photoionized stellar wind.

SXP 1062 is an exceptional case of a young neutron star in a wind-fed high-mass X-ray binary associated with a supernova remnant. A unique combination of measured spin period, its derivative, luminosity and young age makes this source a key probe for the physics of accretion and neutron star evolution. Theoretical models proposed to explain the properties of SXP 1062 shall be tested with new data.

\mu\ Columbae is a prototypical weak-wind O-star for which we have obtained a high-resolution X-ray spectrum with the Chandra LETG/ACIS-S instrument and a low resolution spectrum with Suzaku. This allows us, for the first time, to investigate the role of X-rays on the wind structure in a bona fide weak-wind system and to determine whether there actually is a massive, hot wind. The X-ray emission measure indicates that the outflow is an order of magnitude greater than that derived from UV lines and is commensurate with the nominal wind-luminosity relationship for O-stars. Therefore, the ``weak-wind problem''---identified from cool wind UV/optical spectra---is largely resolved by accounting for the hot wind seen in X-rays. From X-ray line profiles, Doppler shifts, and relative strengths, we find that this weak-wind star is typical of other late O dwarfs. The X-ray spectra do not suggest a magnetically confined plasma---the spectrum is soft and lines are broadened; Suzaku spectra confirm the lack of emission above 2 keV. Nor do the relative line shifts and widths suggest any wind decoupling by ions. The He-like triplets indicate that the bulk of the X-ray emission is formed rather close to the star, within 5 stellar radii. Our results challenge the idea that some OB stars are ``weak-wind'' stars that deviate from the standard wind-luminosity relationship. The wind is not weak, but it is hot and its bulk is only detectable in X-rays.

The planetary nebula (PN) A30 is believed to have undergone a very late thermal pulse resulting in the ejection of knots of hydrogen-poor material. Using HST images we have detected the angular expansion of these knots and derived an age of 850+280-150 yr. To investigate the spectral and spatial properties of the soft X-ray emission detected by ROSAT, we have obtained Chandra and XMM-Newton observations of A30. The X-ray emission from A30 can be separated into two components: a point-source at the central star and diffuse emission associated with the hydrogen-poor knots and the cloverleaf structure inside the nebular shell. To help us assess the role of the current stellar wind in powering this X-ray emission, we have determined the stellar parameters of the central star of A 30 using a non-LTE model fit to its optical and UV spectrum. The spatial distribution and spectral properties of the diffuse X-ray emission is suggestive that it is generated by the post-born-again and present fast stellar winds interacting with the hydrogen-poor ejecta of the born-again event. This emission can be attributed to shock-heated plasma, as the hydrogen-poor knots are ablated by the stellar winds, under which circumstances the efficient mass-loading of the present fast stellar wind raises its density and damps its velocity to produce the observed diffuse soft X-rays. Charge transfer reactions between the ions of the stellar winds and material of the born-again ejecta has also been considered as a possible mechanism for the production of diffuse X-ray emission, and upper limits on the expected X-ray production by this mechanism have been derived. The origin of the X-ray emission from the central star of A 30 is puzzling: shocks in the present fast stellar wind and photospheric emission can be ruled out, while the development of a new, compact hot bubble confining the fast stellar wind seems implausible.

Sk 183 is the visually-brightest star in the N90 nebula, a young star-forming region in the Wing of the Small Magellanic Cloud (SMC). We present new optical spectroscopy from the Very Large Telescope which reveals Sk 183 to be one of the most massive O-type stars in the SMC. Classified as an O3-type dwarf on the basis of its nitrogen spectrum, the star also displays broadened He I absorption which suggests a later type. We propose that Sk 183 has a composite spectrum and that it is similar to another star in the SMC, MPG 324. This brings the number of rare O2- and O3-type stars known in the whole of the SMC to a mere four. We estimate physical parameters for Sk 183 from analysis of its spectrum. For a single-star model, we estimate an effective temperature of 46+/-2 kK, a low mass-loss rate of ~10^-7 Msun yr^-1, and a spectroscopic mass of 46^+9_-8 Msun (for an adopted distance modulus of 18.7 mag to the young population in the SMC Wing). An illustrative binary model requires a slightly hotter temperature (~47.5 kK) for the primary component. In either scenario, Sk 183 is the earliest-type star known in N90 and will therefore be the dominant source of hydrogen-ionising photons. This suggests Sk 183 is the primary influence on the star formation along the inner edge of the nebula.

We present the first high-resolution X-ray spectrum of a putatively single Wolf-Rayet star. 400 ks observations of WR 6 by the XMM-Newton-telescope resulted in a superb quality high-resolution X-ray spectrum. Spectral analysis reveals that the X-rays originate far out in the stellar wind, more than 30 stellar radii from the photosphere, and thus outside the wind acceleration zone where the line-driving instability could create shocks. The X-ray emitting plasma reaches temperatures up to 50\,MK, and is embedded within the un-shocked, "cool" stellar wind as revealed by characteristic spectral signatures. We detect a fluorescent Fe line at approx 6.4 keV. The presence of fluorescence is consistent with a two-component medium, where the cool wind is permeated with the hot X-ray emitting plasma. The wind must have a very porous structure to allow the observed amount of X-rays to escape. We find that neither the line-driving instability nor any alternative binary scenario can explain the data. We suggest a scenario where X-rays are produced when the fast wind rams into slow "sticky clumps" that resist acceleration. Our new data show that the X-rays in single WR-star are generated by some special mechanism different from the one operating in the O-star winds.

The Quintuplet, one of three massive stellar clusters in the Galactic center, is located about 30pc in projection from Sagittarius A*. Based on near-infrared K-band spectra we determine temperatures and luminosities for all stars in our sample and construct the Herztsprung-Russell diagram. We find two distinct groups: early-type OB stars and late-type KM stars, well separated from each other. By comparison with Geneva stellar evolution models we derive initial masses exceeding 8 solar masses for the OB stars, that are located along an isochrone corresponding to a cluster age of about 4 million years. In addition, we derive number ratios (e. g. N_WR/N_O) and compare them with predictions of population synthesis models. We find that an instantaneous burst of star formation at about 3.3 to 3.6\,Myr ago is the most likely scenario to form the Quintuplet cluster. The late-type stars in the sample are red giant branch (RGB) stars or red supergiants (RSGs) according to their spectral signatures. It is discussed if they could physically belong to the Quintuplet cluster. Furthermore, we apply a mass-luminosity relation to construct the initial mass function (IMF) of the cluster. We find indications for a slightly top-heavy IMF.

The true mass-loss rates from massive stars are important for many branches of astrophysics. For the correct modeling of the resonance lines, which are among the key diagnostics of stellar mass-loss, the stellar wind clumping turned out to be very important. In order to incorporate clumping into radiative transfer calculation, 3-D models are required. Various properties of the clumps may have strong impact on the resonance line formation and, therefore, on the determination of empirical mass-loss rates. We incorporate the 3-D nature of the stellar wind clumping into radiative transfer calculations and investigate how different model parameters influence the resonance line formation. We develop a full 3-D Monte Carlo radiative transfer code for inhomogeneous expanding stellar winds. The number density of clumps follows the mass conservation. For the first time, realistic 3-D models that describe the dense as well as the tenuous wind components are used to model the formation of resonance lines in a clumped stellar wind. At the same time, non-monotonic velocity fields are accounted for. The 3-D density and velocity wind inhomogeneities show very strong impact on the resonance line formation. The different parameters describing the clumping and the velocity field results in different line strengths and profiles. We present a set of representative models for various sets of model parameters and investigate how the resonance lines are affected. Our 3-D models show that the line opacity is reduced for larger clump separation and for more shallow velocity gradients within the clumps. Our new model demonstrates that to obtain empirically correct mass-loss rates from the UV resonance lines, the wind clumping and its 3-D nature must be taken into account.

We study the effects of optically thick clumps, non-void inter-clump medium, variation of the onset of clumping, and velocity dispersion inside clumps on the formation of resonance lines. For this purpose we developed a full 3-D Monte Carlo Radiative Transfer (MCRT) code that is able to handle 3-D shapes of clumps and arbitrary 3-D velocity fields. The method we developed allows us to take into account contributions from density and velocity wind inhomogeneities to the total opacity very precisely. The first comparison with observation shows that 3-D density and velocity wind inhomogeneities have a very strong influence on the resonance line formation, and that they have to be accounted for in order to obtain reliable mass-loss rate determinations.

The clumping of massive star winds is an established paradigm confirmed by multiple lines of evidence and supported by stellar wind theory. The purpose of this paper is to bridge the gap between detailed models of inhomogeneous stellar winds in single stars and the phenomenological description of donor winds in supergiant high-mass X-ray binaries (HMXBs). We use results from time-dependent hydrodynamical models of the instability in the line-driven wind of a massive supergiant star to derive the time-dependent accretion rate onto a compact object in the Bondi-Hoyle-Lyttleton approximation. The strong density and velocity fluctuations in the wind result in strong variability of the synthetic X-ray light curves. The model predicts a large scale X-ray variability, up to eight orders of magnitude, on relatively short timescales. The apparent lack of evidence for such strong variability in the observed HMXBs indicates that the details of accretion process act to reduce the variability due to the stellar wind velocity and density jumps. We, also, study the absorption of X-rays in the clumped stellar wind by means of a 2-D stochastic wind model and find that absorption of X-rays changes strongly at different orbital phases. Furthermore, we address the photoionization in the clumped wind, and show that the degree of ionization is affected by the wind clumping. A correction factor for the photoionization parameter is derived. It is shown that the photoionization parameter is reduced by a factor Xi compared to the smooth wind models with the same mass-loss rate, where Xi is the wind inhomogeneity parameter. We conclude that wind clumping must also be taken into account when comparing the observed and model spectra of the photoionized stellar wind.

We report on a new Be/X-ray pulsar binary located in the Wing of the Small Magellanic Cloud (SMC). The strong pulsed X-ray source was discovered with the Chandra and XMM-Newton X-ray observatories. The X-ray pulse period of 1062 s is consistently determined from both Chandra and XMM-Newton observations, revealing one of the slowest rotating X-ray pulsars known in the SMC. The optical counterpart of the X-ray source is the emission-line star 2dFS 3831. Its B0-0.5(III)e+ spectral type is determined from VLT-FLAMES and 2dF optical spectroscopy, establishing the system as a Be/X-ray binary (Be-XRB). The hard X-ray spectrum is well fitted by a power-law with additional thermal and blackbody components, the latter reminiscent of persistent Be-XRBs. This system is the first evidence of a recent supernova in the low density surroundings of NGC 602. We detect a shell nebula around 2dFS 3831 in H-alpha and [O III] images and conclude that it is most likely a supernova remnant. If it is linked to the supernova explosion that created this new X-ray pulsar, its kinematic age of (2-4)x10^4 yr provides a constraint on the age of the pulsar.

We investigate the connections between the magnetic fields and the X-ray emission from massive stars. Our study shows that the X-ray properties of known strongly magnetic stars are diverse: while some comply to the predictions of the magnetically confined wind model, others do not. We conclude that strong, hard, and variable X-ray emission may be a sufficient attribute of magnetic massive stars, but it is not a necessary one. We address the general properties of X-ray emission from "normal" massive stars, especially the long standing mystery about the correlations between the parameters of X-ray emission and fundamental stellar properties. The recent development in stellar structure modeling shows that small scale surface magnetic fields may be common. We suggest a "hybrid" scenario which could explain the X-ray emission from massive stars by a combination of magnetic mechanisms on the surface and shocks in the stellar wind. The magnetic mechanisms and the wind shocks are triggered by convective motions in sub-photospheric layers. This scenario opens the door for a natural explanation of the well established correlation between bolometric and X-ray luminosities.

We report on the status of our spectropolarimetric observations of massive stars. During the last years, we have discovered magnetic fields in many objects of the upper main sequence, including Be stars, beta Cephei and Slowly Pulsating B stars, and a dozen O stars. Since the effects of those magnetic fields have been found to be substantial by recent models, we are looking into their impact on stellar rotation, pulsation, stellar winds, and chemical abundances. Accurate studies of the age, environment, and kinematic characteristics of the magnetic stars are also promising to give us new insight into the origin of the magnetic fields. Furthermore, longer time series of magnetic field measurements allow us to observe the temporal variability of the magnetic field and to deduce the stellar rotation period and the magnetic field geometry. Studies of the magnetic field in massive stars are indispensable to understand the conditions controlling the presence of those fields and their implications on the stellar physical parameters and evolution.

We present a comprehensive study of X-ray emission and wind properties of massive magnetic early B-type stars. Dedicated XMM-Newton observations were obtained for three stars xi1 CMa, V2052 Oph, and zeta Cas. We report the first detection of X-ray emission from V2052 Oph and zeta Cas. The observations show that the X-ray spectra of our program stars are quite soft. We compile the complete sample of early B-type stars with detected magnetic fields to date and existing X-ray measurements, in order to study whether the X-ray emission can be used as a general proxy for stellar magnetism. We find that hard and strong X-ray emission does not necessarily correlate with the presence of a magnetic field. We analyze the UV spectra of five non-supergiant B stars with magnetic fields by means of non-LTE iron-blanketed model atmospheres. The latter are calculated with the PoWR code, which treats the photosphere as well as the the wind, and also accounts for X-rays. Our models accurately fit the stellar photospheric spectra in the optical and the UV. The parameters of X-ray emission, temperature and flux are included in the model in accordance with observations. We confirm the earlier findings that the filling factors of X-ray emitting material are very high. Our analysis reveals that the magnetic early type B stars studied here have weak winds. The mass-loss rates are significantly lower than predicted by hydrodynamically consistent models. We find that, although the X-rays strongly affect the ionization structure of the wind, this effect is not sufficient in reducing the total radiative acceleration. When the X-rays are accounted for at the intensity and temperatures observed, there is still sufficient radiative acceleration to drive stronger mass-loss than we empirically infer from the UV spectral lines. (abridged)

X-ray observations of the double-binary OB-star system QZ Car (HD 93206) obtained with the Chandra X-ray Observatory over a period of roughly 2 years are presented. The orbit of systems A (O9.7 I+b2 v, PA = 21 d) and B (O8 III+o9 v, PB = 6 d) are reasonably well sampled by the observations, allowing the origin of the X-ray emission to be examined in detail. The X-ray spectra can be well fitted by an attenuated three temperature thermal plasma model, characterised by cool, moderate, and hot plasma components at kT ~ 0.2, 0.7, and 2 keV, respectively, and a circumstellar absorption of ~ 0.2 x 10^22 cm-2. Although the hot plasma component could be indicating the presence of wind-wind collision shocks in the system, the model fluxes calculated from spectral fits, with an average value of ~ 7 x 10^-13 erg s-1 cm-2, do not show a clear correlation with the orbits of the two constituent binaries. A semi-analytical model of QZ Car reveals that a stable momentum balance may not be established in either system A or B. Yet, despite this, system B is expected to produce an observed X-ray flux well in excess of the observations. If one considers the wind of the O8 III star to be disrupted by mass transfer the model and observations are in far better agreement, which lends support to the previous suggestion of mass-transfer in the O8 III + o9 v binary. We conclude that the X-ray emission from QZ Car can be reasonably well accounted for by a combination of contributions mainly from the single stars and the mutual wind-wind collision between systems A and B.

We have developed lists of likely B3--A0 stars (called "late B" stars) in the young cluster Trumpler 16. The following criteria were used: location within 3' of Eta Car, an appropriate V and B-V combination, and proper motion (where available). Color and magnitude cuts have been made assuming an E(B-V) =0.55 mag +/- 0.1, which is a good approximation close to the center of Trumpler 16. These lists have been cross-correlated with X-ray sources found in the Chandra Carina Complex Project (CCCP). Previous studies have shown that only very rarely (if at all) do late main sequence B stars produce X-rays. We present evidence that the X-ray detected sources are binaries with low-mass companions, since stars less massive than 1.4 Msun are strong X-ray sources at the age of the cluster. Both the median X-ray energies and X-ray luminosities of these sources are in good agreement with values for typical low-mass coronal X-ray sources. We find that 39% of the late B stars based on a list with proper motions have low-mass companions. Similarly, 32% of a sample without proper motions have low-mass companions. We discuss the X-ray detection completeness. These results on low-mass companions of intermediate mass stars are complementary to spectroscopic and interferometric results, and probe new parameter space of low mass companions at all separations. They do not support a steeply rising distribution of mass ratios to low masses for intermediate-mass (5 Msun) primaries, such as would be found by random pairing from the Initial Mass Function.

The Chandra Carina Complex contains 200 known O- and B type stars. The Chandra survey detected 68 of the 70 O stars and 61 of 127 known B0-B3 stars. We have assembled a publicly available optical/X-ray database to identify OB stars that depart from the canonical Lx/Lbol relation, or whose average X-ray temperatures exceed 1 keV. Among the single O stars with high kT we identify two candidate magnetically confined wind shock sources: Tr16-22, O8.5 V, and LS 1865, O8.5 V((f)). The O4 III(fc) star HD 93250 exhibits strong, hard, variable X-rays, suggesting it may be a massive binary with a period of >30 days. The visual O2 If* binary HD 93129A shows soft 0.6 keV and hard 1.9 keV emission components, suggesting embedded wind shocks close to the O2 If* Aa primary, and colliding wind shocks between Aa and Ab. Of the 11 known O-type spectroscopic binaries, the long orbital-period systems HD 93343, HD 93403 and QZ Car have higher shock temperatures than short-period systems such as HD 93205 and FO 15. Although the X-rays from most B stars may be produced in the coronae of unseen, low-mass pre-main-sequence companions, a dozen B stars with high Lx cannot be explained by a distribution of unseen companions. One of these, SS73 24 in the Treasure Chest cluster, is a new candidate Herbig Be star.

The key empirical property of the X-ray emission from O stars is a strong correlation between the bolometric and X-ray luminosities. In the framework of the Chandra Carina Complex Project, 129 O and B stars have been detected as X-ray sources; 78 of those, all with spectral type earlier than B3, have enough counts for at least a rough X-ray spectral characterization. This leads to an estimate of the Lx/Lbol ratio for an exceptional number of 60 O stars belonging to the same region and triples the number of Carina massive stars studied spectroscopically in X-rays. The derived log(Lx/Lbol) is -7.26 for single objects, with a dispersion of only 0.21dex. Using the properties of hot massive stars listed in the literature, we compare the X-ray luminosities of different types of objects. In the case of O stars, the Lx/Lbol ratios are similar for bright and faint objects, as well as for stars of different luminosity classes or spectral types. Binaries appear only slightly harder and slightly more luminous in X-rays than single objects; the differences are not formally significant (at the 1% level), except for the Lx/Lbol ratio in the medium (1.0--2.5keV) energy band. Weak-wind objects have similar X-ray luminosities but they display slightly softer spectra compared to "normal" O stars with the same bolometric luminosity. Discarding three overluminous objects, we find a very shallow trend of harder emission in brighter objects. The properties of the few B stars bright enough to yield some spectral information appear to be different overall (constant X-ray luminosities, harder spectra), hinting that another mechanism for producing X-rays, besides wind shocks, might be at work. However, it must be stressed that the earliest and X-ray brightest amongst these few detected objects are similar to the latest O stars, suggesting a possibly smooth transition between the two processes.

The Great Nebula in Carina provides an exceptional view into the violent massive star formation and feedback that typifies giant HII regions and starburst galaxies. We have mapped the Carina star-forming complex in X-rays, using archival Chandra data and a mosaic of 20 new 60ks pointings using the Chandra X-ray Observatory's Advanced CCD Imaging Spectrometer, as a testbed for understanding recent and ongoing star formation and to probe Carina's regions of bright diffuse X-ray emission. This study has yielded a catalog of properties of >14,000 X-ray point sources; >9800 of them have multiwavelength counterparts. Using Chandra's unsurpassed X-ray spatial resolution, we have separated these point sources from the extensive, spatially-complex diffuse emission that pervades the region; X-ray properties of this diffuse emission suggest that it traces feedback from Carina's massive stars. In this introductory paper, we motivate the survey design, describe the Chandra observations, and present some simple results, providing a foundation for the 15 papers that follow in this Special Issue and that present detailed catalogs, methods, and science results.

The star zeta Ophiuchi is one of the brightest massive stars in the northern hemisphere and was intensively studied in various wavelength domains. The currently available observational material suggests that certain observed phenomena are related to the presence of a magnetic field. We acquired spectropolarimetric observations of zeta Oph with FORS1 mounted on the 8-m Kueyen telescope of the VLT to investigate if a magnetic field is indeed present in this star. Using all available absorption lines, we detect a mean longitudinal magnetic field B_z_all = 141+-45G, confirming the magnetic nature of this star. We review the X-ray properties of zeta Oph with the aim to understand whether the X-ray emission of zeta Oph is dominated by magnetic or by wind instability processes.

X-ray emission is ubiquitous among massive stars. In the last decade, X-ray observations revolutionized our perception of stellar winds but opened a Pandora's box of urgent problems. X-rays penetrating stellar winds suffer mainly continuum absorption, which greatly simplifies the radiative transfer treatment. The small and large scale structures in stellar winds must be accounted for to understand the X-ray emission from massive stars. The analysis of X-ray spectral lines can help to infer the parameters of wind clumping, which is prerequisite for obtaining empirically correct stellar mass-loss rates. The imprint of large scale structures, such as CIRs and equatorial disks, on the X-ray emission is predicted, and new observations are testing theoretical expectations. The X-ray emission from magnetic stars proves to be more diverse than anticipated from the direct application of the magnetically-confined wind model. Many outstanding questions about X-rays from massive stars will be answered when the models and the observations advance.

Some of the Wolf-Rayet stars are found to have very high bolometric luminosities (more than 1000000 solar). We employ the Potsdam Wolf-Rayet (PoWR) model atmospheres for their spectral analysis, which yields the bolometric corrections. Distance and interstellar reddening also enter the luminosity estimates. Among the Galactic stars, there is a group of very luminous WNL stars (i.e. WR stars of late subtype from nitrogen sequence with hydrogen being depleted in their atmospheres, but not absent). Their distances are often the major source of uncertainty. From K-band spectroscopy we found a very luminous star ($\log L/L_\odot = 6.5$) in the Galactic center region, which we termed the Peony Star because of the form of its surrounding dusty nebula. A similar group of very luminous WNL stars is found in the Large Magellanic Cloud (LMC). In the Small Magellanic Cloud (SMC) the majority of WR stars resides in binary systems. The single WNL stars in the SMC are not very luminous. We conclude that a significant number of very luminous WNL stars exist in the Galaxy and the LMC. With initial masses above 60$M_\odot$, they apparently evolved directly to the WNL stage without a prior excursion to the red side of the HRD. At the low metallicity of the SMC, the binary channel may be dominant for the formation of WR stars.

Based on $K$-band integral-field spectroscopy, we analyze four Wolf-Rayet stars of the nitrogen sequence (WN) found in the inner part of the Quintuplet cluster. All WN stars (WR102d, WR102i, WR102hb, and WR102ea) are of spectral subtype WN9h. One further star, LHO110, is included in the analysis which has been classified as Of/WN? previously but turns out to be most likely a WN9h star as well. The Potsdam Wolf-Rayet (PoWR) models for expanding atmospheres are used to derive the fundamental stellar and wind parameters. The stars turn out to be very luminous, $\log{(L/L_\odot)} > 6.0$, with relatively low stellar temperatures, $T_* \approx$ 25--35\,kK. Their stellar winds contain a significant fraction of hydrogen, up to $X_\mathrm{H} \sim 0.45$ (by mass). We discuss the position of the Galactic center WN stars in the Hertzsprung-Russell diagram and find that they form a distinct group. In this respect, the Quintuplet WN stars are similar to late-type WN stars found in the Arches cluster and elsewhere in the Galaxy. Comparison with stellar evolutionary models reveals that the Quintuplet WN stars should have been initially more massive than 60\,$M_\odot$. They are about 2.1 to 3.6\,Million years old, and might still be central hydrogen burning objects. The analysis of the spectral energy distributions of the program stars results in a mean extinction of $A_K = 3.1 \pm 0.5$\,mag ($A_V = 27 \pm 4$\,mag) towards the Quintuplet cluster.

We obtained relatively high signal-to-noise X-ray spectral data of the early massive star tau Sco (B0.2V) with the Suzaku XIS instrument. This source displays several unusual features that motivated our study: (a) redshifted absorption in UV P Cygni lines to approximately +250 km/s suggestive of infalling gas, (b) unusually hard X-ray emission requiring hot plasma at temperatures in excess of 10 MK whereas most massive stars show relatively soft X-rays at a few MK, and (c) a complex photospheric magnetic field of open and closed field lines. In an attempt to understand the hard component better, X-ray data were obtained at six roughly equally spaced phases within the same epoch of tau Sco's 41 day rotation period. The XIS instrument has three operable detectors: XIS1 is back illuminated with sensitivity down to 0.2 keV; XIS0 and XIS2 are front illuminated with sensivitity only down to 0.4 keV and have overall less effective area than XIS1. The XIS0 and XIS3 detectors show relatively little variability. In contrast, there is a 4sigma detection of a 4% drop in the count rate of the XIS1 detector at one rotational phase. In addition, all three detectors show a 3% increase in count rate at the same phase. The most optimistic prediction of X-ray variability allows for a 40% change in the count rate, particularly near phases where we have pointings. Observed modulations in the X-ray light curve on the rotation cycle is an order of magnitude smaller than this, which places new stringent constraints on future modeling of this interesting magnetic massive star.

We obtained X-ray XMM-Newton observations of the open cluster Berkeley 87 and the massive star-forming region (SFR) ON 2. In addition, archival infrared Spitzer Space Telescope observations were used. It is likely that the SFR ON 2 and Berkeley 87 are at the same distance, 1.23 kpc, and hence are associated. The XMM-Newton observations detected X-rays from massive stars in Berkeley 87 as well as diffuse emission from the SFR ON 2. The two patches of diffuse X-ray emission are encompassed in the shell-like H II region GAL 75.84+0.40 in the northern part of ON 2 and in the ON 2S region in the southern part of ON 2. The diffuse emission from GAL 75.84+0.40 suffers an absorption column equivalent to A_V approx. 28 mag. Its spectrum can be fitted either with a thermal plasma model at T < 30 MK or by an absorbed power-law model with gamma; approx. -2.6. The X-ray luminosity of GAL 75.84+0.40 is L_X approx. 6 10^31 erg/s. The diffuse emission from ON 2S is adjacent to the ultra-compact H II (UCHII) region Cygnus 2N, but does not coincide with it or with any other known UCHII region. It has a luminosity of L_X approx. 4 10^31 erg/s. The spectrum can be fitted with an absorbed power-law model with gamma; approx.-1.4. We adopt the view of Turner and Forbes (1982) that the SFR ON 2 is physically associated with the massive star cluster Berkeley 87 hosting the WO type star WR 142. We discuss different explanations for the apparently diffuse X-ray emission in these SFRs. These include synchrotron radiation, invoked by the co-existence of strongly shocked stellar winds and turbulent magnetic fields in the star-forming complex, cluster wind emission, or an unresolved population of discrete sources.

Hot stars are sources of X-ray emission originating in their winds. Although hydrodynamical simulations that are able to predict this X-ray emission are available, the inclusion of X-rays in stationary wind models is usually based on simplifying approximations. To improve this, we use results from time-dependent hydrodynamical simulations of the line-driven wind instability (seeded by the base perturbation) to derive the analytical approximation of X-ray emission in the stellar wind. We use this approximation in our non-LTE wind models and find that an improved inclusion of X-rays leads to a better agreement between model ionization fractions and those derived from servations. Furthermore, the slope of the L_x-L relation is in better agreement with observations, however the X-ray luminosity is underestimated by a factor of three. We propose a possible solution for this discrepancy.

We report the serendipitous discovery of a ring nebula around a candidate Wolf-Rayet (WR) star, HBHA 4202-22, in Cygnus using the Spitzer Space Telescope archival data. Our spectroscopic follow-up observations confirmed the WR nature of this star (we named it WR138a) and showed that it belongs to the WN8-9h subtype. We thereby add a new example to the known sample of late WN stars with circumstellar nebulae. We analyzed the spectrum of WR138a by using the Potsdam Wolf-Rayet (PoWR) model atmospheres, obtaining a stellar temperature of 40 kK. The stellar wind composition is dominated by helium with 20 per cent of hydrogen. The stellar spectrum is highly reddened and absorbed (E_{B-V} = 2.4 mag, A_V = 7.4 mag). Adopting a stellar luminosity of log L/Lsun = 5.3, the star has a mass-loss rate of 10^{-4.7} Msun yr^{-1}, and resides in a distance of 4.2 kpc. We measured the proper motion for WR138a and found that it is a runaway star with a peculiar velocity of ~50 km/s. Implications of the runaway nature of WR138a for constraining the mass of its progenitor star and understanding the origin of its ring nebula are discussed.

Analyzing the spectra from Wolf-Rayet stars requires adequate non-LTE modeling of their expanding atmosphere. The numerical schemes for solving the radiative transfer in the co-moving frame of reference have been developed by Mihalas and co-workers 30 years ago. The most elaborate codes can cope today with many hundred explicit non-LTE levels or super-levels and account for metal-line blanketing. The limited agreement with observed spectra indicates that the model simplifications are still severe. One approximation that has to be blamed is homogeneity. Stellar-wind clumping on small scales was easily implemented, while "macro-clumping" is still a big challenge. First studies showed that macro-clumping can reduce the strength of predicted P-Cygni line profiles in O-star spectra, and largely affects the X-ray line spectra from stellar winds. The classical model for radiation-driven winds by Castor, Abbot and Klein fails to explain the very dense winds from Wolf-Rayet stars. Only when we solved the detailed non-LTE radiative transfer consistently with the hydrodynamic equations, mass-loss rates above the single-scattering limit have been obtained.

We report the discovery of weak yet hard X-ray emission from the Wolf-Rayet (WR) star WR142 with the XMM-Newton X-ray telescope. Being of spectral subtype WO2, WR142 is a massive star in a very advanced evolutionary stage, short before its explosion as a supernova or gamma-ray burst. This is the first detection of X-ray emission from a WO-type star. We rule out any serendipitous X-ray sources within approx 1" of WR142. WR142 has an X-ray luminosity of L_X=7\times10^{30} erg/s, which constitutes only $\lsim 10^{-8}$ of its bolometric luminosity. The hard X-ray spectrum suggests a plasma temperature of about 100 MK. Commonly, X-ray emission from stellar winds is attributed to embedded shocks due to the intrinsic instability of the radiation driving. From qualitative considerations we conclude that this mechanism cannot account for the hardness of the observed radiation. There are no hints for a binary companion. Therefore the only remaining, albeit speculative explanation must refer to magnetic activity. Possibly related, WR142 seems to rotate extremely fast, as indicated by the unusually round profiles of its optical emission lines. Our detection implies that the wind of WR142 must be relatively transparent to X-rays, which can be due to strong wind ionization, wind clumping, or non-spherical geometry from rapid rotation.

Three very massive clusters are known to reside in the Galactic Center region, the Arches cluster, the Quintuplet cluster and the Central parsec cluster. We obtained spectroscopic observations of the Quintuplet cluster with the Integral Field Spectrograph SINFONI-SPIFFI at the ESO-VLT. The spectral range comprises the near-IR K-band from 1.94 to 2.45 micrometer. The 3D data cubes of the individual fields were flux-calibrated and combined to one contiguous cube, from which the spectra of all detectable point sources were extracted. We present a catalog of 160 stellar sources in the inner part of the Quintuplet cluster.

We report the discovery of variability in the X-ray emission from the Wolf-Rayet type star WR 65. Using archival Chandra data spanning over 5 yr we detect changes of the X-ray flux by a factor of 3 accompanied by changes in the X-ray spectra. We believe that this X-ray emission originates from wind-wind collision in a massive binary system. The observed changes can be explained by the variations in the emission measure of the hot plasma, and by the different absorption column along the binary orbit. The X-ray spectra of WR 65 display prominent emission features at wavelengths corresponding to the lines of strongly ionized Fe, Ca, Ar, S, Si, and Mg. WR 65 is a carbon rich WC9d star that is a persistent dust maker. This is the first investigation of any X-ray spectrum for a star of this spectral type. There are indications that the dust and the complex geometry of the colliding wind region are pivotal in explaining the X-ray properties of WR 65.

We study relatively isolated massive WN-type stars in the Galactic center. The K-band spectra of WR102ka and WR102c are exploited to infer the stellar parameters and to compute synthetic stellar spectra using the Potsdam Wolf-Rayet (PoWR) model atmosphere code. These models are combined with dust-shell models for analyzing the Spitzer IRS spectra of these objects. Archival IR images complement the interpretation. We report that WR102ka and WR102c are among the most luminous stars in the Milky Way. The mid-IR continua for both objects are dominated by dust emission. For the first time we report the presence of dust in the close vicinity of WN stars. Also for the first time, we have detected lines of pure-rotational transitions of molecular hydrogen in a massive-star nebula. A peony-shaped nebula around 102ka is resolved by the Spitzer MIPS camera. We attribute the formation of this nebula to the recent evolutionary history of WR102ka.

X-ray spectroscopy is a sensitive probe of stellar winds. X-rays originate from optically thin shock-heated plasma deep inside the wind and propagate outwards throughout absorbing cool material. Recent analyses of the line ratios from He-like ions in the X-ray spectra of O-stars highlighted problems with this general paradigm: the measured line ratios of highest ions are consistent with the location of the hottest X-ray emitting plasma very close to the base of the wind, perhaps indicating the presence of a corona, while measurements from lower ions conform with the wind-embedded shock model. Generally, to correctly model the emerging X-ray spectra, a detailed knowledge of the cool wind opacities based on stellar atmosphere models is prerequisite. A nearly grey stellar wind opacity for the X-rays is deduced from the analyses of high-resolution X-ray spectra. This indicates that the stellar winds are strongly clumped. Furthermore, the nearly symmetric shape of X-ray emission line profiles can be explained if the wind clumps are radially compressed. In massive binaries the orbital variations of X-ray emission allow to probe the opacity of the stellar wind; results support the picture of strong wind clumping. In high-mass X-ray binaries, the stochastic X-ray variability and the extend of the stellar-wind part photoionized by X-rays provide further strong evidence that stellar winds consist of dense clumps.

We report on observations of the eclipsing and interacting binary beta Lyrae from the Suzaku X-ray telescope. This system involves an early B star embedded in an optically and geometrically thick disk that is siphoning atmospheric gases from a less massive late B II companion. Motivated by an unpublished X-ray spectrum from the Einstein X-ray telescope suggesting unusually hard emission, we obtained time with Suzaku for pointings at three different phases within a single orbit. From the XIS detectors, the softer X-ray emission appears typical of an early-type star. What is surprising is the remarkably unchanging character of this emission, both in luminosity and in spectral shape, despite the highly asymmetric geometry of the system. We see no eclipse effect below 10 keV. The constancy of the soft emission is plausibly related to the wind of the embedded B star and Thomson scattering of X-rays in the system, although it might be due to extended shock structures arising near the accretion disk as a result of the unusually high mass-transfer rate. There is some evidence from the PIN instrument for hard emission in the 10-60 keV range. Follow-up observations with the RXTE satellite will confirm this preliminary detection.

Context: The mass-loss rate is a key parameter of massive stars. Adequate stellar atmosphere models are required for spectral analyses and mass-loss determinations. Present models can only account for the inhomogeneity of stellar winds in the approximation of small-scale structures that are optically thin. This treatment of ``microclumping'' has led to reducing empirical mass-loss rates by factors of two and more. Aims: Stellar wind clumps can be optically thick in spectral lines. We investigate how this ``macroclumping'' impacts on empirical mass-loss rates. Methods: The Potsdam Wolf-Rayet (PoWR) model atmosphere code is generalized in the ``formal integral'' to account for clumps that are not necessarily optically thin. Results: Optically thick clumps reduce the effective opacity. This has a pronounced effect on the emergent spectrum. Our modeling for the O-type supergiant zeta Puppis reveals that the optically thin H-alpha line is not affected by wind porosity, but that the PV resonance doublet becomes significantly weaker when macroclumping is taken into account. The reported discrepancies between resonance-line and recombination-line diagnostics can be resolved entirely with the macroclumping modeling without downward revision of the mass-loss rate. Conclusions: Mass-loss rates inferred from optically thin emission, such as the H-alpha line in O stars, are not influenced by macroclumping. The strength of optically thick lines, however, is reduced because of the porosity effects. Therefore, neglecting the porosity in stellar wind modeling can lead to underestimating empirical mass-loss rates.

Archival X-ray spectra of the four prominent single, non-magnetic O stars Zeta Pup, Zeta Ori, Ksi Per and Zeta Oph, obtained in high resolution with Chandra HETGS/MEG have been studied. The resolved X-ray emission line profiles provide information about the shocked, hot gas which emits the X-radiation, and about the bulk of comparably cool stellar wind material which partly absorbs this radiation. In this paper, we synthesize X-ray line profiles with a model of a clumpy stellar wind. We find that the geometrical shape of the wind inhomogeneities is important: better agreement with the observations can be achieved with radially compressed clumps than with spherical clumps. The parameters of the model, i.e. chemical abundances, stellar radius, mass-loss rate and terminal wind velocity, are taken from existing analyses of UV and optical spectra of the programme stars. On this basis, we also calculate the continuum-absorption coefficient of the cool-wind material, using the Potsdam Wolf-Rayet (PoWR) model atmosphere code. The radial location of X-ray emitting gas is restricted from analysing the fir line ratios of helium-like ions. The only remaining free parameter of our model is the typical distance between the clumps; here, we assume that at any point in the wind there is one clump passing by per one dynamical time-scale of the wind. The total emission in a model line is scaled to the observation. There is a good agreement between synthetic and observed line profiles. We conclude that the X-ray emission line profiles in O stars can be explained by hot plasma embedded in a cool wind which is highly clumped in the form of radially compressed shell fragments.

Radiative transfer in a clumped winds is used to describe X-ray emission line profiles observed in the XMM-Newton RGS spectrum of the OI star Zeta Puppis. It is shown that this X-ray spectrum can be explained as originating from a multi-temperature collisional plasma located in the wind acceleration zone. The X-rays are attenuated in the clumped stellar wind, which gives characteristic profiles to the emergent lines. We specifically study the N VII emission line in the spectrum of Zeta Pup. Long RGS exposures reveal that the N VII line profile is structured. On the basis of our Zeta Pup atmosphere model, we rule out the presence of N VII in the cool wind component. We suggest that the detailed N VII line structure is due to self-absorption in the hot plasma. Wind clumping also affects the transfer of ionizing radiation in high-mass X-ray binaries (HMXBs). We derive analytical formulae for the ionisation parameter in dependence on the parameters of wind clumping.

It is commonly adopted that X-rays from O stars are produced deep inside the stellar wind, and transported outwards through the bulk of the expanding matter which attenuates the radiation and affects the shape of emission line profiles. The ability of Chandra and XMM-Newton to resolve these lines spectroscopically provided a stringent test for the theory of X-ray production. It turned out that none of the existing models was able to reproduce the observations consistently. The major caveat of these models was the underlying assumption of a smooth stellar wind. Motivated by the various observational evidence that the stellar winds are in fact structured, we present a 2-D model of a stochastic, inhomogeneous wind. The X-ray radiative transfer is derived for such media. It is shown that profiles from a clumped wind differ drastically from those predicted by conventional homogeneous models. We review the up-to-date observations of X-ray line profiles from stellar winds and present line fits obtained from the inhomogeneous wind model. The necessity to account for inhomogeneities in calculating the X-ray transport in massive star winds, including for HMXB is highlighted.

The Galactic center harbors some of the most massive star clusters known in the Galaxy: the Arches and the Quintuplet. Based on the Chandra observations of these clusters (PI: Wang) which recently became public, I discuss the X-ray emission from the massive stars in these clusters. Confirming the general trend for Wolf-Rayet (WR) stars being X-ray dim, none of them is detected in the Quintuplet cluster. The most massive star known in the Galaxy, the Pistol star, is also not detected, invoking questions regarding the proposed binary nature of this object. X-ray emission in the Arches cluster is dominated by three stellar point sources. All three sources as well as the cluster's diffuse radiation show strong emission at 6.4-6.7 keV, indicating the presence of fluorescenting cool material. The Arches point sources may be identified as colliding wind binaries, albeit other possibilities cannot be ruled out.

The evolution of X-ray emission from young massive star clusters is modeled, taking into account the emission from the stars as well as from the cluster wind. It is shown that the level and character of the soft (0.2-10 keV) X-ray emission change drastically with cluster age and are tightly linked with stellar evolution. Using the modern X-ray observations of massive stars we show that the correlation between bolometric and X-ray luminosity known for single O stars also holds for O+O and O+Wolf-Rayet (WR) binaries. The diffuse emission originates from the cluster wind heated by the kinetic energy of stellar winds and supernova explosions. To model the evolution of the cluster wind, the mass and energy yields from a population synthesis are used as input to a hydrodynamic model. It is shown that in a very young clusters the emission from the cluster wind is low. When the cluster evolves, WR stars are formed. Their strong stellar winds power an increasing X-ray emission of the cluster wind. Subsequent supernova explosions pump the level of diffuse emission even higher. Clusters at this evolutionary stage may have no X-ray bright stellar point sources, but a relatively high level of diffuse emission. A supernova remnant may become a dominant X-ray source, but only for a short time interval of a few thousand years. We retrieve and analyse Chandra and XMM-Newton observations of six massive star clusters located in the Large Magellanic Cloud. Our model reproduces the observed diffuse and point-source emission from these LMC clusters, as well as from the Galactic clusters Arches, Quintuplet and NGC 3603.

It is commonly adopted that X-rays from O stars are produced deep inside the stellar wind, and transported outwards through the bulk of the expanding matter which attenuates the radiation and affects the shape of emission line profiles. None of the existing models was able to fit the observationed lines consistently. The possible caveat of these models was the assumption of a smooth stellar wind. Motivated by the various evidences that the stellar winds are in fact structured, we present a 2-D numerical model of a stochastic, inhomogeneous wind. Small parcels of hot, X-ray emitting gas are permeated with cool, absorbing wind material which is compressed into thin shell fragments. Wind fragmentation alters the radiative transfer drastically, compared to homogeneous models of the same mass-loss rate. X-rays produced deep inside the wind, which would be totally absorbed in a homogeneous flow, can effectively escape from a fragmented wind. The wind absorption becomes wavelength independent if the individual fragments are optically thick. The X-ray line profiles are flat-topped in the blue part and decline steeply in the red part for the winds with short acceleration zone. For the winds where the acceleration extends over significant distances, the lines can appear nearly symmetric and only slightly blueshifted, in contrast to the skewed, triangular line profiles typically obtained from homogeneous wind models of high optical depth. We show that profiles from a fragmented wind model can reproduce the observed line profiles from zeta Orionis.

The carbon-rich WC5 star WR 114 was not detected during a 15.9 ksec XMM-Newton observation, implying an upper limit to the X-ray luminosity of Lx < 2.5*10^{30} erg/s and to the X-ray to bolometric luminosity ratio of Lx/Lbol < 4*10^{-9}. This confirms indications from earlier less sensitive measurements that there has been no convincing X-ray detection of any single WC star. This lack of detections is reinforced by XMM-Newton and Chandra observations of WC stars. Thus the conclusion has to be drawn that the stars with radiatively-driven stellar winds of this particular class are insignificant X-ray sources. We attribute this to photoelectronic absorption by the stellar wind. The high opacity of the metal-rich and dense winds from WC stars puts the radius of optical depth unity at hundreds or thousands of stellar radii for much of the X-ray band. We believe that the essential absence of hot plasma so far out in the wind exacerbated by the large distances and correspondingly high ISM column densities makes the WC stars too faint to be detectable with current technology. The result also applies to many WC stars in binary systems, of which only about 20% are identified X-ray sources, presumably due to colliding winds.

We discuss X-ray line formation in dense O star winds. A random distribution of wind shocks is assumed to emit X-rays that are partially absorbed by cooler wind gas. The cool gas resides in highly compressed fragments oriented perpendicular to the radial flow direction. For fully opaque fragments, we find that the blueshifted part of X-ray line profiles remains flat-topped even after severe wind attenuation, whereas the red part shows a steep decline. These box-type, blueshifted profiles resemble recent Chandra observations of the O3 star zeta Pup. For partially transparent fragments, the emission lines become similar to those from a homogeneous wind.

We describe X-ray production in the atmospheres of hot, early-type stars in the framework of a ``stochastic shock model''. The extended envelope of a star is assumed to possess numerous X-ray emitting ``hot'' zones that are produced by shocks and embedded in the ambient ``cold'' medium in dynamical equilibrium. It is shown that the apparent lack of X-ray variability on short (order of hours) timescales do not contradict a shock model for X-ray production. The character of the X-ray variability is found to depend on the frequency with which hot zones are generated, the cool wind opacity to X-rays, and the wind flow parameters, such as mass loss rate and terminal speed.