We report the discovery of a handful of optical hydrogen-poor knots in the central part of an extended infrared nebula centred on the [WO1] star WR 72, obtained by spectroscopic and imaging observations with the Southern African Large Telescope (SALT). Wide-field Infrared Survey Explorer (WISE) images of the nebula show that it is composed of an extended almost circular halo (of $\approx6$ arcmin or $\approx2.4$ pc in diameter) and an elongated and apparently bipolar inner shell (of a factor of six smaller size), within which the knots are concentrated. Our findings indicate that WR 72 is a new member of the rare group of hydrogen-poor planetary nebulae, which may be explained through a very late thermal pulse of a post-AGB star, or by a merger of two white dwarfs.

Wolf-Rayet star's winds can be so dense and so optically thick that the photosphere appears in the highly supersonic part of the outflow, veiling the underlying subsonic part of the star, and leaving the initial acceleration of the wind inaccessible to observations. We investigate the conditions and the structure of the subsonic part of the outflow of Galactic WR stars, in particular of the WNE subclass; our focus is on the conditions at the sonic point. We compute 1D hydrodynamic stellar structure models for massive helium stars adopting outer boundaries at the sonic point. We find that the outflows of our models are accelerated to supersonic velocities by the radiative force from opacity bumps either at temperatures of the order of 200kK by the Fe opacity bump or of the order of 50kK by the HeII opacity bump. For a given mass-loss rate, the conditions in the subsonic part of the outflow are independent from the detailed physical conditions in the supersonic part. The close proximity to the Eddington limit at the sonic point allows us to construct a Sonic HR diagram, relating the sonic point temperature to the L/M ratio and the stellar mass-loss rate, thereby constraining the sonic point conditions, the subsonic structure, and the stellar wind mass-loss rates from observations. The minimum mass-loss rate necessary to have the flow accelerated to supersonic velocities by the Fe opacity bump is derived. A comparison of the observed parameters of Galactic WNE stars to this minimum mass-loss rate indicates that their winds are launched to supersonic velocities by the radiation pressure arising from the Fe-bump. Conversely, models which do not show transonic flows from the Fe opacity bump form inflated envelopes. We derive an analytic criterion for the appearance of envelope inflation in the subphotospheric layers.

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.

Theoretically, rotation-induced chemical mixing in massive stars has far reaching evolutionary consequences, affecting the sequence of morphological phases, lifetimes, nucleosynthesis, and supernova characteristics. Using a sample of 72 presumably single O-type giants to supergiants observed in the context of the VLT-FLAMES Tarantula Survey (VFTS), we aim to investigate rotational mixing in evolved core-hydrogen burning stars initially more massive than $15\,M_\odot$ by analysing their surface nitrogen abundances. Using stellar and wind properties derived in a previous VFTS study, we constrained the nitrogen abundance by fitting the equivalent widths of relatively strong lines that are sensitive to changes in the abundance of this element. Given the quality of the data, we constrained the nitrogen abundance in 38 cases; for 34 stars only upper limits could be derived, which includes almost all stars rotating at $v_\mathrm{e}\sin i >200\,\mathrm{km s^{-1}}$. We analysed the nitrogen abundance as a function of projected rotation rate $v_\mathrm{e}\sin i$ and confronted it with predictions of rotational mixing. The upper limits on the nitrogen abundance of the rapidly rotating stars are not in apparent violation with theoretical expectations. However, we found a group of N-enhanced slowly-spinning stars that is not in accordance with predictions of rotational mixing in single stars. Among O-type stars with (rotation-corrected) gravities less than $\log\,g_c = 3.75$ this group constitutes 30$-$40 percent of the population. We found a correlation between nitrogen and helium abundance which is consistent with expectations, suggesting that, whatever the mechanism that brings N to the surface, it displays CNO-processed material.

The progenitors of many Type II supernovae have been observationally identified but the search for Type Ibc supernova (SN Ibc) progenitors has thus far been unsuccessful, despite the expectation that they are luminous Wolf-Rayet (WR) stars. We investigate how the evolution of massive helium stars affects their visual appearances, and discuss the implications for the detectability of SN Ibc progenitors. Massive WR stars that rapidly lose their helium envelopes through stellar-wind mass-loss end their lives when their effective temperatures -- related to their hydrostatic surfaces -- exceed about 150kK.At their pre-supernova stage, their surface properties resemble those of hot Galactic WR stars of WO sub-type. These are visually faint with narrow-band visual magnitudes Mv = -1.5 ~ -2.5, despite their high bolometric luminosities (log L/Lsun = 5.6 ~ 5.7), compared to the bulk of Galactic WR stars (Mv < -4). In contrast, relatively low-mass helium stars that retain a thick helium envelope appear fairly bright in optical bands, depending on the final masses and the history of the envelope expansion during the late evolutionary stages. We conclude that SNe Ibc observations have so far not provided strong constraints on progenitor bolometric luminosities and masses, even with the deepest searches. We also argue that Ic progenitors are more challenging to identify than Ib progenitors in any optical images.

Some studies have claimed the existence of a stellar upper-mass limit of 150 Msun. A factor that is often overlooked concerns the issue that there might be a significant difference between the present-day and the initial mass of the most massive stars - as a result of mass loss. The upper-mass limit may be substantially higher, possibly exceeding 200 Msun. The issue of the upper mass-limit will however remain uncertain as long as there is only limited quantitative knowledge of mass loss in close proximity to the Eddington (= Gamma) limit. For this reason, we present mass-loss predictions from Monte Carlo radiative transfer models for very massive stars up to 300 Msun. Using our new dynamical approach, we find an upturn or "kink" in the mass-loss versus Gamma dependence, at the point where our model winds become optically thick. These are the first mass-loss predictions where the transition from optically thin O-star winds to optically thick Wolf-Rayet winds has been resolved.

Long gamma-ray bursts involve the most powerful cosmic explosions since the Big Bang. Whilst it has been established that GRBs are related to the death throes of massive stars, the identification of their progenitors has proved challenging. Theory suggests that rotating Wolf-Rayet stars are the best candidates, but their strong stellar winds shroud their surfaces, preventing a direct measurement of their rotation. Fortunately, linear spectropolarimetry may be used to probe the flattening of their winds due to stellar spin. Spectropolarimetry surveys show that an 80% majority of WR stars have spherically symmetric winds and are thus rotating slowly, yet a small 20% minority display a spectropolarimetric signature indicative of rotation. Here we find a highly significant correlation between WR objects that carry the signature of stellar rotation and the subset of WR stars with ejecta nebulae that have only recently transitioned from a red sugergiant or luminous blue variable phase. As these youthful WR stars have yet to spin-down due to mass loss, they are the best candidate GRB progenitors identified to date. When we take recently published WR ejecta nebula numbers we find that five out of the six line-effect WR stars are surrounded by ejecta nebulae. The statistics imply that the null hypothesis of no correlation between line-effect WR stars and ejecta nebulae can be rejected at the 0.0004% level. Given that four line-effect and WR ejecta nebula have spectroscopically been confirmed to contain nucleosynthetic products, we argue that the correlation is both statistically significant and physically convincing. The implication is that we have identified a WR sub-population that fulfills the necessary criteria for making GRBs. Finally, we discuss the potential of identifying GRB progenitors via spectropolarimetry with extremely large telescopes.

We report magnetic and spectroscopic observations and modeling of the Of?p star HD 148937 within the context of the MiMeS LP at the CFHT. Thirty-two high signal-to-noise ratio circularly polarised (Stokes V) spectra and 13 unpolarised (Stokes I) spectra of HD 148937 were acquired in 2009 and 2010. A definite detection of a Stokes V Zeeman signature is obtained in the grand mean of all observations (in both LSD mean profiles and individual spectral lines). The longitudinal magnetic field inferred from the Stokes V LSD profiles is consistently negative, in contrast to the essentially zero field strength measured from the diagnostic null profiles. A period search of equivalent width measurements confirms the previously-reported 7.03 d variability period. The variation of equivalent widths is not strictly periodic: we present evidence for evolution of the amount or distribution of circumstellar plasma. Interpreting the 7.03 d period as the stellar rotational period within the context of the ORM, we have phased the equivalent widths and longitudinal field measurements. The longitudinal field measurements show a weak sinusoidal variation of constant sign, with extrema out of phase with the H{\alpha} variation by about 0.25 cycles. The inferred magnetic configuration confirms the suggestion of Naz\'e et al (2010), who proposed that the weaker variability of HD 148937 as compared to other members of this class is a consequence of the stellar geometry. Based on the derived magnetic properties and published wind characteristics, we find a wind magnetic confinement parameter \eta\ast \simeq 20 and rotation parameter W = 0.12, supporting a picture in which the Halpha emission and other line variability have their origin in an oblique, rigidly rotating magnetospheric structure resulting from a magnetically channeled wind. (Abridged.)

Some studies have claimed a universal stellar upper-mass limit of 150 Msun. A factor that is often overlooked is that there might be a difference between the current and initial masses of the most massive stars, as a result of mass loss. We present Monte Carlo mass-loss predictions for very massive stars in the range 40-300 Msun, with large luminosities and Eddington factors Gamma. Using our new dynamical approach, we find an upturn in the mass-loss vs. Gamma dependence, at the point where the winds become optically thick. This coincides with the location where wind efficiency numbers surpass the single-scattering limit of Eta = 1, reaching values up to Eta = 2.5. Our modelling suggests a transition from common O-type winds to Wolf-Rayet characteristics at the point where the winds become optically thick. This transitional behaviour is also revealed with respect to the wind acceleration parameter beta, which starts at values below 1 for the optically thin O-stars, and naturally reaches values as high as 1.5-2 for the optically thick Wolf-Rayet models. An additional finding concerns the transition in spectral morphology of the Of and WN characteristic He II line at 4686 Angstrom. When we express our mass-loss predictions as a function of the electron scattering Gamma_e (=L/M) only, we obtain a mass-loss Gamma dependence that is consistent with a previously reported power-law Mdot propto Gamma^5 (Vink 2006) that was based on our semi-empirical modelling approach. When we express Mdot in terms of both Gamma and stellar mass, we find Mdot propto M^0.8 Gamma^4.8 for our high Gamma models. Finally, we confirm that the Gamma-effect on the mass-loss predictions is much stronger than that of an increased helium abundance, calling for a fundamental revision in the way mass loss is incorporated in evolutionary models of the most massive stars.

The VLT-FLAMES Tarantula Survey (VFTS) is an ESO Large Programme that has obtained multi-epoch optical spectroscopy of over 800 massive stars in the 30 Doradus region of the Large Magellanic Cloud (LMC). Here we introduce our scientific motivations and give an overview of the survey targets, including optical and near-infrared photometry and comprehensive details of the data reduction. One of the principal objectives was to detect massive binary systems via variations in their radial velocities, thus shaping the multi-epoch observing strategy. Spectral classifications are given for the massive emission-line stars observed by the survey, including the discovery of a new Wolf-Rayet star (VFTS 682, classified as WN5h), 2' to the northeast of R136. To illustrate the diversity of objects encompassed by the survey, we investigate the spectral properties of sixteen targets identified by Gruendl & Chu from Spitzer photometry as candidate young stellar objects or stars with notable mid-infrared excesses. Detailed spectral classification and quantitative analysis of the O- and B-type stars in the VFTS sample, paying particular attention to the effects of rotational mixing and binarity, will be presented in a series of future articles to address fundamental questions in both stellar and cluster evolution.

The Tarantula survey is an ESO Large Programme which has obtained multi-epochs spectroscopy of over 800 massive stars in the 30 Dor region in the Large Magelanic Cloud. Here we briefly describe the main drivers of the survey and the observational material derived.

The location of B supergiants in the Hertzsprung-Russell diagram (HRD) represents a long-standing problem in massive star evolution. Here we propose their nature may be revealed utilising their rotational properties, and we highlight a steep drop in massive star rotation rates at an effective temperature of 22000 K. We discuss two potential explanations for it. On the one hand, the feature might be due to the end of the main sequence, which could potentially constrain the core overshooting parameter. On the other hand, the feature might be the result of enhanced mass loss at the predicted location of the bi-stability jump. We term this effect "bi-stability breaking" and discuss its potential consequences for the evolution of massive stars.

The Tarantula Survey is an ambitious ESO Large Programme that has obtained multi-epoch spectroscopy of over 1,000 massive stars in the 30 Doradus region of the Large Magellanic Cloud. Here we introduce the scientific motivations of the survey and give an overview of the observational sample. Ultimately, quantitative analysis of every star, paying particular attention to the effects of rotational mixing and binarity, will be used to address fundamental questions in both stellar and cluster evolution.

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

The mass loss from Wolf-Rayet (WR) stars is of fundamental importance for the final fate of massive stars and their chemical yields. Its Z-dependence is discussed in relation to the formation of long-duration Gamma Ray Bursts (GRBs) and the yields from early stellar generations. However, the mechanism of formation of WR-type stellar winds is still under debate. We present the first fully self-consistent atmosphere/wind models for late-type WN stars. We investigate the mechanisms leading to their strong mass loss, and examine the dependence on stellar parameters, in particular on the metallicity Z. We identify WNL stars as very massive stars close to the Eddington limit, potentially still in the phase of central H-burning. Due to their high L/M ratios, these stars develop optically thick, radiatively driven winds. These winds show qualitatively different properties than the thin winds of OB stars. The resultant mass loss depends strongly on Z, but also on the Eddington factor, and the stellar temperature. We combine our results in a parametrized mass loss recipe for WNL stars. According to our present model computations, stars close to the Eddington limit tend to form strong WR-type winds, even at very low Z. Our models thus predict an efficient mass loss mechanism for low metallicity stars. For extremely metal-poor stars, we find that the self-enrichment with primary nitrogen can drive WR-type mass loss. These first WN stars might play an important role in the enrichment of the early ISM with freshly produced nitrogen.

CONTEXT: Very massive stars pass through the Wolf-Rayet (WR) stage before they finally explode. Details of their evolution have not yet been safely established, and their physics are not well understood. Their spectral analysis requires adequate model atmospheres, which have been developed step by step during the past decades and account in their recent version for line blanketing by the millions of lines from iron and iron-group elements. However, only very few WN stars have been re-analyzed by means of line-blanketed models yet. AIMS: The quantitative spectral analysis of a large sample of Galactic WN stars with the most advanced generation of model atmospheres should provide an empirical basis for various studies about the origin, evolution, and physics of the Wolf-Rayet stars and their powerful winds. METHODS: We analyze a large sample of Galactic WN stars by means of the Potsdam Wolf-Rayet (PoWR) model atmospheres, which account for iron line blanketing and clumping. The results are compared with a synthetic population, generated from the Geneva tracks for massive star evolution. RESULTS: We obtain a homogeneous set of stellar and atmospheric parameters for the Galactic WN stars, partly revising earlier results. CONCLUSIONS: Comparing the results of our spectral analyses of the Galactic WN stars with the predictions of the Geneva evolutionary calculations, we conclude that there is rough qualitative agreement. However, the quantitative discrepancies are still severe, and there is no preference for the tracks that account for the effects of rotation. It seems that the evolution of massive stars is still not satisfactorily understood.