The mass of a star is the most fundamental parameter for its structure, evolution, and final fate. It is particularly important for any kind of stellar archaeology and characterization of exoplanets. There exists a variety of methods in astronomy to estimate or determine it. In this review we present a significant number of such methods, beginning with the most direct and model-independent approach using detached eclipsing binaries. We then move to more indirect and model-dependent methods, such as the quite commonly used isochrone or stellar track fitting. The arrival of quantitative asteroseismology has opened a completely new approach to determine stellar masses and to complement and improve the accuracy of other methods. We include methods for different evolutionary stages, from the pre-main sequence to evolved (super)giants and final remnants. For all methods uncertainties and restrictions will be discussed. We provide lists of altogether more than 200 benchmark stars with relative mass accuracies between $[0.3,2]\%$ for the covered mass range of $M\in [0.1,16]\,M_\odot$, $75\%$ of which are stars burning hydrogen in their core and the other $25\%$ covering all other evolved stages. We close with a recommendation how to combine various methods to arrive at a "mass-ladder" for stars.

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.

Mass loss through stellar winds plays a dominant role in the evolution of massive stars. In particular the mass-loss rates of very massive stars (VMSs, $> 100\,M_{\odot}$) are highly uncertain. Such stars display Wolf-Rayet spectral morphologies (WNh) whilst on the main-sequence. Metal-poor VMSs are progenitors of gamma-ray bursts and pair instability supernovae. In this study we extended the widely used stellar wind theory by Castor, Abbott & Klein from the optically thin (O star) to the optically thick main-sequence (WNh) wind regime. In particular we modify the mass-loss rate formula in a way that we are able to explain the empirical mass-loss dependence on the Eddington parameter ($\Gamma_{\rm e}$). The new mass-loss recipe is suitable for incorporation into current stellar evolution models for massive and very massive stars. It makes verifiable predictions, namely how the mass-loss rate scales with metallicity and at which Eddington parameter the transition from optically thin O star to optically thick WNh star winds occurs. In the case of the star cluster R136 in the Large Magellanic Cloud we find in the optically thin wind regime $\dot{M} \propto \Gamma_{\rm e}^{3}$ while in the optically thick wind regime $\dot{M} \propto 1/ (1 - \Gamma_{\rm e})^{3.5}$. The transition from optically thin to optically thick winds occurs at $\Gamma_{\rm e, trans} \approx 0.47$. The transition mass-loss rate is $\log \dot{M}~(M_{\odot} \mathrm{yr}^{-1}) \approx -4.76 \pm 0.18$, which is in line with the prediction by Vink & Gr\"afener assuming a volume filling factor of $f_{\rm V} = 0.23_{-0.15}^{+0.40}$.

The process of massive star ($M\geq8~M_\odot$) formation is still poorly understood. Observations of massive young stellar objects (MYSOs) are challenging due to their rarity, short formation timescale, large distances, and high circumstellar extinction. Here, we present the results of a spectroscopic analysis of a population of MYSOs in the Large Magellanic Cloud (LMC). We took advantage of the spectral resolution and wavelength coverage of X-shooter (300-2500 nm), mounted on the European Southern Observatory Very Large Telescope, to detect characteristic spectral features in a dozen MYSO candidates near 30 Doradus, the largest starburst region in the Local Group hosting the most massive stars known. The X-shooter spectra are strongly contaminated by nebular emission. We used a scaling method to subtract the nebular contamination from our objects. We detect H$\alpha,\beta$, [O I] 630.0 nm, Ca II infrared triplet, [Fe II] 1643.5 nm, fluorescent Fe II 1687.8 nm, H$_2$ 2121.8 nm, Br$\gamma$, and CO bandhead emission in the spectra of multiple candidates. This leads to the spectroscopic confirmation of 10 candidates as bona fide MYSOs. We compare our observations with photometric observations from the literature and find all MYSOs to have a strong near-infrared excess. We compute lower limits to the brightness and luminosity of the MYSO candidates, confirming the near-infrared excess and the massive nature of the objects. No clear correlation is seen between the Br$\gamma$ luminosity and metallicity. Combining our sample with other LMC samples results in a combined detection rate of disk features such as fluorescent Fe II and CO bandheads which is consistent with the Galactic rate (40\%). Most of our MYSOs show outflow features.

We report the first APOGEE metallicities and alpha-element abundances measured for 3600 red giant stars spanning a large radial range of both the Large (LMC) and Small Magellanic Clouds (SMC), the largest Milky Way 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 Clouds (MCs) stars detected to date. The LMC [alpha/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 [alpha/Fe] with decreasing [Fe/H], and constrain the "alpha-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 Milky Way (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.

Cosmic History has witnessed the lives and deaths of multiple generations of massive stars, all of them invigorating their host galaxies with ionizing photons, kinetic energy, fresh material and stellar-mass black holes. Ubiquitous engines as they are, Astrophysics needs a good understanding of their formation, evolution, properties and yields throughout the history of the Universe, and with decreasing metal content mimicking the environment at the earliest epochs. Ultimately, a physical model that could be extrapolated to zero metallicity would enable tackling long-standing questions such as "What did the First, very massive stars of the Universe look like?" or "What was their role in the re-ionization of the Universe?". Yet, most our knowledge of metal-poor massive stars is drawn from one single point in metallicity. Massive stars in the Small Magellanic Cloud (SMC, $\sim 1/5 Z_{\odot}$) currently serve as templates for low-metallicity objects in the early Universe, even though significant differences with respect to massive stars with poorer metal content have been reported. This White Paper summarizes the current knowledge on extremely (sub-SMC) metal poor massive stars, highlighting the most outstanding open questions and the need to supersede the SMC as standard. A new paradigm can be built from nearby extremely metal-poor galaxies that make a new metallicity ladder, but massive stars in these galaxies are out of reach to current observational facilities. Such task would require an L-size mission, consisting of a 10m-class space telescope operating in the optical and the ultraviolet ranges. Alternatively, we propose that ESA unites efforts with NASA to make the LUVOIR mission concept a reality, thus continuing the successful partnership that made Hubble Space Telescope one of the greatest observatories of all time.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

We introduce the 4-metre Multi-Object Spectroscopic Telescope (4MOST), a new high-multiplex, wide-field spectroscopic survey facility under development for the four-metre-class Visible and Infrared Survey Telescope for Astronomy (VISTA) at Paranal. Its key specifications are: a large field of view (FoV) of 4.2 square degrees and a high multiplex capability, with 1624 fibres feeding two low-resolution spectrographs ($R = \lambda/\Delta\lambda \sim 6500$), and 812 fibres transferring light to the high-resolution spectrograph ($R \sim 20\,000$). After a description of the instrument and its expected performance, a short overview is given of its operational scheme and planned 4MOST Consortium science; these aspects are covered in more detail in other articles in this edition of The Messenger. Finally, the processes, schedules, and policies concerning the selection of ESO Community Surveys are presented, commencing with a singular opportunity to submit Letters of Intent for Public Surveys during the first five years of 4MOST operations.

The One Thousand and One Magellanic Fields (1001MC) survey aims to measure the kinematics and elemental abundances of many different stellar populations that sample the history of formation and interaction of the Magellanic Clouds. The survey will collect spectra of about half a million stars with $G < 19.5$ magnitudes (Vega) distributed over an area of about 1000 square degrees and will provide an invaluable dataset for a wide range of scientific applications.

The 4MOST instrument is a multi-object spectrograph that will address Galactic and extragalactic science cases simultaneously by observing targets from a large number of different surveys within each science exposure. This parallel mode of operation and the survey nature of 4MOST require some distinct 4MOST-specific operational features within the overall operations model of ESO. The main feature is that the 4MOST Consortium will deliver, not only the instrument, but also contractual services to the user community, which is why 4MOST is also described as a facility. This white paper concentrates on information particularly useful to answering the forthcoming Call for Letters of Intent.

The Maunakea Spectroscopic Explorer (MSE) is a planned 11.25-m aperture facility with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. A rebirth of the 3.6m Canada-France-Hawaii Telescope on Maunakea, MSE will use 4332 fibers operating at three different resolving powers (R ~ 2500, 6000, 40000) across a wavelength range of 0.36-1.8mum, with dynamical fiber positioning that allows fibers to match the exposure times of individual objects. MSE will enable spectroscopic surveys with unprecedented scale and sensitivity by collecting millions of spectra per year down to limiting magnitudes of g ~ 20-24 mag, with a nominal velocity precision of ~100 m/s in high-resolution mode. This white paper describes science cases for stellar astrophysics and exoplanet science using MSE, including the discovery and atmospheric characterization of exoplanets and substellar objects, stellar physics with star clusters, asteroseismology of solar-like oscillators and opacity-driven pulsators, studies of stellar rotation, activity, and multiplicity, as well as the chemical characterization of AGB and extremely metal-poor stars.

Here we confirm Melnick 34, an X-ray bright star in the 30 Doradus region of the Large Magellanic Cloud, as an SB2 binary comprising WN5h+WN5h components. We present orbital solutions using 26 epochs of VLT/UVES spectra and 22 epochs of archival Gemini/GMOS spectra. Radial-velocity monitoring and automated template fitting methods both reveal a similar high eccentricity system with a mass ratio close to unity, and an orbital period in agreement with the 155.1 +/- 1 day X-ray light curve period previously derived by Pollock et al. Our favoured solution derived an eccentricity of 0.68 +/- 0.02 and mass ratio of 0.92 +/- 0.07, giving minimum masses of Ma_sin^{3}(i) = 65 +/- 7 Msun and Mb_sin^{3}(i) = 60 +/- 7 Msun. Spectral modelling using WN5h templates with CMFGEN reveals temperatures of T ~53 kK for each component and luminosities of log(La/Lsun) = 6.43 +/- 0.08 and log(Lb/Lsun) = 6.37 +/- 0.08, from which BONNSAI evolutionary modelling gives masses of Ma = 139 (+21,-18) Msun and Mb = 127 (+17,-17) Msun and ages of ~0.6 Myrs. Spectroscopic and dynamic masses would agree if Mk34 has an inclination of i ~50{\deg}, making Mk34 the most massive binary known and an excellent candidate for investigating the properties of colliding wind binaries. Within 2-3 Myrs, both components of Mk34 are expected to evolve to stellar mass black holes which, assuming the binary system survives, would make Mk34 a potential binary black hole merger progenitor and gravitational wave source.

How very massive stars form is still an open question in astrophysics. VFTS682 is among the most massive stars known, with an inferred initial mass of $\sim$150 $M_\odot$ . It is located in 30 Doradus at a projected distance of 29 pc from the central cluster R136. Its apparent isolation led to two hypotheses: either it formed in relative isolation or it was ejected dynamically from the cluster. We investigate the kinematics of VFTS682 as obtained by Gaia and Hubble Space Telescope astrometry. We derive a projected velocity relative to the cluster of $38 \pm 17 \mathrm{km \ s^{-1}}$ (1$\sigma$ confidence interval). Although the error bars are substantial, two independent measures suggest that VFTS682 is a runaway ejected from the central cluster. This hypothesis is further supported by a variety of circumstantial clues. The central cluster is known to harbor other stars more massive than 150 $M_\odot$ of similar spectral type and recent astrometric studies on VFTS16 and VFTS72 provide direct evidence that the cluster can eject some of its most massive members, in agreement with theoretical predictions. If future data confirm the runaway nature, this would make VFTS682 the most massive runaway star known to date.

Eta Car is one of the most intriguing luminous blue variables in the Galaxy. Observations and models at different wavelengths suggest a central binary with a 5.54 yr period residing in its core. 2D and 3D radiative transfer and hydrodynamic simulations predict a primary with a dense and slow stellar wind that interacts with the faster and lower density wind of the secondary. The wind-wind collision scenario suggests that the secondary's wind penetrates the primary's wind creating a low-density cavity in it, with dense walls where the two winds interact. We aim to trace the inner ~5-50 au structure of Eta Car's wind-wind interaction, as seen through BrG and, for the first time, through the He I 2s-2p line. We have used spectro-interferometric observations with GRAVITY at the VLTI. Our modeling of the continuum allows us to estimate its FWHM angular size close to 2 mas and an elongation ratio of 1.06 +/- 0.05 over a PA = 130 +/- 20 deg. Our CMFGEN modeling helped us to confirm that the role of the secondary should be taken into account to properly reproduce the observed BrG and He I lines. Chromatic images across BrG reveal a southeast arc-like feature, possibly associated to the hot post-shocked winds flowing along the cavity wall. The images of He I 2s-2p served to constrain the 20 mas structure of the line-emitting region. The observed morphology of He I suggests that the secondary is responsible for the ionized material that produces the line profile. Both the BrG and the He I 2s-2p maps are consistent with previous hydrodynamical models of the colliding wind scenario. Future dedicated simulations together with an extensive interferometric campaign are necessary to refine our constraints on the wind and stellar parameters of the binary, which finally will help us predict the evolutionary path of Eta Car.

Farr and Mandel reanalyse our data, finding initial-mass-function slopes for high mass stars in 30 Doradus that agree with our results. However, their reanalysis appears to underpredict the observed number of massive stars. Their technique results in more precise slopes than in our work, strengthening our conclusion that there is an excess of massive stars above $30\,\mathrm{M}_\odot$ in 30 Doradus.

The 30 Doradus (30 Dor) nebula in the Large Magellanic Cloud (LMC) is the brightest HII region in the Local Group and a prototype starburst similar to those found in high redshift galaxies. It is thus a stepping stone to understand the complex formation processes of stars in starburst regions across the Universe. Here, we have studied the formation history of massive stars in 30 Dor using masses and ages derived for 452 mainly OB stars from the spectroscopic VLT-FLAMES Tarantula Survey (VFTS). We find that stars of all ages and masses are scattered throughout 30 Dor. This is remarkable because it implies that massive stars either moved large distances or formed independently over the whole field of view in relative isolation. We find that both channels contribute to the 30 Dor massive star population. Massive star formation rapidly accelerated about 8 Myr ago, first forming stars in the field before giving birth to the stellar populations in NGC 2060 and NGC 2070. The R136 star cluster in NGC 2070 formed last and, since then, about 1 Myr ago, star formation seems to be diminished with some continuing in the surroundings of R136. Massive stars within a projected distance of 8 pc of R136 are not coeval but show an age range of up to 6 Myr. Our mass distributions are well populated up to $200\,\mathrm{M}_\odot$. The inferred IMF is shallower than a Salpeter-like IMF and appears to be the same across 30 Dor. By comparing our sample of stars to stellar models in the Hertzsprung-Russell diagram, we find evidence for missing physics in the models above $\log L/\mathrm{L}_\odot=6$ that is likely connected to enhanced wind mass loss for stars approaching the Eddington limit. [abridged]

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.

The Tarantula region in the Large Magellanic Cloud contains the richest population of spatially resolved massive O-type stars known so far. This unmatched sample offers an opportunity to test models describing their main-sequence evolution and mass-loss properties. Using ground-based optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS), we aim to determine stellar, photospheric and wind properties of 72 presumably single O-type giants, bright giants and supergiants and to confront them with predictions of stellar evolution and of line-driven mass-loss theories. We apply an automated method for quantitative spectroscopic analysis of O stars combining the non-LTE stellar atmosphere model {{\sc fastwind}} with the genetic fitting algorithm {{\sc pikaia}} to determine the following stellar properties: effective temperature, surface gravity, mass-loss rate, helium abundance, and projected rotational velocity. We present empirical effective temperature versus spectral subtype calibrations at LMC-metallicity for giants and supergiants. In the spectroscopic and classical Hertzsprung-Russell diagrams, our sample O stars are found to occupy the region predicted to be the core hydrogen-burning phase by Brott et al. (2011) and K\"{o}hler et al. (2015). Except for five stars, the helium abundance of our sample stars is in agreement with the initial LMC composition. The aforementioned five stars present moderate projected rotational velocities (i.e., $v_{\mathrm{e}}\,\sin\,i\,<\,200\,\mathrm{km\,s^{-1}}$) and hence do not agree with current predictions of rotational mixing in main-sequence stars. Adopting theoretical results for the wind velocity law, we find modified wind momenta for LMC stars that are $\sim$0.3 dex higher than earlier results. [Due to the limitation of characters, the abstract appearing here is slightly shorter than that in the PDF file.]

The 30 Doradus star-forming region in the Large Magellanic Cloud is a nearby analogue of large star-formation events in the distant Universe. We determine the recent formation history and the initial mass function (IMF) of massive stars in 30 Doradus based on spectroscopic observations of 247 stars more massive than 15 solar masses ($\mathrm{M}_\odot$). The main episode of massive star formation started about $8\,\mathrm{Myr}$ ago and the star-formation rate seems to have declined in the last $1\,\mathrm{Myr}$. The IMF is densely sampled up to $200\,\mathrm{M}_\odot$ and contains $32\pm12\%$ more stars above $30\,\mathrm{M}_\odot$ than predicted by a standard Salpeter IMF. In the mass range $15-200\,\mathrm{M}_\odot$, the IMF power-law exponent is $1.90^{+0.37}_{-0.26}$, shallower than the Salpeter value of 2.35.

We present a number of notable results from the VLT-FLAMES Tarantula Survey (VFTS), an ESO Large Program during which we obtained multi-epoch medium-resolution optical spectroscopy of a very large sample of over 800 massive stars in the 30 Doradus region of the Large Magellanic Cloud (LMC). This unprecedented data-set has enabled us to address some key questions regarding atmospheres and winds, as well as the evolution of (very) massive stars. Here we focus on O-type runaways, the width of the main sequence, and the mass-loss rates for (very) massive stars. We also provide indications for the presence of a top-heavy initial mass function (IMF) in 30 Dor.

Characterization of the dynamics of massive star systems and the astrophysical properties of the interacting components are a prerequisite for understanding their formation and evolution. Optical interferometry at milliarcsecond resolution is a key observing technique for resolving high-mass multiple compact systems. Here we report on VLTI/GRAVITY, Magellan/FIRE, and MPG2.2m/FEROS observations of the late-O/early-B type system HD 93206 A, which is a member of the massive cluster Collinder 228 in the Carina nebula complex. With a total mass of about 90 M_sun, it is one of the most compact massive-quadruple systems known. In addition to measuring the separation and position angle of the outer binary Aa - Ac, we observe BrG and HeI variability in phase with the orbital motion of the two inner binaries. From the differential phases analysis, we conclude that the BrG emission arises from the interaction regions within the components of the individual binaries, which is consistent with previous models for the X-ray emission of the system based on wind-wind interaction. With an average 3-sigma deviation of ~15 deg, we establish an upper limit of p ~ 0.157 mas (0.35 AU) for the size of the BrG line-emitting region. Future interferometric observations with GRAVITY using the 8m UTs will allow us to constrain the line-emitting regions down to angular sizes of 20 uas (0.05 AU at the distance of the Carina nebula).

Red giant stars are perhaps the most important type of stars for Galactic and extra-galactic archaeology: they are luminous, occur in all stellar populations, and their surface temperatures allow precise abundance determinations for many different chemical elements. Yet, the full star formation and enrichment history of a galaxy can be traced directly only if two key observables can be determined for large stellar samples - age and chemical composition. While spectroscopy is a powerful method to analyse the detailed abundances of stars, stellar ages are the "missing link in the chain", since they are not a direct observable. However, spectroscopy should be able to estimate stellar masses, which for red giants directly infer ages provided their chemical composition is known. Here we establish a new empirical relation between the shape of the hydrogen line in the observed spectra of red giants and stellar mass determined from asteroseismology. The relation allows to determine stellar masses and ages with the accuracy of 10-15%. The method can be used with confidence for stars in the following range of stellar parameters: 4000 < Teff < 5000 K, 0.5 < log g < 3.5, -2.0 < [Fe/H] < 0.3, and luminosities log L/LSun < 2.5. Our analysis provides observational evidence that the Halpha spectral characteristics of red giant stars are tightly correlated with their mass and therefore their age. We also show that the method samples well all stellar populations with ages above 1 Gyr. Targeting bright giants, the method allows to obtain simultaneous age and chemical abundance information far deeper than would be possible with asteroseismology, extending the possible survey volume to remote regions of the Milky Way and even to neighbouring galaxies like Andromeda or the Magellanic Clouds already with present instrumentation, like VLT and Keck facilities.

The envelopes of stars near the Eddington limit are prone to various instabilities. A high Eddington factor in connection with the Fe opacity peak leads to convective instability, and a corresponding envelope inflation may induce pulsational instability. Here, we investigate the occurrence and consequences of both instabilities in models of Wolf-Rayet stars. We determine the convective velocities in the sub-surface convective zones to estimate the amplitude of the turbulent velocity at the base of the wind that potentially leads to the formation of small-scale wind structures, as observed in several WR stars. We also investigate the effect of mass loss on the pulsations of our models. We approximated solar metallicity WR stars by models of mass-losing helium stars, and we characterized the properties of convection in the envelope adopting the standard MLT. Our results show the occurrence of sub-surface convective regions in all studied models. Small surface velocity amplitudes are predicted for models with masses below 10Msun. For models with M>10Msun, the surface velocity amplitudes are of the order of 10km/s. Moreover we find the occurrence of pulsations for stars in the mass range 9-14Msun, while mass loss appears to stabilize the more massive WR stars. We confront our results with observationally derived line variabilities of 17 WN stars. The data suggest variability to occur for stars above 10Msun, which is increasing linearly with mass above this value, in agreement with our results. We further find some of our models to be unstable to radial pulsations, and predict local magnetic fields of the order of hundreds of Gauss in WR stars more massive than 10Msun. Our study relates the surface velocity fluctuations induced by sub-surface convection to the formation of clumping in the inner part of the wind. From this mechanism, we expect a stronger variability in more massive WR stars.

Massive stars form on different scales ranging from large, dispersed OB associations to compact, dense starburst clusters. The complex structure of regions of massive star formation, and the involved short timescales provide a challenge for our understanding of their birth and early evolution. As one of the most massive and luminous star-forming region in our Galaxy, W49 is the ideal place to study the formation of the most massive stars. By classifying the massive young stars deeply embedded into the molecular cloud of W49, we aim to investigate and trace the star formation history of this region. We analyse near-infrared $K$-band spectroscopic observations of W49 from LBT/LUCI combined with $JHK$ images obtained with NTT/SOFI and LBT/LUCI. Based on $JHK$-band photometry and K-band spectroscopy the massive stars are placed in a Hertzsprung Russell diagram. By comparison with evolutionary models, their age and hence the star formation history of W49 can be investigated. Fourteen O type stars as well as two young stellar objects (YSOs) are identified by our spectroscopic survey. Eleven O-stars are main sequence stars with subtypes ranging from O3 to O9.5, with masses ranging from $\sim 20$ $M_{\odot}$ to $\sim 120$ $M_{\odot}$. Three of the O-stars show strong wind features, and are considered to be Of-type supergiants with masses beyond 100 $M_{\odot}$ . The two YSOs show CO emission, indicative for the presence of circumstellar disks in the central region of the massive cluster. The age of the cluster is estimated as $\sim1.5$ Myr, with star formation still ongoing in different parts of the region. The ionising photons from the central massive stars have not yet cleared the molecular cocoon surrounding the cluster. W49 is comparable to extragalactic star-forming regions and provides us with an unique possibility to study a starburst in detail.

We report the results of spectrophotometric observations of the massive star MN18 revealed via discovery of a bipolar nebula around it with the Spitzer Space Telescope. Using the optical spectrum obtained with the Southern African Large Telescope, we classify this star as B1 Ia. The evolved status of MN18 is supported by the detection of nitrogen overabundance in the nebula, which implies that it is composed of processed material ejected by the star. We analysed the spectrum of MN18 by using the code CMFGEN, obtaining a stellar effective temperature of \approx 21 kK. The star is highly reddened, E(B-V)\approx 2 mag. Adopting an absolute visual magnitude of M_V=-6.8\pm0.5 (typical of B1 supergiants), MN18 has a luminosity of log L/Lsun \approx 5.42\pm0.30, a mass-loss rate of \approx (2.8-4.5)\times10^{-7} Msun/yr, and resides at a distance of \approx 5.6^{+1.5} _{-1.2} kpc. We discuss the origin of the nebula around MN18 and compare it with similar nebulae produced by other blue supergiants in the Galaxy (Sher 25, HD 168625, [SBW2007] 1) and the Large Magellanic Cloud (Sk-69 202). The nitrogen abundances in these nebulae imply that blue supergiants can produce them from the main sequence stage up to the pre-supernova stage. We also present a K-band spectrum of the candidate luminous blue variable MN56 (encircled by a ring-like nebula) and report the discovery of an OB star at \approx 17 arcsec from MN18. The possible membership of MN18 and the OB star of the star cluster Lynga 3 is discussed.

The evolution of massive stars even on the main sequence is not yet well understood. Due to the steep mass-luminosity relation, massive main sequence stars become very 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 investigate their observational consequences. We use the grids of detailed stellar models by Brott et al. (2011) and Koehler et al. (2015), to investigate the envelope properties of core hydrogen burning massive stars. We find that at the stellar surface, the Eddington limit is almost never reached, even for stars up to 500 Msun. When an appropriate Eddington limit is defined locally in the stellar envelope, most stars more massive than 40 Msun 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 Lsun i.e. stars above 40 Msun 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.

We report the discovery of a circular mid-infrared shell around the emission-line star Wray 16-137 using archival data of the Spitzer Space Telescope. Follow-up optical spectroscopy of Wray 16-137 with the Southern African Large Telescope revealed a rich emission spectrum typical of the classical luminous blue variables (LBVs) like P Cygni. Subsequent spectroscopic and photometric observations showed drastic changes in the spectrum and brightness during the last three years, meaning that Wray 16-137 currently undergoes an S Dor-like outburst. Namely, we found that the star has brightened by \approx 1 mag in the V and I_c bands, while its spectrum became dominated by Fe ii lines. Taken together, our observations unambiguously show that Wray 16-137 is a new member of the family of Galactic bona fide LBVs.

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 Hii 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 to denser WNh 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 stellar and wind parameters. For the first time, we observationally resolve the transition between optically thin O and optically thick WNh star winds. Our results suggest the existence of a kink between both mass-loss regimes, in agreement with recent MC simulations. For the optically thick regime, we confirm the steep dependence on the Eddington factor from previous theoretical and observational studies. The transition occurs on the MS near a luminosity of 10^6.1Lsun, or a mass of 80...90Msun. Above this limit, we find that - even when accounting for moderate wind clumping (with f = 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. 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 commonly used extinction laws of Cardelli et al. (1989) have limitations that, among other issues, hamper the determination of the effective temperatures of O and early B stars from optical+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+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+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.

We compile the first comprehensive census of hot luminous stars in the 30 Doradus (30 Dor) star forming region of the LMC. The census extends to a radius of 10arcmin (150pc) from the central cluster, R136. Stars were selected photometrically and combined with the latest spectral types. 1145 candidate hot luminous stars were identified of which >700 were considered genuine early type stars that contribute to feedback. We assess the spectroscopic completeness to be 85% in outer regions (>5pc) but fall to 35% in the vicinity of R136, giving a total of 500 hot luminous stars with spectroscopy. Stellar calibrations and models were used to obtain their physical parameters before integrated values were compared to global observations and the population synthesis code, Starburst99. The 31 W-R and Of/WN stars made large contributions to the total ionising and wind luminosities of ~40% and ~50%, respectively. Stars with Minit>100Msun 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 underestimated the ionising and wind luminosities of R136 by factors of ~2 and ~9, respectively. The census inferred a SFR of 0.073+/-0.04Msun/yr for 30 Dor, typically higher than results from popular SFR calibrations. However, it remained consistent with a far-UV luminosity tracer and a combined Halpha and mid-infrared tracer, but only after correcting for Halpha extinction. The global ionising output exceeded measurements from the associated gas and dust, suggesting ~6(+55/-6)% of ionising photons escape the region. When studying very luminous star forming regions, it is therefore essential to include the most massive stars to ensure a reliable energy budget. If 30 Dor is typical of other large star forming regions, estimates of the SFR will be underpredicted if this escape fraction is not accounted for.(abridged)

The development of infrared observational facilities has revealed a number of massive stars in obscured environments throughout the Milky Way and beyond. The determination of their stellar and wind properties from infrared diagnostics is thus required to take full advantage of the wealth of observations available in the near and mid infrared. However, the task is challenging. This session addressed some of the problems encountered and showed the limitations and successes of infrared studies of massive stars.

Whether massive stars can occasionally form in relative isolation or if they require a large cluster of lower-mass stars around them is a key test in the differentiation of star formation theories as well as how the initial mass function of stars is sampled. Previous attempts to find O-type stars that formed in isolation were hindered by the possibility that such stars are merely runaways from clusters, i.e., their current isolation does not reflect their birth conditions. We introduce a new method to find O-type stars that are not affected by such a degeneracy. Using the VLT-FLAMES Tarantula Survey and additional high resolution imaging we have identified stars that satisfy the following constraints: 1) they are O-type stars that are not detected to be part of a binary system based on RV time series analysis; 2) they are designated spectral type O7 or earlier ; 3) their velocities are within 1\sigma of the mean of OB-type stars in the 30 Doradus region, i.e. they are not runaways along our line-of-sight; 4) the projected surface density of stars does not increase within 3 pc towards the O-star (no evidence for clusters); 5) their sight lines are associated with gaseous and/or dusty filaments in the ISM, and 6) if a second candidate is found in the direction of the same filament with which the target is associated, both are required to have similar velocities. With these criteria, we have identified 15 stars in the 30 Doradus region, which are strong candidates for being high-mass stars that have formed in isolation. Additionally, we employed extensive MC stellar cluster simulations to confirm that our results rule out the presence of clusters around the candidates. Eleven of these are classified as Vz stars, possibly associated with the zero-age main sequence. We include a newly discovered W-R star as a candidate, although it does not meet all of the above criteria.

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.

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 (RV ~4.7), and obtain a surprisingly high luminosity log(L/Lsun) = 6.5 \pm 0.2, corresponding to a present-day mass of ~150Msun. The high effective temperature of 52.2 \pm 2.5kK 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 Wolf-Rayet 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.

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.

We report the discovery that R139 in 30 Doradus is a massive spectroscopic binary system. Multi-epoch optical spectroscopy of R139 was obtained as part of the VLT-FLAMES Tarantula Survey, revealing a double-lined system. The two components are of similar spectral types; the primary exhibits strong C III 4650 emission and is classified as an O6.5 Iafc supergiant, while the secondary is an O6 Iaf supergiant. The radial-velocity variations indicate a highly eccentric orbit with a period of 153.9 days. Photometry obtained with the Faulkes Telescope South shows no evidence for significant variability within an 18 month period. The orbital solution yields lower mass limits for the components of M1sin^3 i = 78 \pm 8 Msun and M2sin^3 i = 66 \pm 7 Msun. As R139 appears to be the most massive binary system known to contain two evolved Of supergiants, it will provide an excellent test for atmospheric and evolutionary models.

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.

We perform a photometric analysis on deep, seeing-limited near-IR VLT images of the region NGC 346/N66 and a nearby control field of the SMC in order to locate the centres of active high- and intermediate-mass star formation through the identification of near-IR bright objects as candidate stellar sources under formation.We use archived imaging data obtained with the high-resolution camera ISAAC at VLT of NGC 346/N66 and we construct the near-IR colour-magnitude (CMD) and colour-colour diagrams (C-CD) of all detected sources. We investigate the nature of all stellar populations in the observed CMDs, and we identify all stellar sources that show significant near-IR excess emission in the observed C-CD. We select, thus, the best candidates for being young stellar sources. Based on their near-IR colours we select 263 candidate young stellar sources. This sample comprises a variety of objects such as intermediate-mass PMS and Herbig Ae/Be stars and possibly massive YSOs, providing original near-IR colours for them. The spatial distribution of the selected candidate sources shows that they are located along the dusty filamentary structures of N66 seen in mid- and far-IR dust emission and agrees very well with that of previously detected candidate YSOs and PMS stars. Our study provides an original accurate set of near-IR colours for candidate young stellar sources. While this adds significant information on the star formation process in NGC 346/N66, it is still necessary to reveal the exact nature of these objects, a process that requires the combination of multi-wavelength data, so that complete Spectral Energy Distributions can be constructed for individual sources. This would be a quite important follow-up study to that presented here.