The mass, origin and evolutionary stage of the binary system LB-1 has been the subject of intense debate, following the claim that it hosts an $\sim$70$M_{\odot}$ black hole, in stark contrast with the expectations for stellar remnants in the Milky Way. We conducted a high-resolution, phase-resolved spectroscopic study of the near-infrared Paschen lines in this system, using the 3.5-m telescope at Calar Alto Observatory. We find that Pa$\beta$ and Pa$\gamma$ (after proper subtraction of the stellar absorption component) are well fitted with a standard double-peaked model, typical of disk emission. We measured the velocity shifts of the red and blue peaks at 28 orbital phases: the line center has an orbital motion in perfect antiphase with the stellar motion, and the radial velocity amplitude ranges from 8 to 13 km/s for different choices of lines and profile modelling. We interpret this curve as proof that the disk is tracing the orbital motion of the primary, ruling out the circumbinary disk and the hierarchical triple scenarios. The phase-averaged peak-to-peak half-separation (proxy for the projected rotational velocity of the outer disk) is $\sim$70 km s$^{-1}$, larger than the stellar orbital velocity and also inconsistent with a circumbinary disk. From those results, we infer a primary mass 4--8 times higher than the secondary mass. Moreover, we show that the ratio of the blue and red peaks (V/R intensity ratio) has a sinusoidal behaviour in phase with the secondary star, which can be interpreted as the effect of external irradiation by the secondary star on the outer disk. Finally, we briefly discuss our findings in the context of alternative scenarios recently proposed for LB-1. Definitive tests between alternative solutions will require further astrometric data from $Gaia$.

From an analysis of absorption lines in the high resolution spectra we have derived the radial velocities, stellar parameters (Teff, gravity, wind-strength parameter logQ and projected rotational velocity) and abundances (C, N, O, and Si ) of IRAS 17460-3114, IRAS 18131-3008, IRAS 19336-0400, LSE 45 and LSE 163. Abundances are found to be solar, except for a low Si abundance in IRAS 19336-0400 and a mild CNO pattern in LSE 163, that rotates at an unusual large rotational velocity for its spectral classification. Combining the stellar parameters information with Gaia DR2 data we are able to derive absolute magnitudes, radii and luminosities and clarify the possible post-AGB nature of the objects. IRAS 17460-3114 and IRAS 18131-3008 are found to be massive OB stars, whereas IRAS 19336-0400 is found to be a post-AGB star, already showing nebular lines in the spectrum. However, we could not confirm the nature of LSE 45 and LSE 163 as post-AGB stars, although their parameters are much more inconsistent with those of massive stars. In both cases, we find a discrepancy between the spectroscopic mass and that derived from the predictions of post-AGB evolutionary tracks. In addition, LSE 45 lacks nebular lines, that are present in IRAS 19336-0400 at a similar temperature. In the case of LSE 163 the rotational velocity (259+/-15 km/s) would be extremely large for a star evolving to CSPN. The combination of this rotational velocity, the high Galactic latitude, slightly large radial velocity and mild CNO enhancement suggests a history of binary interaction.

Uncertainties in stellar structure and evolution theory are largest for stars undergoing core convection on the main sequence. A powerful way to calibrate the free parameters used in the theory of stellar interiors is asteroseismology, which provides direct measurements of angular momentum and element transport. We report the detection and classification of new variable O and B stars using high-precision short-cadence (2-min) photometric observations assembled by the Transiting Exoplanet Survey Satellite (TESS). In our sample of 154 O and B stars, we detect a high percentage (90%) of variability. Among these we find 23 multiperiodic pulsators, 6 eclipsing binaries, 21 rotational variables, and 25 stars with stochastic low-frequency variability. Several additional variables overlap between these categories. Our study of O and B stars not only demonstrates the high data quality achieved by TESS for optimal studies of the variability of the most massive stars in the Universe, but also represents the first step towards the selection and composition of a large sample of O and B pulsators with high potential for joint asteroseismic and spectroscopic modeling of their interior structure with unprecedented precision.

We present an 80-d long uninterrupted high-cadence K2 light curve of the B1Iab supergiant rho Leo (HD 91316), deduced with the method of halo photometry. This light curve reveals a dominant frequency of $f_{\rmrot}=0.0373$d$^{-1}$ and its harmonics. This dominant frequency corresponds with a rotation period of 26.8d and is subject to amplitude and phase modulation. The K2 photometry additionally reveals multiperiodic low-frequency variability ($<1.5 $d$^{-1}$) and is in full agreement with low-cadence high-resolution spectroscopy assembled during 1800 days. The spectroscopy reveals rotational modulation by a dynamic aspherical wind with an amplitude of about 20km s$^{-1}$ in the H$\alpha$ line, as well as photospheric velocity variations of a few km s$^{-1}$ at frequencies in the range 0.2 to 0.6 d$^{-1}$ in the SiIII 4567\AA\ line. Given the large macroturbulence needed to explain the spectral line broadening of the star, we interpret the detected photospheric velocity as due to travelling super-inertial low-degree large-scale gravity waves with dominant tangential amplitudes and discuss why $\rho$~Leo is an excellent target to study how the observed photospheric variability propagates into the wind.

The O6 Vn star HD 64315 is believed to belong to the star-forming region known as NGC 2467, but previous distance estimates do not support this association. We explore the multiple nature of this star with the aim of determining its distance, and understanding its connection to NGC 2467. A total of 52 high-resolution spectra have been gathered over a decade. We use their analysis, in combination with the photometric data from All Sky Automated Survey and Hipparcos catalogues, to conclude that HD 64315 is composed of at least two spectroscopic binaries, one of which is an eclipsing binary. HD 64315 contains two binary systems, one of which is an eclipsing binary. The two binaries are separated by 0.09 arcsec (or 500 AU) if the most likely distance to the system, around 5 kpc, is considered. The presence of fainter companions is not excluded by current observations. The non-eclipsing binary (HD 64315 AaAb) has a period of 2.70962901+/-0.00000021 d. Its components are hotter than those of the eclipsing binary, and dominate the appearance of the system. The eclipsing binary (HD 64315 BaBb) has a shorter period of 1.0189569+/-0.0000008 d. We derive masses of 14.6+-2.3 M$_\odot$ for both components of the BaBb system. They are almost identical; both stars are overfilling their respective Roche lobes, and share a common envelope in an overcontact configuration. The non-eclipsing binary is a detached system composed of two stars with spectral types around O6 V with minimum masses of 10.8 M$_\odot$ and 10.2 M$_\odot$, and likely masses aprox. 30 M$_\odot$. HD 64315 provides a cautionary tale about high-mass star isolation and multiplicity. Its total mass is likely above 90 M$_\odot$,but it seems to have formed without an accompanying cluster. It contains one the most massive overcontact binaries known, a likely merger progenitor in a very wide multiple system.

Stellar evolution models are most uncertain for evolved massive stars. Asteroseismology based on high-precision uninterrupted space photometry has become a new way to test the outcome of stellar evolution theory and was recently applied to a multitude of stars, but not yet to massive evolved supergiants.Our aim is to detect, analyse and interpret the photospheric and wind variability of the O9.5Iab star HD 188209 from Kepler space photometry and long-term high-resolution spectroscopy. We used Kepler scattered-light photometry obtained by the nominal mission during 1460d to deduce the photometric variability of this O-type supergiant. In addition, we assembled and analysed high-resolution high signal-to-noise spectroscopy taken with four spectrographs during some 1800d to interpret the temporal spectroscopic variability of the star. The variability of this blue supergiant derived from the scattered-light space photometry is in full in agreement with the one found in the ground-based spectroscopy. We find significant low-frequency variability that is consistently detected in all spectral lines of HD 188209. The photospheric variability propagates into the wind, where it has similar frequencies but slightly higher amplitudes. The morphology of the frequency spectra derived from the long-term photometry and spectroscopy points towards a spectrum of travelling waves with frequency values in the range expected for an evolved O-type star. Convectively-driven internal gravity waves excited in the stellar interior offer the most plausible explanation of the detected variability.

Massive star evolution remains only partly constrained. In particular, the exact role of rotation has been questioned by puzzling properties of OB stars in the Magellanic Clouds. Our goal is to study the relation between surface chemical composition and rotational velocity, and to test predictions of evolutionary models including rotation. We have performed a spectroscopic analysis of a sample of fifteen Galactic O7-8 giant stars. This sample is homogeneous in terms of mass, metallicity and evolutionary state. It is made of stars with a wide range of projected rotational velocities. We show that the sample stars are located on the second half of the main sequence, in a relatively narrow mass range (25-40 Msun). Almost all stars with projected rotational velocities above 100 km/s have N/C ratios about ten times the initial value. Below 100 km/s a wide range of N/C values is observed. The relation between N/C and surface gravity is well reproduced by various sets of models. Some evolutionary models including rotation are also able to consistently explain slowly rotating, highly enriched stars. This is due to differential rotation which efficiently transports nucleosynthesis products and allows the surface to rotate slower than the core. In addition, angular momentum removal by winds amplifies surface braking on the main sequence. Comparison of the surface composition of O7-8 giant stars with a sample of B stars with initial masses about four times smaller reveal that chemical enrichment scales with initial mass, as expected from theory. Although evolutionary models that include rotation face difficulties in explaining the chemical properties of O- and B-type stars at low metallicity, some of them can consistently account for the properties of main-sequence Galactic O stars in the mass range 25-40 Msun

Macroturbulence, introduced as a fudge to reproduce the width and shape of stellar absorption lines, reflects gas motions in stellar atmospheres. While in cool stars, it is thought to be caused by convection zones immediately beneath the stellar surface, the origin of macroturbulence in hot stars is still under discussion. Recent works established a correlation between the turbulent-to-total pressure ratio inside the envelope of stellar models and the macroturbulent velocities observed in corresponding Galactic stars. To probe this connection further, we evaluated the turbulent pressure that arises in the envelope convective zones of stellar models in the mass range 1-125 Msun based on the mixing-length theory and computed for metallicities of the Large and Small Magellanic Cloud. We find that the turbulent pressure contributions in models with these metallicities located in the hot high-luminosity part of the Hertzsprung-Russel (HR) diagram is lower than in similar models with solar metallicity, whereas the turbulent pressure in low-metallicity models populating the cool part of the HR-diagram is not reduced. Based on our models, we find that the currently available observations of hot massive stars in the Magellanic Clouds appear to support a connection between macroturbulence and the turbulent pressure in stellar envelopes. Multidimensional simulations of sub-surface convection zones and a larger number of high-quality observations are necessary to test this idea more rigorously.

A significant fraction of massive main-sequence stars show strong, large-scale magnetic fields. The origin of these fields, their lifetimes, and their role in shaping the characteristics and evolution of massive stars are currently not well understood. We compile a catalogue of 389 massive main-sequence stars, 61 of which are magnetic, and derive their fundamental parameters and ages. The two samples contain stars brighter than magnitude 9 in the V band and range in mass between 5 and 100 Msun. We find that the fractional main-sequence age distribution of all considered stars follows what is expected for a magnitude limited sample, while that of magnetic stars shows a clear decrease towards the end of the main sequence. This dearth of old magnetic stars is independent of the choice of adopted stellar evolution tracks, and appears to become more prominent when considering only the most massive stars. We show that the decreasing trend in the distribution is significantly stronger than expected from magnetic flux conservation. We also find that binary rejuvenation and magnetic suppression of core convection are unlikely to be responsible for the observed lack of older magnetic massive stars, and conclude that its most probable cause is the decay of the magnetic field, over a time span longer than the stellar lifetime for the lowest considered masses, and shorter for the highest masses. We then investigate the spin-down ages of the slowly rotating magnetic massive stars and find them to exceed the stellar ages by far in many cases. The high fraction of very slowly rotating magnetic stars thus provides an independent argument for a decay of the magnetic fields.

Since its launch, the X-ray and gamma-ray observatory INTEGRAL satellite has revealed a new class of high-mass X-ray binaries (HMXB) displaying fast flares and hosting supergiant companion stars. Optical and infrared (OIR) observations in a multi-wavelength context are essential to understand the nature and evolution of these newly discovered celestial objects. The goal of this multiwavelength study (from ultraviolet to infrared) is to characterise the properties of IGR J16465-4507, to confirm its HMXB nature and that it hosts a supergiant star. We analysed all OIR, photometric and spectroscopic observations taken on this source, carried out at ESO facilities. Using spectroscopic data, we constrained the spectral type of the companion star between B0.5 and B1 Ib, settling the debate on the true nature of this source. We measured a high rotation velocity of v = 320 +/- 8 km/s from fitting absorption and emission lines in a stellar spectral model. We then built a spectral energy distribution from photometric observations to evaluate the origin of the different components radiating at each energy range. We finally show that, having accurately determined the spectral type of the early-B supergiant in IGR J16465-4507, we firmly support its classification as an intermediate supergiant fast X-ray transient (SFXT).

On the basis of the Galactic O-Star Spectroscopic Survey (GOSSS), a detailed systematic investigation of the O Vz stars is presented. The currently used spectral classification criteria are rediscussed, and the Vz phenomenon is recalibrated through the addition of a quantitative criterion based on the equivalent widths of the He I 4471, He II 4542, and He II 4686 spectral lines. The GOSSS O Vz and O V populations resulting from the newly adopted spectral classification criteria are comparatively analyzed. The locations of the O Vz stars are probed, showing a concentration of the most extreme cases toward the youngest star forming regions. The occurrence of the Vz spectral peculiarity in a solar-metallicity environment, as predicted by the fastwind code, is also investigated, confirming the importance of taking into account several processes for the correct interpretation of the phenomenon.

We investigate the characteristics of two newly discovered short-period, double-lined, massive binary systems, VFTS 450 (O9.7$\;$II--Ib$\,$+$\,$O7::) and VFTS 652 (B1$\;$Ib$\,+\,$O9:$\;$III:). We perform model-atmosphere analyses to characterise the photospheric properties of both members of each binary (denoting the `primary' as the spectroscopically more conspicuous component). Radial velocities and optical photometry are used to estimate the binary-system parameters. We estimate $T_{\rm eff}=27$ kK, $\log{(g)}=2.9$ (cgs) for the VFTS 450 primary spectrum (34kK, 3.6: for the secondary spectrum); and $T_{\rm eff} = 22$kK, $\log{(g)}=2.8$ for the VFTS 652 primary spectrum (35kK, 3.7: for the secondary spectrum). Both primaries show surface nitrogen enrichments (of more than 1 dex for VFTS 652), and probable moderate oxygen depletions relative to reference LMC abundances. We determine orbital periods of 6.89d and 8.59d for VFTS 450 and VFTS 652, respectively, and argue that the primaries must be close to filling their Roche lobes. Supposing this to be the case, we estimate component masses in the range $\sim$20--50M$_\odot$. The secondary spectra are associated with the more massive components, suggesting that both systems are high-mass analogues of classical Algol systems, undergoing case-A mass transfer. Difficulties in reconciling the spectroscopic analyses with the light-curves and with evolutionary considerations suggest that the secondary spectra are contaminated by (or arise in) accretion disks.

The major mass fraction of the envelope of hot luminous stars is radiatively stable. However, the partial ionisation of hydrogen, helium and iron gives rise to extended sub-surface convection zones in all of them. In this work, we investigate the effect of the pressure induced by the turbulent motion in these zones based on the mixing length theory, and search for observable consequences. We find that the turbulent pressure fraction can amount up to ~5% in OB supergiants, and to ~30% in cooler supergiants. The resulting structural changes are, however, not significantly affecting the evolutionary tracks compared to previous calculations. Instead, a comparison of macroturbulent velocities derived from high quality spectra of OB stars with the turbulent pressure fraction obtained in corresponding stellar models reveals a strong correlation of these two quantities. We discuss a possible physical connection, and conclude that turbulent pressure fluctuations may drive high-order oscillations, which - as conjectured earlier - manifest themselves as macroturbulence in the photospheres of hot luminous stars.

Massive stars burn hydrogen through the CNO cycle during most of their evolution. When mixing is efficient, or when mass transfer in binary systems happens, chemically processed material is observed at the surface of O and B stars. ON stars show stronger lines of nitrogen than morphologically normal counterparts. Whether this corresponds to the presence of material processed through the CNO cycle or not is not known. Our goal is to answer this question. We perform a spectroscopic analysis of a sample of ON stars with atmosphere models. We determine the fundamental parameters as well as the He, C, N, and O surface abundances. We also measure the projected rotational velocities. We compare the properties of the ON stars to those of normal O stars. We show that ON stars are usually helium-rich. Their CNO surface abundances are fully consistent with predictions of nucleosynthesis. ON stars are more chemically evolved and rotate - on average - faster than normal O stars. Evolutionary models including rotation cannot account for the extreme enrichment observed among ON main sequence stars. Some ON stars are members of binary systems, but others are single stars as indicated by stable radial velocities. Hence, mass transfer is not a simple explanation for the observed chemical properties. We conclude that ON stars show extreme chemical enrichment at their surface, consistent with nucleosynthesis through the CNO cycle. Its origin is not clear at present.

Within the context of the "B fields in OB stars (BOB)" collaboration, we used the HARPSpol spectropolarimeter to observe the early B-type stars beta CMa (HD44743; B1 II/III) and epsilon CMa (HD52089; B1.5 II). For both stars, we consistently detected the signature of a weak (<30 G in absolute value) longitudinal magnetic field. We determined the physical parameters of both stars and characterise their X-ray spectrum. For beta CMa, our mode identification analysis led to determining a rotation period of 13.6+/-1.2 days and of an inclination angle of the rotation axis of 57.6+/-1.7 degrees, with respect to the line of sight. On the basis of these measurements and assuming a dipolar field geometry, we derived a best fitting obliquity of ~22 degrees and a dipolar magnetic field strength (Bd) of ~100 G (60<Bd<230 G within 1 sigma), below what is typically found for other magnetic massive stars. For epsilon CMa we could only determine a lower limit on the dipolar magnetic field strength of 13 G. For this star, we determine that the rotation period ranges between 1.3 and 24 days. Both stars are expected to have a dynamical magnetosphere. We also conclude that both stars are most likely core hydrogen burning and that they have spent more than 2/3 of their main sequence lifetime. A histogram of the distribution of the dipolar magnetic field strength for the magnetic massive stars known to date does not show the magnetic field "desert" observed instead for intermediate-mass stars. The biases involved in the detection of (weak) magnetic fields in massive stars with the currently available instrumentation and techniques imply that weak fields might be more common than currently observed. Our results show that, if present, even relatively weak magnetic fields are detectable in massive stars and that more observational effort is probably still needed to properly access the magnetic field incidence.

The Fourier Transform method is a popular tool to derive the rotational velocities of stars from their spectral line profiles. However, its domain of validity does not include line-profile variables with time-dependent profiles. We investigate the performance of the method for such cases, by interpreting the line-profile variations of spotted B stars, and of pulsating B tars, as if their spectral lines were caused by uniform surface rotation along with macroturbulence. We perform time-series analysis and harmonic least-squares fitting of various line diagnostics and of the outcome of several implementations of the Fourier Transform method. We find that the projected rotational velocities derived from the Fourier Transform vary appreciably during the pulsation cycle whenever the pulsational and rotational velocity fields are of similar magnitude. The macroturbulent velocities derived while ignoring the pulsations can vary with tens of km/s during the pulsation cycle. The temporal behaviour of the deduced rotational and macroturbulent velocities are in antiphase with each other. The rotational velocity is in phase with the second moment of the line profiles. The application of the Fourier method to stars with considerable pulsational line broadening may lead to an appreciable spread in the values of the rotation velocity, and, by implication, of the deduced value of the macroturbulence. These two quantities should therefore not be derived from single snapshot spectra if the aim is to use them as a solid diagnostic for the evaluation of stellar evolution models of slow to moderate rotators.

Detailed spectral classifications are presented for 352 O-B0 stars in the VLT-FLAMES Tarantula Survey, of which 213 O-type are of sufficient quality for further morphological analysis. Among them, six subcategories of special interest are distinguished. (1) Several new examples of the earliest spectral types O2-O3 have been found. (2) A group of extremely rapidly rotating main-sequence objects has been isolated, including the largest $v\sin i$ values known, the spatial and radial-velocity distributions of which suggest ejection from the two principal ionizing clusters. (3) Several new examples of the evolved, rapidly rotating Onfp class show similar evidence. (4) No fewer than 48 members of the Vz category, hypothesized to be on or near the ZAMS, are found in this sample; in contrast to the rapid rotators, they are strongly concentrated to the ionizing clusters, supporting their interpretation as very young objects, as do their relatively faint absolute magnitudes. (5) A surprisingly large fraction of the main-sequence spectra belong to the recently recognized V((fc)) class, with C III emission lines of similar strength to the usual N III in V((f)) spectra; there are also six objects with very high-quality data but no trace of either mission feature, presenting new challenges to physical interpretations. (6) Five spectra with morphologically enhanced nitrogen lines have been detected. Absolute visual magnitudes have been derived for each star with individual extinction laws, and composite HRDs provide evidence of the multiple generations present in this field. Associations with X-ray sources are noted. Further analyses of this unique dataset underway will provide new insights into the evolution of massive stars and starburst clusters.

Stellar-mass black holes have all been discovered through X-ray emission, which arises from the accretion of gas from their binary companions (this gas is either stripped from low-mass stars or supplied as winds from massive ones). Binary evolution models also predict the existence of black holes accreting from the equatorial envelope of rapidly spinning Be-type stars (stars of the Be type are hot blue irregular variables showing characteristic spectral emission lines of hydrogen). Of the ~80 Be X-ray binaries known in the Galaxy, however, only pulsating neutron stars have been found as companions. A black hole was formally allowed as a solution for the companion to the Be star MWC 656 (also known as HD 215227), although that was based on a single radial velocity curve of the Be star, a mistaken spectral classification and rough estimates of the inclination angle. Here we report observations of an accretion disk line mirroring the orbit of the Be star. This, together with an improved radial velocity curve of the Be star through fitting sharp Fe II profiles from the equatorial disk, and a refined Be classification (to that of a B1.5-B2 III star), reveals a black hole of 3.8 to 6.9 solar masses orbiting MWC 656, the candidate counterpart of the gamma-ray source AGL J2241+4454. The black hole is X-ray quiescent and fed by a radiatively inefficient accretion flow giving a luminosity less than 1.6 x 10-7 times the Eddington luminosity. This implies that Be binaries with black-hole companions are difficult to detect by conventional X-ray surveys.

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

{We aim to detect and interpret photometric and spectroscopic variability of the bright CoRoT B-type supergiant target HD\,46769 ($V=5.79$). We also attempt to detect a magnetic field in the target.} {We analyse a 23-day oversampled CoRoT light curve after detrending, as well as spectroscopic follow-up data, by using standard Fourier analysis and Phase Dispersion Minimization methods. We determine the fundamental parameters of the star, as well as its abundances from the most prominent spectral lines. We perform a Monte Carlo analysis of spectropolarimetric data to obtain an upper limit of the polar magnetic field, assumping a dipole field.} {In the CoRoT data, we detect a dominant period of 4.84\,d with an amplitude of 87\,ppm, and some of its (sub-)multiples. Given the shape of the phase-folded light curve and the absence of binary motion, we interpret the dominant variability in terms of rotational modulation, with a rotation period of 9.69\,d. Subtraction of the rotational modulation signal does not reveal any sign of pulsations. Our results are consistent with the absence of variability in the Hipparcos light curve. The spectroscopy leads to a projected rotational velocity of 72$\pm 2$\,km\,s$^{-1}$ and does not reveal periodic variability nor the need to invoke macroturbulent line broadening. No signature of a magnetic field is detected in our data. A field stronger than $\sim 500$\,G at the poles can be excluded, unless the possible non-detected field were more complex than dipolar.} {The absence of pulsations and of macroturbulence of this evolved B-type supergiant is placed into context of instability computations and of observed variability of evolved B-type stars.}

While collisionally excited lines in HII regions allow one to easily probe the chemical composition of the interstellar medium in galaxies, the possible presence of important temperature fluctuations casts some doubt on the derived abundances. To provide new insights into this question, we have carried out a detailed study of a giant HII region, H 1013, located in the galaxy M101, for which many observational data exist and which has been claimed to harbour temperature fluctuations at a level of t^2 = 0.03-0.06. We have first complemented the already available optical observational datasets with a mid-infrared spectrum obtained with the Spitzer Space Telescope. Combined with optical data, this spectrum provides unprecedented information on the temperature structure of this giant HII region. A preliminary analysis based on empirical temperature diagnostics suggests that temperature fluctuations should be quite weak. We have then performed a detailed modelling using the pyCloudy package based on the photoionization code Cloudy. We have been able to produce photoionization models constrained by the observed Hb surface brightness distribution and by the known properties of the ionizing stellar population than can account for most of the line ratios within their uncertainties. Since the observational constraints are both strong and numerous, this argues against the presence of significant temperature fluctuations in H 1013. The oxygen abundance of our best model is 12 + log O/H = 8.57, as opposed to the values of 8.73 and 8.93 advocated by Esteban et al. (2009) and Bresolin (2007), respectively, based on the significant temperature fluctuations they derived. However, our model is not able to reproduce the intensities of the oxygen recombination lines . This cannot be attributed to observational uncertainties and requires an explanation other than temperature fluctuations.

The Tarantula Nebula (30 Dor) is a spectacular star-forming region in the Large Magellanic Cloud, seen through gas in the Galactic Disc and Halo. Diffuse Interstellar Bands offer a unique probe of the diffuse, cool-warm gas in these regions. The aim is to use DIBs as diagnostics of the local interstellar conditions, whilst at the same time deriving properties of the yet-unknown carriers. Spectra of over 800 early-type stars from the VLT Flames Tarantula Survey (VFTS) were analysed. Maps were created, separately, for the Galactic and LMC absorption in the DIBs at 4428 and 6614 Ang and - in a smaller region near the central cluster R136 - neutral sodium (Na I D); we also measured the DIBs at 5780 and 5797 Ang. The maps show strong 4428 and 6614 Ang DIBs in the quiescent cloud complex to the south of 30 Dor but weak absorption in the harsher environments to the north (bubbles) and near the OB associations. The Na maps show at least five kinematic components in the LMC and a shell-like structure surrounding R136, and small-scale structure in the Milky Way. The strengths of the 4428, 5780, 5797 and 6614 Ang DIBs are correlated, also with Na absorption and visual extinction. The strong 4428 Ang DIB is present already at low Na column density but the 6614, 5780 and 5797 Ang DIBs start to be detectable at subsequently larger Na column densities. The relative strength of the 5780 and 5797 Ang DIBs clearly confirm the Tarantula Nebula and Galactic high-latitude gas to represent a harsh radiation environment. The resilience of the 4428 Ang DIB suggests its carrier is large, compact and neutral. Structure is detected in the distribution of cool-warm gas on scales between one and >100 pc in the LMC and as little as 0.01 pc in the Sun's vicinity. Stellar winds from the central cluster R136 have created an expanding shell; some infalling gas is also detected, reminiscent of a galactic "fountain".

The spectral and spatial characteristics of two special categories of O stars found in the VFTS dataset are presented. One of them comprises very rapid rotators, including several more extreme than any previously known. These objects are distributed around the peripheries of the main 30 Doradus clusters, suggesting a runaway nature for which their radial velocities already provide preliminary supporting evidence. The other category consists of a large number of Vz stars, previously hypothesized on spectroscopic grounds to be on or very near the ZAMS. Their distribution is the inverse of that of the rapid rotators: the Vz are strongly concentrated to the ionizing clusters, plus a newly recognized band of recent and current star formation to the north, which provides strong circumstantial evidence for their extreme youth.

Six O-type stars were observed continuously by the CoRoT satellite during a 34.3-day run. The unprecedented quality of the data allows us to detect even low-amplitude stellar pulsations in some of these stars (HD 46202 and the binaries HD 46149 and Plaskett's star). These cover both opacity-driven modes and solar-like stochastic oscillations, both of importance to the asteroseismological modelling of O stars. Additional effects can be seen in the CoRoT light curves, such as binarity and rotational modulation. Some of the hottest O-type stars (HD 46223, HD 46150 and HD 46966) are dominated by the presence of red-noise: we speculate that this is related to a sub-surface convection zone.

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

The star HD 193322 is a remarkable multiple system of massive stars that lies at the heart of the cluster Collinder 419. Here we report on new spectroscopic observations and radial velocities of the narrow-lined component Ab1 that we use to determine its orbital motion around a close companion Ab2 ($P = 312$ d) and around a distant third star Aa ($P = 35$ y).We have also obtained long baseline interferometry of the target in the $K^\prime$-band with the CHARA Array that we use in two ways. First, we combine published speckle interferometric measurements with CHARA separated fringe packet measurements to improve the visual orbit for the wide Aa,Ab binary. Second, we use measurements of the fringe packet from Aa to calibrate the visibility of the fringes of the Ab1,Ab2 binary, and we analyze these fringe visibilities to determine the visual orbit of the close system. The two most massive stars, Aa and Ab1, have masses of approximately 21 and $23 M_\odot$, respectively, and their spectral line broadening indicates that they represent extremes of fast and slow projected rotational velocity, respectively.

Early B-type stars are invaluable indicators for elemental abundances of their birth environments. In contrast to the surrounding neutral interstellar matter (ISM) and HII regions their chemical composition is unaffected by depletion onto dust grains and by the derivation of different abundances from recombination and collisional lines. In combination with ISM or nebular gas-phase abundances they facilitate the dust-phase composition to be constrained. Precise abundances of C, N, Mg, Ne, Fe in early B-type stars in the Orion star-forming region are determined in order to: a) review previous determinations using a self-consistent quantitative spectral analysis based on modern stellar atmospheres and recently updated model atoms, b) complement results found in Paper I for oxygen and silicon, c) establish an accurate and reliable set of stellar metal abundances to constrain the dust-phase composition of the Orion HII region in Paper II of the series. A detailed, self-consistent spectroscopic study of a sample of 13 narrow-lined B0V-B2V stars in Ori OB1 is performed. High-quality spectra obtained with FIES@NOT are analysed using a non-LTE method and line-profile fitting techniques, validating the approach by comparison with results obtained in Paper I using line-blanketed non-LTE model atmospheres and a curve-of-growth analysis. The two independent analysis strategies give consistent results for basic stellar parameters and abundances of oxygen and silicon. The extended analysis to C, N, Mg, Ne, and Fe finds a high degree of chemical homogeneity, with the 1sigma-scatter adopting values of 0.03--0.07 dex around the mean for the various elements. Present-day abundances from B-type stars in Ori OB1 are compatible at similar precision with cosmic abundance standard values as recently established from early-type stars in the solar neighbourhood and also with the Sun. (abridged)

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

The Tarantula survey is an ESO Large Programme which has obtained multi-epochs spectroscopy of over 800 massive stars in the 30 Dor region in the Large Magelanic Cloud. Here we briefly describe the main drivers of the survey and the observational material derived.

Using the CoRoT space based photometry of the O-type binary HD46149, stellar atmospheric effects related to rotation can be separated from pulsations, because they leave distinct signatures in the light curve. This offers the possibility of characterising and exploiting any pulsations seismologically. Combining high-quality space based photometry, multi-wavelength photometry, spectroscopy and constraints imposed by binarity and cluster membership, the detected pulsations in HD46149 are analyzed and compared with those for a grid of stellar evolutionary models in a proof-of-concept approach. We present evidence of solar-like oscillations in a massive O-type star, and show that the observed frequency range and spacings are compatible with theoretical predictions. Thus, we unlock and confirm the strong potential of this seismically unexplored region in the HR diagram.

The planetary nebula TS 01 (also called PN G 135.9+55.9 or SBS 1150+599A), with its record-holding low oxygen abundance and its double degenerate close binary core (period 3.9 h), is an exceptional object located in the Galactic halo. We have secured observational data in a complete wavelength range in order to pin down the abundances of half a dozen elements in the nebula. The abundances are obtained via detailed photoionization modelling taking into account all the observational constraints (including geometry and aperture effects) using the pseudo-3D photoionization code Cloudy_3D. The spectral energy distribution of the ionizing radiation is taken from appropriate model atmospheres. Both stellar components contribute to the ionization: the ``cool'' one provides the bulk of hydrogen ionization, and the ``hot'' one is responsible for the presence of the most highly charged ions, which explains why previous attempts to model the nebula experienced difficulties. The nebular abundances of C, N, O, and Ne are found to be respectively, 1/3.5, 1/4.2, 1/70, and 1/11 of the Solar value, with uncertainties of a factor 2. Thus the extreme O deficiency of this object is confirmed. The abundances of S and Ar are less than 1/30 of Solar. Standard models of stellar evolution and nucleosynthesis cannot explain the abundance pattern observed in the nebula. To obtain an extreme oxygen deficiency in a star whose progenitor has an initial mass of about 1 msun requires an additional mixing process, which can be induced by stellar rotation and/or by the presence of the close companion. We have computed a stellar model with initial mass of 1 msun, appropriate metallicity, and initial rotation of 100 kms, and find that rotation greatly improves the agreement between the predicted and observed abundances.

The Tarantula Survey is an ambitious ESO Large Programme that has obtained multi-epoch spectroscopy of over 1,000 massive stars in the 30 Doradus region of the Large Magellanic Cloud. Here we introduce the scientific motivations of the survey and give an overview of the observational sample. Ultimately, quantitative analysis of every star, paying particular attention to the effects of rotational mixing and binarity, will be used to address fundamental questions in both stellar and cluster evolution.

(abridged) Non-ionizing stellar continua are a source of photons for continuum pumping in the hydrogen Lyman transitions. In the environments where these transitions are optically thick, deexcitation occurs through higher series lines, so that the flux in these lines has a fluorescent contribution in addition to recombination; in particular, Balmer emissivities are systematically enhanced above case B. The effectiveness of such mechanism in HII regions and the adequacy of photoionization models as a tool to study it are the two main focuses of this work. We find that photoionization models of H II regions illuminated by low-resolution population synthesis models significantly overpredict the fluorescent contribution to the Balmer lines. Conversely, photoionization models in which the non-ionizing part of the continuum is omitted or is not transferred underpredict the fluorescent contribution to the Balmer lines, producing a bias of similar amplitude in the opposite direction. In this paper, we carry out realistic estimations of the fluorescent Balmer intensity and discuss the variations to be expected as the simulated observational setup and the stellar population's parameters are varied. In all the cases explored, we find that fluorescent excitation provides a significant contribution. We also show that differential fluorescent enhancement may produce line-of-sight differences in the Balmer decrement, mimicking interstellar extinction. Fluorescent excitation emerges from our study as a small but important mechanism for the enhancement of Balmer lines, which should be taken into account in the abundance analysis of photoionized regions, particularly in the case of high-precision applications such as the determination of primordial helium.

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

We provide atmospheric parameters and rotational velocities of a large sample (~400) of O- and early B-type stars, analysed in a homogeneous and consistent manner, for use in constraining theoretical models. Comparison of the rotational velocities with evolutionary tracks suggest that the end of core hydrogen burning occurs later than currently predicted. We also show that the large number of the luminous blue supergiants observed in the fields are unlikely to have directly evolved from main-sequence massive O-type stars as neither their low rotational velocities or position on the H-R diagram are predicted. We suggest that blue-loops or mass-transfer binary systems may populate the blue supergiant regime. By comparing the rotational velocity distributions of the Magellanic Cloud stars to a similar Galactic sample we find that (at 3sigma confidence level) massive stars (above 8Msun) in the SMC rotate faster than those in the solar neighbourhood. However there appears to be no significant difference between the rotational velocity distributions in the Galaxy and the LMC. We find that the vsini distributions in the SMC and LMC can modelled with an intrinsic rotational velocity distribution which is a Gaussian peaking at 175km/s (SMC) and 100km/s (LMC). We find that in NGC346 in the SMC, the 10-25Msun main-sequence stars appear to rotate faster than their higher mass counterparts. Recently Yoon et al. (2006) have determined rates of GRBs by modelling rapidly rotating massive star progenitors. Our measured rotational velocity distribution for the 10-25Msun stars is peaked at slightly higher velocities than they assume, supporting the idea that GRBs could come from rapid rotators with initial masses as low as 14Msun at low metallicities. (abridged).

High-resolution spectra for 24 SMC and Galactic B-type supergiants have been analysed to estimate the contributions of both macroturbulence and rotation to the broadening of their metal lines. Two different methodologies are considered, viz. goodness-of-fit comparisons between observed and theoretical line profiles and identifying zeros in the Fourier transforms of the observed profiles. The advantages and limitations of the two methods are briefly discussed with the latter techniques being adopted for estimated projected rotational velocities (\vsini) but the former being used to estimate macroturbulent velocities. Only one SMC supergiant, SK 191, shows a significant degree of rotational broadening (\vsini $\simeq$ 90 \kms). For the remaining targets, the distribution of projected rotational velocities are similar in both our Galactic and SMC samples with larger values being found at earlier spectral types. There is marginal evidence for the projected rotational velocities in the SMC being higher than those in the Galactic targets but any differences are only of the order of 5-10 \kms, whilst evolutionary models predict differences in this effective temperature range of typically 20 to 70 \kms. The combined sample is consistent with a linear variation of projected rotational velocity with effective temperature, which would imply rotational velocities for supergiants of 70 \kms at an effective temperature of 28 000 K (approximately B0 spectral type) decreasing to 32 \kms at 12 000 K (B8 spectral type). For all targets, the macroturbulent broadening would appear to be consistent with a Gaussian distribution (although other distributions cannot be discounted) with an $\frac{1}{e}$ half-width varying from approximately 20 \kms at B8 to 60 \kms at B0 spectral types.

We present the results of a spectroscopic analysis of the Trapezium cluster stars inside the Orion Nebula. The rotational velocities have been obtained using Fourier analysis method, finding agreement with values derived from the usual method, based on linewidth measurements. The rotational velocity derived for theta 1 Ori C through this method is consistent with the variability of some of its spectral features that have a period of 15.42 days. By means of the fit of H, HeI and HeII observed profiles with FASTWIND synthetic profiles, stellar parameters and wind characteristics have been derived. This methodology let us estimate the errors associated with these parameters. It is found that macroturbulence effects have to be included for a good fit to the HeI-II lines in the spectrum of theta 1 Ori C. By means of a very accurate study, oxygen abundances have been derived for the three B0.5V stars theta 1 Ori A, D and theta 2 Ori B. Final abundances are consistent with the nebular gas-phase results presented in Esteban et al. (2004) and are lower than those given by Cunha and Lambert (1994). Our results suggest a lower dust depletion factor of oxygen than previous estimations for the Orion nebula.

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