The apparent lack of massive O-type stars near the zero-age main sequence (at ages < 2 Myr) is a topic widely discussed. Different explanations for this elusive detection have been proposed, but no firm conclusions have been reached yet. We reassess this empirical result benefiting from the high-quality spectroscopic observations of >400 Galactic O-type stars gathered by the IACOB and OWN surveys. We used temperatures and gravities from a iacob-gbat/fastwind spectroscopic analysis to locate our sample in the Kiel and spectroscopic HR diagrams. We evaluated the completeness of our sample of stars, observational biases using information from the Galactic O star catalog (GOSC), systematics of our methodology, and compare with other recent studies using smaller samples of Galactic O-type stars. We base our discussion on the spectroscopic HR diagram to avoid the use of uncertain distances. We performed a detailed study of the young cluster Trumpler-14 as an example of how Gaia cluster distances can help to construct the associated classical HR diagram. The apparent lack of massive O-type stars near the ZAMS with masses between 30 and 70 Msol persist even when spectroscopic results from a large, non-biased sample of stars are used. We do not find correlation between the dearth of stars and observational biases, limitations of our methodology, or the use of spectroscopic HR diagram instead of the classical one. Investigating the efficiency of mass accretion during the formation process we conclude that an adjustment of the accretion rate towards lower values could reconcile the hotter boundary of detected O-type stars and the theoretical birthline. Last, we discuss that the presence of a small sample of O2-O3.5 stars found closer to the ZAMS might be explained taking into account non-standard star evolution (e.g. binary interaction, mergers, or homogeneous evolution).

Context. Massive stars are predicted to excite internal gravity waves (IGWs) by turbulent core convection and from turbulent pressure fluctuations in their near-surface layers. These IGWs are extremely efficient at transporting angular momentum and chemical species within stellar interiors, but they remain largely unconstrained observationally. Aims. We aim to characterise the photometric detection of IGWs across a large number of O and early-B stars in the Hertzsprung-Russell diagram, and explain the ubiquitous detection of stochastic variability in the photospheres of massive stars. Methods. We combined high-precision time-series photometry from the NASA Transiting Exoplanet Survey Satellite with high-resolution ground-based spectroscopy of 70 stars with spectral types O and B to probe the relationship between the photometric signatures of IGWs and parameters such as spectroscopic mass, luminosity, and macroturbulence. Results. A relationship is found between the location of a star in the spectroscopic Hertzsprung-Russell diagram and the amplitudes and frequencies of stochastic photometric variability in the light curves of massive stars. Furthermore, the properties of the stochastic variability are statistically correlated with macroturbulent velocity broadening in the spectral lines of massive stars. Conclusions. The common ensemble morphology for the stochastic low-frequency variability detected in space photometry and its relationship to macroturbulence is strong evidence for IGWs in massive stars, since these types of waves are unique in providing the dominant tangential velocity field required to explain the observed spectroscopy.

Lack of high-precision long-term continuous photometric data for large samples of stars has prevented the large-scale exploration of pulsational variability in the OB star regime. As a result, the candidates for in-depth asteroseismic modelling remained limited to a few tens of dwarfs. The TESS nominal space mission has surveyed the southern sky, yielding continuous data of at least 27 d for hundreds of OB stars. We aim to couple TESS data in the southern sky with spectroscopy to study the variability over mass and evolution. We focus mainly on the presence of coherent pulsation modes that may or may not be present in the theoretical instability domains and unravel all frequency behaviour in the amplitude spectra of the TESS data. We compose a sample of 98 OB-type stars observed by TESS in Sectors 1-13 and with available high-resolution spectroscopy gathered by the IACOB and OWN surveys. We present the short-cadence 2-min light curves of dozens of OB-type stars, that have one or more spectra in the IACOB or OWN database. Based on these light curves and their Lomb-Scargle periodograms we perform variability classification and frequency analysis, and place the stars in the spectroscopic Hertzsprung-Russell diagram to interpret the variability in an evolutionary context. We deduce diverse origins of the variability found in all of the 98 OB stars in the TESS data. Among these we find several new variable stars, including three hybrid pulsators, three eclipsing binaries, high frequency modes in a Be star, and potential heat-driven pulsations in two Oe stars. We identify stars for which future asteroseismic modelling is possible, provided mode identification is achieved. By comparing the position of the variables to theoretical instability strips we discuss the current shortcomings in non-adiabatic pulsation theory, and the distribution of pulsators in the upper Hertzsprung-Russell diagram.

Quantitative spectroscopy is a powerful technique from which we can extract information about the physical properties and surface chemical composition of stars. In this chapter, I guide the reader through the main ideas required to get initiated in the learning process to become an expert in the application of state-of-the-art quantitative spectroscopic techniques to the study of massive OB stars. NB: This chapter is intended to serve to young students as a first approach to a field which has attracted my attention during the last 20 years. I should note that, despite its importance, at present, the number of real experts in the field around the world is limited to less than 50 people, and about one third of them are close to retirement. Hence, I consider that this is a good moment to write a summary text on the subject to serve as guideline for the next generations of students interested in joining the massive star crew. If you are one of them, please, use this chapter as a first working notebook. Do not stop here. Dig also, for further details, into the literature I quote along the text. And, once there, dig even deeper to find all the original sources explaining in more detail the physical and technical concepts that are presently incorporated into our modern (almost) automatized tools.

Being relatively rare, the properties of magnetic O stars are not fully understood. To date fewer than a dozen of these stars have been confirmed, making any inference of their global properties uncertain due to small number statistics. To better understand these objects it is necessary to increase the known sample. The MOBSTER collaboration aims to do this by identifying candidate magnetic O, B, and A stars from the identification of rotational modulation in high-precision photometry from the Transiting Exoplanet Survey Satellite (TESS). Here we discuss the collaboration's efforts to detect rotational modulation in TESS targets to identify candidate magnetic O stars for future spectropolarimetric observations.

LB-1 has recently been proposed to be a binary system at 4 kpc consisting of a B star of 8 Msol and a massive stellar black hole of 70 Msol. This finding challenges our current theories of massive star evolution and formation of BHs at solar metallicity. Our objective is to derive the effective temperature, surface gravity and chemical composition of the B-type component in order to determine its nature and evolutionary status and, indirectly, to constrain the mass of the BH. We use the non-LTE stellar atmosphere code FASTWIND to analyse new and archival high resolution data. We determine (Teff, logg) values of (14000$\pm500$ K, 3.50$\pm0.15$ dex) that, combined with the Gaia parallax, implies a spectroscopic mass, from logg, of $3.2^{+2.1}_{-1.9}$ Msol and an evolutionary mass, assuming single star evolution, of $5.2^{+0.3}_{-0.6}$ Msol. We determine an upper limit of 8 km/s for the projected rotational velocity and derive the surface abundances, finding the star to have a silicon abundance below solar, to be significantly enhanced in nitrogen and iron, and depleted in carbon and magnesium. Complementary evidence derived from a photometric extinction analysis and Gaia yields similar results for Teff and logg and a consistent distance around 2~kpc. We propose that the B star is a slightly evolved main sequence star of 3-5 Msol with surface abundances reminiscent of diffusion in late B/A chemically peculiar stars with low rotational velocities. There is also evidence for CN-processed material in its atmosphere. These conclusions rely critically on the distance inferred from the Gaia parallax. The goodness of fit of the Gaia astrometry also favours a high-inclination orbit. If the orbit is edge-on and the B star has a mass of 3-5 Msol, the mass of the dark companion would be 4-5 Msol, which would be easier to explain with our current stellar evolutionary models.

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

Almost all massive stars explode as supernovae and form a black hole or neutron star. The remnant mass and the impact of the chemical yield on subsequent star formation and galactic evolution strongly depend on the internal physics of the progenitor star, which is currently not well understood. The theoretical uncertainties of stellar interiors accumulate with stellar age, which is particularly pertinent for the blue supergiant phase. Stellar oscillations represent a unique method of probing stellar interiors, yet inference for blue supergiants is hampered by a dearth of observed pulsation modes. Here we report the detection of diverse variability in blue supergiants using the K2 and TESS space missions. The discovery of pulsation modes or an entire spectrum of low-frequency gravity waves in these stars allow us to map the evolution of hot massive stars towards the ends of their lives. Future asteroseismic modelling will provide constraints on ages, core masses, interior mixing, rotation and angular momentum transport. The discovery of variability in blue supergiants is a step towards a data-driven empirical calibration of theoretical evolution models for the most massive stars in the Universe.

[ABRIDGED] AIMS. MONOS is collecting information and studying O-type spectroscopic binaries with delta > -20 deg. In this 1st paper we describe the sample and provide spectral classifications and additional information for objects with previous SB or EB orbits. In the future we will test the validity of previous solutions and calculate new SB orbits. METHODS.The spectra in this paper have 2 sources: GOSSS, which is obtaining blue-violet R~2500 spectroscopy of thousands of massive stars, and LiLiMaRlin, a library of libraries of high-resolution spectroscopy of massive stars obtained from 4 different surveys and additional data from our own observing programs and public archives. We also use lucky images from AstraLux. RESULTS. We present homogeneous spectral classifications for 92 O-type spectroscopic multiple systems and 10 optical companions. We discuss the visual multiplicity of each system with the support of AstraLux images and additional sources. For 11 O-type objects and for 6 B-type objects we present their first GOSSS spectral classifications. For 2 known EBs we detect SB2/SB1 lines for the first time, to which we add a third system already reported by us recently. For 2 previous SB1 systems we detect their SB2 nature for the first time and give their first separate spectral classifications, something we also do for a 3rd object just recently identified as a SB2. We also detect 9 new astrometric companions and provide updated information on several others. For sigma Ori AaAbB we provide spectral classifications for the 3 components with a single observation for the first time thanks to a lucky spectroscopy observation obtained close to the Aa,Ab periastron and for theta^1 Ori CaCb we add it to the class of Galactic Of?p stars, raising the number of its members to 6. Our sample of O-type spectroscopic binaries contains more triple/higher-order systems than double systems.

We present a description of "the Conductor", an automated software package which handles all observations at the first Stellar Observations Network Group (SONG) node telescope at the Teide Observatory on the island of Tenerife. The idea was to provide a complete description on the automated procedures for target selection and execution of observations and to demonstrate how the SONG robotic telescope is being operated. The Conductor is a software package developed in Python and running on a server in Aarhus which makes use of a large set of database tables through which it communicates with the SONG nodes. Based on a number of selection criteria the Conductor is able to identify the optimum target to be observed at any given moment taking into account local weather conditions and technical constraints. The Conductor has made it possible for the Hertzsprung SONG telescope to become a highly sophisticated and efficient robotic telescopic facility without human interaction. It can handle everything from Principal Investigators submitting their proposed targets with specific settings to the data being available for download after the observations has been carried out. At present, and thanks to the availability of the Conductor, the first node of the SONG project can be considered a world leading robotic telescope with respect to needed human interactions, efficiency and flexibility in observing strategy.

Boyajian's star is an apparently normal main sequence F-type star with a very unusual light curve. The dipping activity of the star, discovered during the Kepler mission, presents deep, asymmetric, and aperiodic events. Here we present high resolution spectroscopic follow-up during some dimming events recorded post-Kepler observations, from ground-based telescopes. We analise data from the HERMES, HARPS-N and FIES spectrographs to characterise the stellar atmosphere and to put some constraints on the hypotheses that have appeared in the literature concerning the occulting elements. The star's magnetism, if existing, is not extreme. The spots on the surface, if present, would occupy 0.02% of the area, at most. The chromosphere, irrespective of the epoch of observation, is hotter than the values expected from radiative equilibrium, meaning that the star has some degree of activity. We find no clear evidence of the interstellar medium nor exocoments being responsible for the dimmings of the light curve. However, we detect at 1-2 sigma level, a decrease of the radial velocity of the star during the first dip recorded after the \emph{\emph{Kepler}} observations. We claim the presence of an optically thick object with likely inclined and high impact parameter orbits that produces the observed Rossiter-McLaughlin effect.

Early-type stars are predicted to excite an entire spectrum of internal gravity waves (IGWs) at the interface of their convective cores and radiative envelopes. Numerical simulations of IGWs predict stochastic low-frequency variability in photometric observations, yet the detection of IGWs in early-type stars has been limited by a dearth of high-quality photometric time series. We present observational evidence of stochastic low-frequency variability in the CoRoT photometry of a sample of O, B, A and F stars. The presence of this stochastic low-frequency variability in stars across the upper main-sequence cannot be universally explained as granulation or stellar winds, but its morphology is found to be consistent with predictions from IGW simulations.

Context. Main sequence stars with a convective core are predicted to stochastically excite Internal Gravity Waves (IGWs), which effectively transport angular momentum throughout the stellar interior and explain the observed near-uniform interior rotation rates of intermediate-mass stars. However, there are few detections of IGWs, and fewer still made using photometry, with more detections needed to constrain numerical simulations. Aims. We aim to formalise the detection and characterisation of IGWs in photometric observations of stars born with convective cores (M > 1.5 M$_{\odot}$) and parameterise the low-frequency power excess caused by IGWs. Methods. Using the most recent CoRoT light curves for a sample of O, B, A and F stars, we parameterise the morphology of the flux contribution of IGWs in Fourier space using an MCMC numerical scheme within a Bayesian framework. We compare this to predictions from IGW numerical simulations and investigate how the observed morphology changes as a function of stellar parameters. Results. We demonstrate that a common morphology for the low-frequency power excess is observed in early-type stars observed by CoRoT. Our study shows that a background frequency-dependent source of astrophysical signal is common, which we interpret as IGWs. We provide constraints on the amplitudes of IGWs and the shape of their detected frequency spectrum across a range of mass, which is the first ensemble study of stochastic variability in such a diverse sample of stars. Conclusions. The evidence of a low-frequency power excess across a wide mass range supports the interpretation of IGWs in photometry of O, B, A and F stars. We also discuss the prospects of observing hundreds of massive stars with the Transiting Exoplanet Survey Satellite (TESS) in the near future.

LiLiMaRlin is a library of libraries of massive-star high-resolution optical spectra built by collecting data from [a] our spectroscopic surveys (OWN, IACOB. NoMaDS, and CAF\'E-BEANS) and programs and [b] searches in public archives. The current version has 18 077 spectra of 1665 stars obtained with seven different telescopes (HET 9.2 m, NOT 2.56 m, CAHA 2.2 m, MPG/ESO 2.2 m, OHP 1.93 m, Mercator 1.2 m, and Stella 1.2 m). All the spectra have been filtered to eliminate misidentifications and bad-quality ones, uniformly reprocessed, and placed on a common format. We present applications of this library of libraries to the analysis of spectroscopic binaries (OWN and MONOS, see poster by E. Trigueros P\'aez at this meeting) and the study of the interstellar medium (CollDIBs). We discuss our plans for the future.

Context. Most massive stars are in binary or multiple systems. Several massive stars have been detected as doublelined spectroscopic binaries and among these, the OWN Survey has detected a non-negligible number whose components show very different spectral line broadening (i.e., projected rotational velocities). This fact raises a discussion about the contributing processes, such as angular-momentum transfer and tidal forces. Aims. We seek to constrain the physical and evolutionary status of one of such systems, the O+O binary HD 93343. Methods. We analyzed a series of high-resolution multiepoch optical spectra to determine the orbital parameters, projected rotational velocities, and evolutionary status of the system. Results. HD 93343 is a binary system comprised of two O7.5 Vz stars that each have minimum masses of approximately 22 Mo in a wide and eccentric orbit (e = 0.398$\pm$0.004; P=50.432$\pm$0.001 d). Both stars have very similar stellar parameters, and hence ages. As expected from the qualitative appearance of the combined spectrum of the system, however, these stars have very different projected rotational velocities (~65 and ~325 km/s, respectively). Conclusions. The orbits and stellar parameters obtained for both components seem to indicate that their youth and relative separation is enough to discard the effects of mass transfer and tidal friction. Thus, non-synchronization should be intrinsic to their formation.

Cygnus OB2 is a rich OB association in the Galaxy which has experienced intense star formation in the last 20-25 Myr. Its stellar population shows a correlation between age and Galactic longitude. Exploring the chemical composition of its stellar content we will be able to check the degree of homogeneity of the natal molecular cloud and possible effects of self-enrichment processes. Our aim is to determine silicon and oxygen abundances for a sample of eight early-type slow rotators in Cygnus OB2 in order to check possible inhomogeneities across the whole association and whether there exists a correlation of chemical composition with Galactic longitude. We have performed a spectroscopic analysis of a sample of late O and early B stars with low rotational velocity, which have been chosen so as to cover the whole association area. We have carried out an analysis based on equivalent widths of metal lines, the wings of the H Balmer lines and FASTWIND stellar atmosphere models to determine their stellar fundamental parameters as well as the silicon and oxygen surface abundances. We derive a rather homogeneous distribution of silicon and oxygen abundances across the region, with average values of 12+log(Si/H)=7.53$\pm$0.08 dex and 12+log(O/H)=8.65$\pm$0.12 dex. We find a homogeneous chemical composition in Cygnus OB2 with no clear evidence for significant chemical self-enrichment, despite indications of strong stellar winds and possible supernovae during the history of the region. Comparison with different scenarios of chemical enrichment by stellar winds and supernovae point to star forming efficiencies not significantly above 10%. The degree of homogeneity that we find is consistent with the observed Milky Way oxygen gradient based on HII regions. We also find that the oxygen scatter within Cygnus OB2 is at least of the same order than among HII regions at similar Galactocentric distance.

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.

A previous spectroscopic study identified the very massive O2 III star VFTS 16 in the Tarantula Nebula as a runaway star based on its peculiar line-of-sight velocity. We use the Gaia DR2 catalog to measure the relative proper motion of VFTS 16 and nearby bright stars to test if this star might have been ejected from the central cluster, R136, via dynamical ejection. We find that the position angle and magnitude of the relative proper motion (0.338 +/- 0.046 mas/yr, or approximately 80 +\- 11 km/s) of VFTS 16 are consistent with ejection from R136 approximately 1.5 +/- 0.2 Myr ago, very soon after the cluster was formed. There is some tension with the presumed age of VFTS 16 that, from published stellar parameters, cannot be greater than 0.9 +0.3/-0.2 Myr. Older ages for this star would appear to be prohibited due to the absence of He I lines in its optical spectrum, since this sets a firm lower limit on its effective temperature. The dynamical constraints may imply an unusual evolutionary history for this object, perhaps indicating it is a merger product. Gaia DR2 also confirms that another very massive star in the Tarantula Nebula, VFTS 72 (alias BI253; O2 III-V(n)((f*)), is also a runaway on the basis of its proper motion as measured by Gaia. While its tangential proper motion (0.392 +/-0.062 mas/yr or 93 +/-15 km/s) would be consistent with dynamical ejection from R136 approximately 1 Myr ago, its position angle is discrepant with this direction at the 2$\sigma$ level. From their Gaia DR2 proper motions we conclude that the two ~100 solar mass O2 stars, VFTS 16 and VFTS72, are fast runaway stars, with space velocities of around 100 km/s relative to R136 and the local massive star population. The dynamics of VFTS16 are consistent with it having been ejected from R136, and this star therefore sets a robust lower limit on the age of the central cluster of ~1.3 Myr.

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

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.

CONTEXT.The first Gaia Data Release (DR1) significantly improved the previously available proper motions for the majority of the Tycho-2 stars. AIMS. We want to detect runaway stars using Gaia DR1 proper motions and compare our results with previous searches. METHODS. Runaway O stars and BA supergiants are detected using a 2-D proper-motion method. The sample is selected using Simbad, spectra from our GOSSS project, literature spectral types, and photometry processed using CHORIZOS. RESULTS. We detect 76 runaway stars, 17 (possibly 19) of them with no prior identification as such, with an estimated detection rate of approximately one half of the real runaway fraction. An age effect appears to be present, with objects of spectral subtype B1 and later having travelled for longer distances than runaways of earlier subtypes. We also tentatively propose that the fraction of runaways is lower among BA supergiants that among O stars but further studies using future Gaia data releases are needed to confirm this. The frequency of fast rotators is high among runaway O stars, which indicates that a significant fraction of them (and possibly a majority) is produced in supernova explosions.

CONTEXT: Many massive stars have nearby companions whose presence hamper their characterization through spectroscopy. AIMS: We want to obtain spatially resolved spectroscopy of close massive visual binaries to derive their spectral types. METHODS: We obtain a large number of short long-slit spectroscopic exposures of five close binaries under good seeing conditions, select those with the best characteristics, extract the spectra using multiple-profile fitting, and combine the results to derive spatially separated spectra. RESULTS: We demonstrate the usefulness of Lucky Spectroscopy by presenting the spatially resolved spectra of the components of each system, in two cases with separations of only ~0.3". Those are delta Ori Aa+Ab (resolved in the optical for the first time) and sigma Ori AaAb+B (first time ever resolved). We also spatially resolve 15 Mon AaAb+B, zeta Ori AaAb+B (both previously resolved with GOSSS, the Galactic O-Star Spectroscopic Survey), and eta Ori AaAb+B, a system with two spectroscopic B+B binaries and a fifth visual component. The systems have in common that they are composed of an inner pair of slow rotators orbited by one or more fast rotators, a characteristic that could have consequences for the theories of massive star formation.

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

Many massive stars are part of binary or higher multiplicity systems. The present work focusses on two higher multiplicity systems: HD17505A and HD206267A. Determining the fundamental parameters of the components of the inner binary of these systems is mandatory to quantify the impact of binary or triple interactions on their evolution. We analysed high-resolution optical spectra to determine new orbital solutions of the inner binary systems. After subtracting the spectrum of the tertiary component, a spectral disentangling code was applied to reconstruct the individual spectra of the primary and secondary. We then analysed these spectra with the non-LTE model atmosphere code CMFGEN to establish the stellar parameters and the CNO abundances of these stars. The inner binaries of these systems have eccentric orbits with e ~ 0.13 despite their relatively short orbital periods of 8.6 and 3.7 days for HD17505Aa and HD206267Aa, respectively. Slight modifications of the CNO abundances are found in both components of each system. The components of HD17505Aa are both well inside their Roche lobe, whilst the primary of HD206267Aa nearly fills its Roche lobe around periastron passage. Whilst the rotation of the primary of HD206267Aa is in pseudo-synchronization with the orbital motion, the secondary displays a rotation rate that is higher. The CNO abundances and properties of HD17505Aa can be explained by single star evolutionary models accounting for the effects of rotation, suggesting that this system has not yet experienced binary interaction. The properties of HD206267Aa suggest that some intermittent binary interaction might have taken place during periastron passages, but is apparently not operating anymore.

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.

Despite the important advances in space asteroseismology during the last decade, the early phases of evolution of stars with masses above $\sim$15 M$_{\odot}$ have been only vaguely explored up to now. Our goal is to detect, analyze and interpret variability in the early-B type supergiant HD\,2905 using long-term, ground based, high resolution spectroscopy. We gather a total of 1141 high-resolution spectra covering some 2900 days. We complement these observations with the $Hipparcos$ light curve, which includes 160 data points obtained during a time span of $\sim$1200 days. We investigate spectroscopic variability of up to 12 diagnostic lines by using the zero and first moments of the line profiles. We perform a frequency analysis of both the spectroscopic and photometric dataset. HD\,2905 is a spectroscopic variable with peak-to-peak amplitudes in the zero and first moments of the photospheric lines of up to 15\% and 30 \kms, respectively. The amplitude of the line-profile variability is correlated with the line formation depth in the photosphere and wind. All investigated lines present complex temporal behavior indicative of multi-periodic variability with timescales of a few days to several weeks. The Scargle periodograms of the {\it Hipparcos\} light curve and the first moment of purely photospheric lines reveal a low-frequency amplitude excess and a clear dominant frequency at $\sim$0.37 d$^{-1}$. In the spectroscopy, several additional frequencies are present in the range 0.1 - 0.4 d$^{-1}$. These may be associated with heat-driven gravity modes, convectively-driven gravity waves, or sub-surface convective motions. Additional frequencies are detected below 0.1 d$^{-1}$. In the particular case of H$\alpha$, these are produced by rotational modulation of a non-spherically symmetric stellar wind.

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.

We present C and O abundances in the Magellanic Clouds derived from deep spectra of HII regions. The data have been taken with the Ultraviolet-Visual Echelle Spectrograph at the 8.2-m VLT. The sample comprises 5 HII regions in the Large Magellanic Cloud (LMC) and 4 in the Small Magellanic Cloud (SMC). We measure pure recombination lines (RLs) of CII and OII in all the objects, permitting to derive the abundance discrepancy factors (ADFs) for O^2+, as well as their O/H, C/H and C/O ratios. We compare the ADFs with those of other HII regions in different galaxies. The results suggest a possible metallicity dependence of the ADF for the low-metallicity objects, but more uncertain for high-metallicity objects. We compare nebular and B-type stellar abundances and we find that the stellar abundances agree better with the nebular ones derived from collisionally excited lines (CELs). Comparing these results with other galaxies we observe that stellar abundances seem to agree better with the nebular ones derived from CELs in low-metallicity environments and from RLs in high-metallicity environments. The C/H, O/H and C/O ratios show almost flat radial gradients, in contrast with the spiral galaxies where such gradients are negative. We explore the chemical evolution analysing C/O vs. O/H and comparing with the results of HII regions in other galaxies. The LMC seems to show a similar chemical evolution to the external zones of small spiral galaxies and the SMC behaves as a typical star-forming dwarf galaxy.

We have detected 13 new runaway-star candidates of spectral type O combining the TGAS (Tycho-Gaia Astrometric Solution) proper motions from Gaia Data Release 1 (DR1) and the sample from GOSSS (Galactic O-Star Spectroscopic Survey). We have also combined TGAS and Hipparcos proper motions to check that our technique recovers many of the previously known O-type runaways in the sample.

The Gaia-ESO Survey (GES) is obtaining high quality spectra of $\sim$10$^{5}$ stars in our Galaxy, providing an homogeneous and unique overview of all the main components of the Milky Way, its formation history and the evolution of young, mature and ancient Galactic populations. Our group is in charge of the early-type massive stars that define the youngest population in the survey. In this contribution, we present the results of the quantitative spectroscopic analysis of O-type stars in the Carina Nebula within the Gaia-ESO Survey. For this aim, we have used FASTWIND and CMFGEN stellar atmosphere codes, providing stellar parameters for the current sample (GES data release iDR4).

The origin and evolution of magnetism in OB stars is far from being well understood. With approximately 70 magnetic OB stars known, any new object with unusual characteristics may turn out to be a key piece of the puzzle. We report the detection of an exceptionally strong magnetic field in the He-strong B2IV star CPD-62 2124. Spectropolarimetric FORS2 and HARPSpol observations were analysed by two independent teams and procedures, concluding on a strong longitudinal magnetic field of approximately 5.2 kG. The quantitative characterisation of the stellar atmosphere yields an effective temperature of 23650$\pm$250 K, a surface gravity of 3.95$\pm$0.10 dex and a surface helium fraction of 0.35$\pm$0.02 by number. The metal composition is in agreement with the cosmic abundance standard, except for Mg, Si and S, which are slightly non-solar. The strong and broad ($\sim$300 km/s) disc-like emission displayed by the H$\alpha$ line suggests a centrifugal magnetosphere supported by the strong magnetic field. Our results imply that CPD-62 2124 is an early B-type star hosting one of the strongest magnetic fields discovered to date, and one of the most evolved He-strong stars known, with a fractional main-sequence lifetime of approximately 0.6.

We present first results from the quantitative spectroscopic analysis of 266 Galactic O-type stars targeted by the IACOB and OWN surveys (implying the largest sample of stars of this type analyzed homogeneously). We also evaluate what is the present situation regarding available information about distances, as provided by the Hipparcos and Gaia missions.

Massive OB stars are critical to the ecology of galaxies, and yet our knowledge of OB stars in the Milky Way, fainter than $V \sim 12$, remains patchy. Data from the VST Photometric H$\alpha$ Survey (VPHAS+) permit the construction of the first deep catalogues of blue excess-selected OB stars, without neglecting the stellar field. A total of 14900 candidates with 2MASS cross-matches are blue-selected from a 42 square-degree region in the Galactic Plane, capturing the Carina Arm over the Galactic longitude range $282^{\circ} \lesssim \ell \lesssim 293^{\circ}$. Spectral energy distribution fitting is performed on these candidates' combined VPHAS+ $u,g,r,i$ and 2MASS $J,H,K$ magnitudes. This delivers: effective temperature constraints, statistically separating O from early-B stars; high-quality extinction parameters, $A_0$ and $R_V$ (random errors typically $< 0.1$). The high-confidence O-B2 candidates number 5915 and a further 5170 fit to later B spectral type. Spectroscopy of 276 of the former confirms 97% of them. The fraction of emission line stars among all candidate B stars is 7--8% . Greyer ($R_V > 3.5$) extinction laws are ubiquitous in the region, over the distance range 2.5--3 kpc to $\sim$10~kpc. Near prominent massive clusters, $R_V$ tends to rise, with particularly large and chaotic excursions to $R_V \sim 5$ seen in the Carina Nebula. The data reveal a hitherto unnoticed association of 108 O-B2 stars around the O5If$+$ star LSS 2063 ($\ell = 289.77^{\circ}$, $b = -1.22^{\circ}$). Treating the OB star scale-height as a constant within the thin disk, we find an orderly mean relation between extinction ($A_0$) and distance in the Galactic longitude range, $287.6^{\circ} < \ell < 293.5^{\circ}$, and infer the subtle onset of thin-disk warping. A halo around NGC 3603, roughly a degree in diameter, of $\sim$500 O-B2 stars with $4 < A_0 (\rm{mag}) < 7$ is noted.

Asteroseismology is a powerful tool to access the internal structure of stars. Apart from the important impact of theoretical developments, progress in this field has been commonly associated with the analysis of time-resolved observations. Recently, the so-called macroturbulent broadening has been proposed to be a complementary and less expensive way -- in terms of observational time -- to investigate pulsations in massive stars. We assess to what extent this ubiquitous non-rotational broadening component shaping the line profiles of O stars and B supergiants is a spectroscopic signature of pulsation modes driven by a heat mechanism. We compute stellar main sequence and post-main sequence models from 3 to 70Msun with the ATON stellar evolution code and determine the instability domains for heat-driven modes for degrees l=1-20 using the adiabatic and non-adiabatic codes LOSC and MAD. We use the observational material presented in Sim\'on-D\'iaz et al. (2016) to investigate possible correlations between the single snapshot line-broadening properties of a sample of ~260 O and B-type stars and their location inside/outside the various predicted instability domains. We present an homogeneous prediction for the non-radial instability domains of massive stars for degree l up to 20. We provide a global picture of what to expect from an observational point of view in terms of frequency range of excited modes, and investigate the behavior of the instabilities with stellar evolution and increasing degree of the mode. Furthermore, our pulsational stability analysis, once compared to the empirical results of Sim\'on-D\'iaz et al. (2016), indicates that stellar oscillations originated by a heat mechanism can not explain alone the occurrence of the large non-rotational line-broadening component commonly detected in the O star and B supergiant domain.

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.

During the summer of 2013, a 4-month spectroscopic campaign took place to observe the variabilities in three Wolf-Rayet stars. The spectroscopic data have been analyzed for WR 134 (WN6b), to better understand its behaviour and long-term periodicity, which we interpret as arising from corotating interaction regions (CIRs) in the wind. By analyzing the variability of the He II $\lambda$5411 emission line, the previously identified period was refined to P = 2.255 $\pm$ 0.008 (s.d.) days. The coherency time of the variability, which we associate with the lifetime of the CIRs in the wind, was deduced to be 40 $\pm$ 6 days, or $\sim$ 18 cycles, by cross-correlating the variability patterns as a function of time. When comparing the phased observational grayscale difference images with theoretical grayscales previously calculated from models including CIRs in an optically thin stellar wind, we find that two CIRs were likely present. A separation in longitude of $\Delta \phi \simeq$ 90$^{\circ}$ was determined between the two CIRs and we suggest that the different maximum velocities that they reach indicate that they emerge from different latitudes. We have also been able to detect observational signatures of the CIRs in other spectral lines (C IV $\lambda\lambda$5802,5812 and He I $\lambda$5876). Furthermore, a DAC was found to be present simultaneously with the CIR signatures detected in the He I $\lambda$5876 emission line which is consistent with the proposed geometry of the large-scale structures in the wind. Small-scale structures also show a presence in the wind, simultaneously with the larger scale structures, showing that they do in fact co-exist.

Context. The eclipsing binary GU Mon is located in the star-forming cluster Dolidze 25, which has the lowest metallicity measured in a Milky Way young cluster. Aims. GU Mon has been identified as a short-period eclipsing binary with two early B-type components. We set out to derive its orbital and stellar parameters. Methods. We present a comprehensive analysis, including B and V light curves and 11 high-resolution spectra, to verify the orbital period and determine parameters. We use the stellar atmosphere code FASTWIND to obtain stellar parameters and create templates for cross-correlation. We obtain a model to fit the light and radial-velocity curves using the Wilson-Devinney code iteratively and simultaneously. Results. The two components of GU Mon are identical stars of spectral type B1 V, with the same mass and temperature. The lightcurves are typical of an EW-type binary. The spectroscopic and photometric analyses agree on a period of 0.896640 +- 0.000007 d. We determine a mass of 9.0 +- 0.6 Msun for each component and temperatures of 28 000 +- 2 000 K. Both values are consistent with the spectral type. The two stars are overfilling their respective Roche lobes, sharing a common envelope, and therefore the orbit is synchronised and circularised. Conclusions. The GU Mon system has a fill-out factor above 0.8, containing two dwarf B-type stars on the main sequence. The two stars are in a very advanced stage of interaction, with their extreme physical similarity likely due to the common envelope. The expected evolution of such a system will very probably lead to a merger while still on the main sequence.

This is the third installment of GOSSS, a massive spectroscopic survey of Galactic O stars, based on new homogeneous, high signal-to-noise ratio, R~2500 digital observations selected from the Galactic O-Star Catalog (GOSC). In this paper we present 142 additional stellar systems with O stars from both hemispheres, bringing the total of O-type systems published within the project to 590. Among the new objects there are 20 new O stars. We also identify 11 new double-lined spectroscopic binaries (SB2s), of which 6 are of O+O type and 5 of O+B type, and an additional new tripled-lined spectroscopic binary (SB3) of O+O+B type. We also revise some of the previous GOSSS classifications, present some egregious examples of stars erroneously classified as O-type in the past, introduce the use of luminosity class IV at spectral types O4-O5.5, and adapt the classification scheme to the work of Arias et al. (2016).

Our aim is to find and classify OB stars in Sextans A, to later determine accurate stellar parameters of these blue massive stars in this low metallicity region $(Z \sim 0.1 \rm Z_{\odot})$. Using UBV photometry, the reddening-free index Q and GALEX imaging, we built a list of blue massive star candidates in Sextans A. We obtained low resolution (R $\sim$ 1000) GTC-OSIRIS spectra for a fraction of them and carried out spectral classification. For the confirmed O-stars we derive preliminary stellar parameters. The target selection criteria and observations were successful and have produced the first spectroscopic atlas of OB-type stars in Sextans A. From the whole sample of 18 observed stars, 12 were classified as early OB-types, including 5 O-stars. The radial velocities of all target stars are in agreement with their Sextans A membership, although three of them show significant deviations. We determined the stellar parameters of the O-type stars using the stellar atmosphere code FASTWIND, and revisited the sub-SMC temperature scale. Two of the O-stars are consistent with relatively strong winds and enhanced helium abundances, although results are not conclusive. We discuss the position of the OB stars in the HRD. Initial stellar masses run from slightly below 20 up to 40 solar masses. The target selection method worked well for Sextans A, confirming the procedure developed in Garcia \& Herrero (2013). The stellar temperatures are consistent with findings in other galaxies. Some of the targets deserve follow-up spectroscopy because of indications of a runaway nature, an enhanced helium abundance or a relatively strong wind. We observe a correlation between HI and OB associations similar to the irregular galaxy IC1613, confirming the previous result that the most recent star formation of Sextans A is currently on-going near the rim of the H\,{\sc I} cavity.

Context: Nearly one century after their discovery, the carrier(s) of Diffuse Interstellar Bands is/are still unknown and there are few sightlines studied in detail for a large number of DIBs. Aims: We want to study the ISM sightlines towards LS III +46 11 and LS III +46 12, two early-O-type stellar systems, and LS III +46 11 B, a mid-B-type star. The three targets are located in the stellar cluster Berkeley 90 and have a high extinction. Methods: We use the multi-epoch high-S/N optical spectra presented in paper I (Ma\'iz Apell\'aniz et al. 2015), the extinction results derived there, and additional spectra. Results: We have measured equivalent widths, velocities, and FWHMs for a large number of absorption lines in the rich ISM spectrum in front of Berkeley 90. The absorbing ISM has at least two clouds at different velocities, one with a lower column density (thinner) in the K I lines located away from Berkeley 90 and another one with a higher column density (thicker) associated with the cluster. The first cloud has similar properties for both O-star sightlines but the second one is thicker for LS III +46 11. The comparison between species indicate that the cloud with a higher column density has a denser core, allowing us to classify the DIBs in a sigma-zeta scale, some of them for the first time. The LS III +46 12 sightline also has a high-velocity redshifted component.

The young open cluster Dolidze 25, in the direction of the Galactic Anticentre, has been attributed a very low metallicity, with typical abundances between $-0.5$ and $-0.7$ dex below solar. We intend to derive accurate cluster parameters and accurate stellar abundances for some of its members. We have obtained a large sample of intermediate- and high-resolution spectra for stars in and around Dolidze 25. We used the FASTWIND code to generate stellar atmosphere models to fit the observed spectra. We derive stellar parameters for a large number of OB stars in the area, and abundances of oxygen and silicon for a number of stars with spectral types around B0. We measure low abundances in stars of Dolidze 25. For the three stars with spectral types around B0, we find $0.3$ dex (Si) and $0.5$ dex (O) below the values typical in the solar neighbourhood. These values, even though not as low as those given previously, confirm Dolidze 25 and the surrounding H II region Sh2-284 as the most metal-poor star-forming environment known in the Milky Way. We derive a distance $4.5\pm0.3\:$kpc to the cluster ($r_{\textrm{G}}\approx12.3\:$kpc). The cluster cannot be older than $\sim3\:$Myr, and likely is not much younger. One star in its immediate vicinity, sharing the same distance, has Si and O abundances at most $0.15\:$dex below solar. The low abundances measured in Dolidze 25 are compatible with currently accepted values for the slope of the Galactic metallicity gradient, if we take into account that variations of at least $\pm0.15\:$dex are observed at a given radius. The area traditionally identified as Dolidze 25 is only a small part of a much larger star-forming region that comprises the whole dust shell associated with Sh2-284 and very likely several other smaller H II regions in its vicinity.

The number of magnetic stars detected among massive stars is small; nevertheless, the role played by the magnetic field in stellar evolution cannot be disregarded. Links between line profile variability, enhancements/depletions of surface chemical abundances, and magnetic fields have been identified for low-mass B-stars, but for the O-type domain this is almost unexplored. Based on FORS2 and HARPS spectropolarimetric data, we present the first detection of a magnetic field in HD54879, a single slowly rotating O9.7 V star. Using two independent and different techniques we obtained the firm detection of a surface average longitudinal magnetic field with a maximum amplitude of about 600 G, in modulus. A quantitative spectroscopic analysis of the star with the stellar atmosphere code FASTWIND results in an effective temperature and a surface gravity of 33000$\pm1000$ K and 4.0$\pm0.1$ dex. The abundances of carbon, nitrogen, oxygen, silicon, and magnesium are found to be slightly lower than solar, but compatible within the errors. We investigate line-profile variability in HD54879 by complementing our spectra with spectroscopic data from other recent OB-star surveys. The photospheric lines remain constant in shape between 2009 and 2014, although H$\alpha$ shows a variable emission. The H$\alpha$ emission is too strong for a standard O9.7 V and is probably linked to the magnetic field and the presence of circumstellar material. Its normal chemical composition and the absence of photospheric line profile variations make HD54879 the most strongly magnetic, non-variable single O-star detected to date.

We investigate the multiplicity properties of 408 B-type stars observed in the 30 Doradus region of the Large Magellanic Cloud with multi-epoch spectroscopy from the VLT-FLAMES Tarantula Survey (VFTS). We use a cross-correlation method to estimate relative radial velocities from the helium and metal absorption lines for each of our targets. Objects with significant radial-velocity variations (and with an amplitude larger than 16 km/s) are classified as spectroscopic binaries. We find an observed spectroscopic binary fraction (defined by periods of <10^3.5 d and mass ratios >0.1) for the B-type stars, f_B(obs) = 0.25 +/- 0.02, which appears constant across the field of view, except for the two older clusters (Hodge 301 and SL 639). These two clusters have significantly lower fractions of 0.08 +/- 0.08 and 0.10 +/- 0.09, respectively. Using synthetic populations and a model of our observed epochs and their potential biases, we constrain the intrinsic multiplicity properties of the dwarf and giant (i.e. relatively unevolved) B-type stars in 30 Dor. We obtain a present-day binary fraction f_B(true) = 0.58 +/- 0.11, with a flat period distribution. Within the uncertainties, the multiplicity properties of the B-type stars agree with those for the O stars in 30 Dor from the VFTS.

Context: It appears that most (if not all) massive stars are born in multiple systems. At the same time, the most massive binaries are hard to find due to their low numbers throughout the Galaxy and the implied large distances and extinctions. AIMS: We want to study: [a] LS III +46 11, identified in this paper as a very massive binary; [b] another nearby massive system, LS III +46 12; and [c] the surrounding stellar cluster, Berkeley 90. Methods: Most of the data used in this paper are multi-epoch high-S/N optical spectra though we also use Lucky Imaging and archival photometry. The spectra are reduced with devoted pipelines and processed with our own software, such as a spectroscopic-orbit code, CHORIZOS, and MGB. Results: LS III +46 11 is identified as a new very-early-O-type spectroscopic binary [O3.5 If* + O3.5 If*] and LS III +46 12 as another early O-type system [O4.5 V((f))]. We measure a 97.2-day period for LS III +46 12 and derive minimum masses of 38.80$\pm$0.83 M_Sol and 35.60$\pm$0.77 M_Sol for its two stars. We measure the extinction to both stars, estimate the distance, search for optical companions, and study the surrounding cluster. In doing so, a variable extinction is found as well as discrepant results for the distance. We discuss possible explanations and suggest that LS III +46 12 may be a hidden binary system, where the companion is currently undetected.

The IACOB project is an ambitious long-term project which is contributing to step forward in our knowledge about the physical properties and evolution of Galactic massive stars. The project aims at building a large database of high-resolution, multi-epoch, spectra of Galactic OB stars, and the scientific exploitation of the database using state-of-the-art models and techniques. In this proceeding, we summarize the latest updates of the IACOB spectroscopic database and highlight some of the first scientific results from the IACOB project; we also announce the first data release and the first public version of the iacob-broad tool for the line-broadening characterization of OB-type spectra.

We present spectral classifications for 438 B-type stars observed as part of the VLT-FLAMES Tarantula Survey (VFTS) in the 30 Doradus region of the Large Magellanic Cloud. Radial velocities are provided for 307 apparently single stars, and for 99 targets with radial-velocity variations which are consistent with them being spectroscopic binaries. We investigate the spatial distribution of the radial velocities across the 30 Dor region, and use the results to identify candidate runaway stars. Excluding potential runaways and members of two older clusters in the survey region (SL 639 and Hodge 301), we determine a systemic velocity for 30 Dor of 271.6 +/- 12.2 km/s from 273 presumed single stars. Employing a 3-sigma criterion we identify nine candidate runaway stars (2.9% of the single stars with radial-velocity estimates). The projected rotational velocities of the candidate runaways appear to be significantly different to those of the full B-type sample, with a strong preference for either large (>345 km/s) or small (<65 km/s) rotational velocities. Of the candidate runaways, VFTS 358 (classified B0.5: V) has the largest differential radial velocity (-106.9 +/- 16.2 km/s), and a preliminary atmospheric analysis finds a significantly enriched nitrogen abundance of 12+log(N/H) > ~8.5. Combined with a large rotational velocity (vsini = 345 +/- 22 km/s), this is suggestive of past binary interaction for this star.

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

CAF\'E-BEANS is an on-going survey running on the 2.2 m telescope at Calar Alto. For more than two years, CAF\'E-BEANS has been collecting high-resolution spectra of early-type stars with the aim of detecting and characterising spectroscopic binaries. The main goal of this project is a thorough characterisation of multiplicity in high-mass stars by detecting all spectroscopic and visual binaries in a large sample of Galactic O-type stars, and solving their orbits. Our final objective is eliminating all biases in the high-mass-star IMF created by undetected binaries.

In this poster we present three developments related to the Galactic O-Star Spectroscopic Survey (GOSSS). First, we are making public the first version of MGB, an IDL code that allows the user to compare oberved spectra to a grid of spectroscopic standards to measure spectral types, luminosity classes, rotation indexes, and spectral qualifiers. Second, we present the associated grid of standard stars for the spectral types O2 to O9.7, with several improvements over the original GOSSS grid of Sota et al. (2011). Third, we present a list of egregious classification errors in SIMBAD: stars that are or have been listed there as being of O type but that in reality are late-type stars.

Context. The early-type binary MY Cam belongs to the young open cluster Alicante 1, embedded in Cam OB3. Aims. MY Cam consists of two early-O type main-sequence stars and shows a photometric modulation suggesting an orbital period slightly above one day. We intend to confirm this orbital period and derive orbital and stellar parameters. Methods. Timing analysis of a very exhaustive (4607 points) light curve indicates a period of 1.1754514 +- 0.0000015 d. High- resolution spectra and the cross-correlation technique implemented in the TODCOR program were used to derive radial velocities and obtain the corresponding radial velocity curves for MY Cam. Modelling with the stellar atmosphere code FASTWIND was used to obtain stellar parameters and create templates for cross-correlation. Stellar and orbital parameters were derived using the Wilson-Devinney code, such that a complete solution to the binary system could be described. Results. The determined masses of the primary and secondary stars in MY Cam are 37.7 +- 1.6 and 31.6 +- 1.4 Msol, respectively. The corresponding temperatures, derived from the model atmosphere fit, are 42 000 and 39 000 K, with the more massive component being hotter. Both stars are overfilling their Roche lobes, sharing a common envelope. Conclusions. MY Cam contains the most massive dwarf O-type stars found so far in an eclipsing binary. Both components are still on the main sequence, and probably not far from the zero-age main sequence. The system is a likely merger progenitor, owing to its very short period.

The distribution of stars in the Hertzsprung-Russell diagram narrates their evolutionary history and directly assesses their properties. Placing stars in this diagram however requires the knowledge of their distances and interstellar extinctions, which are often poorly known for Galactic stars. The spectroscopic Hertzsprung-Russell diagram (sHRD) tells similar evolutionary tales, but is independent of distance and extinction measurements. Based on spectroscopically derived effective temperatures and gravities of almost 600 stars, we derive for the first time the observational distribution of Galactic massive stars in the sHRD. While biases and statistical limitations in the data prevent detailed quantitative conclusions at this time, we see several clear qualitative trends. By comparing the observational sHRD with different state-of-the-art stellar evolutionary predictions, we conclude that convective core overshooting may be mass-dependent and, at high mass ($\geq 15\,M_\odot$), stronger than previously thought. Furthermore, we find evidence for an empirical upper limit in the sHRD for stars with $T_{\rm{eff}}$ between 10000 and 32000 K and, a strikingly large number of objects below this line. This over-density may be due to inflation expanding envelopes in massive main-sequence stars near the Eddington limit.

(Abridged) Main sequence massive stars embedded in an HII region should have the same chemical abundances as the surrounding nebular gas+dust. The Cocoon nebula, a close-by Galactic HII region ionized by a narrow line B0.5 V single star (BD+46 3474), is an ideal target to perform a detailed comparison of nebular and stellar abundances in the same Galactic HII region. We investigate the chemical content of O, N and S in the Cocoon nebula from two different points of view: an empirical analysis of the nebular spectrum and a detailed spectroscopic analysis of its ionizing B-type star using state-of-the-art stellar atmosphere modeling. By comparing the stellar and nebular abundances, we aim to indirectly address the long-standing problem of the discrepancy found between abundances obtained from collisionally excited lines (CELs) and optical recombination lines in photoionized nebulae. We collect spatially resolved spectroscopy of the Cocoon nebula and a high resolution optical spectrum of its ionizing star. Standard nebular techniques are used to compute the physical conditions and gaseous abundances of O, N and S. We perform a self-consistent spectroscopic abundance analysis of BD+46 3474 based on the atmosphere code FASTWIND to determine the stellar parameters and Si, O, and N abundances. The Cocoon nebula and its ionizing star, located at a distance of 800+-80 pc, have a very similar chemical composition as the Orion nebula and other B-type stars in the solar vicinity. This result agrees with the high degree of homogeneity of the present-day composition of the solar neighbourhood as derived from the study of the local cold-gas ISM. The comparison of stellar and nebular CELs abundances in the Cocoon nebula indicates that O and N gas+dust nebular values are in better agreement with stellar ones assuming small temperature fluctuations, of the order of those found in the Orion nebula.

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.

Imagine we could do asteroseismology of large samples of OB-type stars by using just one spectrum per target. That would be great! But this is probably a crazy and stupid idea. Or maybe not. Maybe we have the possibility to open a new window to investigate stellar oscillations in massive stars that has been in front of us for many years, but has not attracted very much our attention: the characterization and understanding of the so-called macroturbulent broadening in OB-type stars.

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.

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

The blue supergiant Sher25 is surrounded by an asymmetric, hourglass-shaped circumstellar nebula, which shows similarities to the triple-ring structure seen around SN1987A. From optical spectroscopy over six consecutive nights, we detect periodic radial velocity variations in the stellar spectrum of Sher25 with a peak-to-peak amplitude of ~12 km/s on a timescale of about 6 days, confirming the tentative detec-tion of similar variations by Hendry et al. From consideration of the amplitude and timescale of the signal, coupled with observed line profile variations, we propose that the physical origin of these variations is related to pulsations in the stellar atmosphere, rejecting the previous hypothesis of a massive, short-period binary companion. The radial velocities of two other blue supergiants with similar bipolar nebulae, SBW1 and HD 168625, were also monitored over the course of six nights, but these did not display any significant radial velocity variations.

BD+6073 is the optical counterpart of the X-ray source IGR J00370+6122, a probable accretion-powered X-ray pulsar. The X-ray light curve of this binary system shows clear periodicity at 15.7 d, which has been interpreted as repeated outbursts around the periastron of an eccentric orbit. We obtained high-resolution spectra of BD+6073 at different epochs. We used the FASTWind code to generate a stellar atmosphere model to fit the observed spectrum and obtain physical magnitudes. The synthetic spectrum was used as a template for cross-correlation with the observed spectra to measure radial velocities. The radial velocity curve provided an orbital solution for the system. We have also analysed the RXTE/ASM and Swift/BAT light curves to confirm the stability of the periodicity. BD +6073 is a BN0.7 Ib low-luminosity supergiant located at an approximate distance of 3.1 kpc, in the CasOB4 association. We derive Teff=24000 K and log gc=3.0, and chemical abundances consistent with a moderately high level of evolution. The spectroscopic and evolutionary masses are consistent at the 1 sigma level with a mass of 15 solar masses. The recurrence time of the X-ray flares is the orbital period of the system. The NS is in a high eccentricity (e=0.56) orbit, and the X-ray emission is strongly peaked around orbital phase 0.2, though the observations are consistent with some level of X-ray activity happening at all orbital phases. The X-ray behaviour of IGR J00370+6122 is reminiscent of intermediate SFXTs, though its peak luminosity is rather low. The orbit is somewhat wider than those of classical persistent supergiant X-ray binaries, which, combined with the low luminosity of the mass donor, explains the low X-ray luminosity. IGR J00370+6122 will likely evolve towards a persistent supergiant system, highlighting the evolutionary connection between different classes of wind-accreting X-ray sources.

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

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.

Stellar rotation is an important parameter in the evolution of massive stars. Accurate and reliable measurements of projected rotational velocities in large samples of OB stars are crucial to confront the predictions of stellar evolutionary models with observational constraints. We reassess previous determinations of projected rotational velocities (vsini) in Galactic OB stars using a large, high quality spectroscopic dataset, and a strategy which account for other sources of broadening appart from rotation affecting the diagnostic lines We present a versatile and user friendly IDL tool, based on a combined Fourier Transform (FT) + goodness of fit (GOF) methodology, for the line-broadening characterization in OB-type stars. We use this tool to (a) investigate the impact of macroturbulent and microturbulent broadenings on vsini measurements, and (b) determine vsini in a sample of 200 Galactic OB-type stars, also characterizing the amount of macroturbulent broadening (\vmacro) affecting the line profiles. We present observational evidence illustrating the strengths and limitations of the proposed FT+GOF methodology for the case of OB stars. We confirm previous statements (based on indirect arguments or smaller samples) that the macroturbulent broadening is ubiquitous in the massive star domain. We compare the newly derived vsini with previous determinations not accounting for this extra line-broadening contribution, and show that those cases with vsini< 120 km/s need to be systematically revised downwards by ~25 (+/-20) km/s. We suggest that microturbulence may impose an upper limit below which vsini and \vmacro\ could be incorrectly derived by means of the proposed methodology as presently used, and discuss the implications of this statement on the study of relatively narrow line massive stars.

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.

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 are using five different surveys to compile the largest sample of diffuse interstellar band (DIB) measurements ever collected. GOSSS is obtaining intermediate-resolution blue-violet spectroscopy of ~2500 OB stars, of which 60% have already been observed and processed. The other four surveys have already collected multi-epoch high-resolution optical spectroscopy of 700 OB stars with different telescopes, including the 9 m Hobby-Eberly Telescope in McDonald Observatory. Some of our stars are highly-extinguished targets for which no good-quality optical spectra have ever been published. For all of the targets in our sample we have obtained accurate spectral types, measured non-DIB ISM lines, and compiled information from the literature to calculate the extinction. Here we present the first results of the project, the properties of twenty DIBs in the 4100-5500 {\AA} range. We clearly detect a couple of previously elusive DIBs at 4170 {\AA} and 4591 {\AA}; the latter could have coronene and ovalene cations as carriers.

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

This paper presents a first observational investigation of the faint Of?cp star NGC 1624-2, yielding important new constraints on its spectral and physical characteristics, rotation, magnetic field strength, X-ray emission and magnetospheric properties. Modeling the spectrum and spectral energy distribution, we conclude that NGC 1624-2 is a main sequence star of mass M {\simeq} 30 M{\odot}, and infer an effective temperature of 35 {\pm} 2 kK and log g = 4.0 {\pm} 0.2. Based on an extensive time series of optical spectral observations we report significant variability of a large number of spectral lines, and infer a unique period of 157.99 {\pm} 0.94 d which we interpret as the rotational period of the star. We report the detection of a very strong - 5.35 {\pm} 0.5 kG - longitudinal magnetic field <Bz>, coupled with probable Zeeman splitting of Stokes I profiles of metal lines confirming a surface field modulus <B> of 14 {\pm} 1 kG, consistent with a surface dipole of polar strength >~ 20 kG. This is the largest magnetic field ever detected in an O-type star, and the first report of Zeeman splitting of Stokes I profiles in such an object. We also report the detection of reversed Stokes V profiles associated with weak, high-excitation emission lines of O iii, which we propose may form in the close magnetosphere of the star. We analyze archival Chandra ACIS-I X-ray data, inferring a very hard spectrum with an X-ray efficiency log Lx/Lbol = -6.4, a factor of 4 larger than the canonical value for O-type stars and comparable to that of the young magnetic O-type star {\theta}1 Ori C and other Of?p stars. Finally, we examine the probable magnetospheric properties of the star, reporting in particular very strong magnetic confinement of the stellar wind, with {\eta}* {\simeq} 1.5 {\times} 10^4, and a very large Alfven radius, RAlf = 11.4 R*.

The quantitative study of the physical properties and chemical abundances of large samples of massive blue stars at different metallicities is a powerful tool to understand the nature and evolution of these objects. Their analysis beyond the Milky Way is challenging, nonetheless it is doable and the best way to investigate their behavior in different environments. Fulfilling this task in an objective way requires the implementation of automatic analysis techniques that can perform the analyses systematically, minimizing at the same time any possible bias. As part of the ARAUCARIA project we carry out the first quantitative spectroscopic analysis of a sample of 12 B-type supergiants in the galaxy NGC55 at 1.94 Mpc away. By applying the methodology developed in this work, we derive their stellar parameters, chemical abundances and provide a characterization of the present-day metallicity of their host galaxy. Based on the characteristics of the stellar atmosphere/line formation code FASTWIND, we designed and created a grid of models for the analysis of massive blue supergiant stars. Along with this new grid, we implemented a spectral analysis algorithm. Both tools were specially developed to perform fully consistent quantitative spectroscopic analyses of low spectral resolution of B-type supergiants in a fast and objective way. We present the main characteristics of our FASTWIND model grid and perform a number of tests to investigate the reliability of our methodology. The automatic tool is applied afterward to a sample of 12 B-type supergiant stars in NGC55, deriving the stellar parameters and abundances. The results indicate that our stars are part of a young population evolving towards a red supergiant phase. The derived chemical composition hints to an average metallicity similar to the one of the Large Magellanic Cloud, with no indication of a spatial trend across the galaxy.

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.

The IACOB spectroscopic survey of Galactic OB stars is an ambitious observational project aimed at compiling a large, homogeneous, high-resolution database of optical spectra of massive stars observable from the Northern hemisphere. The quantitative spectroscopic analysis of this database, complemented by the invaluable information provided by Gaia (mainly regarding photometry and distances), will represent a major step forward in our knowledge of the fundamental physical characteristics of Galactic massive stars. In addition, results from this analysis will be of interest for other scientific questions to be investigated using Gaia observations. In this contribution we outline the present status of the IACOB spectroscopic database and indicate briefly some of the synergy links between the IACOB and Gaia scientific projects.

There are several ongoing massive-star (mostly of spectral type O) surveys that are significantly increasing the quality and quantity of the spectroscopic information about these objects. Here we discuss and present results for two of them, GOSSS and NoMaDS. We also discuss recent and future developments on the Galactic O-Star Catalog and announce the upcoming availability of Marxist Ghost Buster, an IDL code that attacks spectral classification (hence the name) by using an interactive comparison with spectral libraries.

Context: The discovery of new massive star clusters and massive stellar populations in previously known clusters in our Galaxy by means of infrared studies has changed our view of the Milky Way from an inactive to an active star-forming machine. Within this scenario, we present a near-infrared spectrophotometric study of the stellar content of the compact H II region Sh2-152. Aims: We aim to determine the distance, extinction, age, and mass of Sh2-152, using for the first time near-infrared stellar classification for several sources in the region. Methods: Using our near-infrared (J, H, and Ks) photometry and the colour-magnitude diagram for the cluster field, we selected 13 bright stars, candidate members of the reddened cluster's main sequence, for H- and K- spectroscopy and spectral classification. This near-infrared information was complemented with an optical spectrum of the ionizing central star to confirm its spectral nature. Results: From the 13 spectroscopically observed stars, 5 were classified as B-type, 3 as G-type, 2 were young stellar objects (YSOs), and 3 remained unclassified (because of the poor data quality). The cluster's extinction varies from A_{Ks}=0.5 to 2.6 magnitudes (A_V=4.5 to 24 magnitudes) and the distance is estimated to be 3.21\pm0.21 kpc. The age of the cluster is younger than 9.4 Myr and the lower limit to the total mass of the cluster is (2.45\pm0.79)\cdot10^3 M_{Sun}. We compare the number of ionizing photons emitted from the OB-type stars with the Lyman continuum photons derived from the radio observations and conclude that both quantities are consistent for the central region of Sh2-152. In contrast, the main ionizing source of the lower region remains unidentified.

We report on the detection of a third massive star component in the \sigma\ Orionis AB system, traditionally considered as a binary system. The system has been monitored by the IACOB spectroscopic survey of Northern Massive stars program, obtaining 23 high-resolution FIES@NOT spectra with a time-span of ~2.5 years. The analysis of the radial velocity curves of the two spectroscopic components observed in the spectra has allowed us obtain the orbital parameters of the system, resulting in a high eccentric orbit (e~0.78) with an orbital period of 143.5 +/- 0.5 d. This result implies the actual presence of three stars in the \sigma\ Orionis AB system when combined with previous results obtained from the study of the astrometric orbit (with an estimated period of ~157 yr).

We have selected the Galactic HII region M43, a close-by apparently spherical nebula ionized by a single star (HD37061) to investigate several topics of recent interest in the field of HII regions and massive stars. In a series of two papers we perform a combined, comprehensive study of the nebula and its ionizing star by using as many observational constraints as possible. We collected for this study a set of high-quality observations, including the optical spectrum of HD37061, along with nebular optical imaging and long-slit spatially resolved spectroscopy. The first part of our study comprises a quantitative spectroscopic analysis of the ionizing star, and the empirical analysis of the nebular images and spectroscopy. We determine the stellar parameters of HD37061 and the total number of ionizing photons emitted by the star. We find observational evidence of the presence of scattered light from the Huygens region (brightest part of the Orion nebula) in the M43 region. We show the importance of an adequate correction of this scattered light from the imagery and spectroscopic observations of M43 for a proper determination of the total nebular H_alpha luminosity, the nebular physical conditions and chemical abundances. We perform a detailed nebular empirical analysis of 9 apertures extracted from a long-slit located to the west of HD37061, obtaining the spatial distribution of the physical conditions and ionic abundances. For three of the analyzed elements (O, S, and N) we could determine total abundances directly from observable ions (no ionization correction factors were needed). The comparison of these abundances with those derived from the spectrum of the Orion nebula indicates the importance of the atomic data and, specially in the case of M\,42, the considered ionization correction factors.

We re-examine the recombination/collisional emission line (RL/CEL) nebular abundance discrepancy problem in the light of recent high-quality abundance determinations in young stars in the Orion star-forming region. We re-evaluate the CEL and RL abundances of several elements in the Orion nebula and estimate the associated uncertainties, taking into account the uncertainties in the ionization correction factors for unseen ions. We estimate the amount of oxygen trapped in dust grains for several scenarios of dust formation. We compare the resulting gas+dust nebular abundances with the stellar abundances of a sample of 13 B-type stars from the Orion star-forming region (Ori\,OB1), analyzed in Papers I and III of this series. We find that the oxygen nebular abundance based on recombination lines agrees much better with the stellar abundances than the one derived from the collisionally excited lines. This result calls for further investigation. If the CEL/RL abundance discrepancy were caused by temperature fluctuations in the nebula, as argued by some authors, the same kind of discrepancy should be seen for the other elements, such as C, N and Ne, which is not what we find in the present study. Another problem is that with the RL abundances, the energy balance of the Orion nebula is not well understood. We make some suggestions concerning the next steps to undertake to solve this problem.

We present the IACOB spectroscopic database, an homogeneous set of high quality, high resolution spectra of Galactic O- and B-type stars obtained with the FIES spectrograph attached to the Nordic Optical Telescope. We also present some results from ongoing projects using the IACOB database.

We present a summary of main results from the studies performed in the series of papers "The chemical composition of the Orion star forming region". We reinvestigate the chemical composition of B-type stars in the Orion OB1 association by means of state-of-the-art stellar atmosphere codes, atomic models and techniques, and compare the resulting abundances with those obtained from the emission line spectra of the Orion nebula (M42), and recent determinations of the Solar chemical composition.

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.

We present the IACOB spectroscopic database, the largest homogeneous database of high-resolution, high signal-to-noise ratio spectra of Northern Galactic OB-type stars compiled up to date. The spectra were obtained with the FIES spectrograph attached to the Nordic Optical Telescope. We briefly summarize the main characeristics and present status of the IACOB, first scientific results, and some future plans for its extension and scientific exploitation.

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.

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.

Recent accurate abundance analyses of B-type main sequence stars in the solar vicinity has shown that abundances derived from these stellar objects are more homogeneous and metal-rich than previously thought. We investigate whether the inhomogeneity of abundances previously found in B-type stars in the Ori OB1 association is real (hence a signature of enrichment of the newly formed stars in an induced star formation scenario) or a consequence of intrinsic errors induced by the use of photometric indices to establish the stellar parameters prior to the abundance analysis. We obtained a new (improved) spectroscopic data set comprising 13 B-type stars in the various Ori OB1 associations, and performed a detailed, self-consistent spectroscopic abundance analysis by means of the modern stellar atmosphere code FASTWIND. We detect systematic errors in the stellar parameters determined previously which affect the derived abundances. Once these errors are accounted for, we find a high degree of homogeneity in the O and Si abundances for stars in the four Ori OB1 subgroups. The derived abundances are in very good agreement with recent determinations in other B-type stars in the solar vicinity. We also compare our results with those obtained for the Sun during the epoch of the "solar crisis", and the Orion nebula.

Several works have recently shown that there is an important extra line-broadening (usually called $macroturbulence$) affecting the spectra of O and B Supergiants that adds to stellar rotation. So far, the only (very recent) physical explanation for the appearance of $macroturbulence$ relates to oscillations. This is a plausible explanation, but no direct evidence confirming its validity has been presented yet. We recently started an observational project to obtain constraints on the time-scales of variability associated to this extra line-broadening and its possible origin. Our observational strategy consists of the study of a well selected group of O and B stars, for which we obtain time series of high-quality spectra. We present some preliminary results from our first campaign with FIES@NOT2.5m.

We present a detailed comparison of the ionizing spectral energy distributions (SEDs) predicted by four modern stellar atmosphere codes, TLUSTY, CMFGEN, WMbasic, and FASTWIND. We consider three sets of stellar parameters representing a late O-type dwarf (O9.5 V), a mid O-type (O7 V) dwarf, and an early O-type dwarf (O5.5 V). We explore two different possibilities for such a comparison, following what we called evolutionary and observational approaches: in the evolutionary approach one compares the SEDs of stars defined by the same values of Teff and logg; in the observational approach the models to be compared do not necessarily have the same Teff and logg, but produce similar H and HeI-II optical lines. We find that there is a better agreement, in terms of Q(H0), the ratio Q(He0)/Q(H0), and the shape of the SEDs predicted by the four codes in the spectral range between 13 and 30 eV, when models are compared following the observational approach. However, even in this case, large differences are found at higher energies. We then discuss how the differences in the SEDs may affect the overall properties of surrounding nebulae in terms of temperature and ionization structure. We find that the effect over the nebular temperature is not larger than 300-350 K. Contrarily, the different SEDs produce significantly different nebular ionization structures. This will lead to important consequences on the establishment of the ionization correction factors that are used in the abundance determination of HII regions, as well as in the characterization of the ionizing stellar population from nebular line ratios.

We present some results of an on-going project aimed at studying a sample of Galactic HII regions ionized by a single massive star to test the predictions of modern generation stellar atmosphere codes in the H Lyman continuum. The observations collected for this study comprise the optical spectra of the corresponding ionizing stars, along with imaging and long-slit spatially resolved nebular observations. The analysis of the stellar spectra allows to obtain the stellar parameters of the ionizing star, while the nebular observations provide constraints on the nebular abundances and gas distribution. All this information is then used to construct tailored photoionization models of the HII regions. The reliability of the stellar ionizing fluxes is hence tested by comparing the photoionization model results with the observations in terms of the spatial variation across the nebula of an appropriate set of nebular line ratios.

We present a comprehensive study that applies the Fourier transform to a sample of O and early B-type stars (either dwarfs, giants, or supergiants) to determine their projected rotational velocities, compare with previous values obtained with other methods, and seek for evidence of extra broadening in the spectral lines The Fourier technique, extensively used in the study of cooler stars, has remained only marginally applied for the case of early-type stars. The comparison of \vsini values obtained through the \ft and \fwhm methods shows that the \fwhm technique must be used with care in the analysis of OB giants and supergiants, and when it is applied to \ion{He}{i} lines. Contrarily, the \ft method appears to be a powerful tool to derive reliable projected rotational velocities, and separate the effect of rotation from other broadening mechanisms present in these stars. The analysis of the sample of OB stars shows that while dwarfs and giants display a broad range of projected rotational velocities, from less than 30 up to 450 \kms, supergiants have in general values close to or below 100 \kms. The analysis has also definitely shown that while the effect of extra broadening is negligible in OB dwarfs, it is clearly present in supergiants. When examining the behavior of the projected rotational velocities with the stellar parameters and across the HR diagram, we conclude, in agreement with previous researchers, that the rotational velocity should decrease when the stars evolve. On the contrary, macroturbulence may be constant, resulting therefore in an increasing importance as compared to rotation when the stars evolve.

We present the first results of a project aimed at the combined study of massive stars and their surrounding nebulae by means of a detailed study of Galactic HII regions ionized by only one massive star. With this, we intend to check the validity of the new generation of massive star model atmosphere codes in terms of ionizing flux distribution. We take into account the effect of the nebular density distribution in our analyses. Various types of stellar and nebular observations have been collected for this purpose.