The intriguing binary LS V +22 25 (LB-1) has drawn much attention following claims of it being a single-lined spectroscopic binary with a 79-day orbit comprising a B-type star and a ~70 Msun black hole. Recent analyses have implied that the visible primary star is a stripped He-rich star. However, the nature of the secondary, which was proposed to be a black hole, a neutron star, or a main sequence star, remains unknown. Based on 26 newly acquired spectroscopic observations, we perform an orbital analysis and spectral disentangling of LB-1 to elucidate the nature of the system. Our analysis reveals that LB-1 contains two components of comparable brightness in the optical. The narrow-lined primary, which we estimate to contribute ~55% in the optical, has spectral properties that suggest that it is a stripped star: it has a small spectroscopic mass (~1 Msun) for a B-type star and it is He- and N-rich. The "hidden" secondary, which contributes about 45% of the optical flux, is a rapidly rotating (vsini ~ 300 km/s) B3 V star with a decretion disk -- a Be star. Hence, LB-1 does not contain a compact object. Instead, it is a rare Be binary system consisting of a stripped star (the former mass donor) and a Be star rotating at near its critical velocity (the former mass accretor). This system is a clear example that binary interactions play a decisive role in the production of rapid stellar rotators and Be stars.

It is now well established that the majority of massive stars reside in multiple systems. However, the effect of multiplicity is not sufficiently understood, resulting in a plethora of uncertainties about the end stages of massive-star evolution. In order to investigate these uncertainties, it is useful to study massive stars just before their demise. Classical Wolf-Rayet stars represent the final end stages of stars at the upper-mass end. The multiplicity fraction of these stars was reported to be ${\sim}0.4$ in the Galaxy but no correction for observational biases has been attempted. The aim of this study is to conduct a homogeneous radial-velocity survey of a magnitude-limited ($V$ $\leq 12$) sample of Galactic Wolf-Rayet stars to derive their bias-corrected multiplicity properties. The present paper focuses on 12 northern Galactic carbon-rich (WC) Wolf-Rayet stars observable with the 1.2m Mercator telescope on the island of La Palma. We homogeneously measured relative radial velocities (RVs) for carbon-rich Wolf-Rayet stars using cross-correlation. Variations in the derived RVs were used to flag binary candidates. We investigated probable orbital configurations and provide a first correction of observational biases through Monte-Carlo simulations. Of the 12 northern Galactic WC stars in our sample, seven show peak-to-peak RV variations larger than 10 km s$^{-1}$, which we adopt as our detection threshold. This results in an observed spectroscopic multiplicity fraction of 0.58 with a binomial error of 0.14. In our campaign, we find a clear lack of short-period (P~$<~\sim$100\,d), indicating that a large number of Galactic WC binaries likely reside in long-period systems. Finally, our simulations show that at the 10% significance level, the intrinsic multiplicity fraction of the Galactic WC population is at least 0.72.

About 20% of all B-type stars are classical Be stars. The Be phenomenon is strongly correlated with rapid rotation, the origin of which remains unclear. It may be rooted in single- or binary-star evolution. In the framework of the binary channel, the initially more massive star transfers mass and angular momentum to the original secondary, which becomes a Be star. The system then evolves into a Be binary with a post-main-sequence companion, which may later be disrupted in a supernova event. Hence, if the binary channel dominates the formation of Be stars, one may expect a strong lack of close Be binaries with main sequence (MS) companions. Through an extensive, star-by-star review of the literature of a magnitude-limited sample of Galactic early-type Be stars, we investigate whether Be binaries with MS companions are known to exist. Our sample is constructed from the BeSS database and cross-matched with all available literature on the individual stars. Out of an initial list of 505 Be stars, we compile a final sample of 287 Galactic Be stars earlier than B1.5 with V<=12 mag. Out of those, 13 objects were reported as Be binaries with known post-MS companions and 11 as binaries with unknown, uncertain or debated companions. We find no confirmed reports of Be binaries with MS companions. For the remaining 263 targets, no significant reports of multiplicity exist in the literature, implying that they are either Be binaries with faint companions, or truly single. The clear lack of reported MS companions to Be stars, which stands in contrast to the high number of detected B+B MS binaries, strongly supports the hypothesis that early-type Be stars are binary interaction products that spun up after mass and angular momentum transfer from a companion star. Taken at face value, our results may suggest that a large majority of the early-type Be stars have formed through binary mass-transfer.

HR 6819 was recently proposed to be a triple system consisting of an inner B-type giant + black hole binary with an orbital period of 40d and an outer Be tertiary. This interpretation is mainly based on two inferences: that the emission attributed to the outer Be star is stationary, and that the inner star, which is used as mass calibrator for the black hole, is a B-type giant. We re-investigate the properties of HR 6819 by spectral disentangling and an atmosphere analysis of the disentangled spectra to search for a possibly simpler alternative explanation for HR 6819. Disentangling implies that the Be component is not a static tertiary, but rather a component of the binary in the 40-d orbit. The inferred radial velocity amplitudes imply an extreme mass ratio of M_2/M_1 = 15 +/- 3. We infer spectroscopic masses of 0.4$^{+0.3}_{-0.1}$ Msun and 6$^{+5}_{-3}$ Msun for the primary and secondary, which agree well with the dynamical masses for an inclination of i = 32 deg. This indicates that the primary might be a stripped star rather than a B-type giant. Evolutionary modelling suggests that a possible progenitor system would be a tight (P_i ~ 2d) B+B binary system that experienced conservative mass transfer. While the observed nitrogen enrichment of the primary conforms with the predictions of the evolutionary models, we find no indications for the predicted He enrichment. We suggest that HR 6819 is a binary system consisting of a stripped B-type primary and a rapidly-rotating Be star that formed from a previous mass-transfer event. In the framework of this interpretation, HR 6819 does not contain a black hole. Interferometry can distinguish between these two scenarios by providing an independent measurement of the separation between the visible components.

The recent gravitational wave measurements have demonstrated the existence of stellar mass black hole binaries. It is essential for our understanding of massive star evolution to identify the contribution of binary evolution to the formation of double black holes. A promising way to progress is investigating the progenitors of double black hole systems and comparing predictions with local massive star samples such as the population in 30 Doradus in the Large Magellanic Cloud (LMC). Methods. To this purpose, we analyse a large grid of detailed binary evolution models at LMC metallicity with initial primary masses between 10 and 40 Msun, and identify which model systems potentially evolve into a binary consisting of a black hole and a massive main sequence star. We then derive the observable properties of such systems, as well as peculiarities of the OB star component. We find that about 3% of the LMC late O and early B stars in binaries are expected to possess a black hole companion, when assuming stars with a final helium core mass above 6.6 M to form black holes. While the vast majority of them may be X-ray quiet, our models suggest that these may be identified in spectroscopic binaries, either by large amplitude radial velocity variations ( > 50 km s ) and simultaneous nitrogen surface enrichment, or through a moderate radial velocity ( > 10 km/s ) and simultaneously rapid rotation of the OB star. The predicted mass ratios are such that main sequence companions could be excluded in most cases. A comparison to the observed OB+WR binaries in the LMC, Be/X-ray binaries, and known massive BH binaries supports our conclusion. We expect spectroscopic observations to be able to test key assumptions in our models, with important implications for massive star evolution in general, and for the formation of double-black hole mergers in particular.

Classical Wolf-Rayet (WR) stars are massive, hydrogen depleted, post main-sequence stars that exhibit emission-line dominated spectra. For a given metallicity Z, stars exceeding a certain initial mass M_single(Z) can reach the WR phase through intrinsic mass-loss (single-star channel). Stars of lower masses can reach the WR phase via binary mass transfer (binary channel). It is commonly assumed that the binary channel dominates the formation of WR stars in environments with low Z such as the SMC and LMC. However, their reported WR binary fractions of 30-40% are comparable to that of the Galaxy. Here, we explain this apparent contradiction by considering the minimum initial mass M_spec(Z) needed for the stripped product to appear as a WR star. We calibrate M_spec(Z) using the lowest-luminosity WR stars in the Clouds and the Galaxy. A range of M_single(Z) values are explored using various evolution codes. We estimate the additional contribution of the binary channel by considering the interval [M_spec(Z), M_single(Z)], which characterises the initial-mass range in which binaries can form additional WR stars. Results: The WR-phenomenon ceases below luminosities of logL = 4.9, 5.25, and 5.6 [Lsun] in the Galaxy, LMC, and SMC, which translates to He-star masses of 7.5, 11, 17 Msun and initial masses of M_spec = 18, 23, 37 Msun. Stripped stars with lower initial masses in the respective galaxies would tend to not appear as WR stars. M_single lies in the range 20-30, 30-60, and > 40 Msun for the Galaxy, LMC, and SMC. We find that that the additional contribution of the binary channel is a non-trivial function of Z that cannot be conclusively claimed to be monotonically increasing with decreasing Z. Hence, one should not a-priori expect that binary interactions become increasingly important in forming WR stars at low Z, or that the WR binary fraction grows with decreasing Z.

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.

Accurate stellar parameters of individual objects in binary systems are essential to constrain the effects of binarity on stellar evolution. These parameters serve as a prerequisite to probing existing and future theoretical evolutionary models. We aim to derive the atmospheric parameters of the 31 SB2s in the TMBM sample. This sample, composed of detached, semi-detached and contact systems with at least one of the components classified as an O star, is an excellent test-bed to study how binarity can impact our knowledge of the evolution of massive stars. 32 epochs of FLAMES/GIRAFFE spectra are analysed using spectral disentangling to construct the individual spectra of 62 components. We apply the CMFGEN atmosphere code to determine their stellar parameters and their He, C and N surface abundances. From these properties, we show that the effects of tides on chemical mixing are limited. Components on longer-period orbits show higher nitrogen enrichment at their surface than those on shorter-period orbits, in contrast to expectations of rotational or tidal mixing, implying that other mechanisms play a role in this process. Components filling their Roche lobe are mass donors. They exhibit higher nitrogen content at their surface and rotate more slowly than their companions. By accreting new material, their companions spin faster and are rejuvenated. Their locations in the N-vsini diagram tend to show that binary products are good candidates to populate the two groups of stars (slowly rotating, nitrogen-enriched and rapidly rotating non-enriched) that cannot be reproduced through single-star population synthesis. This sample is the largest sample of binaries to be studied in such a homogeneous way. The study of these objects gives us strong observational constraints to test theoretical binary evolutionary tracks.

A high fraction of massive stars are found to be binaries but only a few of them are reported as photometrically variable. By studying the populations of SB2 in the 30 Doradus region, we found a subset of them that have photometry from the OGLE project and that display variations in their light curves related to orbital motions. The goal of this study is to determine the dynamical masses and radii of the 26 binary components to investigate the mass-discrepancy problem and to provide an empirical mass-luminosity relation for the LMC. We use the PHOEBE programme to perform a systematic analysis of the OGLE V and I light curves obtained for 13 binary systems in 30 Dor. We adopt Teff, and orbital parameters derived previously to obtain the inclinations of the systems and the parameters of the individual components. Three systems display eclipses in their light curves, while the others only display ellipsoidal variations. We classify two systems as over-contact, five as semi-detached, and four as detached. The two remaining systems have uncertain configurations due to large uncertainties on their inclinations. The fact that systems display ellipsoidal variations has a significant impact on the inclination errors. From the dynamical masses, luminosities, and radii, we provide LMC-based empirical mass-luminosity and mass-radius relations, and we compare them to other relations given for the Galaxy, the LMC, and the SMC. These relations differ for different mass ranges, but do not seem to depend on the metallicity regimes. We also compare the dynamical, spectroscopic, and evolutionary masses of the stars in our sample. While the dynamical and spectroscopic masses agree with each other, the evolutionary masses are systematically higher, at least for stars in semi-detached systems. This suggests that the mass discrepancy can be partly explained by past or ongoing interactions between the stars.

Liu et al. (2019) recently reported the detection of a 68 [+11/-13] solar mass (Msun) black hole (BH) paired with an 8.2 [+0.9/-1.2] Msun B-type sub-giant star in the 78.9-day spectroscopic binary system LB-1. Such a black hole is over twice as massive as any other known stellar-mass black hole with non-compact companions2 and its mass approaches those that result from BH-BH coalescences that are detected by gravitational wave interferometers. Its presence in a solar-like metallicity environment challenges conventional theories of massive binary evolution, stellar winds and core-collapse supernovae, so that more exotic scenarios seem to be needed to explain the existence and properties of LB-1. Here, we show that the observational diagnostics used to derive the BH mass results from the orbital motion of the B-type star, not that of the BH. As a consequence, no evidence for a massive BH remains in the data, therefore solving the existing tension with formation models of such a massive BH at solar metallicity and with theories of massive star evolution in general.

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.

The massive O4.5 V + O5.5 V binary VFTS 352 in the Tarantula nebula is one of the shortest-period and most massive overcontact binaries known. Recent theoretical studies indicate that some of these systems could ultimately lead to the formation of gravitational waves via black hole binary mergers through the chemically homogeneous evolution pathway. By analyzing ultraviolet-optical phase-resolved spectroscopic data, we aim to constrain atmospheric and wind properties that could be later used to confront theoretical predictions from binary evolution. In particular, surface abundances are powerful diagnostics of the evolutionary status, mass transfer and the internal mixing processes. From a set of 32 VLT/FLAMES visual and 8 HST/COS ultraviolet spectra, we used spectral disentangling to separate the primary and secondary components. Using a genetic algorithm wrapped around the NLTE model atmosphere and spectral synthesis code FASTWIND, we perform an 11-parameter optimization to derive the atmospheric and wind parameters of both components, including the surface abundances of He, C, N, O and Si. We find that both components are hotter than expected compared to single-star evolutionary models indicating that additional mixing processes may be at play. However the derived chemical abundances do not show significant indications of mixing when adopting baseline values typical for the system environment.

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.

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.

Context. Metal-poor massive stars are supposed to be progenitors of certain supernovae, gamma-ray bursts and compact object mergers, potentially contributing to the early epochs of the Universe with their strong ionizing radiation. However, they remain mainly theoretical as individual spectroscopic observations of such objects have rarely been carried out below the metallicity of the SMC. Aims. This work aims at exploring what our state-of-the-art theories of stellar evolution combined with those of stellar atmospheres predict about a certain type of metal-poor (0.02 Z$_{\odot}$) hot massive stars, the chemically homogeneously evolving ones, called TWUIN stars. Methods. Synthetic spectra corresponding to a broad range in masses (20-130 M$_{\odot}$) and covering several evolutionary phases from the zero-age main-sequence up to the core helium-burning stage were computed. Results. We find that TWUIN stars show almost no emission lines during most of their {core hydrogen-burning} lifetimes. Most metal lines are completely absent, including nitrogen. During their core helium-burning stage, lines switch to emission and even some metal lines (oxygen and carbon, but still almost no nitrogen) show up. Mass loss and clumping play a significant role in line-formation in later evolutionary phases, particularly during core helium-burning. Most of our spectra are classified as an early O type giant or supergiant, and we find Wolf-Rayet stars of type WO in the core helium-burning phase. Conclusions. An extremely hot, early O type star observed in a low-metallicity galaxy could be the outcome of chemically homogeneous evolution $-$ and therefore the progenitor of a long-duration gamma-ray burst or a type Ic supernova. TWUIN stars may play an important role in reionizing the Universe due to their being hot without showing prominent emission lines during the majority of their lifetimes.

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.

Stochastically-triggered photospheric light variations reaching $\sim$$40$ mmag peak-to-valley amplitudes have been detected in the O8Iaf supergiant V973 Scorpii as the outcome of two months of high-precision time-resolved photometric observations with the BRIght Target Explorer (BRITE) nanosatellites. The amplitude spectrum of the time series photometry exhibits a pronounced broad bump in the low-frequency regime ($\lesssim$$0.9$ d$^{-1}$) where several prominent frequencies are detected. A time-frequency analysis of the observations reveals typical mode lifetimes of the order of $5-10$ days. The overall features of the observed brightness amplitude spectrum of V973 Sco match well with those extrapolated from two-dimensional hydrodynamical simulations of convectively-driven internal gravity waves randomly excited from deep in the convective cores of massive stars. An alternative or additional possible source of excitation from a subsurface convection zone needs to be explored in future theoretical investigations.

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.

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.

From $5.5$ months of dual-band optical photometric monitoring at the $1$ mmag level, BRITE-Constellation has revealed two simultaneous types of variability in the O4I(n)fp star $\zeta$ Puppis: one single periodic non-sinusoidal component superimposed on a stochastic component. The monoperiodic component is the $1.78$ d signal previously detected by Coriolis/SMEI, but this time along with a prominent first harmonic. The shape of this signal changes over time, a behaviour that is incompatible with stellar oscillations but consistent with rotational modulation arising from evolving bright surface inhomogeneities. By means of a constrained non-linear light curve inversion algorithm we mapped the locations of the bright surface spots and traced their evolution. Our simultaneous ground-based multi-site spectroscopic monitoring of the star unveiled cyclical modulation of its He II $\lambda4686$ wind emission line with the $1.78$-day rotation period, showing signatures of Corotating Interaction Regions (CIRs) that turn out to be driven by the bright photospheric spots observed by BRITE. Traces of wind clumps are also observed in the He II $\lambda4686$ line and are correlated with the amplitudes of the stochastic component of the light variations probed by BRITE at the photosphere, suggesting that the BRITE observations additionally unveiled the photospheric drivers of wind clumps in $\zeta$ Pup and that the clumping phenomenon starts at the very base of the wind. The origins of both the bright surface inhomogeneities and the stochastic light variations remain unknown, but a subsurface convective zone might play an important role in the generation of these two types of photospheric variability.

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.

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.

We report on the first multi-color precision light curve of the bright Wolf-Rayet binary $\gamma^2$ Velorum, obtained over six months with the nanosatellites in the BRITE- Constellation fleet. In parallel, we obtained 488 high-resolution optical spectra of the system. In this first report on the datasets, we revise the spectroscopic orbit and report on the bulk properties of the colliding winds. We find a dependence of both the light curve and excess emission properties that scales with the inverse of the binary separation. When analyzing the spectroscopic properties in combination with the photometry, we find that the phase dependence is caused only by excess emission in the lines, and not from a changing continuum. We also detect a narrow, high-velocity absorption component from the He I $\lambda$5876 transition, which appears twice in the orbit. We calculate smoothed-particle hydrodynamical simulations of the colliding winds and can accurately associate the absorption from He I to the leading and trailing arms of the wind shock cone passing tangentially through our line of sight. The simulations also explain the general strength and kinematics of the emission excess observed in wind lines such as C III $\lambda$5696 of the system. These results represent the first in a series of investigations into the winds and properties of $\gamma^2$ Velorum through multi-technique and multi-wavelength observational campaigns.

CONTEXT: For more than two decades, stellar atmosphere codes have been used to derive the stellar and wind parameters of massive stars. Although they have become a powerful tool and sufficiently reproduce the observed spectral appearance, they can hardly be used for more than measuring parameters. One major obstacle is their inconsistency between the calculated radiation field and the wind stratification due to the usage of prescribed mass-loss rates and wind-velocity fields. AIMS: We present the concepts for a new generation of hydrodynamically consistent non-local thermodynamical equilibrium (non-LTE) stellar atmosphere models that allow for detailed studies of radiation-driven stellar winds. As a first demonstration, this new kind of model is applied to a massive O star. METHODS: Based on earlier works, the PoWR code has been extended with the option to consistently solve the hydrodynamic equation together with the statistical equations and the radiative transfer in order to obtain a hydrodynamically consistent atmosphere stratification. In these models, the whole velocity field is iteratively updated together with an adjustment of the mass-loss rate. RESULTS: The concepts for obtaining hydrodynamically consistent models using a comoving-frame radiative transfer are outlined. To provide a useful benchmark, we present a demonstration model, which was motivated to describe the well-studied O4 supergiant Zeta Pup. The obtained stellar and wind parameters are within the current range of literature values. CONCLUSIONS: For the first time, the PoWR code has been used to obtain a hydrodynamically consistent model for a massive O star. This has been achieved by a profound revision of earlier concepts used for Wolf-Rayet stars. The velocity field is shaped by various elements contributing to the radiative acceleration, especially in the outer wind. (...)

Massive binaries (MBs) play a crucial role in the Universe. Knowing the distributions of their orbital parameters (OPs) is important for a wide range of topics, from stellar feedback to binary evolution channels, from the distribution of supernova types to gravitational wave progenitors, yet, no direct measurements exist outside the Milky Way. The Tarantula Massive Binary Monitoring was designed to help fill this gap by obtaining multi-epoch radial velocity monitoring of 102 MBs in the 30 Dor. In this paper, we analyse 32 VLT/FLAMES observations of 93 O- and 7 B-type binaries. We performed a Fourier analysis and obtained orbital solutions for 82 systems: 51 single- and 31 double-lined spectroscopic binaries. Overall, the OPs and binary fraction are remarkably similar across the 30 Dor region and compared to existing Galactic samples (GSs). This indicates that within these domains environmental effects are of second order in shaping the properties of MBs. A small difference is found in the distribution of orbital periods (OrbPs), which is slightly flatter (in log space) in 30 Dor than in the Galaxy, although this may be compatible within error estimates and differences in the fitting methodology. Also, OrbPs in 30 Dor can be as short as 1.1 d; somewhat shorter than seen in GSs. Equal mass binaries q>0.95 in 30 Dor are all found outside NGC 2070 the very young and massive cluster at 30 Dor's core. One outstanding exception however is the fact that earliest spectral types (O2-O7) tend to have shorter OrbPs than latter (O9.2-O9.7). Our results point to a relative universality of the incidence rate of MBs and their OPs in the metallicity range from solar ($Z_{\odot}$) to about $0.5Z_{\odot}$. This provides the first direct constraints on MB properties in massive star-forming galaxies at the Universes peak of star formation at redshifts z~1 to 2, which are estimated to have $Z~0.5Z_{\odot}$.

Wolf-Rayet (WR) stars are the most advanced stage in the evolution of the most massive stars. The strong feedback provided by these objects and their subsequent supernova (SN) explosions are decisive for a variety of astrophysical topics such as the cosmic matter cycle. Consequently, understanding the properties of WR stars and their evolution is indispensable. A crucial but still not well known quantity determining the evolution of WR stars is their mass-loss rate. Since the mass loss is predicted to increase with metallicity, the feedback provided by these objects and their spectral appearance are expected to be a function of the metal content of their host galaxy. This has severe implications for the role of massive stars in general and the exploration of low metallicity environments in particular. Hitherto, the metallicity dependence of WR star winds was not well studied. In this contribution, we review the results from our comprehensive spectral analyses of WR stars in environments of different metallicities, ranging from slightly super-solar to SMC-like metallicities. Based on these studies, we derived empirical relations for the dependence of the WN mass-loss rates on the metallicity and iron abundance, respectively.

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.

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.

WR 148 (HD 197406) is an extreme runaway system considered to be a potential candidate for a short-period (4.3173 d) rare WR + compact object binary. Provided with new high resolution, high signal-to-noise spectra from the Keck observatory, we determine the orbital parameters for both the primary WR and the secondary, yielding respective projected orbital velocity amplitudes of $88.1\pm3.8$ km s$^{-1}$ and $79.2\pm3.1$ km s$^{-1}$ and implying a mass ratio of $1.1\pm0.1$. We then apply the shift-and-add technique to disentangle the spectra and obtain spectra compatible with a WN7ha and an O4-6 star. Considering an orbital inclination of $\sim67^\circ$, derived from previous polarimetry observations, the system's total mass would be a mere 2-3 M$_{\odot}$ , an unprecedented result for a putative massive binary system. However, a system comprising a 37 M$_{\odot}$ secondary (typical mass of an O5V star) and a 33 M$_{\odot}$ primary (given the mass ratio) would infer an inclination of $\sim18^\circ$. We therefore reconsider the previous methods of deriving the orbital inclination based on time-dependent polarimetry and photometry. While the polarimetric results are inconclusive requiring better data, the photometric results favour low inclinations. Finally, we compute WR 148's space velocity and retrace the runaway's trajectory back to the Galactic plane (GP). With an ejection velocity of $198\pm27$ km s$^{-1}$ and a travel time of $4.7\pm0.8$ Myr to reach its current location, WR 148 was most likely ejected via dynamical interactions in a young cluster.

We report on interferometric observations with the CHARA Array of two classical Wolf-Rayet stars in suspected binary systems, namely WR 137 and WR 138. In both cases, we resolve the component stars to be separated by a few milliarcseconds. The data were collected in the H-band, and provide a measure of the fractional flux for both stars in each system. We find that the WR star is the dominant H-band light source in both systems ($f_{\rm WR, 137} = 0.59\pm0.04$; $f_{\rm WR, 138} = 0.67\pm0.01$), which is confirmed through both comparisons with estimated fundamental parameters for WR stars and O dwarfs, as well as through spectral modeling of each system. Our spectral modeling also provides fundamental parameters for the stars and winds in these systems. The results on WR 138 provide evidence that it is a binary system which may have gone through a previous mass-transfer episode to create the WR star. The separation and position of the stars in the WR 137 system together with previous results from the IOTA interferometer provides evidence that the binary is seen nearly edge-on. The possible edge-on orbit of WR 137 aligns well with the dust production site imaged by the Hubble Space Telescope during a previous periastron passage, showing that the dust production may be concentrated in the orbital plane.

During the summer of 2013, a 4-month spectroscopic campaign took place to observe the variabilities in three Wolf-Rayet stars. The spectroscopic data have been analyzed for WR 134 (WN6b), to better understand its behaviour and long-term periodicity, which we interpret as arising from corotating interaction regions (CIRs) in the wind. By analyzing the variability of the He II $\lambda$5411 emission line, the previously identified period was refined to P = 2.255 $\pm$ 0.008 (s.d.) days. The coherency time of the variability, which we associate with the lifetime of the CIRs in the wind, was deduced to be 40 $\pm$ 6 days, or $\sim$ 18 cycles, by cross-correlating the variability patterns as a function of time. When comparing the phased observational grayscale difference images with theoretical grayscales previously calculated from models including CIRs in an optically thin stellar wind, we find that two CIRs were likely present. A separation in longitude of $\Delta \phi \simeq$ 90$^{\circ}$ was determined between the two CIRs and we suggest that the different maximum velocities that they reach indicate that they emerge from different latitudes. We have also been able to detect observational signatures of the CIRs in other spectral lines (C IV $\lambda\lambda$5802,5812 and He I $\lambda$5876). Furthermore, a DAC was found to be present simultaneously with the CIR signatures detected in the He I $\lambda$5876 emission line which is consistent with the proposed geometry of the large-scale structures in the wind. Small-scale structures also show a presence in the wind, simultaneously with the larger scale structures, showing that they do in fact co-exist.

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.

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.

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.

Wolf-Rayet (WR) stars have a severe impact on their environments owing to their strong ionizing radiation fields and powerful stellar winds. Since these winds are considered to be driven by radiation pressure, it is theoretically expected that the degree of the wind mass-loss depends on the initial metallicity of WR stars. Following our comprehensive studies of WR stars in the Milky Way, M31, and the LMC, we derive stellar parameters and mass-loss rates for all seven putatively single WN stars known in the SMC. Based on these data, we discuss the impact of a low-metallicity environment on the mass loss and evolution of WR stars. The quantitative analysis of the WN stars is performed with the Potsdam Wolf-Rayet (PoWR) model atmosphere code. The physical properties of our program stars are obtained from fitting synthetic spectra to multi-band observations. In all SMC WN stars, a considerable surface hydrogen abundance is detectable. The majority of these objects have stellar temperatures exceeding 75 kK, while their luminosities range from 10^5.5 to 10^6.1 Lsun. The WN stars in the SMC exhibit on average lower mass-loss rates and weaker winds than their counterparts in the Milky Way, M31, and the LMC. By comparing the mass-loss rates derived for WN stars in different Local Group galaxies, we conclude that a clear dependence of the wind mass-loss on the initial metallicity is evident, supporting the current paradigm that WR winds are driven by radiation. A metallicity effect on the evolution of massive stars is obvious from the HRD positions of the SMC WN stars at high temperatures and high luminosities. Standard evolution tracks are not able to reproduce these parameters and the observed surface hydrogen abundances. Homogeneous evolution might provide a better explanation for their evolutionary past.

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.

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.

CONTEXT: Spectroscopic analysis remains the most common method to derive masses of massive stars, the most fundamental stellar parameter. While binary orbits and stellar pulsations can provide much sharper constraints on the stellar mass, these methods are only rarely applicable to massive stars. Unfortunately, spectroscopic masses of massive stars heavily depend on the detailed physics of model atmospheres. AIMS: We demonstrate the impact of a consistent treatment of the radiative pressure on inferred gravities and spectroscopic masses of massive stars. Specifically, we investigate the contribution of line and continuum transitions to the photospheric radiative pressure. We further explore the effect of model parameters, e.g., abundances, on the deduced spectroscopic mass. Lastly, we compare our results with the plane-parallel TLUSTY code, commonly used for the analysis of massive stars with photospheric spectra. METHODS: We calculate a small set of O-star models with the Potsdam Wolf-Rayet (PoWR) code using different approaches for the quasi-hydrostatic part. These models allow us to quantify the effect of accounting for the radiative pressure consistently. We further use PoWR models to show how the Doppler widths of line profiles and abundances of elements such as iron affect the radiative pressure, and, as a consequence, the derived spectroscopic masses. RESULTS: Our study implies that errors on the order of a factor of two in the inferred spectroscopic mass are to be expected when neglecting the contribution of line and continuum transitions to the radiative acceleration in the photosphere. Usage of implausible microturbulent velocities, or the neglect of important opacity sources such as Fe, may result in errors of approximately 50% in the spectroscopic mass. A comparison with TLUSTY model atmospheres reveals a very good agreement with PoWR at the limit of low mass-loss rates.

We report the first-ever discovery of a Wolf-Rayet (WR) star in the Large Magellanic Cloud via detection of a circular shell with the Spitzer Space Telescope. Follow-up observations with Gemini-South resolved the central star of the shell into two components separated from each other by approx 2 arcsec (or approx 0.5 pc in projection). One of these components turns out to be a WN3 star with H and He lines both in emission and absorption (we named it BAT99 3a using the numbering system based on extending the Breysacher et al. catalogue). Spectroscopy of the second component showed that it is a B0 V star. Subsequent spectroscopic observations of BAT99 3a with the du Pont 2.5-m telescope and the Southern African Large Telescope revealed that it is a close, eccentric binary system, and that the absorption lines are associated with an O companion star. We analyzed the spectrum of the binary system using the non-LTE Potsdam Wolf-Rayet (PoWR) code, confirming that the WR component is a very hot (approx 90 kK) WN star. For this star, we derived a luminosity of log L/Lsun =5.45 and a mass-loss rate of 10^{-5.8} Msun/yr, and found that the stellar wind composition is dominated by helium with 20 per cent of hydrogen. Spectroscopy of the shell revealed an He iii region centred on BAT99 3a and having the same angular radius (approx 15 arcsec) as the shell. We thereby add a new example to a rare class of high-excitation nebulae photoionized by WR stars. Analysis of the nebular spectrum showed that the shell is composed of unprocessed material, implying that the shell was swept-up from the local interstellar medium. We discuss the physical relationship between the newly identified massive stars and their possible membership of a previously unrecognized star cluster.

Massive Wolf-Rayet stars are recognized today to be in a very common, but short, evolutionary phase of massive stars. While our understanding of Wolf-Rayet stars has increased dramatically over the past decades, it remains unclear whether rapid rotators are among them. There are various indications that rapidly rotating Wolf-Rayet stars should exist. Unfortunately, due to their expanding atmospheres, rotational velocities of Wolf-Rayet stars are very difficult to measure. However, recently observed spectra of several Wolf-Rayet stars reveal peculiarly broad and round emission lines. Could these spectra imply rapid rotation? In this work, we model the effects of rotation on the atmospheres of Wolf-Rayet stars. We further investigate whether the peculiar spectra of five Wolf-Rayet stars may be explained with the help of stellar rotation, infer appropriate rotation parameters, and discuss the implications of our results. We make use of the Potsdam Wolf-Rayet (PoWR) non-LTE model atmosphere code. Since the observed spectra of Wolf-Rayet stars are mainly formed in their expanding atmospheres, rotation must be accounted for with a 3D integration scheme of the formal integral. For this purpose, we assume a rotational velocity field consisting of an inner co-rotating domain and an outer domain, where the angular momentum is conserved. We find that rotation can reproduce the unique spectra analyzed here. However, the inferred rotational velocities at the stellar surface are large (~200 km/s), and the inferred co-rotation radii (~10 stellar radii) suggest the existence of very strong photospheric magnetic fields (~20 kG).