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.

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.

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.

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.

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.

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.

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

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.

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.

We present spectral classifications from optical spectroscopy of 263 massive stars in the north-eastern region of the Large Magellanic Cloud. The observed two-degree field includes the massive 30 Doradus star-forming region, the environs of SN1987A, and a number of star-forming complexes to the south of 30 Dor. These are the first classifications for the majority (203) of the stars and include eleven double-lined spectroscopic binaries. The sample also includes the first examples of early OC-type spectra (AAOmega 30 Dor 248 and 280), distinguished by the weakness of their nitrogen spectra and by C IV 4658 emission. We propose that these stars have relatively unprocessed CNO abundances compared to morphologically normal O-type stars, indicative of an earlier evolutionary phase. From analysis of observations obtained on two consecutive nights, we present radial-velocity estimates for 233 stars, finding one apparent single-lined binary and nine (>3sigma) outliers compared to the systemic velocity; the latter objects could be runaway stars or large-amplitude binary systems and further spectroscopy is required to investigate their nature.

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.

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.

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

Context: Comprehensive studies of Wolf-Rayet stars were performed in the past for the Galactic and the LMC population. The results revealed significant differences, but also unexpected similarities between the WR populations of these different galaxies. Analyzing the WR stars in M31 will extend our understanding of these objects in different galactic environments. Aims: The present study aims at the late-type WN stars in M31. The stellar and wind parameters will tell about the formation of WR stars in other galaxies with different metallicity and star formation histories. The obtained parameters will provide constraints to the evolution of massive stars in the environment of M31. Methods: We used the latest version of the Potsdam Wolf-Rayet model atmosphere code to analyze the stars via fitting optical spectra and photometric data. To account for the relatively low temperatures of the late WN10 and WN11 subtypes, our WN models have been extended into this temperature regime. Results: Stellar and atmospheric parameters are derived for all known late-type WN stars in M31 with available spectra. All of these stars still have hydrogen in their outer envelopes, some of them up to 50% by mass. The stars are located on the cool side of the zero age main sequence in the Hertzsprung-Russell diagram, while their luminosities range from $10^5$ to $10^6$ Lsun. It is remarkable that no star exceeds $10^6$ Lsun. Conclusions: If formed via single-star evolution, the late-type WN stars in M31 stem from an initial mass range between 20 and 60 Msun. From the very late-type WN9-11 stars, only one star is located in the S Doradus instability strip. We do not find any late-type WN stars with the high luminosities known in the Milky Way.

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.

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.