The evolution and fate of very massive stars (VMS) is tightly connected to their mass-loss properties. Their initial and final masses differ significantly as a result of mass loss. VMS have strong stellar winds and extremely high ionising fluxes, which are thought to be critical sources of both mechanical and radiative feedback in giant H  regions. However, how VMS mass-loss properties change during stellar evolution is poorly understood. In the framework of the VLT-Flames Tarantula Survey (VFTS), we explore the mass-loss transition region from optically thin O star winds to denser WNh Wolf-Rayet star winds, thereby testing theoretical predictions. To this purpose we select 62 O, Of, Of/WN, and WNh stars, an unprecedented sample of stars with the highest masses and luminosities known. We perform a spectral analysis of optical VFTS as well as near-infrared VLT/SINFONI data using the non-LTE radiative transfer code CMFGEN to obtain both stellar and wind parameters. For the first time, we observationally resolve the transition between optically thin O star winds and optically thick hydrogen-rich WNh Wolf-Rayet winds. Our results suggest the existence of a "kink" between both mass-loss regimes, in agreement with recent Monte Carlo simulations. For the optically thick regime, we confirm the steep dependence on the classical Eddington factor Γ e from previous theoretical and observational studies. The transition occurs on the main sequence near a luminosity of 10 6.1 L , or a mass of 80 . . . 90 M . Above this limit, we find that -even when accounting for moderate wind clumping (with f v = 0.1) -wind mass-loss rates are enhanced with respect to standard prescriptions currently adopted in stellar evolution calculations. We also show that this results in substantial helium surface enrichment. Finally, based on our spectroscopic analyses, we are able to provide the most accurate ionising fluxes for VMS known to date, confirming the pivotal role of VMS in ionising and shaping their environments.

Context. The VLT-FLAMES Tarantula Survey has an extensive view of the copious number of massive stars in the 30 Doradus (30 Dor) star forming region of the Large Magellanic Cloud. These stars play a crucial role in our understanding of the stellar feedback in more distant, unresolved star forming regions. Aims. The first comprehensive census of hot luminous stars in 30 Dor is compiled within a 10 arcmin (150 pc) radius of its central cluster, R136. We investigate the stellar content and spectroscopic completeness of the early type stars. Estimates were made for both the integrated ionising luminosity and stellar wind luminosity. These values were used to re-assess the star formation rate (SFR) of the region and determine the ionising photon escape fraction. Methods. Stars were selected photometrically and combined with the latest spectral classifications. Spectral types were estimated for stars lacking spectroscopy and corrections were made for binary systems, where possible. Stellar calibrations were applied to obtain their physical parameters and wind properties. Their integrated properties were then compared to global observations from ultraviolet (UV) to far-infrared (FIR) imaging as well as the population synthesis code, Starburst99. Results. Our census identified 1145 candidate hot luminous stars within 150 pc of R136 of which >700 were considered to be genuine early type stars and contribute to feedback. We assess the survey to be spectroscopically complete to 85% in the outer regions (>5 pc) but only 35% complete in the region of the R136 cluster, giving a total of 500 hot luminous stars in the census which had spectroscopy. Only 31 were found to be Wolf-Rayet (W-R) or Of/WN stars, but their contribution to the integrated ionising luminosity and wind luminosity was ∼40% and ∼50%, respectively. Similarly, stars with M init > 100 M (mostly H-rich WN stars) also showed high contributions to the global feedback, ∼25% in both cases. Such massive stars are not accounted for by the current Starburst99 code, which was found to underestimate the integrated ionising luminosity of R136 by a factor ∼2 and the wind luminosity by a factor ∼9. The census inferred a SFR for 30 Dor of 0.073 ± 0.04 M yr −1 . This was generally higher than that obtained from some popular SFR calibrations but still showed good consistency with the far-UV luminosity tracer as well as the combined Hα and mid-infrared tracer, but only after correcting for Hα extinction. The global ionising output was also found to exceed that measured from the associated gas and dust, suggesting that ∼6 +55 −6 % of the ionising photons escape the region. Conclusions. When studying the most luminous star forming regions, it is essential to include their most massive stars if one is to determine a reliable energy budget. Photon leakage becomes more likely after including their large contributions to the ionising output. If 30 Dor is typical of other massive star forming regions, estimates of the SFR will be underpredicted if this escape fraction is not a...

Context. Massive stars play a vital role in the Universe. However, their evolution even on the main sequence is not yet well understood. Aims. Due to the steep mass-luminosity relation, massive main sequence stars become extremely luminous. This brings their envelopes very close to the Eddington limit. We are analysing stellar evolutionary models in which the Eddington limit is reached and exceeded, and explore the rich diversity of physical phenomena which take place in their envelopes, and we investigate their observational consequences. Methods. We use the published grids of detailed stellar models by Brott et al. (2011) andKöhler et al. (2015), computed with a state-of-the-art one-dimensional hydrodynamic stellar evolution code using LMC composition, to investigate the envelope properties of core hydrogen burning massive stars. Results. We find that at the stellar surface, the Eddington limit is almost never reached, even for stars up to 500 M ⊙ . When we define an appropriate Eddington limit locally in the stellar envelope, we can show that most stars more massive than ∼ 40 M ⊙ actually exceed this limit, in particular in the partial ionization zones of iron, helium or hydrogen. While most models adjust their structure such that the local Eddington limit is exceeded at most by a few per cent, our most extreme models do so by a factor of more than seven. We find that the local violation of the Eddington limit has severe consequences for the envelope structure, as it leads to envelope inflation, convection, density inversions and possibly to pulsations. We find that all models with luminosities higher than 4 × 10 5 L ⊙ , i.e. stars above ∼ 40 M ⊙ show inflation, with a radius increase of up to a factor of about 40. We find that the hot edge of the S Dor variability region coincides with a line beyond which our models are inflated by more than a factor of two, indicating a possible connection between S Dor variability and inflation. Furthermore, our coolest models show highly inflated envelopes with masses of up to several solar masses, and appear to be candidates to produce major LBV eruptions. Conclusions. Our models show that the Eddington limit is expected to be reached in all stars above ∼ 40 M ⊙ in the LMC, and by even lower mass stars in the Galaxy, or in close binaries or rapid rotators. While our results do not support the idea of a direct super-Eddington wind driven by continuum photons, the consequences of the Eddington limit in the form of inflation, pulsations and possibly eruptions may well give rise to a significant enhancement of the the time averaged mass loss rate.

We present an optical analysis of 55 members of R136, the central cluster in the Tarantula Nebula of the Large Magellanic Cloud. Our sample was observed with STIS aboard the Hubble Space Telescope, is complete down to about 40 M⊙, and includes 7 very massive stars with masses over 100 M⊙. We performed a spectroscopic analysis to derive their physical properties. Using evolutionary models we find that the initial mass function (IMF) of massive stars in R136 is suggestive of being top-heavy with a power-law exponent γ ≈ 2 ± 0.3, but steeper exponents cannot be excluded. The age of R136 lies between 1 and 2 Myr with a median age of around 1.6 Myr. Stars more luminous than log L/L⊙ = 6.3 are helium enriched and their evolution is dominated by mass loss, but rotational mixing or some other form of mixing could be still required to explain the helium composition at the surface. Stars more massive than 40 M⊙ have larger spectroscopic than evolutionary masses. The slope of the wind-luminosity relation assuming unclumped stellar winds is 2.41 ± 0.13 which is steeper than usually obtained (∼1.8). The ionising (log Q0 [ph/s] = 51.4) and mechanical (log LSW [erg/s] = 39.1) output of R136 is dominated by the most massive stars (>100 M⊙). R136 contributes around a quarter of the ionising flux and around a fifth of the mechanical feedback to the overall budget of the Tarantula Nebula. For a census of massive stars of the Tarantula Nebula region we combined our results with the VLT-FLAMES Tarantula Survey plus other spectroscopic studies. We observe a lack of evolved Wolf-Rayet stars and luminous blue and red supergiants.

VFTS 682 is located in an active star-forming region, at a projected distance of 29 pc from the young massive cluster R136 in the Tarantula Nebula of the Large Magellanic Cloud. It was previously reported as a candidate young stellar object, and more recently spectroscopically revealed as a hydrogen-rich Wolf-Rayet (WN5h) star. Our aim is to obtain the stellar properties, such as its intrinsic luminosity, and to investigate the origin of VFTS 682. To this purpose, we model optical spectra from the VLT-FLAMES Tarantula Survey with the non-LTE stellar atmosphere code cmfgen, as well as the spectral energy distribution from complementary optical and infrared photometry. We find the extinction properties to be highly peculiar (R V ∼ 4.7), and obtain a surprisingly high luminosity log(L/L ) = 6.5 ± 0.2, corresponding to a present-day mass of ∼150 M . The high effective temperature of 52.2 ± 2.5 kK might be explained by chemically homogeneous evolution -suggested to be the key process in the path towards long gamma-ray bursts. Lightcurves of the object show variability at the 10% level on a timescale of years. Such changes are unprecedented for classical WolfRayet stars, and are more reminiscent of Luminous Blue Variables. Finally, we discuss two possibilities for the origin of VFTS 682: (i) the star either formed in situ, which would have profound implications for the formation mechanism of massive stars, or (ii) VFTS 682 is a slow runaway star that originated from the dense cluster R136, which would make it the most massive runaway known to date.

We report the first APOGEE metallicities and α-element abundances measured for 3600 red giant stars spanning a large radial range of both the Large (LMC) and Small Magellanic Clouds, the largest Milky Way (MW) dwarf galaxies. Our sample is an order of magnitude larger than that of previous studies and extends to much larger radial distances. These are the first results presented that make use of the newly installed southern APOGEE instrument on the du Pont telescope at Las Campanas Observatory. Our unbiased sample of the LMC spans a large range in metallicity, from [Fe/H] = −0.2 to very metal-poor stars with [Fe/H] ≈ −2.5, the most metal-poor Magellanic Cloud (MC) stars detected to date. The LMC [α/Fe]–[Fe/H] distribution is very flat over a large metallicity range but rises by ∼0.1 dex at −1.0 < [Fe/H] ≲ −0.5. We interpret this as a sign of the known recent increase in MC star formation activity and are able to reproduce the pattern with a chemical evolution model that includes a recent “starburst.” At the metal-poor end, we capture the increase of [α/Fe] with decreasing [Fe/H] and constrain the “α-knee” to [Fe/H] ≲ −2.2 in both MCs, implying a low star formation efficiency of ∼0.01 Gyr−1. The MC knees are more metal-poor than those of less massive MW dwarf galaxies such as Fornax, Sculptor, or Sagittarius. One possible interpretation is that the MCs formed in a lower-density environment than the MW, a hypothesis that is consistent with the paradigm that the MCs fell into the MW’s gravitational potential only recently.

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.

Context. The commonly used extinction laws of Cardelli et al. (1989, ApJ, 345, 245) have limitations that, among other issues, hamper the determination of the effective temperatures of O and early B stars from optical and near-infrared (NIR) photometry. Aims. We aim to develop a new family of extinction laws for 30 Doradus, check their general applicability within that region and elsewhere, and apply them to test the feasibility of using optical and NIR photometry to determine the effective temperature of OB stars. Methods. We use spectroscopy and NIR photometry from the VLT-FLAMES Tarantula Survey and optical photometry from HST/WFC3 of 30 Doradus and we analyze them with the software code CHORIZOS using different assumptions, such as the family of extinction laws. Results. We derive a new family of optical and NIR extinction laws for 30 Doradus and confirm its applicability to extinguished Galactic O-type systems. We conclude that by using the new extinction laws it is possible to measure the effective temperatures of OB stars with moderate uncertainties and only a small bias, at least up to E(4405 − 5495) ∼ 1.5 mag.