The VLT-FLAMES Tarantula Survey has observed hundreds of O-type stars in the 30 Doradus region of the Large Magellanic Cloud (LMC). We study the properties of 105 apparently single O-type dwarfs. To determine stellar and wind parameters, we used the IACOB-GBAT package, an automatic procedure based on a large grid of atmospheric models calculated with the FASTWIND code. In addition to classical techniques, we applied the Bayesian BONNSAI tool to estimate evolutionary masses. We provide a new calibration of effective temperature vs. spectral type for O-type dwarfs in the LMC, based on our homogeneous analysis of the largest sample of such objects to date and including all spectral subtypes. Good agreement with previous results is found, although the sampling at the earliest subtypes could be improved. Rotation rates and helium abundances are studied in an evolutionary context. We find that most of the rapid rotators (vsini higher than 300 km/s ) in our sample have masses below 25 MSun and intermediate rotation-corrected gravities (log gc between 3.9 and 4.1). Such rapid rotators are scarce at higher gravities (i.e. younger ages) and absent at lower gravities (larger ages). This is not expected from theoretical evolutionary models, and does not appear to be due to a selection bias in our sample. We compare the estimated evolutionary and spectroscopic masses, finding a trend that the former is higher for masses below 20 MSun. This can be explained as a consequence of limiting our sample to the O-type stars, and we see no compelling evidence for a systematic mass discrepancy. For most of the stars in the sample we were unable to estimate the wind-strength parameter (hence mass-loss rates) reliably, particularly for objects with luminosity lower than logL/LSun about 5.1. Ultraviolet spectroscopy is needed to undertake a detailed investigation of the wind properties of these dwarfs.

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.

We introduce a HST/STIS stellar census of R136a, the central ionizing star cluster of 30 Doradus. We present low resolution far-ultraviolet STIS/MAMA spectroscopy of R136 using 17 contiguous 52x0.2 arcsec slits which together provide complete coverage of the central 0.85 parsec (3.4 arcsec). We provide spectral types of 90% of the 57 sources brighter than m_F555W = 16.0 mag within a radius of 0.5 parsec of R136a1, plus 8 additional nearby sources including R136b (O4\,If/WN8). We measure wind velocities for 52 early-type stars from CIV 1548-51, including 16 O2-3 stars. For the first time we spectroscopically classify all Weigelt & Baier members of R136a, which comprise three WN5 stars (a1-a3), two O supergiants (a5-a6) and three early O dwarfs (a4, a7, a8). A complete Hertzsprung-Russell diagram for the most massive O stars in R136 is provided, from which we obtain a cluster age of 1.5+0.3_-0.7 Myr. In addition, we discuss the integrated ultraviolet spectrum of R136, and highlight the central role played by the most luminous stars in producing the prominent HeII 1640 emission line. This emission is totally dominated by very massive stars with initial masses above ~100 Msun. The presence of strong HeII 1640 emission in the integrated light of very young star clusters (e.g A1 in NGC 3125) favours an initial mass function extending well beyond a conventional upper limit of 100 Msun. We include montages of ultraviolet spectroscopy for LMC O stars in the Appendix. Future studies in this series will focus on optical STIS/CCD medium resolution observations.

Aims: We study the evolution of stellar rotation and wind properties for low-mass main-sequence stars. Our aim is to use rotational evolution models to constrain the mass loss rates in stellar winds and to predict how their properties evolve with time on the main-sequence. Methods: We construct a rotational evolution model that is driven by observed rotational distributions of young stellar clusters. Fitting the free parameters in our model allows us to predict how wind mass loss rate depends on stellar mass, radius, and rotation. We couple the results to the wind model developed in Paper I of this series to predict how wind properties evolve on the main-sequence. Results: We estimate that wind mass loss rate scales with stellar parameters as $\dot{M}_\star \propto R_\star^2 \Omega_\star^{1.33} M_\star^{-3.36}$. We estimate that at young ages, the solar wind likely had a mass loss rate that is an order of magnitude higher than that of the current solar wind. This leads to the wind having a higher density at younger ages; however, the magnitude of this change depends strongly on how we scale wind temperature. Due to the spread in rotation rates, young stars show a large range of wind properties at a given age. This spread in wind properties disappears as the stars age. Conclusions: There is a large uncertainty in our knowledge of the evolution of stellar winds on the main-sequence, due both to our lack of knowledge of stellar winds and the large spread in rotation rates at young ages. Given the sensitivity of planetary atmospheres to stellar wind and radiation conditions, these uncertainties can be significant for our understanding of the evolution of planetary environments.

Aims: We develop a method for estimating the properties of stellar winds for low-mass main-sequence stars between masses of 0.4 and 1.1 solar masses at a range of distances from the star. Methods: We use 1D thermal pressure driven hydrodynamic wind models run using the Versatile Advection Code. Using in situ measurements of the solar wind, we produce models for the slow and fast components of the solar wind. We consider two radically different methods for scaling the base temperature of the wind to other stars: in Model A, we assume that wind temperatures are fundamentally linked to coronal temperatures, and in Model B, we assume that the sound speed at the base of the wind is a fixed fraction of the escape velocity. In Paper II of this series, we use observationally constrained rotational evolution models to derive wind mass loss rates. Results: Our model for the solar wind provides an excellent description of the real solar wind far from the solar surface, but is unrealistic within the solar corona. We run a grid of 1200 wind models to derive relations for the wind properties as a function of stellar mass, radius, and wind temperature. Using these results, we explore how wind properties depend on stellar mass and rotation. Conclusions: Based on our two assumptions about the scaling of the wind temperature, we argue that there is still significant uncertainty in how these properties should be determined. Resolution of this uncertainty will probably require both the application of solar wind physics to other stars and detailed observational constraints on the properties of stellar winds. In the final section of this paper, we give step by step instructions for how to apply our results to calculate the stellar wind conditions far from the stellar surface.

Powerful telescopes equipped with multi-fibre or integral field spectrographs combined with detailed models of stellar atmospheres and automated fitting techniques allow for the analysis of large number of stars. These datasets contain a wealth of information that require new analysis techniques to bridge the gap between observations and stellar evolution models. To that end, we develop BONNSAI (BONN Stellar Astrophysics Interface), a Bayesian statistical method, that is capable of comparing all available observables simultaneously to stellar models while taking observed uncertainties and prior knowledge such as initial mass functions and distributions of stellar rotational velocities into account. BONNSAI can be used to (1) determine probability distributions of fundamental stellar parameters such as initial masses and stellar ages from complex datasets, (2) predict stellar parameters that were not yet observationally determined and (3) test stellar models to further advance our understanding of stellar evolution. An important aspect of BONNSAI is that it singles out stars that cannot be reproduced by stellar models through $\chi^{2}$ hypothesis tests and posterior predictive checks. BONNSAI can be used with any set of stellar models and currently supports massive main-sequence single star models of Milky Way and Large and Small Magellanic Cloud composition. We apply our new method to mock stars to demonstrate its functionality and capabilities. In a first application, we use BONNSAI to test the stellar models of Brott et al. (2011a) by comparing the stellar ages inferred for the primary and secondary stars of eclipsing Milky Way binaries. Ages are determined from dynamical masses and radii that are known to better than 3%. We find that the stellar models reproduce the Milky Way binaries well. BONNSAI is available through a web-interface at http://www.astro.uni-bonn.de/stars/bonnsai.

AIMS. We test the hypothesis of O Vz stars (characterized by having HeII4686 stronger in absorption than other He lines in their blue-violet spectra) being at a younger evolutionary stage than are normal O-type dwarfs. METHODS. We have performed a quantitative spectroscopic analysis of a sample of 38 O Vz and 46 O V stars, identified by the VLT-FLAMES Tarantula Survey in the 30 Doradus region of the Large Magellanic Cloud (LMC). We obtained the stellar and wind parameters of both samples using the FASTWIND stellar atmosphere code and the IACOB-GBAT grid-based automatic tool. In the framework of a differential study, we compared the physical and evolutionary properties of both samples, regarding Teff, logg, logQ and logL. We also investigated the predictions of the FASTWIND code about the O Vz phenomenon. RESULTS. We find a differential distribution of objects in terms of effective temperature, with O Vz stars dominant at intermediate values. The O Vz stars in 30 Doradus tend to be younger and less luminous, and they have weaker winds than the O V stars, but we also find examples with ages of 2-4 Myr and with luminosities and winds that are similar to those of normal O dwarfs. Moreover, the O Vz stars do not appear to have higher gravities than the O V stars. In addition to effective temperature and wind strength, our FASTWIND predictions indicate how important it is to take other stellar parameters (gravity and projected rotational velocity) into account for correctly interpreting the O Vz phenomenon. CONCLUSIONS. In general, the O Vz stars appear to be on or very close to the ZAMS, but there are some examples where the Vz classification does not necessarily imply extreme youth. In particular, the presence of O Vz stars in our sample at more evolved phases than expected is likely a consequence of modest O-star winds owing to the low-metallicity environment of the LMC.

Aims: Projected rotational velocities (\vsini) have been estimated for 334 targets in the VLT-FLAMES Tarantula survey that do not manifest significant radial velocity variations and are not supergiants. They have spectral types from approximately O9.5 to B3. The estimates have been analysed to infer the underlying rotational velocity distribution, which is critical for understanding the evolution of massive stars. Methods: Projected rotational velocities were deduced from the Fourier transforms of spectral lines, with upper limits also being obtained from profile fitting. For the narrower lined stars, metal and non-diffuse helium lines were adopted, and for the broader lined stars, both non-diffuse and diffuse helium lines; the estimates obtained using the different sets of lines are in good agreement. The uncertainty in the mean estimates is typically 4% for most targets. The iterative deconvolution procedure of Lucy has been used to deduce the probability density distribution of the rotational velocities. Results: Projected rotational velocities range up to approximately 450 \kms and show a bi-modal structure. This is also present in the inferred rotational velocity distribution with 25% of the sample having $0\leq$\ve$\leq$100\,\kms and the high velocity component having \ve$\sim 250$\,\kms. There is no evidence from the spatial and radial velocity distributions of the two components that they represent either field and cluster populations or different episodes of star formation. Be-type stars have also been identified. Conclusions: The bi-modal rotational velocity distribution in our sample resembles that found for late-B and early-A type stars. While magnetic braking appears to be a possible mechanism for producing the low-velocity component, we can not rule out alternative explanations.

Rotational mixing in massive main sequence stars is predicted to monotonically increase their surface nitrogen abundance with time. We use this effect to design a method for constraining the age and the inclination angle of massive main sequence stars, given their observed luminosity, effective temperature, projected rotational velocity and surface nitrogen abundance. This method relies on stellar evolution models for different metallicities, masses and rotation rates. We use the population synthesis code STARMAKER to show the range of applicability of our method. We apply this method to 79 early B-type main sequence stars near the LMC clusters NGC 2004 and N 11 and the SMC clusters NGC 330 and NGC 346. From all stars within the sample, 17 were found to be suitable for an age analysis. For ten of them, which are rapidly rotating stars without a strong nitrogen enhancement, it has been previously concluded that they did not evolve as rotationally mixed single stars. This is confirmed by our analysis, which flags the age of these objects as highly discrepant with their isochrone ages. For the other seven stars, their nitrogen and isochrone ages are found to agree within error bars, what validates our method. Constraints on the inclination angle have been derived for the other 62 stars,with the implication that the nitrogen abundances of the SMC stars, for which mostly only upper limits are known, fall on average significantly below those limits. Nitrogen chronology is found to be a new useful tool for testing stellar evolution and to constrain fundamental properties of massive main sequence stars. A web version of this tool is provided.

Rotation has a number of important effects on the evolution of stars. Apart from structural changes because of the centrifugal force, turbulent mixing and meridional circulation caused by rotation can dramatically affect a star's chemical evolution. This leads to changes in the surface temperature and luminosity as well as modifying its lifetime. Observationally rotation decreases the surface gravity, causes enhanced mass loss and leads to surface abundance anomalies of various chemical isotopes. The replication of these physical effects with simple stellar evolution models is very difficult and has resulted in the use of numerous different formulations to describe the physics. Using stellar evolution calculations based on several physical models we discuss the features of the resulting simulated stellar populations which can help to distinguish between the models.

We present a spectroscopic analysis of an extremely rapidly rotating late O-type star, VFTS102, observed during a spectroscopic survey of 30 Doradus. VFTS102 has a projected rotational velocity larger than 500\kms\ and probably as large as 600\kms; as such it would appear to be the most rapidly rotating massive star currently identified. Its radial velocity differs by 40\kms\ from the mean for 30 Doradus, suggesting that it is a runaway. VFTS102 lies 12 pcs from the X-ray pulsar PSR J0537-6910 in the tail of its X-ray diffuse emission. We suggest that these objects originated from a binary system with the rotational and radial velocities of VFTS102 resulting from mass transfer from the progenitor of PSR J0537-691 and the supernova explosion respectively.

This is the second paper in a series aiming at the analysis of nitrogen abundances in O-type stars, to enable further constraints on the early evolution of massive stars. We investigate the NIV lambda4058 emission line formation, provide nitrogen abundances for a substantial O-star sample in the LMC, and compare our findings with recent predictions from stellar evolutionary models. Stellar and wind parameters are determined by line profile fitting of hydrogen, helium and nitrogen lines, based on synthetic spectra calculated by FASTWIND. We derive nitrogen abundances for 20 O- and 5 B-stars, by analyzing all nitrogen lines present in the available optical spectra. The dominating process responsible for emission at NIV lambda4058 in O-stars is the strong depopulation of the lower level of the transition, which increases as a function of Mdot. Unlike the NIII triplet emission, resonance lines do not play a role for typical mass-loss rates and below. The bulk of our sample O-stars seems to be strongly nitrogen-enriched, and a clear correlation of nitrogen and helium enrichment is found. By comparing the nitrogen abundances as a function of vsini ('Hunter-plot') with tailored evolutionary calculations, we identify a considerable number of highly enriched objects at low rotation. Due to the low initial abundance, the detection of strong Nitrogen enrichment in the bulk of O-stars indicates that efficient mixing takes place already during the very early phases of stellar evolution of LMC O-stars.

Aims. We compare the predictions of evolutionary models for early-type stars with atmospheric parameters, projected rotational velocities and nitrogen abundances estimated for a sample of Be-type stars. Our targets are located in 4 fields centred on the Large Magellanic Cloud cluster: NGC 2004 and the N 11 region as well as the Small Magellanic Cloud clusters: NGC 330 and NGC 346. Methods. Atmospheric parameters and photospheric abundances have been determined using the non-LTE atmosphere code tlusty. Effective temperature estimates were deduced using three different methodologies depending on the spectral features observed; in general they were found to yield consistent estimates. Gravities were deduced from Balmer line profiles and microturbulences from the Si iii spectrum. Additionally the contributions of continuum emission from circumstellar discs were estimated. Given its importance in constraining stellar evolutionary models, nitrogen abundances (or upper limits) were deduced for all the stars analysed. Results. Our nitrogen abundances are inconsistent with those predicted for targets spending most of their main sequence life rotating near to the critical velocity. This is consistent with the results we obtain from modelling the inferred rotational velocity distribution of our sample and of other investigators. We consider a number of possibilities to explain the nitrogen abundances and rotational velocities of our Be-type sample.

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 present a dense grid of evolutionary tracks and isochrones of rotating massive main-sequence stars. We provide three grids with different initial compositions tailored to compare with early OB stars in the Small and Large Magellanic Clouds and in the Galaxy. Each grid covers masses ranging from 5 to 60 Msun and initial rotation rates between 0 and about 600 km/s. To calibrate our models we used the results of the VLT-FLAMES Survey of Massive Stars. We determine the amount of convective overshooting by using the observed drop in rotation rates for stars with surface gravities log g < 3.2 to determine the width of the main sequence. We calibrate the efficiency of rotationally induced mixing using the nitrogen abundance determinations for B stars in the Large Magellanic cloud. We describe and provide evolutionary tracks and the evolution of the central and surface abundances. In particular, we discuss the occurrence of quasi-chemically homogeneous evolution, i.e. the severe effects of efficient mixing of the stellar interior found for the most massive fast rotators. We provide a detailed set of isochrones for rotating stars. Rotation as an initial parameter leads to a degeneracy between the age and the mass of massive main sequence stars if determined from its observed location in the Hertzsprung-Russell diagram. We show that the consideration of surface abundances can resolve this degeneracy.

Rotational mixing in massive stars is a widely applied concept, with far reaching consequences for stellar evolution. Nitrogen surface abundances for a large and homogeneous sample of massive B-type stars in the LMC were obtained by the VLT-FLAMES Survey of Massive Stars. This sample is the first covering a broad range of projected stellar rotational velocities, with a large enough sample of high quality data to allow for a statistically significant analysis. We use the sample to provide the first rigorous test of the theory of rotational mixing in massive stars. We calculated a grid of stellar evolution models, using the FLAMES sample to calibrate some of the uncertain mixing processes. We developed a new population-synthesis code, which uses this grid to simulate a large population of stars with masses, ages and rotational velocity distributions consistent with those from the FLAMES sample. The synthesized population is then filtered by the selection effects in the observed sample, to enable a direct comparison between the empirical results and theoretical predictions. Our simulations reproduce the fraction of stars without significant nitrogen enrichment. The predicted number of rapid rotators with enhanced nitrogen is about twice as large as found observationally. Furthermore, a group of stars consisting of slowly rotating, nitrogen-enriched objects and another consisting of rapidly rotating un-enriched objects can not be reproduced by our single-star population synthesis. Additional physical processes appear to be required to understand the population of massive main-sequence stars from the FLAMES sample.We discuss the possible role of binary stars and magnetic fields in the interpretation of our results. We find that the population of slowly rotating nitrogen-enriched stars is unlikely produced via mass transfer and subsequent tidal spin-down in close binary systems

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.

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.

We highlight the role of the light elements (Li, Be, B) in the evolution of massive single and binary stars, which is largely restricted to a diagnostic value, and foremost so for the element boron. However, we show that the boron surface abundance in massive early type stars contains key information about their foregoing evolution which is not obtainable otherwise. In particular, it allows to constrain internal mixing processes and potential previous mass transfer event for binary stars (even if the companion has disappeared). It may also help solving the mystery of the slowly rotating nitrogen-rich massive main sequence 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.

We study the convection zones in the outer envelope of hot massive stars which are caused by opacity peaks associated with iron and helium ionization. We determine the occurrence and properties of these convection zones as function of the stellar parameters. We then confront our results with observations of OB stars. A stellar evolution code is used to compute a grid of massive star models at different metallicities. In these models, the mixing length theory is used to characterize the envelope convection zones. We find the iron convection zone (FeCZ) to be more prominent for lower surface gravity, higher luminosity and higher initial metallicity. It is absent for luminosities below about $10^{3.2}\Lsun$, $10^{3.9}\Lsun$, and $10^{4.2}$\Lsun$ for the Galaxy, LMC and SMC, respectively. We map the strength of the FeCZ on the Hertzsprung-Russell diagram for three metallicities, and compare this with the occurrence of observational phenomena in O stars: microturbulence, non-radial pulsations, wind clumping, and line profile variability. The confirmation of all three trends for the FeCZ as function of stellar parameters by empirical microturbulent velocities argues for a physical connection between sub-photospheric convective motions and small scale stochastic velocities in the photosphere of O- and B-type stars. We further suggest that clumping in the inner parts of the winds of OB stars could be caused by the same mechanism, and that magnetic fields produced in the FeCZ could appear at the surface of OB stars as diagnosed by discrete absorption components in ultraviolet absorption lines.

Models of rotating single stars can successfully account for a wide variety of observed stellar phenomena, such as the surface enhancements of N and He. However, recent observations have questioned the idea that rotational mixing is the main process responsible for the surface enhancements, emphasizing the need for a strong and conclusive test. We investigate the consequences of rotational mixing for massive main-sequence stars in short-period binaries. In these systems the tides spin up the stars to rapid rotation. We use a state-of-the-art stellar evolution code including the effect of rotational mixing, tides, and magnetic fields. We discuss the surface abundances expected in massive close binaries (M1~20 solar masses) and we propose using such systems to test the concept of rotational mixing. As these short-period binaries often show eclipses, their parameters can be determined with high accuracy, allowing for a direct comparison with binary evolution models. In more massive close systems (M1~50 solar masses, Porb<~2 days) we find that helium is efficiently mixed throughout the envelope. The star remains blue and compact during the main-sequence phase. It stays within its Roche lobe while it gradually becomes a helium star. It is the less massive star, in which the effects of rotational mixing are less pronounced, which fills its Roche lobe first. We propose that this evolution path provides an alternative channel for the formation of tight Wolf-Rayet binaries with a main-sequence companion and might explain massive black hole binaries such as the intriguing system M33 X-7.

We have previously analysed the spectra of 135 early B-type stars in the LMC and found several groups of stars that have chemical compositions that conflict with the theory of rotational mixing. Here we extend this study to Galactic and SMC metallicities with the analysis of ~50 Galactic and ~100 SMC early B-type stars with rotational velocities up to ~300km/s. The surface nitrogen abundances are utilised as a probe of the mixing process. In the SMC, we find a population of slowly rotating nitrogen-rich stars amongst the early B type core-hydrogen burning stars, similar to the LMC. In the Galactic sample we find no significant enrichment amongst the core hydrogen-burning stars, which appears to be in contrast with the expectation from both rotating single-star and close binary evolution models. However, only a small number of the rapidly rotating stars have evolved enough to produce a significant nitrogen enrichment, and these may be analogous to the non-enriched rapid rotators previously found in the LMC sample. Finally, in each metallicity regime, a population of highly enriched supergiants is observed, which cannot be the immediate descendants of core-hydrogen burning stars. Their abundances are, however, compatible with them having gone through a previous red supergiant phase. Together, these observations paint a complex picture of the nitrogen enrichment in massive main sequence and supergiant stellar atmospheres, where age and binarity cause crucial effects. Whether rotational mixing is required to understand our results remains an open question at this time, but could be answered by identifying the true binary fraction in those groups of stars that do not agree with single-star evolutionary models (abridged).

We have used VLT FLAMES data to constrain the uncertain physics of rotational mixing in stellar evolution models. We have simulated a population of single stars and find two groups of observed stars that cannot be explained: (1) a group of fast rotating stars which do not show evidence for rotational mixing and (2) a group of slow rotators with strong N enrichment. Binary effects and fossil magnetic fields may be considered to explain those two groups. We suggest that the element boron could be used to distinguish between rotational mixing and the binary scenario. Our single star population simulations quantify the expected amount of boron in fast and slow rotators and allow a comparison with measured nitrogen and boron abundances in B-stars.

High resolution synthetic stellar libraries are of fundamental importance for the preparation of the Gaia Mission. We present new sets of spectral stellar libraries covering two spectral ranges: 300 --1100 nm at 0.1 nm resolution, and 840 -- 890 nm at 0.001 nm resolution. These libraries span a large range in atmospheric parameters, from super-metal-rich to very metal-poor (-5.0 $<$[Fe/H]$<$+1.0), from cool to hot (\teff=3000--50000 K) stars, including peculiar abundance variations. The spectral resolution, spectral type coverage and number of models represent a substantial improvement over previous libraries used in population synthesis models and in atmospheric analysis.

We present results from the first extensive study of convection zones in the envelopes of hot massive stars, which are caused by opacity peaks associated with iron and helium ionization. These convective regions can be located very close to the stellar surface. Recent observations of microturbulence in massive stars from the VLT-Flames survey are in good agreement with our predictions concerning the occurrence and the strength of sub-surface convection in hot stars. We argue further that convection close to the surface may trigger clumping at the base of the stellar wind of massive stars.

Rotational mixing is a very important but uncertain process in the evolution of massive stars. We propose to use close binaries to test its efficiency. Based on rotating single stellar models we predict nitrogen surface enhancements for tidally locked binaries. Furthermore we demonstrate the possibility of a new evolutionary scenario for very massive (M > 40 solar mass) close (P < 3 days) binaries: Case M, in which mixing is so efficient that the stars evolve quasi-chemically homogeneously, stay compact and avoid any Roche-lobe overflow, leading to very close (double) WR binaries.

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.

We review the role of rotation in massive close binary systems. Rotation has been advocated as an essential ingredient in massive single star models. However, rotation clearly is most important in massive binaries where one star accretes matter from a close companion, as the resulting spin-up drives the accretor towards critical rotation. Here, we explore our understanding of this process, and its observable consequences. When accounting for these consequences, the question remains whether rotational effects in massive single stars are still needed to explain the observations.

Rotation has become an important element in evolutionary models of massive stars, specifically via the prediction of rotational mixing. Here, we study a sample of stars, including rapid rotators, to constrain such models and use nitrogen enrichments as a probe of the mixing process. Chemical compositions (C, N, O, Mg and Si) have been estimated for 135 early B-type stars in the Large Magellanic Cloud with projected rotational velocities up to ~300km/s using a non-LTE TLUSTY model atmosphere grid. Evolutionary models, including rotational mixing, have been generated attempting to reproduce these observations by adjusting the overshooting and rotational mixing parameters and produce reasonable agreement with 60% of our core hydrogen burning sample. We find (excluding known binaries) a significant population of highly nitrogen enriched intrinsic slow rotators vsini less than 50km/s incompatible with our models ~20% of the sample). Furthermore, while we find fast rotators with enrichments in agreement with the models, the observation of evolved (log g less than 3.7dex) fast rotators that are relatively unenriched (a further ~20% of the sample) challenges the concept of rotational mixing. We also find that 70% of our blue supergiant sample cannot have evolved directly from the hydrogen burning main-sequence. We are left with a picture where invoking binarity and perhaps fossil magnetic fields are required to understand the surface properties of a population of massive main sequence stars.

The VLT-Flames Survey for Massive Stars (Evans05,Evans06) provides recise measurements of rotational velocities and nitrogen surface abundances of massive stars in the Magellanic Clouds. Specifically, for the first time, such abundances have been estimated for stars with significant rotational velocities. This extraordinary data set gives us the unique possibility to calibrate rotationally and magnetically induced mixing processes. Therefore, we have computed a grid of stellar evolution models varying in mass, initial rotational velocity and chemical composition. In our models we find that although magnetic fields generated by the Spruit-Taylor dynamo are essential to understand the internal angular momentum transport (and hence the rotational behavior), the corresponding chemical mixing must be neglected to reproduce the observations. Further we show that for low metallicities detailed initial abundances are of prime importance, as solar-scaled abundances may result in significant calibration errors.

(Abridged) We investigate the effects of metallicity on the broad-band photometric colors of late-type giants, and make a comparison of synthetic colors with observed photometric properties of late-type giants over a wide range of effective temperatures (T_eff=3500-4800 K) and gravities (log g=0.0-2.5), at [M/H]=-1.0 and -2.0. The influence of metallicity on the synthetic photometric colors is generally small at effective temperatures above \~3800 K, but the effects grow larger at lower T_eff, due to the changing efficiency of molecule formation which reduces molecular opacities at lower [M/H]. To make a detailed comparison of the synthetic and observed photometric colors of late type giants in the T_eff--color and color--color planes, we derive a set of new T_eff--log g--color relations based on synthetic photometric colors, at [M/H]=-0.5, -1.0, -1.5, and -2.0. While differences between the new T_eff--color relations and those available from the literature are typically well within ~100 K, effective temperatures predicted by the scales based on synthetic colors tend to be slightly higher than those resulting from the T_eff--color relations based on observations, with the offsets up to ~100 K. This is clearly seen both at [M/H]=-1.0 and -2.0, especially in the T_eff--(B-V) and T_eff--(V-K) planes. The consistency between T_eff--log g--color scales based on synthetic colors calculated with different stellar atmosphere codes is very good, with typical differences being well within \Delta T_eff~70 K at [M/H]=-1.0 and \Delta T_eff~40 K at [M/H]=-2.0.

To assess the current status in the theoretical modeling of the spectral properties of late-type giants, we provide a comparison of synthetic photometric colors of late-type giants (calculated with PHOENIX, MARCS and ATLAS model atmospheres) with observations, at [M/H]=0.0 and -2.0. Overall, there is a good agreement between observed and synthetic colors, and synthetic colors and published Teff-color relations, both at [M/H]=0.0 and -2.0. Deviations from the observed trends in Teff-color planes are generally within \pm 150K (or less) in the effective temperature range Teff=3500-4800K. Synthetic colors calculated with different stellar atmosphere models typically agree to ~100K, within a large range of effective temperatures and gravities. Some discrepancies are seen in the Teff-(B-V) plane below Teff~3800K at [M/H]=0.0, due to difficulties in reproducing the 'turn-off' to the bluer colors which is seen in the observed data at Teff~3600K. Note that at [M/H]=-2.0 effective temperatures given by the scale of Alonso et al. (1999) are generally lower than those resulting from other Teff-color relations based both on observed and synthetic colors.

We present new synthetic broad-band photometric colors for late-type giants based on synthetic spectra calculated with the PHOENIX model atmosphere code. The grid covers effective temperatures T_eff=3000...5000K, gravities log g=-0.5...+3.5, and metallicities [M/H]=+0.5...-4.0. We show that individual broad-band photometric colors are strongly affected by model parameters such as molecular opacities, gravity, microturbulent velocity, and stellar mass. Our exploratory 3D modeling of a prototypical late-type giant shows that convection has a noticeable effect on the photometric colors too, as it alters significantly both the vertical and horizontal thermal structures in the outer atmosphere. The differences between colors calculated with full 3D hydrodynamical and 1D model atmospheres are significant (e.g., \Delta(V-K)~0.2 mag), translating into offsets in effective temperature of up to ~70K. For a sample of 74 late-type giants in the Solar neighborhood, with interferometric effective temperatures and broad-band photometry available in the literature, we compare observed colors with a new PHOENIX grid of synthetic photometric colors, as well as with photometric colors calculated with the MARCS and ATLAS model atmosphere codes. (abridged)

We present the results of a set of model atmospheres and synthetic spectra computed with the PHOENIX code. The models cover a range of effective temperatures 2700 K <= T_eff <= 10 000 K, gravities (-0.5 <= log(g)<= 5.5) and metalicities (-4.0 <= [Z/H] <= +0.5). In addition, variations of alpha elements are considered for each metalicity. The models are computed with a homogeneous set of input data in order to allow for direct relative comparison between the models. For example, all models use a mixing length of l/H_p = 2.0. We provide synthetic spectra with a resolution of 0.2 nm from the UV to the infrared for all models. We give a brief overview of the input physics and show illustrative results. All synthetic spectra are available via ftp: ftp.hs.uni-hamburg.de/pub/outgoing/phoenix/GAIA

Late type giants (RGB/AGB stars) will be important tracers of the Galactic morphology and evolution in the framework of Gaia, as they are intrinsically bright and thus can probe distant stellar populations or those obscured by interstellar extinction. A realistic representation of their atmospheres and spectra with stellar atmosphere models is thus of crucial importance, both for the design and optimization of Gaia instruments, as well as the interpretation of provided astrophysical data. Our analysis of synthetic photometric colors of late-type giants based on PHOENIX, MARCS and ATLAS model atmospheres indicates a general agreement between the current theoretical predictions and observations in the framework of stationary 1-D model atmospheres. Presently available models allow temperature determinations of RGB/AGB stars to an accuracy of ~\pm100 K. In an exploratory study we try to quantify possible residual systematic effects due to the approximations made in 1-D models using full 3-D hydrodynamical models. We find that differences in broad-band photometric colors calculated with 1-D and 3-D models are significant, translating to the offsets in effective temperature of up to ~70 K. Clearly, full 3-D hydrodynamical models will help to alleviate such ambiguities in current theoretical modeling. Additionally, they will allow to study new phenomena, to open qualitatively new windows for stellar astrophysics in the Gaia-era.