Characterization of the dynamics of massive star systems and the astrophysical properties of the interacting components are a prerequisite for understanding their formation and evolution. Optical interferometry at milliarcsecond resolution is a key observing technique for resolving high-mass multiple compact systems. Here we report on VLTI/GRAVITY, Magellan/FIRE, and MPG2.2m/FEROS observations of the late-O/early-B type system HD 93206 A, which is a member of the massive cluster Collinder 228 in the Carina nebula complex. With a total mass of about 90 M_sun, it is one of the most compact massive-quadruple systems known. In addition to measuring the separation and position angle of the outer binary Aa - Ac, we observe BrG and HeI variability in phase with the orbital motion of the two inner binaries. From the differential phases analysis, we conclude that the BrG emission arises from the interaction regions within the components of the individual binaries, which is consistent with previous models for the X-ray emission of the system based on wind-wind interaction. With an average 3-sigma deviation of ~15 deg, we establish an upper limit of p ~ 0.157 mas (0.35 AU) for the size of the BrG line-emitting region. Future interferometric observations with GRAVITY using the 8m UTs will allow us to constrain the line-emitting regions down to angular sizes of 20 uas (0.05 AU at the distance of the Carina nebula).

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 present observations from the Gaia-ESO Survey in the lines of H$\alpha$, [N II], [S II] and He I of nebular emission in the central part of the Carina Nebula. We investigate the properties of the two already known kinematic components (approaching and receding, respectively), which account for the bulk of emission. Moreover, we investigate the features of the much less known low-intensity high-velocity (absolute RV $>$50 km/s) gas emission. We show that gas giving rise to H$\alpha$ and He I emission is dynamically well correlated, but not identical, to gas seen through forbidden-line emission. Gas temperatures are derived from line-width ratios, and densities from [S II] doublet ratios. The spatial variation of N ionization is also studied, and found to differ between the approaching and receding components. The main result is that the bulk of the emission lines in the central part of Carina arises from several distinct shell-like expanding regions, the most evident found around $\eta$ Car, the Trumpler 14 core, and the star WR25. Such "shells" are non-spherical, and show distortions probably caused by collisions with other shells or colder, higher-density gas. Part of them is also obscured by foreground dust lanes, while only very little dust is found in their interior. Preferential directions, parallel to the dark dust lanes, are found in the shell geometries and physical properties, probably related to strong density gradients in the studied region. We also find evidence that the ionizing flux emerging from $\eta$ Car and the surrounding Homunculus nebula varies with polar angle. The high-velocity components in the wings of H$\alpha$ are found to arise from expanding dust reflecting the $\eta$ Car spectrum.

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

Optical spectroscopic monitoring has been conducted of two O stars in the Small and one in the Large Magellanic Cloud, the spectral characteristics of which place them in the Of?p category, which has been established in the Galaxy to consist of oblique magnetic rotators. All of these Magellanic stars show systematic spectral variations typical of the Of?p class, further strengthening their magnetic candidacy to the point of virtual certainty. The spectral variations are related to photometric variations derived from OGLE data by Naze et al. (2015) in a parallel study, which yields rotational periods for two of them. Now circular spectropolarimetry is required to measure their fields, and ultraviolet spectroscopy to further characterize their low-metallicity, magnetically confined winds, in support of hydrodynamical analyses.

We report on both high-precision photometry from the MOST space telescope and ground-based spectroscopy of the triple system delta Ori A consisting of a binary O9.5II+early-B (Aa1 and Aa2) with P = 5.7d, and a more distant tertiary (O9 IV P > 400 yrs). This data was collected in concert with X-ray spectroscopy from the Chandra X-ray Observatory. Thanks to continuous coverage for 3 weeks, the MOST light curve reveals clear eclipses between Aa1 and Aa2 for the first time in non-phased data. From the spectroscopy we have a well constrained radial velocity curve of Aa1. While we are unable to recover radial velocity variations of the secondary star, we are able to constrain several fundamental parameters of this system and determine an approximate mass of the primary using apsidal motion. We also detected second order modulations at 12 separate frequencies with spacings indicative of tidally influenced oscillations. These spacings have never been seen in a massive binary, making this system one of only a handful of such binaries which show evidence for tidally induced pulsations.

The B fields in OB stars (BOB) survey is an ESO large programme collecting spectropolarimetric observations for a large number of early-type stars in order to study the occurrence rate, properties, and ultimately the origin of magnetic fields in massive stars. As of July 2014, a total of 98 objects were observed over 20 nights with FORS2 and HARPSpol. Our preliminary results indicate that the fraction of magnetic OB stars with an organised, detectable field is low. This conclusion, now independently reached by two different surveys, has profound implications for any theoretical model attempting to explain the field formation in these objects. We discuss in this contribution some important issues addressed by our observations (e.g., the lower bound of the field strength) and the discovery of some remarkable objects.

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.

The Photometry Data Model (PhotDM) standard describes photometry filters, photometric systems, magnitude systems, zero points and its interrelation with the other IVOA data models through a simple data model. Particular attention is given necessarily to optical photometry where specifications of magnitude systems and photometric zero points are required to convert photometric measurements into physical flux density units.

We compile the first comprehensive census of hot luminous stars in the 30 Doradus (30 Dor) star forming region of the LMC. The census extends to a radius of 10arcmin (150pc) from the central cluster, R136. Stars were selected photometrically and combined with the latest spectral types. 1145 candidate hot luminous stars were identified of which >700 were considered genuine early type stars that contribute to feedback. We assess the spectroscopic completeness to be 85% in outer regions (>5pc) but fall to 35% in the vicinity of R136, giving a total of 500 hot luminous stars with spectroscopy. Stellar calibrations and models were used to obtain their physical parameters before integrated values were compared to global observations and the population synthesis code, Starburst99. The 31 W-R and Of/WN stars made large contributions to the total ionising and wind luminosities of ~40% and ~50%, respectively. Stars with Minit>100Msun also showed high contributions to the global feedback, ~25% in both cases. Such massive stars are not accounted for by the current Starburst99 code, which underestimated the ionising and wind luminosities of R136 by factors of ~2 and ~9, respectively. The census inferred a SFR of 0.073+/-0.04Msun/yr for 30 Dor, typically higher than results from popular SFR calibrations. However, it remained consistent with a far-UV luminosity tracer and a combined Halpha and mid-infrared tracer, but only after correcting for Halpha extinction. The global ionising output exceeded measurements from the associated gas and dust, suggesting ~6(+55/-6)% of ionising photons escape the region. When studying very luminous star forming regions, it is therefore essential to include the most massive stars to ensure a reliable energy budget. If 30 Dor is typical of other large star forming regions, estimates of the SFR will be underpredicted if this escape fraction is not accounted for.(abridged)

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

Aims. We analyze the multiplicity properties of the massive O-type star population. With 360 O-type stars, this is the largest homogeneous sample of massive stars analyzed to date. Methods. We use multi-epoch spectroscopy and variability analysis to identify spectroscopic binaries. We also use a Monte-Carlo method to correct for observational biases. Results. We observe a spectroscopic binary fraction of 0.35\pm0.03, which corresponds to the fraction of objects displaying statistically significant radial velocity variations with an amplitude of at least 20km/s. We compute the intrinsic binary fraction to be 0.51\pm0.04. We adopt power-laws to describe the intrinsic period and mass-ratio distributions: f_P ~ (log P)^\pi\ (with 0.15 < log P < 3.5) and f_q ~ q^\kappa\ with 0.1 < q=M_2/M_1 < 1.0. The power-law indexes that best reproduce the observed quantities are \pi = -0.45 +/- 0.30 and \kappa = -1.0\pm0.4. The obtained period distribution thus favours shorter period systems compared to an Oepik law. The mass ratio distribution is slightly skewed towards low mass ratio systems but remains incompatible with a random sampling of a classical mass function. The binary fraction seems mostly uniform across the field of view and independent of the spectral types and luminosity classes. The binary fraction in the outer region of the field of view (r > 7.8', i.e. approx117 pc) and among the O9.7 I/II objects are however significantly lower than expected from statistical fluctuations. Conclusions. Using simple evolutionary considerations, we estimate that over 50% of the current O star population in 30 Dor will exchange mass with its companion within a binary system. This shows that binary interaction is greatly affecting the evolution and fate of massive stars, and must be taken into account to correctly interpret unresolved populations of massive stars.

Based on an analysis of data obtained with the Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST) we report the identification of two distinct stellar populations in the core of the giant HII region 30Doradus in the Large Magellanic Cloud. The most compact and richest component coincides with the center of R136 and is ~1 Myr younger than a second more diffuse clump, located ~5.4 pc toward the northeast. We note that published spectral types of massive stars in these two clumps lend support to the proposed age difference. The morphology and age difference between the two sub-clusters suggests that an ongoing merger may occurring within the core of 30Doradus. This finding is consistent with the predictions of models of hierarchical fragmentation of turbulent giant molecular clouds, according to which star clusters would be the final products of merging smaller sub-structures.

With new data from the Galactic O-Star Spectroscopic Survey, we confirm and expand the ONn category of late-O, nitrogen-enriched (N), rapidly rotating (n) giants. In particular, we have discovered two "clones" (HD 102415 and HD 117490) of one of the most rapidly rotating O stars previously known (HD 191423, "Howarth's Star"). We compare the locations of these objects in the theoretical HR Diagram to those of slowly rotating ON dwarfs and supergiants. All ON giants known to date are rapid rotators, whereas no ON dwarf or supergiant is; but all ON stars are small fractions of their respective spectral-type/luminosity-class/rotational subcategories. The ONn giants, displaying both substantial processed material and high rotation at an intermediate evolutionary stage, may provide significant information about the development of those properties. They may have preserved high initial rotational velocities or been spun up by TAMS core contraction; but alternatively and perhaps more likely, they may be products of binary mass transfer. At least some of them are also runaway stars.

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

On the basis of an extensive new spectroscopic survey of Galactic O stars, we introduce the Ofc category, which consists of normal spectra with C III \lambda\lambda4647-4650-4652 emission lines of comparable intensity to those of the Of defining lines N III \lambda\lambda4634-4640-4642. The former feature is strongly peaked to spectral type O5, at all luminosity classes, but preferentially in some associations or clusters and not others. The relationships of this phenomenon to the selective C III \lambda5696 emission throughout the normal Of domain, and to the peculiar, variable Of?p category, for which strong C III \lambda\lambda4647-4650-4652 emission is a defining characteristic, are discussed. Magnetic fields have recently been detected on two members of the latter category. We also present two new extreme Of?p stars, NGC 1624-2 and CPD -28^{\circ}2561, bringing the number known in the Galaxy to five. Modeling of the behavior of these spectral features can be expected to better define the physical parameters of both normal and peculiar objects, as well as the atomic physics involved.

We present the analysis of high-resolution optical spectroscopic observations of the zero-age main-sequence O star Herschel 36 spanning six years. This star is definitely a multiple system, with at least three components detected in its spectrum. Based on our radial-velocity (RV) study, we propose a picture of a close massive binary and a more distant companion, most probably in wide orbit about each other. The orbital solution for the binary, whose components we identify as O9 V and B0.5 V, is characterized by a period of 1.5415 +/- 0.0006 days. With a spectral type O7.5 V, the third body is the most luminous component of the system and also presents RV variations with a period close to 498 days. Some possible hypotheses to explain the variability are briefly addressed and further observations are suggested.

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

We present the first high-spatial resolution near-infrared (NIR) imaging of NGC 604, obtained with the NICMOS camera aboard the Hubble Space Telescope (HST). These NICMOS broadband images reveal new NIR point sources, clusters, and diffuse structures. We found an excellent spatial correlation between the 8.4 GHz radio continuum and the 2.2mu-m nebular emission. Moreover, massive young stellar object candidates appear aligned with these radio peaks, reinforcing the idea that those areas are star-forming regions. Three different scaled OB associations are recognized in the NICMOS images. The brightest NIR sources in our images have properties that suggest that they are red supergiant stars, of which one of them was previously known. This preliminary analysis of the NICMOS images shows the complexity of the stellar content of the NGC 604 nebula.

We present far-UV spectroscopy obtained with HST for 48 blue objects in the core of 47 Tuc. Based on their position in a FUV-optical colour-magnitude diagram, these were expected to include cataclysmic variables (CVs), blue stragglers (BSs), white dwarfs (WDs) and other exotic objects. For a subset of these sources, we also construct FUV-NIR SEDs. Based on our analysis of this extensive data set, we report the following main results. (1) We spectroscopically confirm 3 previously known or suspected CVs via the detection of emission lines and find new evidence for dwarf nova eruptions in two of these. (2) Only one other source in our spectroscopic sample exhibits marginal evidence for line emission, but predicted and observed CV numbers still agree to within a factor of about 2-3. (3) We have discovered a hot (T_eff = 8700 K), low-mass (M = 0.05 M_sun) secondary star in a previously known 0.8 day binary system. This exotic object is probably the remnant of a subgiant that has been stripped of its envelope and may represent the ``smoking gun'' of a recent dynamical encounter. (4) We have found a Helium WD, the second to be optically detected in 47 Tuc, and the first outside a millisecond-pulsar system. (5) We have discovered a BS-WD binary system, the first known in any globular cluster. (6) We have found two additional candidate WD binary systems with putative main sequence and subgiant companions. (7) We estimate the WD binary fraction in the core of 47 Tuc to be 15 +17/-9 (stat) +8/-7 (sys). (8) One BS in our sample may exceed twice the cluster turn-off mass, but the uncertainties are large. Taken as a whole, our study illustrates the wide range of stellar exotica that are lurking in the cores of GCs, most of which are likely to have undergone significant dynamical encounters. [abridged]

We are conducting a spectroscopic monitoring of O and WN -type stars in our Galaxy. In this work, we summarize some of our first results related to the search for radial-velocity variations indicative of orbital motion.

In this study, we test the theoretically predicted mass-loss behaviour as a function of stellar effective temperature across the so-called `bi-stability' jump (BSJ). We gathered radio observations of 30 OB supergiants (O8-B3). We derived the radio mass-loss rates and wind efficiencies, and compared our results with Halpha mass-loss rates and predictions based on radiation-driven wind models. he wind efficiency shows the possible presence of a local maximum around an effective temperature of 21~000 K -- in qualitative agreement with predictions. We also find that the radio mass-loss rates show good agreement with empirical Halpha rates. However, the empirical mass-loss rates are larger than the predicted rates from radiation-driven wind theory for objects above the BSJ temperature, whilst they are smaller for the rest. A new wind momenta-luminosity relation for O8-B0 stars has been derived.

A detailed study of the Hourglass Nebula in the M8 star forming region is presented. The study is mainly based on recent subarcsec-resolution JHKs images taken at Las Campanas Observatory and complemented with archival HST images and longslit spectroscopy retrieved from the ESO Archive Facility. Using the new numerical code CHORIZOS, we estimate the distance to the earliest stars in the region to be 1.25 kpc. Infrared photometry of all the sources detected in the field is given. From analysis of the JHKs colour-colour diagrams, we find that an important fraction of these sources exhibit significant infrared excess. These objects are candidates to be low- and intermediate-mass pre-main sequence stars. Based on HST observations, the spatial distribution of gas, dust and stars in the region is analyzed. The morphological analysis of these images also reveals a rich variety of structures related to star formation (proplyds, jets, bow shocks), similar to those observed in M16 and M42, along with the detection of the first four Herbig-Haro objects in the region. Furthermore, a longslit spectrum obtained with NTT confirms the identification of one of them (HH 870) in the core of the Hourglass nebula, providing the first direct evidence of active star formation by accretion in M8.