The VLT-FLAMES Survey of Massive Stars was an ESO Large Programme to understand rotational mixing and stellar mass-loss in different metallicity environments, in order to better constrain massive star evolution. We gathered high-quality spectra of over 800 stars in the Galaxy and in the Magellanic Clouds. A sample of this size is unprecedented, enabled by the first high-resolution, wide-field, multi-object spectrograph on an 8-m telescope. We developed spectral analysis techniques that, in combination with non-LTE, line-blanketed model atmospheres, were used to quantitatively characterise every star. The large sample, combined with the theoretical developments, has produced exciting new insights into the evolution of the most massive stars.

[Abridged] We present a comprehensive study of the metallicity dependence of the mass-loss rates in stationary stellar winds of hot massive stars. Assuming a power-law dependence of mass loss on metallicity, Mdot \propto Z^{m}, and adopting a theoretical relation between the terminal velocity and metallicity, v_inf \propto Z^{0.13} (Leitherer et al.), we find m = 0.83 +/- 0.16 for non-clumped outflows from an analysis of the wind momentum luminosity relation (WLR) for stars more luminous than 10^{5.2} Lsun. Within the errors, this result agrees with the prediction of m = 0.69 +/- 0.10 from Vink et al. For the high luminosity stars we find the mass loss rates to be greater than the predictions, implying wind clumping factors in their line-forming regions of ~4. For lower luminosity stars, the winds are so weak that their strengths cannot be reliably derived from optical lines, and one must rely on analysis of UV lines. In the low-luminosity domain the Galactic WLR is found to be much steeper than expected from theory, leading to a discrepancy between UV mass-loss rates and the predictions by a factor 100 at luminosities of L ~ 10^{4.75} Lsun, the origin of which is unknown. We emphasize that even if the current mass-loss rates of hot luminous stars are overestimated as a result of wind clumping, the degree of clumping would likely to be independent of metallicity, so the scalings derived in this study are expected to remain correct.

[Abridged] We have studied the optical spectra of 28 O- and early B-type stars in the Large Magellanic Cloud, 22 of which are associated with the young star-forming region N11. Stellar parameters are determined using an automated fitting method, combining the stellar atmosphere code FASTWIND with the genetic-algorithm optimisation routine PIKAIA. Results for stars in the LH9 and LH10 associations of N11 are consistent with a sequential star formation scenario, in which activity in LH9 triggered the formation of LH10. Our sample contains four stars of spectral type O2, of which the hottest is found to be ~49-54 kK (cf. ~45-46 kK for O3 stars). The masses of helium-enriched dwarfs and giants are systematically lower than those implied by non-rotating evolutionary tracks. We interpret this as evidence for efficient rotationally-enhanced mixing, leading to the surfacing of primary helium and to an increase of the stellar luminosity. This result is consistent with findings for SMC stars by Mokiem et al. For bright giants and supergiants no such mass-discrepancy is found, implying that these stars follow tracks of modestly (or non-)rotating objects. Stellar mass-loss properties were found to be intermediate to those found in massive stars in the Galaxy and the SMC, and comparisons with theoretical predictions at LMC metallicity yielded good agreement over the luminosity range of our targets, i.e. 5.0 < log L/L(sun) < 6.1.

We have studied the optical spectra of a sample of 31 O- and early B-type stars in the Small Magellanic Cloud, 21 of which are associated with the young massive cluster NGC 346. Stellar parameters are determined using an automated fitting method. Comparison with predictions of stellar evolution that account for stellar rotation does not result in a unique age, though most stars are best represented by an age of 1-3 Myr. The present day v_sini distribution of the 21 dwarf stars in our sample is consistent with an underlying rotational velocity (v_r) distribution that can be characterised by a mean velocity of about 160-190 km/s and an effective half width of 100-150 km/s. The v_r distribution must include a small percentage of slowly rotating stars. If predictions of the time evolution of the equatorial velocity for massive stars within the environment of the SMC are correct, the young age of the cluster implies that this underlying distribution is representative for the initial rotational velocity distribution. The location in the Hertzsprung-Russell diagram of the stars showing helium enrichment is in qualitative agreement with evolutionary tracks accounting for rotation, but not for those ignoring v_r. The mass loss rates of the SMC objects having luminosities of log L/L_sun > 5.4 are in excellent agreement with predictions. However, for lower luminosity stars the winds are too weak to determine M_dot accurately from the optical spectrum. Two of three spectroscopically classified Vz stars from our sample are located close to the theoretical zero age main sequence, as expected.

We present the first automated fitting method for the quantitative spectroscopy of O- and early B-type stars with stellar winds. The method combines the non-LTE stellar atmosphere code FASTWIND from Puls et al. (2005) with the genetic algorithm based optimizing routine PIKAIA from Charbonneau (1995), allowing for a homogeneous analysis of upcoming large samples of early-type stars (e.g. Evans et al. 2005). In this first implementation we use continuum normalized optical hydrogen and helium lines to determine photospheric and wind parameters. We have assigned weights to these lines accounting for line blends with species not taken into account, lacking physics, and/or possible or potential problems in the model atmosphere code. We find the method to be robust, fast, and accurate. Using our method we analysed seven O-type stars in the young cluster Cyg OB2 and five other Galactic stars with high rotational velocities and/or low mass loss rates (including 10lac, zeta Oph, and tau Sco) that have been studied in detail with a previous version of FASTWIND. The fits are found to have a quality that is comparable or even better than produced by the classical ``by eye'' method. We define errorbars on the model parameters based on the maximum variations of these parameters in the models that cluster around the global optimum. Using this concept, for the investigated dataset we are able to recover mass-loss rates down to ~6d-8 Msun/yr to within an error of a factor of two. Comparison of our derived spectroscopic masses with those derived from stellar evolutionary models are in very good agreement, i.e. based on the limited sample that we have studied we do not find indications for a mass discrepancy.

In this paper we have analyzed 25 Galactic O and early B-stars by means of H and K band spectroscopy, with the primary goal to investigate to what extent a lone near-IR spectroscopy is able to recover stellar and wind parameters derived in the optical. Most of the spectra have been taken with SUBARU IRCS, at a resolution of 12,000, and with a very high S/N (200 or better). In order to synthesize the strategic H/He lines, we have used our recent, line-blanketed version of FASTWIND. First we investigated the predicted behaviour of the strategic lines. In contradiction to what one expects from the optical in the O-star regime, almost all photospheric H/HeI/HeII H/K band lines become stronger if the gravity decreases. Concerning H and HeII, this finding is related to the behaviour of Stark broadening as a function of electron density, which in the line cores is different for members of lower (optical) and higher (IR) series. Regarding HeI, the predicted behaviour is due to some subtle NLTE effects resulting in a stronger overpopulation of the lower level when the gravity decreases. For most of our objects, we obtained good fits, except for the line cores of Br_gamma in early O-stars with significant mass-loss, where this discrepancy might be an indirect effect of clumping. After having derived the stellar and wind parameters from the IR, we have compared them to results from previous optical analyses. Overall, the IR results coincide in most cases with the optical ones within the typical errors usually quoted for the corresponding parameters, i.e, an uncertainty in Teff of 5%, in log g of 0.1 dex and in Mdot of 0.2 dex, with lower errors at higher wind densities. Outliers above the 1-sigma level where found in four cases with respect to log g and in two cases for Mdot.

We introduce a new survey of massive stars in the Galaxy and the Magellanic Clouds using the Fibre Large Array Multi-Element Spectrograph (FLAMES) instrument at the Very Large Telescope (VLT). Here we present observations of 269 Galactic stars with the FLAMES-Giraffe Spectrograph (R ~ 25,000), in fields centered on the open clusters NGC 3293, NGC 4755 and NGC 6611. These data are supplemented by a further 50 targets observed with the Fibre-Fed Extended Range Optical Spectrograph (FEROS, R = 48,000). Following a description of our scientific motivations and target selection criteria, the data reduction methods are described; of critical importance the FLAMES reduction pipeline is found to yield spectra that are in excellent agreement with less automated methods. Spectral classifications and radial velocity measurements are presented for each star, with particular attention paid to morphological peculiarities and evidence of binarity. These observations represent a significant increase in the known spectral content of NGC 3293 and NGC 4755, and will serve as standards against which our subsequent FLAMES observations in the Magellanic Clouds will be compared.

We present new or improved methods for calculating NLTE, line-blanketed model atmospheres for hot stars with winds (spectral types A to O), with particular emphasis on a fast performance. These methods have been implemented into a previous, more simple version of the model atmosphere code FASTWIND (Santolaya-Rey et al.1997) and allow to spectroscopically analyze rather large samples of massive stars in a reasonable time-scale, using state-of-the-art physics. We describe our (partly approximate) approach to solve the equations of statistical equilibrium for those elements which are primarily responsible for line-blocking and blanketing, as well as an approximate treatment of the line-blocking itself, which is based on a simple statistical approach using suitable means for line opacities and emissivities. Furthermore, we comment on our implementation of a consistent temperature structure. In the second part, we concentrate on a detailed comparison with results from those two codes which have been used in alternative spectroscopical investigations, namely CMFGEN (Hillier & Miller 1998) and WM-Basic (Pauldrach et al. 2001). All three codes predict almost identical temperature structures and fluxes for lambda > 400 A, whereas at lower wavelengths a number of discrepancies are found. Optical H/He lines as synthesized by FASTWIND are compared with results from CMFGEN, obtaining a remarkable coincidence, except for the HeI singlets in the temperature range between 36,000 to 41,000 K for dwarfs and between 31,000 to 35,000 K for supergiants, where CMFGEN predicts much weaker lines. Consequences due to these discrepancies are discussed.

We use a grid of 30 line-blanketed unified stellar photosphere and wind models for O-type stars; computed with the code CMFGEN in order to evaluate its potential in the near-infrared spectral domain. The grid includes dwarfs, giants and supergiants. We analyse the equivalent width behaviour of the 20 strongest lines of hydrogen and helium in spectral windows that can be observed using ground-based instrumentation and compare the results with observations. Our main findings are that: i) HeI/HeII line ratios in the J, H and K bands correlate well with the optical ratio employed in spectral classification, and can therefore be used to determine the spectral type; ii) in supergiant stars the transition from the stellar photosphere to the wind follows a shallower density gradient than the standard approach followed in our models, which can be mimicked by adopting a lower gravity in our prescription of the density stratification. iii) the Brackett gamma line poses a number of peculiar problems which partly might be related to wind clumping, and iv) the Brackett alpha line is an excellent mass-loss indicator. For the first and last item we provide quantitative calibrations.

We present a systematic study of the effect of metallicity on the stellar spectral energy distribution (SED) of O main sequence (dwarf) stars, focussing on the hydrogen and helium ionizing continua, and on the optical and near-IR lines used for spectral classification. The spectra are based on non-LTE line blanketed atmosphere models with stellar winds calculated using the CMFGEN code of Hillier & Miller (1998). We draw the following conclusions. First, we find that the total number of Lyman photons emitted is almost independent of line blanketing effects and metallicity for a given effective temperature. This is because the flux that is blocked by the forest of metal lines at wavelengths shorter than 600 Angstrom is redistributed mainly within the Lyman continuum. Second, the spectral type, as defined by the ratio of the equivalent widths of HeI 4471 Angstrom and HeII 4542 Angstrom, is shown to depend noticeably on the microturbulent velocity in the atmosphere, on metallicity and, within the luminosity class of dwarfs, on gravity. Third, we confirm the decrease in the effective temperature for a given spectral type due to the inclusion of line blanketing recently found by e.g. Martins et al. (2002). Finally, we find that the SED below ~450 Angstrom is highly dependent on metallicity. This is reflected in the behaviour of nebular fine-structure line ratios such as [NeIII]/[NeII] 15.5/12.8 and [ArIII]/[ArII] 9.0/7.0 micron. This dependence complicates the use of these nebular ratios as diagnostic tools for the effective temperature determination of the ionizing stars in HII regions and for age dating of starburst regions in galaxies.