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.

Massive stars play an important role in the mechanical and chemical evolution of galaxies. Understanding the internal processes of these stars is vital to our understanding of their evolution and eventual end products. Deformations from spherical geometry are common for massive stars; however, the tools that are currently available for the study of these systems are almost exclusively one-dimensional. We present a new spectroscopic analysis tool tailored for massive stars that deviate from spherical symmetry. This code (entitled SPAMMS) is a spectroscopic patch model that takes the three-dimensional surface geometry of the system into account to produce spectral profiles at given phases and orientations. In using the Wilson-Devinney-like code PHOEBE in combination with the non local thermodynamic equilibrium (NLTE) radiative transfer code FASTWIND, we created a three-dimensional mesh that represents the surface geometry of our system and we assigned FASTWIND emergent intensity line profiles to each mesh triangle, which take the local parameters such as temperature, surface gravity, and radius into account. These line profiles were then integrated across the visible surface, where their flux contribution and radial velocity are taken into account, thus returning a final line profile for the visible surface of the system at a given phase. We demonstrate that SPAMMS can accurately reproduce the morphology of observed spectral line profiles for overcontact systems. Additionally, we show how line profiles of rapidly-rotating single stars differ when taking rotational distortion into account, and the effects that these can have on the determined parameters. Finally, we demonstrate the code's ability to reproduce the Rossiter-Mclaughlin and Struve-Sahade effects.

Despite their rarity, massive stars dominate the ecology of galaxies via their strong, radiatively-driven winds throughout their lives and as supernovae in their deaths. However, their evolution and subsequent impact on their environment can be significantly affected by the presence of a magnetic field. While recent studies indicate that about 7% of OB stars in the Milky Way host strong, stable, organised (fossil) magnetic fields at their surfaces, little is known about the fields of very massive stars, nor the magnetic properties of stars outside our Galaxy. We aim to continue searching for strong magnetic fields in a diverse set of massive and very massive stars (VMS) in the Large and Small Magellanic Clouds (LMC/SMC), and we evaluate the overall capability of FORS2 to usefully search for and detect stellar magnetic fields in extra-galactic environments. We have obtained FORS2 spectropolarimetry of a sample of 41 stars, which principally consist of spectral types B, O, Of/WN, WNh, and classical WR stars in the LMC and SMC. Four of our targets are Of?p stars; one of them was just recently discovered. Each spectrum was analysed to infer the longitudinal magnetic field. No magnetic fields were formally detected in our study, although Bayesian statistical considerations suggest that the Of?p star SMC159-2 is magnetic with a dipolar field of the order of 2.4 to 4.4kG. In addition, our first constraints of magnetic fields in VMS provide interesting insights into the formation of the most massive stars in the Universe.

We combine the NLTE spectral analysis of the detached O-type eclipsing binary OGLE-LMC-ECL-06782 with the analysis of the radial velocity curve and light curve to measure an independent distance to the LMC. In our spectral analysis we study composite spectra of the system at quadrature and use the information from radial velocity and light curve about stellar gravities, radii and component flux ratio to derive effective temperature, reddening, extinction and intrinsic surface brightness. We obtain a distance modulus to the LMC of m - M = 18.53 +/- 0.04 mag. This value is 0.05 mag larger than the precision distance obtained recently from the analysis of a large sample of detached, long period late spectral type eclipsing binaries but agrees within the margin of the uncertainties. We also determine the surface brightnesses of the system components and find good agreement with the published surface brightness color relationship. A comparison of the observed stellar parameters with the prediction of stellar evolution based on the MESA stellar evolution code shows reasonable agreement, but requires a reduction of the internal angular momentum transport to match the observed rotational velocities.

Analysis of late O-type stars observed in the Large Magellanic Cloud (LMC) by the VLT-FLAMES Tarantula Survey (VFTS) revealed a discrepancy between the physical properties estimated from model-atmosphere analysis and those expected from their morphological classifications. Here we revisit the analysis of 32 of these puzzling objects using new hydrogen-helium-silicon FASTWIND models and a different fitting approach to re-evaluate their physical properties. Our new analysis confirms that these stars indeed have properties that are typical of late O-type dwarfs. We also present the first estimates of silicon abundances for O-type stars in the 30 Dor clusters NGC 2060 and NGC 2070, with a weighted mean abundance for our sample of 7.05 +/- 0.03. Our values are about 0.20 dex lower than those previously derived for B-type stars in the LMC clusters N 11 and NGC 2004 using TLUSTY models. Various possibilities (e.g. differences in the analysis methods, effects of microturbulence, and real differences between stars in different clusters) were considered to account for these results. We also used our grid of FASTWIND models to reassess the impact of using the Galactic classification criteria for late O-type stars in the LMC by scrutinising their sensitivity to different stellar properties. At the cool edge of the O star regime the HeII 4686/HeI 4713 ratio used to assign luminosity class for Galactic stars can mimic giants or bright giants in the LMC, even for objects with high gravities (log_g > 4.0 dex). We argue that this line ratio is not a reliable luminosity diagnostic for late O-type stars in the LMC, and that the SiIV 4989/HeI4026 ratio is more robust for these types.

Context: Knowledge about hot, massive stars is usually inferred from quantitative spectroscopy. To analyse non-spherical phenomena, the existing 1D codes must be extended to higher dimensions, and corresponding tools need to be developed. Aims: We present a 3D radiative transfer code that is capable of calculating continuum and line scattering problems in the winds of hot stars. By considering spherically symmetric test models, we discuss potential error sources, and indicate advantages and disadvantages by comparing different solution methods. Further, we analyse the UV resonance line formation in the winds of rapidly rotating O stars. Methods: We consider both a (simplified) continuum model including scattering and thermal sources, and a UV resonance line transition approximated by a two-level-atom. We applied the short-characteristics (SC) method, using linear or monotonic B\'ezier interpolations, to solve the equation of radiative transfer on a non-uniform Cartesian grid. To calculate scattering dominated problems, our solution method is supplemented by an accelerated $\Lambda$-iteration scheme. Results: For the spherical test models, the mean relative error of the source function is on the $5-20\,\%$ level, depending on the applied interpolation technique and the complexity of the considered model. All calculated line profiles are in excellent agreement with corresponding 1D solutions. The predicted line profiles from fast rotating stars show a distinct behaviour as a function of rotational speed and inclination. This behaviour is tightly coupled to the wind structure and the description of gravity darkening and stellar surface distortion. Conclusions: Our SC methods are ready to be used for quantitative analyses of UV resonance line profiles. When calculating optically thick continua, both SC methods give reliable results, in contrast to the alternative finite-volume method.

[Abridged] Context: Radiation-driven mass loss plays a key role in the life-cycles of massive stars. However, basic predictions of such mass loss still suffer from significant quantitative uncertainties. Aims: We develop new radiation-driven, steady-state wind models for massive stars with hot surfaces, suitable for quantitative predictions of global parameters like mass-loss and wind-momentum rates. Methods: The simulations presented here are based on a self-consistent, iterative grid-solution to the spherically symmetric, steady-state equation of motion, using full NLTE radiative transfer solutions in the co-moving frame to derive the radiative acceleration. We do not rely on any distribution functions or parametrization for computation of the line force responsible for the wind driving. Results: In this first paper, we present models representing two prototypical O-stars in the Galaxy, one with a higher stellar mass M/Msun=59 and luminosity log10(L/Lsun)= 5.87 and one with a lower M/Msun= 27 and log10(L/Lsun)= 5.1. For these simulations, basic predictions for global mass-loss rates and velocity laws are given, and the influence from additional parameters like wind clumping and microturbulent speeds discussed. A key result is that our mass-loss rates are significantly lower than those predicted by the mass-loss recipes normally included in models of massive-star evolution. Conclusions: Our results support previous suggestions that Galactic O-star mass-loss rates may be overestimated in present-day stellar evolution models, and that new rates thus might be needed. Indeed, future papers in this series will incorporate our new models into such simulations of stellar evolution, extending the very first simulations presented here toward larger grids covering a range of metallicities, B supergiants across the bistability jump, and possibly also Wolf-Rayet stars.

In this paper, we investigate the photometric variability of magnetic O-type stars. Such stars possess oblique, predominantly dipolar magnetic fields that confine their winds roughly axisymmetrically about the magnetic equator, thus forming a magnetosphere. We interpret their photometric variability as phase-dependent magnetospheric occultations. For massive star winds dominated by electron scattering opacity in the optical and NIR, we can compute synthetic light curves from simply knowing the magnetosphere's mass density distribution. We exploit the newly-developed Analytical Dynamical Magnetosphere model (ADM) in order to obtain the predicted circumstellar density structures of magnetic O-type stars. The simplicity in our light curve synthesis model allows us to readily conduct a parameter space study. For validation purposes, we first apply our algorithm to HD 191612, the prototypical Of?p star. Next, we attempt to model the photometric variability of the Of?p-type stars identified in the Magellanic Clouds using OGLE photometry. We evaluate the compatibility of the ADM predictions with the observed photometric variations, and discuss the magnetic field properties that are implied by our modelling.

The massive O4.5 V + O5.5 V binary VFTS 352 in the Tarantula nebula is one of the shortest-period and most massive overcontact binaries known. Recent theoretical studies indicate that some of these systems could ultimately lead to the formation of gravitational waves via black hole binary mergers through the chemically homogeneous evolution pathway. By analyzing ultraviolet-optical phase-resolved spectroscopic data, we aim to constrain atmospheric and wind properties that could be later used to confront theoretical predictions from binary evolution. In particular, surface abundances are powerful diagnostics of the evolutionary status, mass transfer and the internal mixing processes. From a set of 32 VLT/FLAMES visual and 8 HST/COS ultraviolet spectra, we used spectral disentangling to separate the primary and secondary components. Using a genetic algorithm wrapped around the NLTE model atmosphere and spectral synthesis code FASTWIND, we perform an 11-parameter optimization to derive the atmospheric and wind parameters of both components, including the surface abundances of He, C, N, O and Si. We find that both components are hotter than expected compared to single-star evolutionary models indicating that additional mixing processes may be at play. However the derived chemical abundances do not show significant indications of mixing when adopting baseline values typical for the system environment.

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.

Context: State of the art quantitative spectroscopy of OB-stars compares synthetic spectra (calculated by means of 1D, spherically symmetric computer codes) with observations. Certain stellar atmospheres, however, show strong deviations from spherical symmetry, and need to be treated in 3D. Aims: We present a newly developed 3D radiative transfer code, tailored to the solution of the radiation field in rapidly expanding stellar atmospheres. We apply our code to the continuum transfer in wind-ablation models, and to the UV resonance line formation in magnetic winds. Methods: We have used a 3D finite-volume method for the solution of the equation of radiative transfer, to study continuum- and line-scattering problems. Convergence has been accelerated by a non-local approximate Lambda-iteration scheme. Particular emphasis has been put on careful (spherically symmetric) test cases. Results: Typical errors of the source functions, when compared to 1D solutions, are of the order of 10-20 %, and increase for optically thick continua. In circumstellar discs, the radiation temperatures in the (optically thin) transition region from wind to disc are quite similar to corresponding values in the wind. For MHD simulations of dynamical magnetospheres, the line profiles, calculated with our 3D code, agree well with previous solutions using a 3D-SEI method. When compared with profiles resulting from the `analytic dynamical magnetosphere' (ADM) model, significant differences become apparent. Conclusions: Due to similar radiation temperatures in the wind and the transition region to the disc, the same line-strength distribution can be applied within radiation hydrodynamic calculations for circumstellar discs in `accreting high-mass stars'. To properly describe the UV line formation in dynamical magnetospheres, the ADM model needs to be further developed, at least in a large part of the outer wind.

The IACOB and OWN surveys are two ambitious complementary observational projects which have made available a large multi-epoch spectroscopic database of optical high resolution spectra of Galactic massive O-type stars. As a first step in the study of the full sample of (more than 350) O stars surveyed by the IACOB/OWN projects, we have performed the quantitative spectroscopic analysis of a subsample of 128 stars included in the modern grid of O-type standards for spectral classification. We use semi-automatized tools to determine the set of spectroscopic parameters that can be obtained from the optical spectrum of O-type stars. We also benefit from the multi-epoch character of the surveys to perform a spectroscopic variability study of the sample, accounting for spectroscopic binarity and variability of the main wind diagnostic lines. We provide a general overview of the stellar and wind parameters of this reference sample, and updated recipes for the SpT\,--\,Teff/log g calibrations for Galactic O-type stars. We evaluate our semi-automatized analysis strategy with $\sim$40 stars from the literature, and find a good agreement. The agreement between the synthetic spectra associated with fastwind best fitting models and the observed spectra is good for most targets, but 46 stars present a particular behavior of the wind diagnostic lines that cannot be reproduced by our grid of spherically symmetric unclumped models. These are potential targets of interest for more detailed investigations of clumpy winds and/or the existence of additional circumstellar components. Last, our variability study has led to the detection of signatures of spectroscopic binarity in 27\% of the stars and small amplitude radial velocity variations in the photospheric lines of another 30\%. Additionally, 31\% of the investigated stars show variability in the wind diagnostic lines.

[Abridged] Clumping in the radiation-driven winds of hot, massive stars affects the derivation of synthetic observables across the electromagnetic spectrum. We implement a formalism for treating wind clumping - in particular the light-leakage effects associated with a medium that is porous in physical and velocity space - into the global (photosphere+wind) NLTE model atmosphere code FASTWIND. We assume a stochastic, two-component wind consisting of a mixture of optically thick and thin clumps embedded in a rarefied inter-clump medium. We account fully for the reductions in opacity associated with porosity in physical and velocity-space, and for the well-known effect that opacities depending on rho^2 are higher in clumpy winds than in smooth ones of equal mass-loss rate. By formulating our method in terms of suitable mean and effective opacities for the clumpy wind, we are able to compute models with the same speed (~15 min. on a modern laptop) as in previous code-generations. Some first, generic results of the new models include: i) Confirming earlier results that velocity-space porosity is critical for analysis of UV wind lines in O-stars; ii) for the optical Halpha line, optically thick clumping effects are small for O-stars, but potentially very important for late B and A-supergiants; iii) spatial porosity is a marginal effect for absorption of high-energy X-rays in O-stars, as long as the mean-free path between clumps are kept at realistic values; iv) porosity is negligible at typical O-star radio-photosphere radii; v) regarding the wind ionization balance, a general trend is that increased rates of recombination in simulations with optically thin clumps lead to overall lower degrees of ionization than in corresponding smooth models, but that this effect now is counteracted by the increased levels of light-leakage associated with porosity in physical and velocity space.

Massive stars play a key role in the evolution of the Universe. Our goal is to compare observed and predicted properties of single Galactic O stars to identify and constrain uncertain physical parameters and processes in stellar evolution and atmosphere models. We used a sample of 53 objects with spectral types from O3 to O9.7. For 30 of these, we determined the main photospheric and wind parameters, using optical spectroscopy and applying the FASTWIND code. For the remaining objects, literature data, obtained by means of the CMFGEN code, were used instead. The properties of our sample were compared to published predictions based on two grids evolution models that include rotationally induced mixing. Within each luminosity class, we find a close correlation of N surface abundance and luminosity, and a stronger N enrichment in more massive and evolved O stars. Additionally, a correlation of the surface nitrogen and helium abundances is observed. The large number of nitrogen-enriched stars above ~30 solar masses argues for rotationally induced mixing as the most likely explanation. However, none of the considered models can match the observed trends correctly, especially in the high mass regime. We confirm mass discrepancy for objects in the low mass O-star regime. We conclude that the rotationally induced mixing of helium to the stellar surface is too strong in some of the models. We also suggest that present inadequacies of the models to represent the N enrichment in more massive stars with relatively slow rotation might be related to problematic efficiencies of rotational mixing. We are left with a picture in which invoking binarity and magnetic fields is required to achieve a more complete agreement of the observed surface properties of a population of massive main- sequence stars with corresponding evolutionary models.

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.]

All known Galactic Of?p stars have been shown to host strong, organized, magnetic fields. Recently, five Of?p stars have been discovered in the Magellanic Clouds. They posses photometric \citep{Naze} and spectroscopic \citep{Walborn} variability compatible with the Oblique Rotator Model (ORM). However, their magnetic fields have yet to be directly detected. We have developed an algorithm allowing for the synthesis of photometric observables based on the Analytic Dynamical Magnetosphere (ADM) model of \citet{Owocki}. We apply our model to OGLE photometry in order to constrain their magnetic geometries and surface dipole strengths. We predict that the field strengths for some of these candidate extra-Galactic magnetic stars may be within the detection limits of the FORS2 instrument.

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.

Context: Clumping in the radiation-driven winds of hot, massive stars arises naturally due to the strong, intrinsic instability of line-driving (the `LDI'). But LDI wind models have so far mostly been limited to 1D, mainly because of severe computational challenges regarding calculation of the multi-dimensional radiation force. Aims: To simulate and examine the dynamics and multi-dimensional nature of wind structure resulting from the LDI. Methods: We introduce a `pseudo-planar', `box-in-a-wind' method that allows us to efficiently compute the line-force in the radial and lateral directions, and then use this approach to carry out 2D radiation-hydrodynamical simulations of the time-dependent wind. Results: Our 2D simulations show that the LDI first manifests itself by mimicking the typical shell-structure seen in 1D models, but how these shells then quickly break up into complex 2D density and velocity structures, characterized by small-scale density `clumps' embedded in larger regions of fast and rarefied gas. Key results of the simulations are that density-variations in the well-developed wind statistically are quite isotropic and that characteristic length-scales are small; a typical clump size is ~0.01R at 2R, thus resulting also in rather low typical clump-masses ~10^17 g. Overall, our results agree well with the theoretical expectation that the characteristic scale for LDI-generated wind-structure is of order the Sobolev length. We further confirm some earlier results that lateral `filling-in' of radially compressed gas leads to somewhat lower clumping factors in 2D simulations than in comparable 1D models. We conclude by discussing an extension of our method toward rotating LDI wind models that exhibit an intriguing combination of large- and small-scale structure extending down to the wind base.

The determination of chemical abundances constitutes a fundamental requirement for obtaining a complete picture of a star. Particularly in massive stars, CNO abundances are of prime interest, due to the nuclear CNO-cycle and various mixing processes which bring these elements to the surface. We aim at enabling a reliable carbon spectroscopy for our unified NLTE atmosphere code FASTWIND. We develop a new carbon model atom including CII/III/IV/V, and discuss problems related to carbon spectroscopy in O-type stars. We describe different tests to examine the reliability of our implementation, and investigate which mechanisms influence the carbon ionization balance. By comparing with high-resolution spectra from six O-type stars, we check in how far observational constraints can be reproduced by our new carbon line synthesis. Carbon lines are even more sensitive to a variation of temperature, gravity, and mass-loss rate, than hydrogen/helium lines. We are able to reproduce most of the observed lines from our stellar sample, and to estimate those specific carbon abundances which bring the lines from different ions into agreement. For hot dwarfs and supergiants earlier than O7, X-rays from wind-embedded shocks can impact the synthesized line strengths, particularly for CIV, potentially affecting the abundance determination. We have demonstrated our capability to derive realistic carbon abundances by means of FASTWIND, using our recently developed model atom. We found that complex effects can have a strong influence on the carbon ionization balance in hot stars. For a further understanding, the UV range needs to be explored as well. By means of detailed nitrogen and oxygen model atoms available to use, we will be able to perform a complete CNO abundance analysis for larger samples of massive stars, and to provide constraints on corresponding evolutionary models and aspects.

Massive stars, at least $\sim$ 10 times more massive than the Sun, have two key properties that make them the main drivers of evolution of star clusters, galaxies, and the Universe as a whole. On the one hand, the outer layers of massive stars are so hot that they produce most of the ionizing ultraviolet radiation of galaxies; in fact, the first massive stars helped to re-ionize the Universe after its Dark Ages. Another important property of massive stars are the strong stellar winds and outflows they produce. This mass loss, and finally the explosion of a massive star as a supernova or a gamma-ray burst, provide a significant input of mechanical and radiative energy into the interstellar space. These two properties together make massive stars one of the most important cosmic engines: they trigger the star formation and enrich the interstellar medium with heavy elements, that ultimately leads to formation of Earth-like rocky planets and the development of complex life. The study of massive star winds is thus a truly multidisciplinary field and has a wide impact on different areas of astronomy. [...] This detailed review summarises the current knowledge on the theory and observations of winds from massive stars, as well as on observations and accretion processes in wind-fed high mass X-ray binaries. The aim is to combine in the near future all available theoretical diagnostics and observational measurements to achieve a unified picture of massive star winds in isolated objects and in binary systems.

Strong winds from massive stars are a topic of interest to a wide range of astrophysical fields. In High-Mass X-ray Binaries the presence of an accreting compact object on the one side allows to infer wind parameters from studies of the varying properties of the emitted X-rays; but on the other side the accretor's gravity and ionizing radiation can strongly influence the wind flow. Based on a collaborative effort of astronomers both from the stellar wind and the X-ray community, this presentation attempts to review our current state of knowledge and indicate avenues for future progress.

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.

The groundbreaking detection of gravitational waves produced by the inspiralling and coalescence of the black hole (BH) binary GW150914 confirms the existence of "heavy" stellar-mass BHs with masses >25 Msun. Initial modelling of the system by Abbott et al. (2016a) supposes that the formation of black holes with such large masses from the evolution of single massive stars is only feasible if the wind mass-loss rates of the progenitors were greatly reduced relative to the mass-loss rates of massive stars in the Galaxy, concluding that heavy BHs must form in low-metallicity (Z < 0.25-0.5 Zsun) environments. However, strong surface magnetic fields also provide a powerful mechanism for modifying mass loss and rotation of massive stars, independent of environmental metallicity (ud-Doula & Owocki 2002; ud-Doula et al. 2008). In this paper we explore the hypothesis that some heavy BHs, with masses >25 Msun such as those inferred to compose GW150914, could be the natural end-point of evolution of magnetic massive stars in a solar-metallicity environment. Using the MESA code, we developed a new grid of single, non-rotating, solar metallicity evolutionary models for initial ZAMS masses from 40-80 Msun that include, for the first time, the quenching of the mass loss due to a realistic dipolar surface magnetic field. The new models predict TAMS masses that are significantly greater than those from equivalent non-magnetic models, reducing the total mass lost by a strongly magnetized 80 Msun star during its main sequence evolution by 20 Msun. This corresponds approximately to the mass loss reduction expected from an environment with metallicity Z = 1/30 Zsun.

Stellar evolution models of massive stars are very sensitive to the adopted mass-loss scheme. The magnitude and evolution of mass-loss rates significantly affect the main sequence evolution, and the properties of post-main sequence objects, including their rotational velocities. Driven by potential discrepancies between theoretically predicted and observationally derived mass-loss rates in the OB star range, we particularly aim to investigate the response to mass-loss rates that are lower than currently adopted, in parallel with the mass-loss behavior at the "first" bi-stability jump. We perform 1D hydrodynamical model calculations of single $20 - 60 \, M_{\odot}$ Galactic ($Z = 0.014$) stars where the effects of stellar winds are already significant during the main sequence phase. We develop an experimental wind routine to examine the behavior and response of the models under the influence of different mass-loss rates. This observationally guided, simple and flexible wind routine is not a new mass-loss description but a useful tool based on the Wind-momentum Luminosity Relation and other scaling relations, and provides a meaningful base for various tests and comparisons. The main result of this study indicates a dichotomy when accounting for currently debated problems regarding mass-loss rates of hot massive stars. In a fully diffusive approach, and for commonly adopted initial rotational velocities, lower mass-loss rates than theoretically predicted require to invoke an additional source of angular momentum loss (either due to bi-stability braking, or yet unidentified) to brake down surface rotational velocities. [...]

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 aim to provide new empirical clues about macroturbulent spectral line broadening in O- and B-type stars to evaluate its physical origin. We use high-resolution spectra of ~430 stars with spectral types in the range O4-B9 (all luminosity classes). We characterize the line-broadening of adequate diagnostic metal lines using a combined FT and GOF technique. We perform a quantitative spectroscopic analysis of the whole sample using automatic tools coupled with a huge grid of FASTWIND models. We also incorporate quantitative information about line asymmetries to our observational description of the characteristics of the line-profiles, and present a comparison of the shape and type of line-profile variability found in a small sample of O stars and B supergiants with still undefined pulsational properties and B main sequence stars with variable line-profiles. We present a homogeneous and statistically significant overview of the (single snapshot) line-broadening properties of stars in the whole O and B star domain. We find empirical evidence of the existence of various types of non-rotational broadening agents acting in the realm of massive stars. Even though all of them could be quoted and quantified as a macroturbulent broadening from a practical point of view, their physical origin can be different. Contrarily to the early- to late-B dwarfs/giants, which present a mixture of cases in terms of line-profile shape and variability, the whole O-type and B supergiant domain (or, roughly speaking, stars with M_ZAMS > 15 M_sol) is fully dominated by stars with a remarkable non-rotational broadening component and very similar profiles (including type of variability). We provide some examples illustrating how this observational dataset can be used to evaluate scenarios aimed at explaining the existence of sources of non-rotational broadening in massive stars.

We study the impact of stellar winds and supernovae on the multi-phase interstellar medium using three-dimensional hydrodynamical simulations carried out with FLASH. The selected galactic disc region has a size of (500 pc)$^2$ x $\pm$ 5 kpc and a gas surface density of 10 M$_{\odot}$/pc$^2$. The simulations include an external stellar potential and gas self-gravity, radiative cooling and diffuse heating, sink particles representing star clusters, stellar winds from these clusters which combine the winds from indi- vidual massive stars by following their evolution tracks, and subsequent supernova explosions. Dust and gas (self-)shielding is followed to compute the chemical state of the gas with a chemical network. We find that stellar winds can regulate star (cluster) formation. Since the winds suppress the accretion of fresh gas soon after the cluster has formed, they lead to clusters which have lower average masses (10$^2$ - 10$^{4.3}$ M$_{\odot}$) and form on shorter timescales (10$^{-3}$ - 10 Myr). In particular we find an anti-correlation of cluster mass and accretion time scale. Without winds the star clusters easily grow to larger masses for ~5 Myr until the first supernova explodes. Overall the most massive stars provide the most wind energy input, while objects beginning their evolution as B-type stars contribute most of the supernova energy input. A significant outflow from the disk (mass loading $\gtrsim$ 1 at 1 kpc) can be launched by thermal gas pressure if more than 50% of the volume near the disc mid-plane can be heated to T > 3x10$^5$ K. Stellar winds alone cannot create a hot volume-filling phase. The models which are in best agreement with observed star formation rates drive either no outflows or weak outflows.

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.

X-rays/EUV radiation emitted from wind-embedded shocks in hot, massive stars can affect the ionization balance in their outer atmospheres, and can be the mechanism responsible for the production of highly ionized species. To allow for these processes in the context of spectral analysis, we have implemented such emission into our unified, NLTE model atmosphere/spectrum synthesis code FASTWIND. The shock structure and corresponding emission is calculated as a function of user-supplied parameters. We account for a temperature and density stratification inside the post-shock cooling zones, calculated for radiative and adiabatic cooling in the inner and outer wind, respectively. The high-energy absorption of the cool wind is considered by adding important K-shell opacities, and corresponding Auger ionization rates have been included into the NLTE network. We tested and verified our implementation carefully against corresponding results from various alternative model atmosphere codes, and studied the effects from shock emission for important ions from He, C, N, O, Si, and P. Surprisingly, dielectronic recombination turned out to play an essential role for the ionization balance of OIV/OV around Teff = 45,000 K. Finally, we investigated the behavior of the mass absorption coefficient, kappa_nu(r), important in the context of X-ray line formation in massive star winds. In almost all considered cases, direct ionization is of major influence, and Auger ionization significantly affects only NVI and OVI. The approximation of a radially constant kappa_nu is justified for r > 1.2 Rstar and lambda < 18 A, and also for many models at longer wavelengths. To estimate the actual value of this quantity, however, the HeII opacities need to be calculated from detailed NLTE modeling, at least for wavelengths longer than 18 to 20 A, and information on the individual CNO abundances has to be present.

This hexennial report covers the activities of IAU Commission 36 -- 'Theory of Stellar Atmospheres' -- during the years 2009 to 2015, and will be the last report from this Commission, being replaced by Commission C.G5. After outlining the composition of the Organization Committee(s), we list the scientific meetings held between 2009 and 2015 that were of relevance for our Commission members, and comment on the establishment and objectives of the new Commission C.G5 ('Stellar and Planetary Atmospheres') within the re-structuring process of the IAU. In the main part of the report, we briefly review specific contributions and achievements within our research field during the last six years, concentrating on the theoretical aspect, and dividing between late-type and massive star atmospheres. We also provide a more general overview of primary research areas, and finish our report with a collection of useful web links.

We take advantage of more than 10 years of monitoring of the eclipsing HMXB systems LMC X-4, Cen X-3, 4U 1700-377, 4U 1538-522, SMC X-1, IGR J18027-2016, Vela X-1, IGR J17252-3616, XTE J1855-026, and OAO 1657-415 with the ASM on-board RXTE and ISGRI on-board INTEGRAL to update their ephemeris. These results are used to refine previous measurements of the orbital period decay of all sources (where available) and provide the first accurate values of the apsidal advance in Vela X-1 and 4U 1538-522. Updated values for the masses of the neutron stars hosted in the ten HMXBs are also provided, as well as the long-term lightcurves folded on the sources best determined orbital parameters. These lightcurves reveal complex eclipse ingresses and egresses, that are understood mostly as being due to the presence of accretion wakes. The results reported in this paper constitute a database to be used for population and evolutionary studies of HMXBs, as well as theoretical modelling of long-term accretion in wind-fed X-ray binaries.

TLUSTY non-LTE model atmosphere calculations have been used to determine atmospheric parameters and nitrogen (N) abundances for 34 single and 18 binary B-type supergiants (BSGs). The effects of flux contribution from an unseen secondary were considered for the binary sample. We present the first systematic study of the incidence of binarity for a sample of BSGs across the theoretical terminal age main sequence (TAMS). To account for the distribution of effective temperatures of the BSGs it may be necessary to extend the TAMS to lower temperatures. This is consistent with the derived distribution of mass discrepancies, projected rotational velocities (vsini) and N abundances, provided that stars cooler than this temperature are post RSG objects. For the BSGs in the Tarantula and previous FLAMES surveys, most have small vsini. About 10% have larger vsini (>100 km/s) but surprisingly these show little or no N enhancement. All the cooler BSGs have low vsini of <70km/s and high N abundance estimates, implying that either bi-stability braking or evolution on a blue loop may be important. A lack of cool binaries, possibly reflects the small sample size. Single star evolutionary models, which include rotation, can account for the N enhancement in both the single and binary samples. The detailed distribution of N abundances in the single and binary samples may be different, possibly reflecting differences in their evolutionary history. The first comparative study of single and binary BSGs has revealed that the main sequence may be significantly wider than previously assumed, extending to Teff=20000K. Some marginal differences in single and binary atmospheric parameters and abundances have been identified, possibly implying non-standard evolution for some of the sample. This sample as a whole has implications for several aspects of our understanding of the evolution of BSGs. Full abstract in paper

Using own and literature data for a large sample of O stars in the Milky Way, we investigate the correspondence between their spectroscopic and evolutionary masses, and try to put constraints on various parameters that might influence the estimates of these two quantities.

We are now in an era of large spectroscopic surveys of OB-type stars. Quantitative spectroscopic analysis of these modern datasets is enabling us to review the physical properties of blue massive stars with robust samples, not only revisiting mean properties and general trends, but also incorporating information about the effects of second-order parameters. We investigate the spectral type -- effective temperature (SpT - Teff) calibration for O-type dwarfs, and its claimed dependence on metallicity, using statistically-meaningful samples of stars extracted from the IACOB and VFTS surveys. We perform a homogeneous differential spectroscopic analysis of 33 Galactic and 53 LMC O~dwarfs (spanning spectral types of O4 -- O9.7) using the IACOB-GBAT package, a $\chi^2$-fitting algorithm based on a large pre-computed grid of FASTWIND models, and standard techniques for the hydrogen/helium analysis of O-type stars. We compare the estimated effective temperatures and gravities as a function of (internally consistent) spectral classifications. While the general trend is that the temperature of a star increases with earlier spectral types and decreasing metallicity, we show that the large range of gravities found for O-type dwarfs -- spaning up to 0.45-0.50 dex in some spectral bins -- plays a critical role on the dependence of the effective temperature calibrations as a function of spectral type and metallicity. This result warns us about the use of SpT - Teff calibrations for O-dwarfs which ignore the effects of gravity, and highlights the risks of employing calibrations based on small samples. The effects of this scatter in gravities (evolutionary status) for O-type dwarfs should be included in future recipes which employ SpT - Teff calibrations.

We review potential mass-loss mechanisms in the various evolutionary stages of massive stars, from the well-known line-driven winds of O-stars and BA-supergiants to the less-understood winds from Red Supergiants. We discuss optically thick winds from Wolf-Rayet stars and Very Massive Stars, and the hypothesis of porosity-moderated, continuum-driven mass loss from stars formally exceeding the Eddington limit, which might explain the giant outbursts from Luminous Blue Variables. We finish this review with a glance on the impact of rapid rotation, magnetic fields and small-scale inhomogeneities in line-driven winds.

Rotation is a key parameter in the evolution of massive stars, affecting their evolution, chemical yields, ionizing photon budget, and final fate. We determined the projected rotational velocity, $v_e\sin i$, of $\sim$330 O-type objects, i.e. $\sim$210 spectroscopic single stars and $\sim$110 primaries in binary systems, in the Tarantula nebula or 30 Doradus (30\,Dor) region. The observations were taken using VLT/FLAMES and constitute the largest homogeneous dataset of multi-epoch spectroscopy of O-type stars currently available. The most distinctive feature of the $v_e\sin i$ distributions of the presumed-single stars and primaries in 30 Dor is a low-velocity peak at around 100\,$\rm{km s^{-1}}$. Stellar winds are not expected to have spun-down the bulk of the stars significantly since their arrival on the main sequence and therefore the peak in the single star sample is likely to represent the outcome of the formation process. Whereas the spin distribution of presumed-single stars shows a well developed tail of stars rotating more rapidly than 300\,$\rm{km s^{-1}}$, the sample of primaries does not feature such a high-velocity tail. The tail of the presumed-single star distribution is attributed for the most part -- and could potentially be completely due -- to spun-up binary products that appear as single stars or that have merged. This would be consistent with the lack of such post-interaction products in the binary sample, that is expected to be dominated by pre-interaction systems. The peak in this distribution is broader and is shifted toward somewhat higher spin rates compared to the distribution of presumed-single stars. Systems displaying large radial velocity variations, typical for short period systems, appear mostly responsible for these differences.

[Abridged] We develop and benchmark a fast and easy-to-use effective-opacity formalism for line and continuum radiative transfer in an accelerating two-component clumpy medium. The formalism bridges the limits of optically thin and thick clumps, and is here used to i) design a simple vorosity-modified Sobolev with exact integration (vmSEI) method for analyzing UV wind resonance lines in hot, massive stars, and ii) derive simple correction factors to the line force driving the outflows of such stars. We show that (for a given ionization factor) UV resonance doublets may be used to analytically predict the upward corrections in empirically inferred mass-loss rates associated with porosity in velocity space (a.k.a. velocity-porosity, or vorosity), but that severe solution degeneracies exist. For an inter-clump density set to 1 % of the mean density, we for O and B supergiants derive upward empirical mass-loss corrections of typically factors of either ~5 or ~50, depending on which of the two applicable solutions is chosen. Overall, our results indicate this solution dichotomy severely limits the use of UV resonance lines as direct mass-loss indicators of clumped hot stellar winds. We next apply the effective-opacity formalism to the standard CAK theory of line-driven winds. By analytic and numerical hydrodynamics calculations, we show that in cases where vorosity is important at the critical point setting the mass-loss rate, the reduced line-force leads to a lower theoretical mass loss, by a factor scaling with the normalized velocity filling factor fvel. On the other hand, if vorosity is important only above this critical point, the predicted mass loss is not affected, but the wind terminal speed is reduced. This shows that porosity in velocity space can have a significant impact not only on the diagnostics, but also on the dynamics and theory of radiatively driven winds.

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.

Rotation is of key importance for the evolution of hot massive stars, however, the rotational velocities of these stars are difficult to determine. Based on our own data for 31 Galactic O stars and incorporating similar data for 86 OB supergiants from the literature, we aim at investigating the properties of rotational and extra line-broadening as a function of stellar parameters and at testing model predictions about the evolution of stellar rotation. Fundamental stellar parameters were determined by means of the code FASTWIND. Projected rotational and extra broadening velocities originate from a combined Ft + GOF method. Model calculations published previously were used to estimate the initial evolutionary masses. The sample O stars with Minit > 50 Msun rotate with less that 26% of their break-up velocity, and they also lack objects with v sin i < 50 km/s. For the stars with Minit > 35 Msun on the hotter side of the bi-stability jump, the observed and predicted rotational rates agree quite well; for those on the cooler side of the jump, the measured velocities are systematically higher than the predicted ones. In general, the derived extra broadening velocities decrease toward cooler Teff, whilst for later evolutionary phases they appear, at the same v sin i, higher for high-mass stars than for low-mass ones. None of the sample stars shows extra broadening velocities higher than 110 km/s. For the majority of the more massive stars, extra broadening either dominates or is in strong competition with rotation. Conclusions: For OB stars of solar metallicity, extra broadening is important and has to be accounted for in the analysis. When appearing at or close to the zero-age main sequence, most of the single and more massive stars rotate slower than previously thought. Model predictions for the evolution of rotation in hot massive stars may need to be updated.

Context: Rotation plays a key role in the life cycles of stars with masses above 8 Msun. Hence, accurate knowledge of the rotation rates of such massive stars is critical for understanding their properties and for constraining models of their evolution. Aims: This paper investigates the reliability of current methods used to derive projected rotation speeds v sin i from line-broadening signatures in the photospheric spectra of massive stars, focusing on stars that are not rapidly rotating. Methods: We use slowly rotating magnetic O-stars with well-determined rotation periods to test the Fourier transform (FT) and goodness-of-fit (GOF) methods typically used to infer projected rotation rates of massive stars. Results: For our two magnetic test stars with measured rotation periods longer than one year, i.e., with v sin i < 1 km/s, we derive v sin i ~ 40-50 km/s from both the FT and GOF methods. These severe overestimates are most likely caused by an insufficient treatment of the competing broadening mechanisms referred to as microturbulence and macroturbulence. Conclusions: These findings warn us not to rely uncritically on results from current standard techniques to derive projected rotation speeds of massive stars in the presence of significant additional line broadening, at least when v sin i <~ 50 km/s. This may, for example, be crucial for i) determining the statistical distribution of observed rotation rates of massive stars, ii) interpreting the evolutionary status and spin-down histories of rotationally braked B-supergiants, and iii) explaining the deficiency of observed O-stars with spectroscopically inferred v sin i ~ 0 km/s. Further investigations of potential shortcomings of the above techniques are presently under way.

The 30 Doradus (30\,Dor) region in the Large Magellanic Cloud (also known as the Tarantula Nebula) is the nearest massive starburst region, containing the richest sample of massive stars in the Local Group. It is the best possible laboratory to investigate aspects of the formation and evolution of massive stars. Here, we focus on rotation which is a key parameter in the evolution of these objects. We establish the projected rotational velocity, $v_{e}\sin i$, distribution of an unprecedented sample of 216 radial velocity constant ($\rm{\Delta RV\, \leq\, 20 \,km s^{-1}}$) O-type stars in 30\,Dor observed in the framework of the VLT-FLAMES Tarantula Survey (VFTS). The distribution of $v_{e}\sin i$ shows a two-component structure: a peak around 80 $\rm{km s^{-1}}$ and a high-velocity tail extending up to $\sim$600 $\rm{km s^{-1}}$. Around 75% of the sample has 0 $\leq\, v_{e}\sin i \leq$ 200 $\rm{km s^{-1}}$ with the other 25% distributed in the high-velocity tail. The presence of the low-velocity peak is consistent with that found in other studies of late-O and early-B stars. The high-velocity tail is compatible with expectations from binary interaction synthesis models and may be predominantly populated by post-binary interaction, spun-up, objects and mergers. This may have important implications for the nature of progenitors of long-duration gamma ray bursts.

Aims. Using ground based multi-object optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS), we aim to establish the (projected) rotational velocity distribution for a sample of 216 presumably single O-type stars in 30 Doradus (30 Dor). Methods. We measured projected rotational velocities, \vrot, by means of a Fourier transform method and a profile fitting method applied on a set of isolated spectral lines. We also used an iterative deconvolution procedure to infer the probability density, $\rm{P(\veq)}$, of the equatorial rotational velocity, \veq. Results. The distribution of \vrot\ shows a two-component structure: a peak around 80 \kms\ and a high-velocity tail extending up to $\sim$600 \kms. This structure is also present in the inferred distribution $\rm{P(\veq)}$ with around 80% of the sample having 0 $<$ \veq\, $\leq\, 300$ \kms\ and the other 20% distributed in the high-velocity region. Conclusions. Most of the stars in our sample rotate with a rate less than 20%\ of their break-up velocity. For the bulk of the sample, mass-loss in a stellar wind and/or envelope expansion is not efficient enough to significantly spin down these stars within the first few Myr of evolution. The presence of a sizeable population of fast rotators is compatible with recent population synthesis computations that investigate the influence of binary evolution on the rotation rate of massive stars. Despite the fact that we have excluded stars that show significant radial velocity variations, our sample may have remained contaminated by post-interaction binary products. The fact that the high-velocity tail may be preferentially (and perhaps even exclusively), populated by post-binary interaction products, has important implications for the evolutionary origin of systems that produce gamma-ray bursts.

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)

Massive, hot OB-stars show clear evidence of strong macroscopic broadening (in addition to rotation) in their photospheric spectral lines. This paper examines the occurrence of such "macro-turbulence" in slowly rotating O-stars with strong, organised surface magnetic fields. Focusing on the CIV 5811A line, we find evidence for significant macro-turbulent broadening in all stars except NGC1624-2, which also has (by far) the strongest magnetic field. Instead, the very sharp CIV lines in NGC1624-2 are dominated by magnetic Zeeman broadening, from which we estimate a dipolar field of approximately 20 kG. By contrast, magnetic broadening is negligible in the other stars (due to their weaker field strengths, on order 1 kG), and their CIV profiles are typically very broad and similar to corresponding lines observed in non-magnetic O-stars. Quantifying this by an isotropic, Gaussian macro-turbulence, we derive vmac = 2.2 (+- 0.9/2.2) km/s for NGC1624-2, and vmac = 20-65 km/s for the rest of the magnetic sample. We use these observational results to test the hypothesis that the field can stabilise the atmosphere and suppress the generation of macro-turbulence down to stellar layers where the magnetic pressure PB and the gas pressure Pg are comparable. Using a simple grey atmosphere to estimate the temperature T0 at which PB = Pg, we find that T0 > Teff for all investigated magnetic stars, but that T0 reaches the ~160000 K layers associated with the iron opacity-bump in hot stars only for NGC1624-2. This is consistent with the view that the responsible physical mechanism for photospheric O-star macro-turbulence may be stellar gravity-mode oscillations excited by sub-surface convection zones, and suggests that a sufficiently strong magnetic field can suppress such iron-bump generated convection and associated pulsational excitation.

The model atmosphere programs FASTWIND and CMFGEN are both elegantly designed to perform non-LTE analyses of the spectra of hot massive stars, and include sphericity and mass-loss. The two codes differ primarily in their approach towards line blanketing, with CMFGEN treating all of the lines in the co-moving frame and FASTWIND taking an approximate approach which speeds up execution times considerably. Although both have been extensively used to model the spectra of O-type stars, no studies have used the codes to independently model the same spectra of the same stars and compare the derived physical properties. We perform this task on ten O-type stars in the Magellanic Clouds. For the late-type O supergiants, both CMFGEN and FASTWIND have trouble fitting some of the He I lines, and we discuss causes and cures. We find that there is no difference in the average effective temperatures found by the two codes for the stars in our sample, although the dispersion is large, due primarily to the various difficulties each code has with He I. The surface gravities determined using FASTWIND are systematically lower by 0.12 dex compared to CMFGEN, a result we attribute to the better treatment of electron scattering by CMFGEN. This has implications for the interpretation of the origin of the so-called mass discrepancy, as the masses derived by FASTWIND are on average lower than inferred from stellar evolutionary models, while those found by CMFGEN are in better agreement.

Recent studies have claimed the existence of very massive stars (VMS) up to 300 solar masses in the local Universe. As this finding may represent a paradigm shift for the canonical stellar upper-mass limit of 150 Msun, it is timely to discuss the status of the data, as well as the far-reaching implications of such objects. We held a Joint Discussion at the General Assembly in Beijing to discuss (i) the determination of the current masses of the most massive stars, (ii) the formation of VMS, (iii) their mass loss, and (iv) their evolution and final fate. The prime aim was to reach broad consensus between observers and theorists on how to identify and quantify the dominant physical processes.

Context. Results from the theory of radiatively driven winds are incorporated in stellar evolutionary and population synthesis models, and used in our interpretation of the observations of the deep Universe. Yet, the theory has been confirmed only until Small Magellanic Cloud metallicities. Analyses of O-stars at lower metallicities are needed to prove the theory. Aims. We have observed GHV-62024, an O6.5 IIIf star in the low-metallicity galaxy IC1613. According to a previous preliminary analysis this star could challenge the radiatively driven wind theory at low metallicities. Methods. Our observations were obtained with VIMOS at VLT, at R~2000 and were analysed using the latest version of the model atmosphere code FASTWIND, which includes N III Results. We obtain the stellar parameters and conclude that the star follows the average wind momentum-luminosity relationship (WLR) expected for its metallicity, but with a high value for the exponent of the wind velocity law, beta. We suggest that this high value may be reached because GHV-62024 could be a fast rotator seen at a low inclination angle. While the derived beta value does not change by adopting a lower wind terminal velocity, a wrong $V_\infty$ has a clear impact on the position of the star in the WLR diagram. The N and He abundances are very high, consistent with strong CNO mixing that could have been caused by the fast rotation, although we cannot discard a different origin. We find again the well-known mass-discrepancy. Conclusions. We conclude that the star follows the WLR expected for its metallicity. The results are consistent with GHV-62024 being a fast rotator seen close to pole-on, contaminated at the surface with CNO products and with a wind structure altered by the fast rotation without modifying the global WLR. We suggest that this could be a general property of fast rotators.

This is the third 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 provide first theoretical predictions for the NIV4058/NIII4640 emission line ratio in dependence of various parameters, and confront them with results from the analysis of a sample of early LMC/SMC O-stars. Stellar and wind parameters are determined by line profile fitting of H/He/N lines, exploiting the helium and nitrogen ionization balance. Corresponding synthetic spectra are calculated using the NLTE atmospheric code FASTWIND. Though there is a monotonic relationship between the emission line ratio and Teff, all other parameters being equal, theoretical predictions indicate additional dependencies, most notably, on the nitrogen abundance. These basic predictions are confirmed by results from atmospheric code CMFGEN. The effective temperatures for the earliest O-stars, inferred from the nitrogen ionization balance, are partly considerably hotter than indicated by previous studies. Consistent with earlier results, effective temperatures increase from supergiants to dwarfs for all spectral types in the LMC. The relation between observed NIV4058/NIII4640 emission line ratio and Teff, for a given luminosity class, turned out to be quite monotonic for our sample stars, and fairly consistent with our model predictions. The scatter within a spectral sub-type is mainly produced by abundance effects. Our findings suggest that the Walborn et al. (2002) classification scheme is able to provide a meaningful relation between spectral type and Teff, provided that it is possible to discriminate for the luminosity class. This might be difficult to achieve in low-Z environments such as the SMC, owing to rather low wind-strengths. According to our predictions, the major bias of the classification scheme is due to nitrogen content.

We present a brief progress report of our project aimed at the investigation of the so-called macroturbulent broadening affecting the line-profiles of O and early B-type stars, and speculate about the possible use of this spectroscopic feature as a single snap-shot approach to investigate pulsations in massive stars.

We present the IACOB grid-based automatic tool for the quantitative spectroscopic analysis of O-stars. The tool consists of an extensive grid of FASTWIND models, and a variety of programs implemented in IDL to handle the observations, perform the automatic analysis, and visualize the results. The tool provides a fast and objective way to determine the stellar parameters and the associated uncertainties of large samples of O-type stars within a reasonable computational time.

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.

This review describes the evidence for small-scale structure, `clumping', in the radiation line-driven winds of hot, massive stars. In particular, we focus on examining to what extent simulations of the strong instability inherent to line-driving can explain the multitude of observational evidence for wind clumping, as well as on how to properly account for extensive structures in density and velocity when interpreting the various wind diagnostics used to derive mass-loss rates.

This is the first paper in a series dealing with optical Nitrogen spectroscopy of O-type stars, aiming at the analysis of Nitrogen abundances. We implemented a new Nitrogen model atom into the NLTE atmosphere/spectrum synthesis code FASTWIND, and compare the resulting optical NIII lines at 4634/40/42 A with other predictions, mostly from Mihalas & Hummer (1973, ApJ 179, 827,`MH'), and from the alternative code CMFGEN. Using similar model atmospheres as MH (not blanketed and wind-free), we are able to reproduce their results, in particular the triplet emission lines. According to MH, these should be strongly related to dielectronic recombination (DR) and the drain by certain two-electron transitions. However, using realistic, fully line-blanketed atmospheres at solar abundances, the key role of DR controlling these emission features is superseded -- for O-star conditions -- by the strength of the stellar wind and metallicity. In the case of wind-free models, the resulting lower ionizing EUV-fluxes severely suppress the emission. As the mass-loss rate is increased, pumping through the NIII resonance line(s) in the presence of a near-photospheric velocity field results in a net optical triplet line emission. A comparison with results from CMFGEN is mostly satisfactory, except for the range 30 kK < Teff < 35 kK, where CMFGEN triggers the triplet emission at lower Teff than FASTWIND. This effect could be traced down to line overlap effects between the NIII and OIII resonance lines that so far cannot be simulated by FASTWIND. Since the efficiency of DR and `two electron drain' strongly depends on the degree of line-blanketing/-blocking, we predict the emission to become stronger in a metal-poor environment, though lower wind-strengths and Nitrogen abundances might counteract this effect. Weak winded stars should display less triplet emission than stars with `normal' winds.

Context. Mass-loss, occurring through radiation driven supersonic winds, is a key issue throughout the evolution of massive stars. Two outstanding problems are currently challenging the theory of radiation-driven winds: wind clumping and the weak-wind problem. Aims. We seek to obtain accurate mass-loss rates of OB stars at different evolutionary stages to constrain the impact of both problems in our current understanding of massive star winds. Methods. We perform a multi-wavelength quantitative analysis of a sample of ten Galactic OB-stars by means of the atmospheric code CMFGEN, with special emphasis on the L-band window. A detailed investigation is carried out on the potential of Br\alpha\ and Pf\gamma\ as mass-loss and clumping diagnostics. Results. For objects with dense winds, Br\alpha\ samples the intermediate wind while Pf\gamma\ maps the inner one. In combination with other indicators (UV, H\alpha, Br\gamma) these lines enable us to constrain the wind clumping structure and to obtain "true" mass-loss rates. For objects with weak winds, Br\alpha\ emerges as a reliable diagnostic tool to constrain the mass-loss rates. The emission component at the line Doppler-core superimposed on the rather shallow Stark absorption wings reacts very sensitively to mass loss already at very low mass-loss values. On the other hand, the line wings display similar sensitivity to mass loss as H\alpha, the classical optical mass loss diagnostics. Conclusions. Our investigation reveals the great diagnostic potential of L-band spectroscopy to derive clumping properties and mass-loss rates of hot star winds. We are confident that Br\alpha\ will become the primary diagnostic tool to measure very low mass-loss rates with unprecedented accuracy

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

The modern-era spectral classification of O-stars relies on either the Walborn or the Conti-Mathys scheme. Since both of these approaches have been developed using low-quality photographic data, their application to high-quality digital data might not be straightforward and be hampered by problems and complications that have not yet been appreciated. Using high-resolution spectra obtained with the ESO/MPG 2.2\,m telescope in La Silla and following the premises of the Walborn and Conti classification schemes, we determined the spectral types and luminosity classes of 19 Galactic O-type stars and compared them to those attributed by Walborn and Mathys based on low-quality data. Our analysis reveals that the morphological spectral types assigned using high-resolution data are systematically later (by up to 1.5 subtypes) then those attributed by Walborn. By means of line-profile simulations, we show that part of this discrepancy is more likely caused by the combined effect of stellar rotation and high spectral resolution on the depth of helium lines used as spectral type indicators. In addition, we demonstrate that at least for narrow-lined stars the "rotational effect" does not disappear when the high-resolution spectra are degraded to the resolution of the Walborn standards. We also find evidence of a systematic difference between our high-resolution quantitative spectral types and those assigned by Mathys. Rotation and spectral resolution are important third parameters in the spectral classification of O-type stars. To obtain reliable spectral classes within the Walborn approach, the unknown and the standard spectra must be compared at the same resolution and \vsini. Owing to resolution effects, the Conti approach might also need to be updated.

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.

Mass-loss rates currently in use for hot, massive stars have recently been seriously questioned, mainly because of the effects of wind clumping. We investigate the impact of clumping on diagnostic ultraviolet resonance and optical recombination lines. Optically thick clumps, a non-void interclump medium, and a non-monotonic velocity field are all accounted for in a single model. We used 2D and 3D stochastic and radiation-hydrodynamic (RH) wind models, constructed by assembling 1D snapshots in radially independent slices. To compute synthetic spectra, we developed and used detailed radiative transfer codes for both recombination lines (solving the "formal integral") and resonance lines (using a Monte-Carlo approach). In addition, we propose an analytic method to model these lines in clumpy winds, which does not rely on optically thin clumping. Results: Synthetic spectra calculated directly from current RH wind models of the line-driven instability are unable to in parallel reproduce strategic optical and ultraviolet lines for the Galactic O-supergiant LCep. Using our stochastic wind models, we obtain consistent fits essentially by increasing the clumping in the inner wind. A mass-loss rate is derived that is approximately two times lower than predicted by the line-driven wind theory, but much higher than the corresponding rate derived from spectra when assuming optically thin clumps. Our analytic formulation for line formation is used to demonstrate the potential impact of optically thick clumping in weak-winded stars and to confirm recent results that resonance doublets may be used as tracers of wind structure and optically thick clumping. (Abridged)

Small-scale inhomogeneities, or `clumping', in the winds of hot, massive stars are conventionally included in spectral analyses by assuming optically thin clumps. To reconcile investigations of different diagnostics using this microclumping technique, very low mass-loss rates must be invoked for O stars. Recently it has been suggested that by using the microclumping approximation one may actually drastically underestimate the mass-loss rates. Here we demonstrate this, present a new, improved description of clumpy winds, and show how corresponding models, in a combined UV and optical analysis, can alleviate discrepancies between previously derived rates and those predicted by the line-driven wind theory. Furthermore, we show that the structures obtained in time-dependent, radiation-hydrodynamic simulations of the intrinsic line-driven instability of such winds, which are the basis to our current understanding of clumping, in their present-day form seem unable to provide a fully self-consistent, simultaneous fit to both UV and optical lines. The reasons for this are discussed.

In this poster we present the results of our analyses of three early massive stars in IC 1613, whose spectra have been observed with VIMOS and analyzed with CMFGEN and FASTWIND. One of the targets resulted a possible LBV and the other two are Of stars with unexpectedly strong winds. The Of stars seem to be strongly contaminated by CNO products. Our preliminary results may represent a challenge for the theory of stellar atmospheres, but they still have to be confirmed by the analysis of more objects and a more complete coverage of the parameter space.

As part of a long term observational project, we are investigating the macroturbulent broadening in O and B supergiants (Sgs) and its possible connection with spectroscopic variability phenomena and stellar oscillations. We present the first results of our project, namely firm observational evidence for a strong correlation between the extra broadening and photospheric line-profile variations in a sample of 13 Sgs with spectral types ranging from O9.5 to B8.

Massive stars are inherently extreme objects, in terms of radiation, mass loss, rotation, and sometimes also magnetic fields. Concentrating on a (personally biased) subset of processes related to pulsations, rapid rotation and its interplay with mass-loss, and the bi-stability mechanism, we will discuss how active (and normal) OB stars can serve as appropriate laboratories to provide further clues.

The spectra of O and B supergiants are known to be affected by a significant form of extra line broadening (usually referred to as macroturbulence) in addition to that produced by stellar rotation. Recent analyses of high resolution spectra have shown that the interpretation of this line broadening as a consequence of large scale turbulent motions would imply highly supersonic velocity fields in photospheric regions, making this scenario quite improbable. Stellar oscillations have been proposed as a likely alternative explanation. As part of a long term observational project, we are investigating the macroturbulent broadening in O and B supergiants and its possible connection with spectroscopic variability phenomena and stellar oscillations. In this letter, we present the first encouraging results of our project, namely firm observational evidence for a strong correlation between the extra broadening and photospheric line-profile variations in a sample of 13 supergiants with spectral types ranging from O9.5 to B8.

The spectrum of O and B Supergiants is known to be affected by an important extra line-broadening (usually called macroturbulence) that adds to stellar rotation. Recent analysis of high resolution spectra has shown that the interpretation of this line-broadening as a consequence of large-scale turbulent motions would imply highly super-sonic velocity fields, making this scenario quite improbable. Stellar oscillations have been proposed as a likely alternative explanation. We present first encouraging results of an observational project aimed at investigating the $macroturbulent$ broadening in O and B Supergiants, and its possible connection with spectroscopic variability phenomena and stellar oscillations: a) all the studied B Supergiants show line profile variations, quantified by means of the first (< v >) and third velocity (< v^3 >) moments of the lines, b) there is a strong correlation between the peak-to-peak amplitudes of the < v > and < v^3 > variability and the size of the extra-broadening.

Small-scale clumping in the winds of hot, massive stars is conventionally included in spectral analyses by assuming optically thin clumps, a void inter-clump medium, and a smooth velocity field. To reconcile investigations of different diagnostics within such models, a highly clumped wind with very low mass-loss rates needs to be invoked. Particularly, unsaturated UV resonance lines seem to indicate rates an order of magnitude (or even more) lower than previously accepted values. We investigate resonance line formation in inhomogeneous hot star winds with non-monotonic velocity fields by means of 2D stochastic and pseudo-2D radiation-hydrodynamic wind models. A Monte-Carlo radiative transfer code is presented and used to produce synthetic line spectra. Results: The optically thin clumping limit is only valid for very weak lines. For intermediate strong lines, the velocity spans of the clumps are of central importance. Current hydrodynamical models predict spans that are too large to reproduce observed profiles unless a very low mass-loss rate is invoked. By simulating lower spans in 2D stochastic models, the profile strengths become drastically reduced, and are consistent with higher mass-loss rates. To simultaneously meet the constraints from strong lines, the inter-clump medium must be non-void. A first comparison to the observed PV doublet in the O6 supergiant lam Cep confirms that a stochastic 2D model reproduces observations with a mass-loss rate roughly ten times higher than that derived from the same lines but assuming optically thin clumping. Tentatively this may resolve discrepancies between theoretical predictions, evolutionary constraints, and recent derived mass-loss rates, and suggests a re-investigation of the structure predicted by current hydrodynamical models.

We aim to determine the fundamental parameters of a sample of B stars with apparent visual magnitudes below 8 in the field-of-view of the CoRoT space mission, from high-resolution spectroscopy. We developed an automatic procedure for the spectroscopic analysis of B-type stars with winds, based on an extensive grid of FASTWIND model atmospheres. We use the equivalent widths and/or the line profile shapes of continuum normalized hydrogen, helium and silicon line profiles to determine the fundamental properties of these stars in an automated way. After thorough tests, both on synthetic datasets and on very high-quality, high-resolution spectra of B stars for which we already had accurate values of their physical properties from alternative analyses, we applied our method to 66 B-type stars contained in the ground-based archive of the CoRoT space mission. We discuss the statistical properties of the sample and compare them with those predicted by evolutionary models of B stars. Our spectroscopic results provide a valuable starting point for any future seismic modelling of the stars, should they be observed by CoRoT.

We aimed at finding a physical explanation for the occurrence of macroturbulence in the atmospheres of hot massive stars, a phenomenon found in observations since more than a decade but yet unexplained. We computed time series of line profiles for evolved massive stars broadened by rotation and by hundreds of low-amplitude nonradial gravity-mode pulsations which are predicted to be excited for evolved massive stars. In general, line profiles based on macrotubulent broadening can mimic those subject to pulsational broadening. In several cases, though, good fits require macroturbulent velocities that pass the speed of sound for realistic pulsation amplitudes. Moreover, we find that the rotation velocity can be seriously underestimated by using a simple parameter description for macroturbulence rather than an appropriate pulsational model description to fit the line profiles. We conclude that macroturbulence is a likely signature of the collective effect of pulsations. We provide line diagnostics and their typical values to decide whether or not pulsational broadening is present in observed line profiles, as well as a procedure to avoid an inaccurate estimation of the rotation velocity.

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 review recent developments regarding radiation driven mass loss from OB-stars. We first summarize the fundamental theoretical predictions, and then compare these to observational results (including the VLT-FLAMES survey of massive stars). Especially we focus on the mass loss-metallicity dependence and on the so-called bi-stability jump. Subsequently we concentrate on two urgent problems, weak winds and wind clumping, that have been identified from various diagnostics and that challenge our present understanding of radiation driven winds. We discuss the problems of ``measuring'' mass-loss rates from weak winds and the potential of the near infrared, Br-alpha line as a tool to enable a more precise quantification, and comment on physical explanations for mass-loss rates that are much lower than predicted by the standard model. Wind clumping, conventionally interpreted as the consequence of a strong instability inherent to radiative line-driving, has severe implications for the interpretation of observational diagnostics, since derived mass-loss rates are usually overestimated when clumping is present but ignored in the analyses. Simplified techniques to account for clumping indicate overestimates by factors of 2 to 10, or even more. If actually true, these results would have a dramatic impact on the evolution of, and the feedback from, massive stars. We discuss ongoing attempts (including own work) to interprete the corresponding observations in terms of more sophisticated models. By allowing for porosity in density and velocity space, and for a non-void inter-clump medium, such models might require only moderate reductions of mass-loss rates.

Based on an exteded grid of NLTE, line blanketed model atmospheres with stellar winds as calculated by means of FASTWIND, we have investigated the change in the strengths of strategic Helium transitions in the optical as caused by a 0.3 decrease in metallicity with respect to solar abundances. Our calculations predict that only part of the observed increase in Teff, of O-type dwarfs could be explained by metallicity effects on the spectral type indicators, while the rest must be attributed to other reasons (e.g., different stellar structures as a function of metallicity or differences between observed and theoretical wind parameters etc.). In addition, we found that using the He II 4686 line to classify stars in low metallicity environments (Z < 0.3 solar) might artificially increase the number of low luminosity (dwarfs and giants) O-stars, on the expense of the number of O-supergiants.

Since the use of high-resolution high signal-to-noise spectroscopy in the study of massive stars, it became clear that an ad-hoc velocity field at the stellar surface, termed macroturbulence, is needed to bring the observed shape of spectral lines into agreement with observations. We seek a physical explanation of this unknown broadening mechanism. We interprete the missing line broadening in terms of collective pulsational velocity broadening due to non-radial gravity-mode oscillations. We also point out that the rotational velocity can be seriously underestimated whenever the line profiles are fitted assuming a Gaussian macroturbulent velocity rather than an appropriate pulsational velocity expression.

Mass loss is a key process in the evolution of massive stars, and must be understood quantitatively to be successfully included in broader astrophysical applications. In this review, we discuss various aspects of radiation driven mass loss, both from the theoretical and the observational side. We focus on winds from OB-stars, with some excursions to the Luminous Blue Variables, Wolf- Rayet stars, A-supergiants and Central Stars of Planetary Nebulae. After reca- pitulating the 1-D, stationary standard model of line-driven wind, extensions accounting for rotation and magnetic fields are discussed. The relevance of the so-called bi-stability jump is outlined. We summarize diagnostical methods to infer wind properties from observations, and compare the results with theore- tical predictions, featuring the massloss-metallicity dependence. Subsequently, we concentrate on two urgent problems which challenge our present understanding of radiation driven winds: weak winds and wind- clumping. We discuss problems of measuring mass-loss rates from weak winds and the potential of NIR- spectroscopy. Wind-clumping has severe implications for the interpretation of observational diagnostics, as derived mass-loss rates can be overestimated by factors of 2 to 10 if clumping is ignored, and we describe ongoing attempts to allow for more uniform results. We point out that independent arguments from stellar evolution favor a moderate reduction of present- day mass-loss rates. We also consider larger scale wind structure, interpreted in terms of co-rotating interacting regions, and complete this review with a discussion of recent progress on the X-ray line emission from massive stars, highlighting as to how far the analysis of such X-ray line emission can give further clues regarding an adequate description of wind clumping. (Abridged abstract)

In order to better determine the physical properties of hot, massive stars as a function of metallicity, we obtained very high SNR optical spectra of 26 O and early B stars in the Magellanic Clouds. These allow accurate modeling even in cases where the He I 4471 line has an equivalent width of only a few tens of mA. The spectra were modeled with FASTWIND, with good fits obtained for 18 stars; the remainder show signatures of being binaries. We include stars in common to recent studies to investigate possible systematic differences. The "automatic" FASTWIND modeling method of Mokiem and collaborators produced temperatures 1100 K hotter on the average, presumably due to the different emphasis given to various temperature-sensitive lines. More significant, however, is that the automatic method always produced some "best" answer, even for stars we identify as composite (binaries). The temperatures found by the TLUSTY/CMFGEN modeling of Bouret, Heap, and collaborators yielded temperatures 1000 K cooler than ours, on average. Significant outliers were due either to real differences in the data (some of the Bouret/Heap data were contaminated by moonlight continua) or the fact we could detect the HeI line needed to better constrain the temperature. Our new data agrees well with the effective temperature scale we presented previously. We confirm that the "Of" emission-lines do not track luminosity classes in the exact same manner as in Milky Way stars. We revisit the the issue of the "mass discrepancy", finding that some of the stars in our sample do have spectroscopic masses that are significantly smaller than those derived from stellar evolutionary models. We do not find that the size of the mass discrepancy is simply related to either effective temperature or surface gravity.

We report on multi-wavelength observations, ranging from the X-ray to radio wave bands, of the IPN-localized gamma-ray burst GRB 070125. Spectroscopic observations reveal the presence of absorption lines due to O I, Si II, and C IV, implying a likely redshift of z = 1.547. The well-sampled light curves, in particular from 0.5 to 4 days after the burst, suggest a jet break at 3.7 days, corresponding to a jet opening angle of ~7.0 degrees, and implying an intrinsic GRB energy in the 1 - 10,000 keV band of around E = (6.3 - 6.9)x 10^(51) erg (based on the fluences measured by the gamma-ray detectors of the IPN network). GRB 070125 is among the brightest afterglows observed to date. The spectral energy distribution implies a host extinction of Av < 0.9 mag. Two rebrightening episodes are observed, one with excellent time coverage, showing an increase in flux of 56% in ~8000 seconds. The evolution of the afterglow light curve is achromatic at all times. Late-time observations of the afterglow do not show evidence for emission from an underlying host galaxy or supernova. Any host galaxy would be subluminous, consistent with current GRB host-galaxy samples. Evidence for strong Mg II absorption features is not found, which is perhaps surprising in view of the relatively high redshift of this burst and the high likelihood for such features along GRB-selected lines of sight.

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 apply the NLTE atmosphere code FASTWIND to perform a spectroscopic study of a small sample of Galactic B-supergiants from B0 to B9. By means of the resulting data and incorporating additional datasets from alternative studies, we investigate the properties of OB-supergiants and compare our findings with theoretical predictions. As a result we find that due to the combined effects of line- and wind-blanketing, the temperature scale of Galactic B-supergiants needs to be revised downwards, by 10 to 20 percent, the latter value being appropriate for stronger winds. In fair accordance with recent results, our sample furthermore indicates a gradual decrease in wind terminal velocities over the bi-stability region, where the limits of this region are located at lower temperatures than the predicted ones. Introducing a distance-independent quantity Q' related to wind-strength, we also show that this quantity is a well defined, monotonically increasing function of Teff outside this region. Inside and from hot to cool, the mass loss rate changes by a factor (in between 0.4 and 2.5) which is (much) smaller than the predicted factor of 5. All this indicates that the decrease in wind terminal velocity over the bi-stability region is not over-compensated by an increase of mass loss rate, as frequently argued (provided the wind-clumping properties on both sides of this region do not differ substantially).

[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.

We study the influence of clumping on the predicted wind structure of O-type stars. For this purpose we artificially include clumping into our stationary wind models. When the clumps are assumed to be optically thin, the radiative line force increases compared to corresponding unclumped models, with a similar effect on either the mass-loss rate or the terminal velocity (depending on the onset of clumping). Optically thick clumps, alternatively, might be able to decrease the radiative force.

Recent studies of massive O-type stars present clear evidences of inhomogeneous and clumped winds. O-type (H-rich) central stars of planetary nebulae (CSPNs) are in some ways the low mass-low luminosity analogous of those massive stars. In this contribution, we present preliminary results of our on-going multi-wavelength (FUV, UV and optical) study of the winds of Galactic CSPNs. Particular emphasis will be given to the clumping factors derived by means of optical lines (Halpha and HeII4686) and "classic" FUV (and UV) lines.

We present extensive optical spectroscopy of the early-type magnetic star HD 191612 (O6.5f?pe-O8fp). The Balmer and HeI lines show strongly variable emission which is highly reproducible on a well-determined 538-d period. Metal lines and HeII absorptions (including many selective emission lines but excluding He II 4686A emission) are essentially constant in line strength, but are variable in velocity, establishing a double-lined binary orbit with P(orb) = 1542d, e=0.45. We conduct a model-atmosphere analysis of the primary, and find that the system is consistent with a O8: giant with a B1: main-sequence secondary. Since the periodic 538-d changes are unrelated to orbital motion, rotational modulation of a magnetically constrained plasma is strongly favoured as the most likely underlying `clock'. An upper limit on the equatorial rotation is consistent with this hypothesis, but is too weak to provide a strong constraint.

[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.

Optical high-resolution spectra of five central stars of planetary nebulae (CSPN) in the Galactic Bulge have been obtained with Keck/HIRES in order to derive their parameters. Since the distance of the objects is quite well known, such a method has the advantage that stellar luminosities and masses can in principle be determined without relying on theoretical relations between both quantities. By alternatively combining the results of our spectroscopic investigation with evolutionary tracks, we obtain so-called spectroscopic distances, which can be compared with the known (average) distance of the Bulge-CSPN. This offers the possibility to test the validity of model atmospheres and present date post-AGB evolution. We analyze optical H/He profiles of five Galactic Bulge CSPN (plus one comparison object) by means of profile fitting based on state of the art non-LTE modeling tools, to constrain their basic atmospheric parameters (Teff, log g, helium abundance and wind strength). Masses and other stellar radius dependent quantities are obtained from both the known distances and from evolutionary tracks, and the results from both approaches are compared. The major result of the present investigation is that the derived spectroscopic distances depend crucially on the applied reddening law. An "average extinction law" leads to a distance of 10.7 +- 1.2 kpc, which is considerably larger than the Galactic Center distance of 8 kpc. However, we find a remarkable internal agreement of the individual spectroscopic distances of our sample objects. Due to the uncertain reddening correction, the analysis presented here cannot yet be regarded as a consistency check for our method, and a rigorous test of the CSPN evolution theory becomes only possible if this problem has been solved.

We have studied a sample of 28 periodically variable B-type supergiants selected from the HIPPARCOS mission and 12 comparison stars covering the whole B-type spectral range. Our goal is to test if their variability is compatible with opacity-driven non-radial oscillations. We have used the NLTE atmosphere code FASTWIND to derive the atmospheric and wind parameters of the complete sample through line profile fitting. We applied the method to selected H, He and Si line profiles, measured with the high resolution CES spectrograph attached to the ESO CAT telescope in La Silla, Chile. From the location of the stars in the (log Teff, log g) diagram, we suggest that variability of our sample supergiants is indeed due to the gravity modes resulting from the opacity mechanism. We find nine of the comparison stars to be periodically variable as well, and suggest them to be new alpha Cyg variables. We find marginal evidence of a correlation between the amplitude of the photometric variability and the wind density. We investigate the Wind Momentum Luminosity Relation for the whole range of B spectral type supergiants, and find that the later types (> B5) perfectly follow the relation for A supergiants. Additionally, we provide a new spectral type - Teff calibration for B supergiants. Our results imply the possibility to probe internal structure models of massive stars of spectral type B through seismic tuning of gravity modes.

In the last few years our knowledge of the physics of massive stars has improved tremendously. However, further investigations are still needed, especially regarding accurate calibrations of their fundamental parameters. To this end, we have constructed a comprehensive grid of NLTE model atmospheres and corresponding synthetic spectra in the massive star domain. The grid covers the complete B type spectral range, extended to late O on the hot side and early A on the cool side, from supergiants to dwarfs and from weak stellar winds to very strong ones. It has been calculated with the latest version of the FASTWIND code. The analysis of an extensive sample of OB stars in the framework of the COROT space mission will lead to accurate calibrations of effective temperatures, gravities, mass loss rates etc. This paper contains a detailed description of the grid, which has been baptised as BSTAR06 and which will be available for further research in the near future.

The progress over the last years in modelling the atmospheres and winds of PN central stars is reviewed. We discuss the effect of the inclusion of the blanketing by millions of metal lines in NLTE on the diagnostics of photospheric and stellar wind lines, which can be used to determine stellar parameters such as effective temperature, gravity, radius, mass loss rate and distance. We also refer to recent work on the winds of massive O-type stars, which indicates that their winds are possibly inhomogeneous and clumped. We investigate implications from this work on the spectral diagnostics of PN central stars and introduce a method to determine wind clumping factors from the relative strengths of Halpha and HeII 4686. Based on new results we discuss the wind properties of CSPN.

We discuss recent evidence that currently accepted mass-loss rates may need to be revised downwards, as a consequence of previously neglected ``clumping'' of the wind. New results on the radial stratification of the corresponding clumping factors are summarized. We investigate the influence of clumping on the ionization equilibrium of phosphorus, which is of major relevance when deriving constraints on the clumping factors from an analysis of the FUV PV resonance lines.

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.

Recent calculations and analyses of O star spectra have revealed discrepancies between theory and observations, and between different theoretical calculations, for the strength of optical HeI singlet transitions.We investigate the source of these discrepancies. Using a non-LTE radiative transfer code we have undertaken detailed test calculations for a range of O star properties. Our principal test model has parameters similar to those of the O9V star, 10 Lac. We show that the discrepancies arise from uncertainties in the radiation field in the HeI resonance transition near 584Angs. The radiation field at 584Angs. is influenced by model assumptions, such as the treatment of line-blanketing and the adopted turbulent velocity, and by the FeIV atomic data. It isshown that two FeIV transitions near 584Angs can have a substantial influence on the strength of the HeI singlet transitions. Because of the difficulty of modeling the HeI singlet lines, particularly in stars with solar metalicity, the HeI triplet lines should be preferred in spectral analyses. These lines are much less sensitive to model assumptions.

Recent results strongly challenge the canonical picture of massive star winds: various evidence indicates that currently accepted mass-loss rates, Mdot, may need to be revised downwards significantly. This is because the most commonly used mass-loss diagnostics are affected by ``clumping'' (small-scale density inhomogeneities), influencing our interpretation of observed spectra and fluxes. Such downward revisions would have dramatic consequences for the evolution of, and feedback from, massive stars, and thus robust determinations of the clumping properties and mass-loss rates are urgently needed. Here, we present a first attempt to constrain the radial stratification of the so-called clumping factor. To this end, we have analyzed a sample of 19 Galactic O-type supergiants/giants, by combining data for Halpha, IR, mm and radio fluxes, and using appropriate analysis methods. Clumping has been included into our analysis in the ``conventional'' way, by assuming the inter-clump matter to be void. Because (almost) all our diagnostics depends on the square of density, we cannot derive absolute clumping factors, but only factors normalized to a certain minimum. This minimum was usually found to be located in the outermost, radio-emitting region, i.e., the radio mass-loss rates are the lowest ones, compared to Mdot derived from Halpha and the IR. The radio rates agree well with those predicted by theory, but are only upper limits, due to unknown clumping in the outer wind. Our most important result concerns a (physical) difference between denser and thinner winds: for denser winds, the innermost region is more strongly clumped than the outermost one (with a normalized clumping factor of 4.1+/-1.4), whereas thinner winds have similar clumping properties in the inner and outer regions.

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.

We present intermediate resolution (R ~ 8,000 - 12,000) high signal-to-noise H- and K-band spectroscopy of a sample of 37 optically visible stars, ranging in spectral type from O3 to B3 and representing most luminosity classes. Spectra of this quality can be used to constrain the temperature, luminosity and general wind properties of OB stars, when used in conjunction with sophisticated atmospheric model codes. Most important is the need for moderately high resolutions (R > 5000) and very high signal-to-noise (S/N > 150) spectra for a meaningful profile analysis. When using near-infrared spectra for a classification system, moderately high signal-to-noise (S/N ~ 100) is still required, though the resolution can be relaxed to just a thousand or two. In the appendix we provide a set of very high quality near-infrared spectra of Brackett lines in six early-A dwarfs. These can be used to aid in the modeling and removal of such lines when early-A dwarfs are used for telluric spectroscopic standards.

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 investigate the line-profile variability (lpv) of H-alpha for a large sample of O-type supergiants. We used the Temporal Variance Spectrum (TVS) analysis, developed by Fullerton et al 1996 and modified by us to take into account the effects of wind emission. By means of a comparative analysis we put a number of constraints on the properties of the variability as a function of stellar and wind parameters. The results of our analysis show that all the stars in the sample show evidence of significant lpv in H-alpha, mostly dominated by processes in the wind. The variations occur between zero and 0.3 v_inf (i.e., below ~1.5 R_star), in good agreement with the results from similar studies. A comparison between observations and line-profile simulations indicates that for stars with intermediate wind densities the H-alpha variability can be explained by simple models, consisting of coherent or broken shells (blobs) uniformly distributed over the wind volume, with an intrinsic scatter in the maximum density contrast of about a factor of two. For stars at lower and higher wind densities, on the other hand, we found certain inconsistencies between observations and our predictions, most importantly concerning the mean amplitude and the symmetry properties of the TVS. This disagreement might be explained with the presence of coherent large-scale structures, partly confined in a volume close to the star. Interpreted in terms of a variable mass-loss rate, the observed variations of H-alpha indicate changes of 4% with respect to the mean value of M_dot for stars with stronger winds and of 16% for stars with weaker winds. The effect of these variations on the corresponding wind momenta is rather insignificant (<0.16 dex), increasing the local scatter without affecting the Wind Momentum Luminosity Relationship.

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 analyze the optical and UV spectra of an additional sample of 20 Magellanic Cloud O stars, and draw conclusions from the complete sample of 40 stars. We find (1) The SMC O3-7 dwarfs are about 4000 K cooler than their Galactic counterparts; this is in the sense expected from the decreased signficiance of line-blanketing and wind-blanketing at lower metallicities. The difference decreases with later types, becoming negligible by B0 V. A similar difference is found for the supergiants. (2) The wind momentum of these stars scales with luminosity and metallicty in the ways predicted by radiatively-driven wind theory. (3) A comparison of the masses derived from spectroscopy with those derived from stellar evolutionary theory shows a significant discrepancy for stars hotter than 45000, although good agreement is found for cooler stars. (4) For the hottest O stars (O2-3.5) neither the NIII/NIV ratio, nor even the HeI/HeII ratio, does a good job of predicting the effective temperature by itself. Instead, a full analysis is needed to derive physical parameters. Thus there are O3.5V stars which are as hot or hotter than stars classified as O2V. (5) The two stars with the most discordant radial velocities in our sample happen to be O3 "field stars". This provides the first compelling observational evidence that the "field" O stars in the Magellanic Clouds may be runaway OB stars, ejected from their birth place.

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 have obtained HST and ground-based observations of a sample of 20 O-type stars in the LMC and SMC, including six of the hottest massive stars known (subtypes O2-3) in the R136 cluster. In general, these data include (a) the HST UV spectra in order to measure the terminal velocities of the stellar winds, (b) high signal-to-noise, blue-optical data where the primary temperature- and gravity-sensitive photospheric lines are found, and (c) nebular-free H-alpha profiles, which provide the mass-loss rates. The line-blanketed non-LTE atmosphere code FASTWIND was then used to determine the physical parameters of this sample of stars. We find good agreement between the synthetic line profiles for the hydrogen, He I, and He II lines in the majority of the stars we analyzed; the three exceptions show evidence of being incipiently resolved spectroscopic binaries or otherwise spectral composites. One such system is apparently an O3 V+O3 V eclipsing binary, and a follow-up radial velocity study is planned to obtain Keplerian masses. Although we did not use them to constrain the fits, good agreement is also found for the He I $\lambda 3187$ and He II $\lambda 3203$ lines in the near-UV, which we plan to exploit in future studies. Our effective temperatures are compared to those recently obtained by Repolust, Puls & Herrero for a sample of Galactic stars using the same techniques. We find that the Magellanic Cloud sample is 3,000-4,000$^\circ$K hotter than their Galactic counterparts for the early through mid-O's. These higher temperatures are the consequence of a decreased importance of wind emission, wind blanketing, and metal-line blanketing at lower metallicities.

Spectroscopic analyses of 7 SMC B-type supergiants and 1 giant have been undertaken using high resolution optical data obtained on the VLT with UVES. FASTWIND, a non-LTE, spherical,line-blanketed model atmosphere code was used to derive atmospheric and wind parameters of these stars as well as their absolute abundances. The implications of these results for stellar evolution and line driven wind theory are discussed.

The spectra of two early B-type supergiant stars in the Sculptor spiral galaxy NGC 300 are analysed by means of non-LTE line blanketed unified model atmospheres, aimed at determining their chemical composition and the fundamental stellar and wind parameters. For the first time a detailed chemical abundance pattern (He, C, N, O, Mg and Si) is obtained for a B-type supergiant beyond the Local Group. The derived stellar properties are consistent with those of other Local Group B-type supergiants of similar types and metallicities. One of the stars shows a near solar metallicity while the other one resembles more a SMC B supergiant. The effects of the lower metallicity can be detected in the derived wind momentum.

We have analyzed six OB supergiants and one giant covering spectral types from O3 to B1 in the Galactic OB association Cyg OB2 by means of an updated version of FASTWIND (SAntolaya-Rey, Puls, Herrero, 1997, A&A 323, 348) that includes an approximate treatment of metal line blocking and blanketing. This large coverage in spectral type allows us to derive a new temperature scale for Galactic O supergiants that is lower than the one obtained by using pure H--He models, either plane-parallel and hydrostatic or spherical with mass-loss. The lower temperatures are thus a combined effect of line blanketing and the large mass-loss rates. In some cases, the newly derived effective temperature is reduced by up to 8000 K. Changes are larger for earlier stars with large mass--loss rates. As a consequence, luminosities are modified as well, which results in a lower number of emerging ionizing photons and reduces the mass discrepancy. Although there are still significant differences between spectroscopic and evolutionary masses, we do not find any obvious systematic pattern of those differences. We derive mass--loss rates and the corresponding wind momentum--luminosity relation for the analyzed stars. Although consistent with previous results by Puls et al. (1996, A&A 305, 171) for Galactic stars, our relation is better defined due to a reduction of errors related to stellar distances and points to a possible separation between extreme Of stars (Of+, Of*) and stars with more moderate morphologies. However this finding is only tentative, as the statistics are still scarce.

We analyze STIS spectra in the 1150-1700 Angstrom wavelength range obtained for six early B supergiants in the neighboring galaxy M31. Because of their likely high (nearly solar) abundance, these stars were originally chosen to be directly comparable to their Galactic counterparts, and represent a much-needed addition to our current sample of B-type supergiants, in our efforts to study the dependence of the Wind Momentum-Luminosity Relationship on spectral type and metallicity. As a first step to determine wind momenta we fit the P-Cygni profiles of the resonance lines of N V, Si IV and C IV with standard methods, and derive terminal velocities for all of the STIS targets. From these lines we also derive ionic stellar wind column densities. Our results are compared with those obtained previously in Galactic supergiants, and confirm earlier claims of `normal' wind line intensities and terminal velocities in M31. For half of the sample we find evidence for an enhanced maximum turbulent velocity when compared to Galactic counterparts.

We present terminal velocities of M33 B-supergiants, obtained from STIS HST spectra as part of our programme to investigate the Wind Momentum -- Luminosity Relationship (WLR) in the Local Group. Terminal velocities are derived from their NV, CIV, and SiIV resonance lines in UV spectra. Comparing with IUE spectra of Galactic B-supergiants we found evidence of low metallicity in three of our objects. %%The terminal velocities are found to be consistent with the %%values found for Galactic stars, except for B-133. For this The terminal velocities are consistent with the corresponding values of Galactic stars, except for B-133. For this star we find a very large vterm and a red SiIV component deeper than the blue one, that might be an indication of binarity. The average ratio between terminal and turbulent wind velocities is 0.25, well above the value found for Galactic stars.

We present observations and analyses of seven Galactic O stars of type O6 and earlier. Analyses are carried out using NLTE plane-parallel, hydrostatic models as well as NLTE spherical models with mass-loss. With detailed calculations for the former and simulations for the latter, it is shown that the flux blocking due to UV metal lines is important for these objects and the way the mechanism operates is explained. We find that the plane-parallel, hydrostatic unblanketed model atmospheres have increasing difficulties in fitting the early-type spectra of massive stars, and for 50000 K and above a fit seems to be impossible. We then perform an analysis using spherical models with mass-loss. Gravities should be increased by 0.1-0.25 dex, reducing, but not solving, the mass discrepancy. We show that spectroscopic masses are in better agreement with the theory of radiatively driven winds than evolutionary masses are. A helium abundance larger than solar is obtained for most objects. Some additional effects (partly related to present approximations) that have an influence in our analyses are studied. Some stars of our sample have such high mass-loss rates that the derivation of gravities from the wings of Balmer lines, in particular Hgamma, becomes doubtful. From the point of view of individual stars, we have analysed some of the most massive and luminous stars in the Milky Way. According to our analysis, three of them (Cyg OB2-7, HD 15570 and HD15558) have particularly large initial masses, close to or in excess of 100 Msun. HD 5689 could have been erroneously assigned to Cas OB7 and might be a runaway star.

The Balmer lines of four A Ia - supergiants (spectral type A0 to A3) and fourteen B Ia and Ib - supergiants (spectral type B0 to B3) in the solar neighbourhood are analyzed by means of NLTE unified model atmospheres to determine the properties of their stellar winds, in particular their wind momenta. As in previous work for O-stars (Puls et al.) a tight relationship between stellar wind momentum and luminosity (`WLR'') is found. However, the WLR varies as function of spectral type. Wind momenta are strongest for O-supergiants, then decrease from early B (B0 and B1) to mid B (B1.5 to B3) spectral types and become stronger again for A-supergiants. The slope of the WLR appears to be steeper for A- and mid B-supergiants than for O-supergiants. The spectral type dependence is interpreted as an effect of ionization changing the effective number and the line strength distribution function of spectral lines absorbing photon momentum around the stellar flux maximum. This interpretation needs to be confirmed by theoretical calculations for radiation driven winds. A-supergiants in M31 observed with HIRES at the Keck telescope have wind momenta compatible with their galactic counterparts. The potential of the WLR as a new, independent extragalactic distance indicator is discussed. It is concluded that with ten to twenty objects, photometry with HST and medium resolution spectroscopy with 8m-telescopes from the ground distance moduli can be obtained with an accuracy of about 0.1m out to the Virgo and Fornax clusters of galaxies.