The intriguing binary LS V +22 25 (LB-1) has drawn much attention following claims of it being a single-lined spectroscopic binary with a 79-day orbit comprising a B-type star and a ~70 Msun black hole. Recent analyses have implied that the visible primary star is a stripped He-rich star. However, the nature of the secondary, which was proposed to be a black hole, a neutron star, or a main sequence star, remains unknown. Based on 26 newly acquired spectroscopic observations, we perform an orbital analysis and spectral disentangling of LB-1 to elucidate the nature of the system. Our analysis reveals that LB-1 contains two components of comparable brightness in the optical. The narrow-lined primary, which we estimate to contribute ~55% in the optical, has spectral properties that suggest that it is a stripped star: it has a small spectroscopic mass (~1 Msun) for a B-type star and it is He- and N-rich. The "hidden" secondary, which contributes about 45% of the optical flux, is a rapidly rotating (vsini ~ 300 km/s) B3 V star with a decretion disk -- a Be star. Hence, LB-1 does not contain a compact object. Instead, it is a rare Be binary system consisting of a stripped star (the former mass donor) and a Be star rotating at near its critical velocity (the former mass accretor). This system is a clear example that binary interactions play a decisive role in the production of rapid stellar rotators and Be stars.

The extreme ultraviolet (EUV) spectra of distant star-forming (SF) regions cannot be probed directly using either ground- or space-based telescopes due to the high cross-section for interaction of EUV photons with the interstellar medium. The mm/submm recombination lines of H and He, which can be observed from the ground can serve as a reliable probe of the EUV. These lines are weekly impacted by dust absorption, their fluxes are independent of metallicity, and depend linearly on the number of EUV continuum photons. We present a study based on ALMA observations of three Galactic ultra-compact HII regions and the starburst region Sgr B2(M), in which we reconstruct the key parameters of the EUV spectra of these sources using mm recombination lines of HI, HeI and HeII. We find that the EUV spectra between 13.6 and 54.4eV of the objects studied here have similar frequency dependence: L_nu~nu^{-4.5 \pm 0.4}. We compare the inferred values of the EUV spectral slopes with the values expected for a single stellar evolution model (Starburst99) and the Binary Population and Spectral Synthesis code (BPASS v2.1). We find that the observed spectral slope gamma~4.5 \pm 0.4 differs from the model predictions. For an instantaneous starburst model with age <5Myrs the observed spectral slopes are consistently steeper than predictions of either of the models. For an instantaneous starburst with age >5Myrs the observed spectral slope is in-between the values predicted by the two models. This may imply that the fraction of interacting binaries in HII regions is substantially lower than assumed in BPASS v2.1. The technique demonstrated here allows one to deduce the EUV spectra of SF regions providing critical insight into photon production rates at lambda>912A and can serve as calibration to starburst synthesis models, improving our understanding of SF and the properties of ionizing flux during reionization.

It is now well established that the majority of massive stars reside in multiple systems. However, the effect of multiplicity is not sufficiently understood, resulting in a plethora of uncertainties about the end stages of massive-star evolution. In order to investigate these uncertainties, it is useful to study massive stars just before their demise. Classical Wolf-Rayet stars represent the final end stages of stars at the upper-mass end. The multiplicity fraction of these stars was reported to be ${\sim}0.4$ in the Galaxy but no correction for observational biases has been attempted. The aim of this study is to conduct a homogeneous radial-velocity survey of a magnitude-limited ($V$ $\leq 12$) sample of Galactic Wolf-Rayet stars to derive their bias-corrected multiplicity properties. The present paper focuses on 12 northern Galactic carbon-rich (WC) Wolf-Rayet stars observable with the 1.2m Mercator telescope on the island of La Palma. We homogeneously measured relative radial velocities (RVs) for carbon-rich Wolf-Rayet stars using cross-correlation. Variations in the derived RVs were used to flag binary candidates. We investigated probable orbital configurations and provide a first correction of observational biases through Monte-Carlo simulations. Of the 12 northern Galactic WC stars in our sample, seven show peak-to-peak RV variations larger than 10 km s$^{-1}$, which we adopt as our detection threshold. This results in an observed spectroscopic multiplicity fraction of 0.58 with a binomial error of 0.14. In our campaign, we find a clear lack of short-period (P~$<~\sim$100\,d), indicating that a large number of Galactic WC binaries likely reside in long-period systems. Finally, our simulations show that at the 10% significance level, the intrinsic multiplicity fraction of the Galactic WC population is at least 0.72.

About 20% of all B-type stars are classical Be stars. The Be phenomenon is strongly correlated with rapid rotation, the origin of which remains unclear. It may be rooted in single- or binary-star evolution. In the framework of the binary channel, the initially more massive star transfers mass and angular momentum to the original secondary, which becomes a Be star. The system then evolves into a Be binary with a post-main-sequence companion, which may later be disrupted in a supernova event. Hence, if the binary channel dominates the formation of Be stars, one may expect a strong lack of close Be binaries with main sequence (MS) companions. Through an extensive, star-by-star review of the literature of a magnitude-limited sample of Galactic early-type Be stars, we investigate whether Be binaries with MS companions are known to exist. Our sample is constructed from the BeSS database and cross-matched with all available literature on the individual stars. Out of an initial list of 505 Be stars, we compile a final sample of 287 Galactic Be stars earlier than B1.5 with V<=12 mag. Out of those, 13 objects were reported as Be binaries with known post-MS companions and 11 as binaries with unknown, uncertain or debated companions. We find no confirmed reports of Be binaries with MS companions. For the remaining 263 targets, no significant reports of multiplicity exist in the literature, implying that they are either Be binaries with faint companions, or truly single. The clear lack of reported MS companions to Be stars, which stands in contrast to the high number of detected B+B MS binaries, strongly supports the hypothesis that early-type Be stars are binary interaction products that spun up after mass and angular momentum transfer from a companion star. Taken at face value, our results may suggest that a large majority of the early-type Be stars have formed through binary mass-transfer.

HR 6819 was recently proposed to be a triple system consisting of an inner B-type giant + black hole binary with an orbital period of 40d and an outer Be tertiary. This interpretation is mainly based on two inferences: that the emission attributed to the outer Be star is stationary, and that the inner star, which is used as mass calibrator for the black hole, is a B-type giant. We re-investigate the properties of HR 6819 by spectral disentangling and an atmosphere analysis of the disentangled spectra to search for a possibly simpler alternative explanation for HR 6819. Disentangling implies that the Be component is not a static tertiary, but rather a component of the binary in the 40-d orbit. The inferred radial velocity amplitudes imply an extreme mass ratio of M_2/M_1 = 15 +/- 3. We infer spectroscopic masses of 0.4$^{+0.3}_{-0.1}$ Msun and 6$^{+5}_{-3}$ Msun for the primary and secondary, which agree well with the dynamical masses for an inclination of i = 32 deg. This indicates that the primary might be a stripped star rather than a B-type giant. Evolutionary modelling suggests that a possible progenitor system would be a tight (P_i ~ 2d) B+B binary system that experienced conservative mass transfer. While the observed nitrogen enrichment of the primary conforms with the predictions of the evolutionary models, we find no indications for the predicted He enrichment. We suggest that HR 6819 is a binary system consisting of a stripped B-type primary and a rapidly-rotating Be star that formed from a previous mass-transfer event. In the framework of this interpretation, HR 6819 does not contain a black hole. Interferometry can distinguish between these two scenarios by providing an independent measurement of the separation between the visible components.

Massive stars like company. However, low-mass companions have remained extremely difficult to detect at angular separations ($\rho$) smaller than 1" (~1000-3000 au considering typical distance to nearby massive stars) given the large brightness contrast between the companion and the central star. Constraints on the low-mass end of the companions mass-function for massive stars are however needed, for example to help distinguishing between various scenarios for the formation of massive stars. To obtain statistically significant constraint on the presence of low-mass companions beyond the typical detection limit of current surveys ($\Delta \mathrm{mag} \lesssim 5$ at $\rho \lesssim 1$"), we initiated a survey of O and Wolf-Rayet stars in the Carina region using the SPHERE coronagraphic instrument on the VLT. In this first paper, we aim to introduce the survey, to present the methodology and to demonstrate the capability of SPHERE for massive stars using the multiple system QZ~Car. High-contrast imaging techniques, such as angular- and spectral-differential imaging techniques as well as PSF-fitting, were applied to detect and measure the relative flux of companions in each spectral channel of the instrument. We detected 19 sources around the QZ~Car system with detection limits of 9~mag at $\rho > 200$~mas for IFS and as faint as 13~mag at $\rho$ > 1".8 for IRDIS (corresponding to sub-solar masses for potential companions). All but two are reported here for the first time. Based on this proof of concept, we showed that VLT/SPHERE allows us to reach the sub-solar mass regime of the companion mass function. This paves the way for this type of observation with a large sample of massive stars to provide novel constraints on the multiplicity of massive stars in a region of the parameter space that has remained inaccessible so far.

We conducted an observational campaign towards one of the most massive and luminous colliding wind binaries in the Galaxy, HD~93129A, close to its periastron passage in 2018. During this time the source was predicted to be in its maximum of high-energy emission. Here we present our data analysis from the X-ray satellites \textit{Chandra} and \textit{NuSTAR} and the $\gamma$-ray satellite \textit{AGILE}. High-energy emission coincident with HD~93129A was detected in the X-ray band up to $\sim$18~keV, whereas in the $\gamma$-ray band only upper limits were obtained. We interpret the derived fluxes using a non-thermal radiative model for the wind-collision region. We establish a conservative upper limit for the fraction of the wind kinetic power that is converted into relativistic electron acceleration, $f_\mathrm{NT,e} < 0.02$. In addition, we set a lower limit for the magnetic field in the wind-collision region as $B_\mathrm{WCR} > 0.3$~G. We also argue a putative interpretation of the emission from which we estimate $f_\mathrm{NT,e} \approx 0.006$ and $B_\mathrm{WCR} \approx 0.5$~G. We conclude that multi-wavelength, dedicated observing campaigns during carefully selected epochs are a powerful tool for characterising the relativistic particle content and magnetic field intensity in colliding wind binaries.

The recent gravitational wave measurements have demonstrated the existence of stellar mass black hole binaries. It is essential for our understanding of massive star evolution to identify the contribution of binary evolution to the formation of double black holes. A promising way to progress is investigating the progenitors of double black hole systems and comparing predictions with local massive star samples such as the population in 30 Doradus in the Large Magellanic Cloud (LMC). Methods. To this purpose, we analyse a large grid of detailed binary evolution models at LMC metallicity with initial primary masses between 10 and 40 Msun, and identify which model systems potentially evolve into a binary consisting of a black hole and a massive main sequence star. We then derive the observable properties of such systems, as well as peculiarities of the OB star component. We find that about 3% of the LMC late O and early B stars in binaries are expected to possess a black hole companion, when assuming stars with a final helium core mass above 6.6 M to form black holes. While the vast majority of them may be X-ray quiet, our models suggest that these may be identified in spectroscopic binaries, either by large amplitude radial velocity variations ( > 50 km s ) and simultaneous nitrogen surface enrichment, or through a moderate radial velocity ( > 10 km/s ) and simultaneously rapid rotation of the OB star. The predicted mass ratios are such that main sequence companions could be excluded in most cases. A comparison to the observed OB+WR binaries in the LMC, Be/X-ray binaries, and known massive BH binaries supports our conclusion. We expect spectroscopic observations to be able to test key assumptions in our models, with important implications for massive star evolution in general, and for the formation of double-black hole mergers in particular.

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

The process of massive star ($M\geq8~M_\odot$) formation is still poorly understood. Observations of massive young stellar objects (MYSOs) are challenging due to their rarity, short formation timescale, large distances, and high circumstellar extinction. Here, we present the results of a spectroscopic analysis of a population of MYSOs in the Large Magellanic Cloud (LMC). We took advantage of the spectral resolution and wavelength coverage of X-shooter (300-2500 nm), mounted on the European Southern Observatory Very Large Telescope, to detect characteristic spectral features in a dozen MYSO candidates near 30 Doradus, the largest starburst region in the Local Group hosting the most massive stars known. The X-shooter spectra are strongly contaminated by nebular emission. We used a scaling method to subtract the nebular contamination from our objects. We detect H$\alpha,\beta$, [O I] 630.0 nm, Ca II infrared triplet, [Fe II] 1643.5 nm, fluorescent Fe II 1687.8 nm, H$_2$ 2121.8 nm, Br$\gamma$, and CO bandhead emission in the spectra of multiple candidates. This leads to the spectroscopic confirmation of 10 candidates as bona fide MYSOs. We compare our observations with photometric observations from the literature and find all MYSOs to have a strong near-infrared excess. We compute lower limits to the brightness and luminosity of the MYSO candidates, confirming the near-infrared excess and the massive nature of the objects. No clear correlation is seen between the Br$\gamma$ luminosity and metallicity. Combining our sample with other LMC samples results in a combined detection rate of disk features such as fluorescent Fe II and CO bandheads which is consistent with the Galactic rate (40\%). Most of our MYSOs show outflow features.

Classical Wolf-Rayet (WR) stars are massive, hydrogen depleted, post main-sequence stars that exhibit emission-line dominated spectra. For a given metallicity Z, stars exceeding a certain initial mass M_single(Z) can reach the WR phase through intrinsic mass-loss (single-star channel). Stars of lower masses can reach the WR phase via binary mass transfer (binary channel). It is commonly assumed that the binary channel dominates the formation of WR stars in environments with low Z such as the SMC and LMC. However, their reported WR binary fractions of 30-40% are comparable to that of the Galaxy. Here, we explain this apparent contradiction by considering the minimum initial mass M_spec(Z) needed for the stripped product to appear as a WR star. We calibrate M_spec(Z) using the lowest-luminosity WR stars in the Clouds and the Galaxy. A range of M_single(Z) values are explored using various evolution codes. We estimate the additional contribution of the binary channel by considering the interval [M_spec(Z), M_single(Z)], which characterises the initial-mass range in which binaries can form additional WR stars. Results: The WR-phenomenon ceases below luminosities of logL = 4.9, 5.25, and 5.6 [Lsun] in the Galaxy, LMC, and SMC, which translates to He-star masses of 7.5, 11, 17 Msun and initial masses of M_spec = 18, 23, 37 Msun. Stripped stars with lower initial masses in the respective galaxies would tend to not appear as WR stars. M_single lies in the range 20-30, 30-60, and > 40 Msun for the Galaxy, LMC, and SMC. We find that that the additional contribution of the binary channel is a non-trivial function of Z that cannot be conclusively claimed to be monotonically increasing with decreasing Z. Hence, one should not a-priori expect that binary interactions become increasingly important in forming WR stars at low Z, or that the WR binary fraction grows with decreasing Z.

The multiplicity properties of massive stars are one of the important outstanding issues in stellar evolution. Quantifying the binary statistics of all evolutionary phases is essential to paint a complete picture of how and when massive stars interact with their companions, and to determine the consequences of these interactions. We investigate the multiplicity of an almost complete census of red supergiant stars (RSGs) in NGC 330, a young massive cluster in the SMC. Using a combination of multi-epoch HARPS and MUSE spectroscopy, we estimate radial velocities and assess the kinematic and multiplicity properties of 15 RSGs in NGC 330. Radial velocities are estimated to better than +/-100 m/s for the HARPS data. The line-of-sight velocity dispersion for the cluster is estimated as 3.20 +0.69-0.52 km/s. When virial equilibrium is assumed, the dynamical mass of the cluster is log (M{dyn} /M{sun}) = 5.20+/-0.17, in good agreement with previous upper limits. We detect significant radial velocity variability in our multi-epoch observations and distinguish between variations caused by atmospheric activity and those caused by binarity. The binary fraction of NGC 330 RSGs is estimated by comparisons with simulated observations of systems with a range of input binary fractions. In this way, we account for observational biases and estimate the intrinsic binary fraction for RSGs in NGC 330 as f{RSG} = 0.3+/-0.1 for orbital periods in the range 2.3< log P [days] <4.3, with q>0.1. Using the distribution of the luminosities of the RSG population, we estimate the age of NGC 330 to be 45+/-5 Myr and estimate a red straggler fraction of 50%. We estimate the binary fraction of RSGs in NGC 330 and conclude that it appears to be lower than that of main-sequence massive stars, which is expected because interactions between an RSG and a companion are assumed to effectively strip the RSG envelope.

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

A high fraction of massive stars are found to be binaries but only a few of them are reported as photometrically variable. By studying the populations of SB2 in the 30 Doradus region, we found a subset of them that have photometry from the OGLE project and that display variations in their light curves related to orbital motions. The goal of this study is to determine the dynamical masses and radii of the 26 binary components to investigate the mass-discrepancy problem and to provide an empirical mass-luminosity relation for the LMC. We use the PHOEBE programme to perform a systematic analysis of the OGLE V and I light curves obtained for 13 binary systems in 30 Dor. We adopt Teff, and orbital parameters derived previously to obtain the inclinations of the systems and the parameters of the individual components. Three systems display eclipses in their light curves, while the others only display ellipsoidal variations. We classify two systems as over-contact, five as semi-detached, and four as detached. The two remaining systems have uncertain configurations due to large uncertainties on their inclinations. The fact that systems display ellipsoidal variations has a significant impact on the inclination errors. From the dynamical masses, luminosities, and radii, we provide LMC-based empirical mass-luminosity and mass-radius relations, and we compare them to other relations given for the Galaxy, the LMC, and the SMC. These relations differ for different mass ranges, but do not seem to depend on the metallicity regimes. We also compare the dynamical, spectroscopic, and evolutionary masses of the stars in our sample. While the dynamical and spectroscopic masses agree with each other, the evolutionary masses are systematically higher, at least for stars in semi-detached systems. This suggests that the mass discrepancy can be partly explained by past or ongoing interactions between the stars.

Accurate stellar parameters of individual objects in binary systems are essential to constrain the effects of binarity on stellar evolution. These parameters serve as a prerequisite to probing existing and future theoretical evolutionary models. We aim to derive the atmospheric parameters of the 31 SB2s in the TMBM sample. This sample, composed of detached, semi-detached and contact systems with at least one of the components classified as an O star, is an excellent test-bed to study how binarity can impact our knowledge of the evolution of massive stars. 32 epochs of FLAMES/GIRAFFE spectra are analysed using spectral disentangling to construct the individual spectra of 62 components. We apply the CMFGEN atmosphere code to determine their stellar parameters and their He, C and N surface abundances. From these properties, we show that the effects of tides on chemical mixing are limited. Components on longer-period orbits show higher nitrogen enrichment at their surface than those on shorter-period orbits, in contrast to expectations of rotational or tidal mixing, implying that other mechanisms play a role in this process. Components filling their Roche lobe are mass donors. They exhibit higher nitrogen content at their surface and rotate more slowly than their companions. By accreting new material, their companions spin faster and are rejuvenated. Their locations in the N-vsini diagram tend to show that binary products are good candidates to populate the two groups of stars (slowly rotating, nitrogen-enriched and rapidly rotating non-enriched) that cannot be reproduced through single-star population synthesis. This sample is the largest sample of binaries to be studied in such a homogeneous way. The study of these objects gives us strong observational constraints to test theoretical binary evolutionary tracks.

Liu et al. (2019) recently reported the detection of a 68 [+11/-13] solar mass (Msun) black hole (BH) paired with an 8.2 [+0.9/-1.2] Msun B-type sub-giant star in the 78.9-day spectroscopic binary system LB-1. Such a black hole is over twice as massive as any other known stellar-mass black hole with non-compact companions2 and its mass approaches those that result from BH-BH coalescences that are detected by gravitational wave interferometers. Its presence in a solar-like metallicity environment challenges conventional theories of massive binary evolution, stellar winds and core-collapse supernovae, so that more exotic scenarios seem to be needed to explain the existence and properties of LB-1. Here, we show that the observational diagnostics used to derive the BH mass results from the orbital motion of the B-type star, not that of the BH. As a consequence, no evidence for a massive BH remains in the data, therefore solving the existing tension with formation models of such a massive BH at solar metallicity and with theories of massive star evolution in general.

Rapid rotation in massive stars imposes a latitudinal variation in the mass loss from radiatively driven winds that can lead to enhanced mass loss at the poles (with little angular momentum loss) and/or equator (with maximal angular momentum loss). Here we present an examination of the stellar wind lines of the two O-type stars with the fastest known equatorial velocities, VFTS 102 ($V\sin i = 610 \pm 30$ km/s; O9:Vnnne+) and VFTS 285 ($V\sin i = 609 \pm 29$ km/s; O7.5 Vnnn) in the Large Magellanic Cloud. Ultraviolet spectra of both stars were obtained with the Hubble Space Telescope Cosmic Origins Spectrograph. The spectrum of VFTS 285 displays a fast outflow in N V and a much slower wind in Si IV, and we argue that there is a two-wind regime in which mass loss is strong at the poles (fast and tenuous wind) but dominant at the equator (slow and dense winds). These ions and wind lines are not present in the spectrum of the cooler star VFTS 102, but the double-peaked H$\alpha$ emission in its spectrum implies equatorial mass loss into a circumstellar disk. The results suggest that in the fastest rotating O-stars, most mass is lost as an equatorial outflow, promoting angular momentum loss that contributes to a spin down over time.

Context: The identification and characterisation of populations of young massive stars in (giant) HII regions provides important constraints on i) the formation process of massive stars and their early feedback on the environment, and ii) the initial conditions for population synthesis models predicting the evolution of ensembles of stars. Aims: We identify and characterise the stellar populations of the following young giant HII regions: M8, G333.6-0.2, and NGC6357. Methods: We acquired H- and K-band spectra of around 200 stars using The K-band KMOS on the ESO Very Large Telescope. The targets for M8 and NGC6357 were selected from the MYStIX project, which combines X-ray observations with near-infrared and mid-infrared data. For G333.6-0.2, the sample selection is based on the near-infrared colours combined with X-ray data. We introduce an automatic spectral classification method in order to obtain temperatures and luminosities for the observed stars. We analyse the stellar populations using their photometric, astrometric, and spectroscopic properties and compared the position of the stars in the Hertzprung-Russell diagram with stellar evolution models to constrain their ages and mass ranges. Results: We confirm the presence of candidate ionising sources in the three regions and report new ones, including the first spectroscopically identified O stars in G333.6-0.2. In M8 and NGC6357, two populations are identified: i) OB main-sequence stars ($M > 5~\rm{M_{\odot}}$) and ii) pre-main sequence stars ($M\approx0.5-5~\rm{M_{\odot}}$). The ages of the clusters are $\sim$1-3~Myr, $< 3$~Myr, and $\sim$0.5-3~Myr for M8, G333.6-0.2, and NGC6357, respectively. We show that MYStIX selected targets have $>$ 90\% probability of being members of the HII region, whereas a selection based on near infrared (NIR) colours leads to a membership probability of only $\sim$70\%.

PolStar is an Explorer-class far ultraviolet (FUV) spectropolarimetry mission designed to target massive stars and their environments. PolStar will take advantage of resonance lines only available in the FUV to measure for the first time the magnetic and wind environment around massive stars to constrain models of rotation and mass loss.

A majority of massive stars are part of binary systems, a large fraction of which will inevitably interact during their lives. Binary-interaction products (BiPs), i.e. stars affected by such interaction, are expected to be commonly present in stellar populations. BiPs are thus a crucial ingredient in the understanding of stellar evolution. We aim to identify and characterize a statistically significant sample of BiPs by studying clusters of 10-40 Myr, an age at which binary population models predict the abundance of BiPs to be highest. One example of such a cluster is NGC 330 in the Small Magellanic Cloud. Using MUSE WFM-AO observations of NGC 330, we resolve the dense cluster core for the first time and are able to extract spectra of its entire massive star population. We develop an automated spectral classification scheme based on the equivalent widths of spectral lines in the red part of the spectrum. We characterize the massive star content of the core of NGC 330 which contains more than 200 B stars, 2 O stars, 6 A-type supergiants and 11 red supergiants. We find a lower limit on the Be star fraction of $32 \pm 3\%$ in the whole sample. It increases to at least $46 \pm 10\%$ when only considering stars brighter than V=17 mag. We estimate an age of the cluster core between 35 and 40 Myr and a total cluster mass of $88^{+17}_{-18} \times 10^3 M_{\odot}$. We find that the population in the cluster core is different than the population in the outskirts: while the stellar content in the core appears to be older than the stars in the outskirts, the Be star fraction and the observed binary fraction are significantly higher. Furthermore, we detect several BiP candidates that will be subject of future studies.

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.

Massive Wolf-Rayet (WR) stars dominate the radiative and mechanical energy budget of galaxies and probe a critical phase in the evolution of massive stars prior to core-collapse. It is not known whether core He-burning WR stars (classical WR, cWR) form predominantly through wind-stripping (w-WR) or binary stripping (b-WR). With spectroscopy of WR binaries so-far largely avoided due to its complexity, our study focuses on the 44 WR binaries / binary candidates of the Large Magellanic Cloud (LMC, metallicity Z~0.5 Zsun), identified on the basis of radial velocity variations, composite spectra, or high X-ray luminosities. Relying on a diverse spectroscopic database, we aim to derive the physical and orbital parameters of our targets, confronting evolution models of evolved massive stars at sub-solar metallicity, and constraining the impact of binary interaction in forming them. Spectroscopy is performed using the Potsdam Wolf-Rayet (PoWR) code and cross-correlation techniques. Disentanglement is performed using the code Spectangular or the shift-and-add algorithm. Evolutionary status is interpreted using the Binary Population and Spectral Synthesis (BPASS) code, exploring binary interaction and chemically-homogeneous evolution. No obvious dichotomy in the locations of apparently-single and binary WN stars on the Hertzsprung-Russell diagram is apparent. According to commonly used stellar evolution models (BPASS, Geneva), most apparently-single WN stars could not have formed as single stars, implying that they were stripped by an undetected companion. Otherwise, it must follow that pre-WR mass-loss/mixing (e.g., during the red supergiant phase) are strongly underestimated in standard stellar evolution models.

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

Anticipating the kinematic constraints from the Gaia mission, we perform an extensive numerical study of the evolution of massive binary systems to predict the peculiar velocities that stars obtain when their companion collapses and disrupts the system. Our aim is to (1) identify which predictions are robust against model uncertainties and assess their implications, (2) investigate which physical processes leave a clear imprint and may therefore be constrained observationally and (3) provide a suite of publicly available model predictions. We find that $22_{-8}^{+26}$% of all massive binary systems merge prior to the first core collapse in the system. Of the remainder, $86_{-9}^{+11}$% become unbound because of the core-collapse. Remarkably, this rarely produce runaway stars (i.e., stars with velocities above 30 km/s). These are outnumbered by more than an order of magnitude by slower unbound companions, or "walkaway stars". This is a robust outcome of our simulations and is due to the reversal of the mass ratio prior to the explosion and widening of the orbit, as we show analytically and numerically. We estimate a $10^{+5}_{-8}$% of massive stars to be walkaways and only $0.5^{+1.0}_{-0.4}$% to be runaways, nearly all of which have accreted mass from their companion. Our findings are consistent with earlier studies, however the low runaway fraction we find is in tension with observed fractions 10%. If Gaia confirms these high fractions of massive runaway stars resulting from binaries, it would imply that we are currently missing physics in the binary models. Finally, we show that high end of the mass distributions of runaway stars is very sensitive to the assumed black hole natal kicks and propose this as a potentially stringent test for the explosion mechanism. We discuss companions remaining bound which can evolve into X-ray and gravitational wave sources.

The incidence of multiplicity in cool, luminous massive stars is relatively unknown compared to their hotter counterparts. Here we present radial velocity (RV) measurements and investigate the multiplicity properties of red supergiants (RSGs) in the 30~Doradus region of the Large Magellanic Cloud. We provide absolute RV measurements for our sample and estimate line-of-sight velocities for the Hodge 301 and SL 639 clusters, which agree well with those of hot stars in the same clusters. By combining results for the RSGs with those for nearby B-type stars, we estimate systemic velocities and velocity dispersions for the two clusters, obtaining estimates for their dynamical masses of $\log (M_{\rm dyn}/M_{\odot})=$3.8$\pm$0.3 for Hodge 301, and an upper limit of $\log (M_{\rm dyn}/M_{\odot})$<3.1$\pm$0.8 for SL 639, assuming Virial equilibrium. Analysis of the multi-epoch data reveals one RV-variable, potential binary candidate (VFTS744), which is likely a semi-regular variable asymptotic giant branch star. We estimate an upper limit on the observed binary fraction for our sample of 0.3, where we are sensitive to maximum periods for individual objects in the range of 1 to 10 000 days and mass-ratios above 0.3 depending on the data quality. From simulations of the RV measurements from binary systems given the current data we conclude that systems within the parameter range q>0.3, $\log$P[days]<3.5, would be detected by our variability criteria, at the 90% confidence level. The intrinsic binary fraction, accounting for observational biases, is estimated using simulations of binary systems with an empirically defined distribution of parameters where orbital periods are uniformly distributed in the 3.3<$\log$P[days]<4.3 range. A range of intrinsic binary fractions are considered; a binary fraction of 0.3 is found to best reproduce the observed data. [Abridged]

We estimate physical parameters for the late-type massive stars observed as part of the VLT-FLAMES Tarantula Survey (VFTS) in the 30 Doradus region of the Large Magellanic Cloud (LMC). The observational sample comprises 20 candidate red supergiants (RSGs) which are the reddest (($B-V$) $>$ 1 mag) and brightest ($V$ $<$ 16 mag) objects in the VFTS. We use optical and near-IR photometry to estimate their temperatures and luminosities, and introduce the luminosity-age diagram to estimate their ages. We derive physical parameters for our targets, including temperatures from a new calibration of $(J-K_{\rm s})_{0}$ colour for luminous cool stars in the LMC, luminosities from their $J$-band magnitudes (thence radii), and ages from comparisons with state-of-the-art evolutionary models. We show that interstellar extinction is a significant factor for our targets, highlighting the need to take it into account in analysis of the physical parameters of RSGs. We find that some of the candidate RSGs could be massive AGB stars. The apparent ages of the RSGs in the Hodge 301 and SL 639 clusters show a significant spread (12-24 Myr). We also apply our approach to the RSG population of the relatively nearby NGC 2100 cluster, finding a similarly large spread. We argue that the effects of mass-transfer in binaries may lead to more massive and luminous RSGs (which we call `red stragglers') than expected from single-star evolution, and that the true cluster ages correspond to the upper limit of the estimated RSG ages. In this way, the RSGs can serve as a new and potentially reliable age tracer in young star clusters. The corresponding analysis yields ages of 24$^{+5}_{-3}$ Myr for Hodge 301, 22$^{+6}_{-5}$ Myr for SL 639, and 23$^{+4}_{-2}$ Myr for NGC 2100.

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]

CubeSats are routinely used for low-cost photometry from space. Space-borne spectroscopy, however, is still the exclusive domain of much larger platforms. Key astrophysical questions in e.g. stellar physics and exoplanet research require uninterrupted spectral monitoring from space over weeks or months. Such monitoring of individual sources is unfortunately not affordable with these large platforms. With CUBESPEC we plan to offer the astronomical community a low-cost CubeSat solution for near-UV/optical/near-IR spectroscopy that enables this type of observations. CUBESPEC is a generic spectrograph that can be configured with minimal hardware changes to deliver both low resolution (R=100) with very large spectral coverage (200-1000nm), as well as high resolution (R=30,000) over a selected wavelength range. It is built around an off-axis Cassegrain telescope and a slit spectrograph with configurable dispersion elements. CUBESPEC will use a compact attitude determination and control system for coarse pointing of the entire spacecraft, supplemented with a fine-guidance system using a fast steering mirror to center the source on the spectrograph slit and to cancel out satellite jitter. An extremely compact optical design allows us to house this instrument in a 6U CubeSat with a volume of only 10x20x30cm$^{3}$, while preserving a maximized entrance pupil of ca. 9x19cm$^{2}$. In this contribution, we give an overview of the CUBESPEC project, discuss its most relevant science cases, and present the design of the instrument.

SMC AB 6 is the shortest-period (6.5d) Wolf-Rayet (WR) binary in the Small Magellanic Cloud, and is therefore crucial for the study of binary interaction and formation of WR stars at low metallicity. The WR component in AB 6 was previously found to be very luminous (logL=6.3[Lsun]) compared to its reported orbital mass (8Msun), placing it significantly above the Eddington limit. Through spectroscopy and orbital analysis of newly acquired optical data taken with UVES, we aim to understand the peculiar results reported for this system and explore its evolutionary history. Results: We find that AB 6 contains at least four stars. The 6.5d period WR binary comprises the WR primary (WN3:h, star A) and a rather rapidly rotating early O-type companion (O5.5 V, star B). Static N and He lines suggest the presence of an emission line star (O5.5 I(f), star C). Finally, narrow absorption lines portraying a long-term radial velocity variation show the existence of a fourth star (O7.5 V, star D). Star D appears to form a second 140d period binary together with a fifth stellar member, which is a B-type dwarf or a black hole. It is not clear that these additional components are bound to the WR binary. The WR star is found to be less luminous than previously thought (logL = 5.9[Lsun]) and, adopting 41Msun for star B, more massive (18Msun). Correspondingly, the WR star does not exceed the Eddington limit. We derive the initial masses of 60 and 40Msun for stars A and B and an age of 3.9 Myr for the system. The WR binary likely experienced nonconservative mass transfer in the past supported by the relatively rapid rotation of star B. Conclusion: Our study shows that AB 6 is a multiple -- probably quintuple -- system. This finding resolves the previously reported puzzle of the WR primary exceeding the Eddington limit and suggests that the WR star exchanged mass with its companion in the past.

We present a catalog of relative proper motions for 368,787 stars in the 30 Doradus region of the Large Magellanic Cloud (LMC), based on a dedicated two-epoch survey with the Hubble Space Telescope (HST) and supplemented with proper motions from our pilot archival study. We demonstrate that a relatively short epoch difference of 3 years is sufficient to reach a $\sim$0.1 mas yr$^{-1}$ level of precision or better. A number of stars have relative proper motions exceeding a 3-sigma error threshold, representing a mixture of Milky Way denizens and 17 potential LMC runaway stars. Based upon 183 VFTS OB-stars with the best proper motions, we conclude that none of them move faster than $\sim$0.3 mas yr$^{-1}$ in each coordinate -- equivalent to $\sim$70 km s$^{-1}$. Among the remaining 351 VFTS stars with less accurate proper motions, only one candidate OB runaway can be identified. We rule out any OB star in our sample moving at a tangential velocity exceeding $\sim$120 km s$^{-1}$. The most significant result of this study is finding 10 stars over wide range of masses, which appear to be ejected from the massive star cluster R136 in the tangential plane to angular distances from $35^{\prime\prime}$ out to $407^{\prime\prime}$, equivalent to 8-98 pc. The tangential velocities of these runaways appear to be correlated with apparent magnitude, indicating a possible dependence on the stellar mass.

Previous analyses of the spectra of OB-type stars in the Magellanic Clouds have identified targets with low projected rotational velocities and relatively high nitrogen abundances. The VFTS obtained spectroscopy for B-type 434 stars. We have estimated atmospheric parameters and nitrogen abundances using TLUSTY model atmospheres for 54 B-type targets that appear to be single, have projected rotational velocities <80 \kms and were not classified as supergiants. In addition, nitrogen abundances for 34 similar stars from a previous survey have been re-evaluated. Approximately 75-80% of the targets have nitrogen enhancements of less than 0.3 dex, consistent with them having experienced only small amounts of mixing. However, stars with low projected rotational velocities <40 \kms and significant nitrogen enrichments are found in both our samples and simulations imply that these cannot all be rapidly rotating objects observed near pole-on. Adopting an enhancement threshold of 0.6 dex, we observed five stars, yet stellar evolution models with rotation predict only 1.25$\pm$1.11 for random stellar viewing inclinations. The excess of such objects is estimated to be 20-30% of all stars with current rotational velocities of less than 40 \kms and 2-4% of the total single B-type sample. These estimates constitute lower limits for stars that appear inconsistent with current grids of stellar evolutionary models. Including targets with smaller nitrogen enhancements of greater than 0.2 dex implies larger percentages of targets that are inconsistent with current evolutionary models, viz. 70% of the stars with rotational velocities less than 40 \kms and 6-8% of the total single stellar population. We consider possible explanations of which the most promising would appear to be breaking due to magnetic fields or stellar mergers with subsequent magnetic braking.

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

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.

We present a number of notable results from the VLT-FLAMES Tarantula Survey (VFTS), an ESO Large Program during which we obtained multi-epoch medium-resolution optical spectroscopy of a very large sample of over 800 massive stars in the 30 Doradus region of the Large Magellanic Cloud (LMC). This unprecedented data-set has enabled us to address some key questions regarding atmospheres and winds, as well as the evolution of (very) massive stars. Here we focus on O-type runaways, the width of the main sequence, and the mass-loss rates for (very) massive stars. We also provide indications for the presence of a top-heavy initial mass function (IMF) in 30 Dor.

Aims. HD166734 is an eccentric eclipsing binary system composed of two supergiant O-type stars, orbiting with a 34.5-day period. In this rare configuration for such stars, the two objects mainly evolve independently, following single-star evolution so far. This system provides a chance to study the individual parameters of two supergiant massive stars and to derive their real masses. Methods. An intensive monitoring was dedicated to HD166734.We analyzed mid- and high-resolution optical spectra to constrain the orbital parameters of this system. We also studied its light curve for the first time, obtained in the VRI filters. Finally, we disentangled the spectra of the two stars and modeled them with the CMFGEN atmosphere code in order to determine the individual physical parameters. Results. HD166734 is a O7.5If+O9I(f) binary. We confirm its orbital period but we revise the other orbital parameters. In comparison to what we found in the literature, the system is more eccentric and, now, the hottest and the most luminous component is also the most massive one. The light curve exhibits only one eclipse and its analysis indicates an inclination of 63.0{\deg} $\pm$ 2.7{\deg}. The photometric analysis provides us with a good estimation of the luminosities of the stars, and therefore their exact positions in the Hertzsprung-Russell diagram. The evolutionary and the spectroscopic masses show good agreement with the dynamical masses of 39.5 Msun for the primary and 33.5 Msun for the secondary, within the uncertainties. The two components are both enriched in helium and in nitrogen and depleted in carbon. In addition, the primary also shows a depletion in oxygen. Their surface abundances are however not different from those derived from single supergiant stars, yielding, for both components, an evolution similar to that of single stars.

Massive stars like company. Here, we provide a brief overview of progresses made over the last 5 years by a number of medium and large surveys. These results provide new insights on the observed and intrinsic multiplicity properties of main sequence massive stars and on the initial conditions for their future evolution. They also bring new interesting constraints on the outcome of the massive star formation process.

The formation process of massive stars is still poorly understood. Massive young stellar objects (mYSOs) are deeply embedded in their parental clouds, they are rare and thus typically distant, and their reddened spectra usually preclude the determination of their photospheric parameters. M17 is one of the best studied HII regions in the sky, is relatively nearby, and hosts a young stellar population. With X-shooter on the ESO Very Large Telescope we have obtained optical to near-infrared spectra of candidate mYSOs, identified by Hanson et al. (1997), and a few OB stars in this region. The large wavelength coverage enables a detailed spectroscopic analysis of their photospheres and circumstellar disks. We confirm the pre-main sequence (PMS) nature of six of the stars and characterise the O stars. The PMS stars have radii consistent with being contracting towards the main sequence and are surrounded by a remnant accretion disk. The observed infrared excess and the (double-peaked) emission lines provide the opportunity to measure structured velocity profiles in the disks. We compare the observed properties of this unique sample of young massive stars with evolutionary tracks of massive protostars by Hosokawa & Omukai (2009), and propose that these mYSOs near the western edge of the HII region are on their way to become main-sequence stars ($\sim 6 - 20$ $M_{\odot}$) after having undergone high mass-accretion rates (${\dot{M}_{\rm acc}} \sim 10^{-4} - 10^{-3}$ $M_{\odot}$ $\rm yr^{-1}$). Their spin distribution upon arrival at the zero age main sequence (ZAMS) is consistent with that observed for young B stars, assuming conservation of angular momentum and homologous contraction.

A spectroscopic analysis has been undertaken for the B-type multiple systems (excluding those with supergiant primaries) in the VLT-FLAMES Tarantula Survey (VFTS). Projected rotational velocities, $v$sin$i$, for the primaries have been estimated using a Fourier Transform technique and confirmed by fitting rotationally broadened profiles. A subset of 33 systems with $v$sin$i$$\leq$ 80 km/s have been analysed using a TLUSTY grid of model atmospheres to estimate stellar parameters and surface abundances for the primaries. The effects of a potential flux contribution from an unseen secondary have also been considered. For 20 targets it was possible to reliably estimate their effective temperatures (T$_{eff}$) but for the other 13 objects it was only possible to provide a constraint of 20,000$\leq$T$_{eff}$$\leq$26,000 K -- the other parameters estimated for these targets will be consequently less reliable. The estimated stellar properties are compared with evolutionary models and are generally consistent with their membership of 30 Doradus, while the nature of the secondaries of 3 SB2 system is discussed. A comparison with a sample of single stars with $v$sin$i$$\leq$ 80 km/s obtained from the VFTS and analysed with the same techniques implies that the atmospheric parameters and nitrogen abundances of the two samples are similar. However, the binary sample may have a lack of primaries with significant nitrogen enhancements, which would be consistent with them having low rotational velocities and having effectively evolved as single stars without significant rotational mixing. This result, which may be actually a consequence of the limitations of the pathfinder investigation presented in this paper, should be considered as a motivation for spectroscopic abundance analysis of large samples of binary stars, with high quality observational data.

Massive binaries (MBs) play a crucial role in the Universe. Knowing the distributions of their orbital parameters (OPs) is important for a wide range of topics, from stellar feedback to binary evolution channels, from the distribution of supernova types to gravitational wave progenitors, yet, no direct measurements exist outside the Milky Way. The Tarantula Massive Binary Monitoring was designed to help fill this gap by obtaining multi-epoch radial velocity monitoring of 102 MBs in the 30 Dor. In this paper, we analyse 32 VLT/FLAMES observations of 93 O- and 7 B-type binaries. We performed a Fourier analysis and obtained orbital solutions for 82 systems: 51 single- and 31 double-lined spectroscopic binaries. Overall, the OPs and binary fraction are remarkably similar across the 30 Dor region and compared to existing Galactic samples (GSs). This indicates that within these domains environmental effects are of second order in shaping the properties of MBs. A small difference is found in the distribution of orbital periods (OrbPs), which is slightly flatter (in log space) in 30 Dor than in the Galaxy, although this may be compatible within error estimates and differences in the fitting methodology. Also, OrbPs in 30 Dor can be as short as 1.1 d; somewhat shorter than seen in GSs. Equal mass binaries q>0.95 in 30 Dor are all found outside NGC 2070 the very young and massive cluster at 30 Dor's core. One outstanding exception however is the fact that earliest spectral types (O2-O7) tend to have shorter OrbPs than latter (O9.2-O9.7). Our results point to a relative universality of the incidence rate of MBs and their OPs in the metallicity range from solar ($Z_{\odot}$) to about $0.5Z_{\odot}$. This provides the first direct constraints on MB properties in massive star-forming galaxies at the Universes peak of star formation at redshifts z~1 to 2, which are estimated to have $Z~0.5Z_{\odot}$.

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

The formation of massive stars remains poorly understood and little is known about their birth multiplicity properties. Here, we investigate the strikingly low radial-velocity dispersion measured for a sample of 11 massive pre- and near-main-sequence stars (sigma_rv = 5.6 +/- 0.2 km/s) in the young massive star forming region M17 to obtain first constraints on the multiplicity properties of young massive stellar objects. Methods: We compute the RV dispersion of synthetic populations of massive stars for various multiplicity properties and we compare the simulated sigma_rv distributions to the observed value. We specifically investigate two scenarios: a low binary fraction and a dearth of short-period binary systems. Results: Simulated populations with low binary fractions (f_bin = 0.12_{-0.09}^{+0.16}) or with truncated period distributions (P_cutoff > 9 months) are able to reproduce the low sigma_rv observed within their 68%-confidence intervals. Parent populations with f_bin > 0.42 or P_cutoff < 47 d can however be rejected at the 5%-significance level. Both constraints are contrast with the high binary fraction and plethora of short-period systems found in few Myr-old, OB-type populations. To explain the difference, the first scenario requires a variation of the outcome of the massive star formation process. In the the second scenario, compact binaries must form later on, and the cut-off period may be related to physical length-scales representative of the bloated pre-main-sequence stellar radii or of their accretion disks. Conclusions: If the obtained constraints are representative of the overall properties of massive young stellar objects, our results may provide support to a formation process in which binaries are initially formed at larger separations, then harden or migrate to produce the typical (untruncated) power-law period distribution observed in few Myr-old OB binaries.

We present the science cases and technological discussions that came from the workshop entitled "Finding the UV-Visible Path Forward" held at NASA GSFC June 25-26, 2015. The material presented outlines the compelling science that can be enabled by a next generation space-based observatory dedicated for UV-visible science, the technologies that are available to include in that observatory design, and the range of possible alternative launch approaches that could also enable some of the science. The recommendations to the Cosmic Origins Program Analysis Group from the workshop attendees on possible future development directions are outlined.

The "B fields in OB stars" (BOB) collaboration is based on an ESO Large Programme, to study the occurrence rate, properties, and ultimately the origin of magnetic fields in massive stars. In the framework of this programme, we carried out low-resolution spectropolarimetric observations of a large sample of massive stars using FORS2 installed at the ESO VLT 8-m telescope. We determined the magnetic field values with two completely independent reduction and analysis pipelines. Our in-depth study of the magnetic field measurements shows that differences between our two pipelines are usually well within 3sigma errors. From the 32 observations of 28 OB stars, we were able to monitor the magnetic fields in CPD-57 3509 and HD164492C, confirm the magnetic field in HD54879, and detect a magnetic field in CPD-62 2124. We obtain a magnetic field detection rate of 6+-3% for the full sample of 69 OB stars observed with FORS2 within the BOB programme. For the pre-selected objects with a v sin i below 60 km/s, we obtain a magnetic field detection rate of 5+-5%. We also discuss X-ray properties and multiplicity of the objects in our FORS2 sample with respect to the magnetic field detections.

We present the first SB2 orbital solution and disentanglement of the massive Wolf-Rayet binary R145 (P = 159d) located in the Large Magellanic Cloud. The primary was claimed to have a stellar mass greater than 300Msun, making it a candidate for the most massive star known. While the primary is a known late type, H-rich Wolf-Rayet star (WN6h), the secondary could not be so far unambiguously detected. Using moderate resolution spectra, we are able to derive accurate radial velocities for both components. By performing simultaneous orbital and polarimetric analyses, we derive the complete set of orbital parameters, including the inclination. The spectra are disentangled and spectroscopically analyzed, and an analysis of the wind-wind collision zone is conducted. The disentangled spectra and our models are consistent with a WN6h type for the primary, and suggest that the secondary is an O3.5 If*/WN7 type star. We derive a high eccentricity of e = 0.78 and minimum masses of M1 sin^3 i ~ M2 sin^3 i ~ 13 +- 2 Msun, with q = M2 / M1 = 1.01 +- 0.07. An analysis of emission excess stemming from a wind-wind collision yields a similar inclination to that obtained from polarimetry (i = 39 +- 6deg). Our analysis thus implies M1 = 53^{+40}_{-20} and M2 = 54^{+40}_{-20} Msun, excluding M1 > 300Msun. A detailed comparison with evolution tracks calculated for single and binary stars, as well as the high eccentricity, suggest that the components of the system underwent quasi-homogeneous evolution and avoided mass-transfer. This scenario would suggest current masses of ~ 80 Msun and initial masses of Mi,1 ~ 105 and Mi,2 ~ 90Msun, consistent with the upper limits of our derived orbital masses, and would imply an age of ~2.2 Myr.

We present a new empirical prescription for the mass-loss rates of carbon and oxygen sequence Wolf-Rayet stars as a function of their luminosity, surface chemical composition, and initial metallicity. The new prescription is based on results of detailed spectral analyses of WC and WO stars, and improves the often applied Nugis & Lamers (2000) relation. We find that the mass-loss rates of WC and WO stars (with $X=0$ and $Y < 0.98$) can be expressed as $\log{\dot{M}} = -9.20 + 0.85\log{(L/L_{\odot})} + 0.44\log{Y} + 0.25\log{(Z_{\mathrm{Fe}}/Z_{\mathrm{Fe}, \odot})}$. This relation is based on mass-loss determinations that assume a volume-filling factor of 0.1, but the prescription can easily be scaled to account for other volume-filling factors. The residual of the fit is $\sigma = 0.06$ dex. We investigated whether the relation can also describe the mass loss of hydrogen-free WN stars and showed that it can when an adjustement of the metallicty dependence ($\log{\dot{M}} \propto 1.3\log{(Z_{\mathrm{Fe}}/Z_{\mathrm{Fe}, \odot})}$) is applied. Compared to Nugis & Lamers (2000), $\dot{M}$ is less sensitive to the luminosity and the surface abundance, implying a stronger mass loss of massive stars in their late stages of evolution. The modest metallicity dependence implies that if WC or WO stars are formed in metal deficient environments, their mass-loss rates are higher than currently anticipated. These effects may result in a larger number of type Ic supernovae and less black holes to be formed, and may favour the production of superluminous type Ic supernovae through interaction with C and O rich circumstellar material or the dense stellar wind.

We have used (a) HST ACS imaging and STIS spectroscopy, (b) ground-based PIONIER/VLT long-baseline interferometry, and (c) ground-based spectroscopy from different instruments to study the orbit of the extreme multiple system HD 93 129 Aa,Ab, which is composed of (at least) two very massive stars in a long-period orbit with e>0.92 that will pass through periastron in 2017/2018. In several ways, the system is an eta Car precursor. Around the time of periastron passage the two very strong winds will collide and generate an outburst of non-thermal hard X-ray emission without precedent in an O+O binary since astronomers have been able to observe above Earth's atmosphere. A coordinated multiwavelength monitoring in the next two years will enable a breakthrough understanding of the wind interactions in such extreme close encounters. Furthermore, we have found evidence that HD 93 129 Aa may be a binary system itself. In that case, we could witness a three-body interaction that may yield a runaway star or a stellar collision close to or shortly after the periastron passage. Either of those outcomes would be unprecedented, as they are predicted to be low-frequency events in the Milky Way.

Our long term aim is to derive model-independent stellar masses and distances for long period massive binaries by combining apparent astrometric orbit with double-lined radial velocity amplitudes (SB2). We follow-up ten O+O binaries with AMBER, PIONIER and GRAVITY at the VLTI. Here, we report about 130 astrometric observations over the last 7 years. We combine this dataset with distance estimates to compute the total mass of the systems. We also compute preliminary individual component masses for the five systems with available SB2 radial velocities. Nine over the ten binaries have their three dimensional orbit well constrained. Four of them are known colliding wind, non-thermal radio emitters, and thus constitute valuable targets for future high angular resolution radio imaging. Two binaries break the correlation between period and eccentricity tentatively observed in previous studies. It suggests either that massive star formation produce a wide range of systems, or that several binary formation mechanisms are at play. Finally, we found that the use of existing SB2 radial velocity amplitudes can lead to unrealistic masses and distances. If not understood, the biases in radial velocity amplitudes will represent an intrinsic limitation for estimating dynamical masses from SB2+interferometry or SB2+Gaia. Nevertheless, our results can be combined with future Gaia astrometry to measure the dynamical masses and distances of the individual components with an accuracy of 5 to 15\%, completely independently of the radial velocities.

Theoretically, rotation-induced chemical mixing in massive stars has far reaching evolutionary consequences, affecting the sequence of morphological phases, lifetimes, nucleosynthesis, and supernova characteristics. Using a sample of 72 presumably single O-type giants to supergiants observed in the context of the VLT-FLAMES Tarantula Survey (VFTS), we aim to investigate rotational mixing in evolved core-hydrogen burning stars initially more massive than $15\,M_\odot$ by analysing their surface nitrogen abundances. Using stellar and wind properties derived in a previous VFTS study, we constrained the nitrogen abundance by fitting the equivalent widths of relatively strong lines that are sensitive to changes in the abundance of this element. Given the quality of the data, we constrained the nitrogen abundance in 38 cases; for 34 stars only upper limits could be derived, which includes almost all stars rotating at $v_\mathrm{e}\sin i >200\,\mathrm{km s^{-1}}$. We analysed the nitrogen abundance as a function of projected rotation rate $v_\mathrm{e}\sin i$ and confronted it with predictions of rotational mixing. The upper limits on the nitrogen abundance of the rapidly rotating stars are not in apparent violation with theoretical expectations. However, we found a group of N-enhanced slowly-spinning stars that is not in accordance with predictions of rotational mixing in single stars. Among O-type stars with (rotation-corrected) gravities less than $\log\,g_c = 3.75$ this group constitutes 30$-$40 percent of the population. We found a correlation between nitrogen and helium abundance which is consistent with expectations, suggesting that, whatever the mechanism that brings N to the surface, it displays CNO-processed material.

We present medium-resolution, near-ultraviolet VLT/FLAMES observations of the star USNO-A0600-15865535. We adapt a standard method of stellar typing to our measurement of the shape of the Balmer epsilon absorption line to demonstrates that USNO-A0600-15865535 is a blue horizontal branch star, residing in the lower stellar halo at a distance of 4.4 kpc from the Sun. We measure the H & K lines of singly-ionized calcium and find two isolated velocity components, one originating in the disk, and one associated with high-velocity cloud complex WD. This detection demonstrated that complex WD is closer than ~4.4 kpc and is the first distance constraint on the +100 km/s Galactic complex of clouds. We find that Complex WD is not in corotation with the Galactic disk as has been assumed for decades. We examine a number of scenarios, and find that the most likely is that Complex WD was ejected from the solar neighborhood and is only a few kpc from the Sun.

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.

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.

Context. Luminous blue variables (LBVs) are rare massive stars with very high luminosity. They are characterized by strong photo-metric and spectroscopic variability related to transient eruptions. The mechanisms at the origin of these eruptions is not well known. In addition, their formation is still problematic and the presence of a companion could help to explain how they form. Aims. This article presents a study of seven LBVs (about 20% of the known Galactic population), some Wolf-Rayet stars, and massive binaries. We probe the environments that surround these massive stars with near-, mid-, and far-infrared images, investigating potential nebula/shells and the companion stars. Methods. To investigate large spatial scales, we used seeing-limited and near diffraction-limited adaptive optics images to obtain a differential diagnostic on the presence of circumstellar matter and to determine their extent. From those images, we also looked for the presence of binary companions on a wide orbit. Once a companion was detected, its gravitational binding to the central star was tested. Tests include the chance projection probability, the proper motion estimates with multi-epoch observations, flux ratio, and star separations. Results. We find that two out of seven of LBVs may have a wide orbit companion. Most of the LBVs display a large circumstellar envelope or several shells. In particular, HD168625, known for its rings, possesses several shells with possibly a large cold shell at the edge of which the rings are formed. For the first time, we have directly imaged the companion of LBV stars.

We present and describe the astro-photometric catalog of more than 800,000 sources found in the Hubble Tarantula Treasury Project (HTTP). HTTP is a Hubble Space Telescope (HST) Treasury program designed to image the entire 30 Doradus region down to the sub-solar (~0.5 solar masses) mass regime using the Wide Field Camera 3 (WFC3) and the Advanced Camera for Surveys (ACS). We observed 30 Doradus in the near ultraviolet (F275W, F336W), optical (F555W, F658N, F775W), and near infrared (F110W, F160W) wavelengths. The stellar photometry was measured using point-spread function (PSF) fitting across all the bands simultaneously. The relative astrometric accuracy of the catalog is 0.4 mas. The astro-photometric catalog, results from artificial star experiments and the mosaics for all the filters are available for download. Color-magnitude diagrams are presented showing the spatial distributions and ages of stars within 30 Dor as well as in the surrounding fields. HTTP provides the first rich and statistically significant sample of intermediate and low mass pre-main sequence candidates and allows us to trace how star formation has been developing through the region. The depth and high spatial resolution of our analysis highlight the dual role of stellar feedback in quenching and triggering star formation on the giant HII region scale. Our results are consistent with stellar sub-clustering in a partially filled gaseous nebula that is offset towards our side of the Large Magellanic Cloud.

In anticipation of the Gaia astrometric mission, a sample of spectroscopic binaries is being observed since 2010 with the Sophie spectrograph at the Haute--Provence Observatory. Our aim is to derive the orbital elements of double-lined spectroscopic binaries (SB2s) with an accuracy sufficient to finally obtain the masses of the components with relative errors as small as 1 % when combined with Gaia astrometric measurements. In order to validate the masses derived from Gaia, interferometric observations are obtained for three SB2s in our sample with F-K components: HIP 14157, HIP 20601 and HIP 117186. The masses of the six stellar components are derived. Due to its edge-on orientation, HIP 14157 is probably an eclipsing binary. We note that almost all the derived masses are a few percent larger than the expectations from the standard spectral-type-mass calibration and mass-luminosity relation. Our calculation also leads to accurate parallaxes for the three binaries, and the Hipparcos parallaxes are confirmed.

We present multi-epoch spectroscopic observations of the massive binary system WR21a, which include the January 2011 periastron passage. Our spectra reveal multiple SB2 lines and facilitate an accurate determination of the orbit and the spectral types of the components. We obtain minimum masses of $64.4\pm4.8 \ M_{\odot}$ and $36.3\pm1.7 \ M_{\odot}$ for the two components of WR21a. Using disentangled spectra of the individual components, we derive spectral types of O3/WN5ha and O3Vz~((f*)) for the primary and secondary, respectively. Using the spectral type of the secondary as an indication for its mass, we estimate an orbital inclination of $i=58.8\pm2.5^{\mathrm{o}}$ and absolute masses of $103.6\pm10.2 \ M_{\odot}$ and $58.3\pm3.7 \ M_{\odot}$, in agreement with the luminosity of the system. The spectral types of the WR21a components indicate that the stars are very young (1$-$2 Myr), similar to the age of the nearby Westerlund 2 cluster. We use evolutionary tracks to determine the mass-luminosity relation for the total system mass. We find that for a distance of 8 kpc and an age of 1.5 Myr, the derived absolute masses are in good agreement with those from evolutionary predictions.

The contact phase expected to precede the coalescence of two massive stars is poorly characterized due to the paucity of observational constraints. Here we report on the discovery of VFTS 352, an O-type binary in the 30 Doradus region, as the most massive and earliest spectral type overcontact system known to date. We derived the 3D geometry of the system, its orbital period $P_{\rm orb}=1.1241452(4)$ d, components' effective temperatures -- $T_1=42\,540\pm280$ K and $T_2=41\,120\pm290$ K -- and dynamical masses -- $M_1=28.63\pm0.30 M_{\odot}$ and $M_2 = 28.85\pm0.30 M_{\odot}$. Compared to single-star evolutionary models, the VFTS 352 components are too hot for their dynamical masses by about 2700 and 1100 K, respectively. These results can be explained naturally as a result of enhanced mixing, theoretically predicted to occur in very short-period tidally-locked systems. The VFTS 352 components are two of the best candidates identified so far to undergo this so-called chemically homogeneous evolution. The future of VFTS 352 is uncertain. If the two stars merge, a very rapidly rotating star will be produced. Instead, if the stars continue to evolve homogeneously and keep shrinking within their Roche Lobes, coalescence can be avoided. In this case, tides may counteract the spin down by winds such that the VFTS 352 components may, at the end of their life, fulfill the requirements for long gamma-ray burst (GRB) progenitors in the collapsar scenario. Independently of whether the VFTS 352 components become GRB progenitors, this scenario makes VFTS 352 interesting as a progenitor of a black hole binary, hence as a potential gravitational wave source through black hole-black hole merger.

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.

Extremely broad emission wings at H$\beta$ and H$\alpha$ have been found in VFTS data for five very luminous BA supergiants in or near 30 Doradus in the Large Magellanic Cloud. The profiles of both lines are extremely asymmetrical, which we have found to be caused by very broad diffuse interstellar bands (DIBs) in the longward wing of H$\beta$ and the shortward wing of H$\alpha$. These DIBs are well known to interstellar but not to many stellar specialists, so that the asymmetries may be mistaken for intrinsic features. The broad emission wings are generally ascribed to electron scattering, although we note difficulties for that interpretation in some objects. Such profiles are known in some Galactic hyper/supergiants and are also seen in both active and quiescent Luminous Blue Variables. No prior or current LBV activity is known in these 30 Dor stars, although a generic relationship to LBVs is not excluded; subject to further observational and theoretical investigation, it is possible that these very luminous supergiants are approaching the LBV stage for the first time. Their locations in the HRD and presumed evolutionary tracks are consistent with that possibility. The available evidence for spectroscopic variations of these objects is reviewed, while recent photometric monitoring does not reveal variability. A search for circumstellar nebulae has been conducted, with an indeterminate result for one of them.

Within the context of the collaboration "B fields in OB stars (BOB)", we used the FORS2 low-resolution spectropolarimeter to search for a magnetic field in 50 massive stars, including two reference magnetic massive stars. Because of the many controversies of magnetic field detections obtained with the FORS instruments, we derived the magnetic field values with two completely independent reduction and analysis pipelines. We compare and discuss the results obtained from the two pipelines. We obtained a general good agreement, indicating that most of the discrepancies on magnetic field detections reported in the literature are caused by the interpretation of the significance of the results (i.e., 3-4 sigma detections considered as genuine, or not), instead of by significant differences in the derived magnetic field values. By combining our results with past FORS1 measurements of HD46328, we improve the estimate of the stellar rotation period, obtaining P = 2.17950+/-0.00009 days. For HD125823, our FORS2 measurements do not fit the available magnetic field model, based on magnetic field values obtained 30 years ago. We repeatedly detect a magnetic field for the O9.7V star HD54879, the HD164492C massive binary, and the He-rich star CPD -57 3509. We obtain a magnetic field detection rate of 6+/-4%, while by considering only the apparently slow rotators we derive a detection rate of 8+/-5%, both comparable with what was previously reported by other similar surveys. We are left with the intriguing result that, although the large majority of magnetic massive stars is rotating slowly, our detection rate is not a strong function of the stellar rotational velocity.

We present measurements of positions and relative proper motions in the 30 Doradus region of the Large Magellanic Cloud (LMC). We detail the construction of a single-epoch astrometric reference frame, based on specially-designed observations obtained with the two main imaging instruments ACS/WFC and WFC3/UVIS onboard the Hubble Space Telescope (HST). Internal comparisons indicate a sub milli-arc-second (mas) precision in the positions and the presence of semi-periodic systematics with a mean amplitude of ~0.8 mas. We combined these observations with numerous archival images taken with WFPC2 and spanning 17 years. The precision of the resulting proper motions for well-measured stars around the massive cluster R 136 can be as good as ~20 microarcsec/yr, although the true accuracy of proper motions is generally lower due to the residual systematic errors. The observed proper-motion dispersion for our highest-quality measurements is ~0.1 mas/yr. Our catalog of positions and proper motions contains 86,590 stars down to V~25 and over a total area of ~70 square arcmin. We examined the proper motions of 105 relatively bright stars and identified a total of 6 candidate runaway stars. We are able to tentatively confirm the runaway status of star VFTS 285, consistent with the findings from line-of-sight velocities, and to show that this star has likely been ejected from R 136. This study demonstrates that with HST it is now possible to reliably measure proper motions of individual stars in the nearest dwarf galaxies such as the LMC.

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.

The initial distribution of spin rates of massive stars is a fingerprint of their elusive formation process. It also sets a key initial condition for stellar evolution and is thus an important ingredient in stellar population synthesis. So far, most studies have focused on single stars. Most O stars are however found in multiple systems. By establishing the spin-rate distribution of a sizeable sample of O-type spectroscopic binaries and by comparing the distributions of binary sub-populations with one another as well as with that of presumed single stars in the same region, we aim to constrain the initial spin distribution of O stars in binaries, and to identify signatures of the physical mechanisms that affect the evolution of the massive stars spin rates. We use ground-based optical spectroscopy obtained in the framework of the VLT-FLAMES Tarantula Survey (VFTS) to establish the projected equatorial rotational velocities (\vrot) for components of 114 spectroscopic binaries in 30 Doradus. The \vrot\ values are derived from the full-width at half-maximum (FWHM) of a set of spectral lines, using a FWHM vs. \vrot\ calibration that we derive based on previous line analysis methods applied to single O-type stars in the VFTS sample. The overall \vrot\ distribution of the primary stars resembles that of single O-type stars in the VFTS, featuring a low-velocity peak (at $\vrot < 200$ kms) and a shoulder at intermediate velocities ($200 < \vrot < 300$ kms). The distributions of binaries and single stars however differ in two ways. First, the main peak at $\vrot \sim$100 kms is broader and slightly shifted toward higher spin rates in the binary distribution compared to that of the presumed-single stars. Second, the \vrot distribution of primaries lacks a significant population of stars spinning faster than 300 kms while such a population is clearly present in the single star sample.

Context. Oxygen sequence Wolf-Rayet (WO) stars represent a very rare stage in the evolution of massive stars. Their spectra show strong emission lines of helium-burning products, in particular highly ionized carbon and oxygen. The properties of WO stars can be used to provide unique constraints on the (post-)helium burning evolution of massive stars, as well as their remaining lifetime and the expected properties of their supernovae. Aims. We aim to homogeneously analyse the currently known presumed-single WO stars to obtain the key stellar and outflow properties and to constrain their evolutionary state. Methods. We use the line-blanketed non-local thermal equilibrium atmosphere code cmfgen to model the X-Shooter spectra of the WO stars and deduce the atmospheric parameters. We calculate dedicated evolutionary models to determine the evolutionary state of the stars. Results. The WO stars have extremely high temperatures that range from 150 kK to 210 kK, and very low surface helium mass fractions that range from 44% down to 14%. Their properties can be reproduced by evolutionary models with helium zero-age main sequence masses of $M_{\mathrm{He, ini}} = 15-25 M_{\odot}$ that exhibit a fairly strong (on the order of a few times $10^{-5} M_{\odot} \mathrm{yr}^{-1}$), homogeneous ($f_\mathrm{c} > 0.3$) stellar wind. Conclusions. WO stars represent the final evolutionary stage of stars with estimated initial masses of $M_{\mathrm{ini}} = 40-60 M_{\odot}$. They are post core-helium burning and predicted to explode as type Ic supernovae within a few thousand years.

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.

This White Paper presents the scientific motivations for a multi-object spectrograph (MOS) on the European Extremely Large Telescope (E-ELT). The MOS case draws on all fields of contemporary astronomy, from extra-solar planets, to the study of the halo of the Milky Way and its satellites, and from resolved stellar populations in nearby galaxies out to observations of the earliest 'first-light' structures in the partially-reionised Universe. The material presented here results from thorough discussions within the community over the past four years, building on the past competitive studies to agree a common strategy toward realising a MOS capability on the E-ELT. The cases have been distilled to a set of common requirements which will be used to define the MOSAIC instrument, entailing two observational modes ('high multiplex' and 'high definition'). When combined with the unprecedented sensitivity of the E-ELT, MOSAIC will be the world's leading MOS facility. In analysing the requirements we also identify a high-multiplex MOS for the longer-term plans for the E-ELT, with an even greater multiplex (>1000 targets) to enable studies of large-scale structures in the high-redshift Universe. Following the green light for the construction of the E-ELT the MOS community, structured through the MOSAIC consortium, is eager to realise a MOS on the E-ELT as soon as possible. We argue that several of the most compelling cases for ELT science, in highly competitive areas of modern astronomy, demand such a capability. For example, MOS observations in the early stages of E-ELT operations will be essential for follow-up of sources identified by the James Webb Space Telescope (JWST). In particular, multi-object adaptive optics and accurate sky subtraction with fibres have both recently been demonstrated on sky, making fast-track development of MOSAIC feasible.

The stellar population of the 30 Doradus star-forming region in the Large Magellanic Cloud contains a subset of apparently single, rapidly rotating O-type stars. The physical processes leading to the formation of this cohort are currently uncertain. One member of this group, the late O-type star VFTS 399, is found to be unexpectedly X-ray bright for its bolometric luminosity - in this study we aim to determine its physical nature and the cause of this behaviour. We find VFTS 399 to be an aperiodic photometric variable with an apparent near-IR excess. Its optical spectrum demonstrates complex emission profiles in the lower Balmer series and select HeI lines - taken together these suggest an OeBe classification. The highly variable X-ray luminosity is too great to be produced by a single star, while the hard, non-thermal nature suggests the presence of an accreting relativistic companion. Finally, the detection of periodic modulation of the X-ray lightcurve is most naturally explained under the assumption that the accretor is a neutron star. VFTS 399 appears to be the first high-mass X-ray binary identified within 30 Dor. Comparison of the current properties of VFTS 399 to binary-evolution models suggests a progenitor mass in excess of 25Msun for the putative neutron star, which may host a magnetic field comparable in strength to those of magnetars. VFTS 399 is now the second member of the cohort of rapidly rotating `single' O-type stars in 30 Dor to show evidence of binary interaction resulting in spin-up, suggesting that this may be a viable evolutionary pathway for the formation of a subset of this stellar population.

Context: Absorption by molecules in the Earth's atmosphere strongly affects ground-based astronomical observations. The resulting absorption line strength and shape depend on the highly variable physical state of the atmosphere, i.e. pressure, temperature, and mixing ratio of the different molecules involved. Usually, supplementary observations of so-called telluric standard stars (TSS) are needed to correct for this effect, which is expensive in terms of telescope time. We have developed the software package molecfit to provide synthetic transmission spectra based on parameters obtained by fitting narrow ranges of the observed spectra of scientific objects. These spectra are calculated by means of the radiative transfer code LBLRTM and an atmospheric model. In this way, the telluric absorption correction for suitable objects can be performed without any additional calibration observations of TSS. Aims: We evaluate the quality of the telluric absorption correction using molecfit with a set of archival ESO-VLT X-Shooter visible and near-infrared spectra. Methods: Thanks to the wavelength coverage from the U to the K band, X-Shooter is well suited to investigate the quality of the telluric absorption correction with respect to the observing conditions, the instrumental set-up, input parameters of the code, the signal-to-noise of the input spectrum, and the atmospheric profiles. These investigations are based on two figures of merit, I_off and I_res, that describe the systematic offsets and the remaining small-scale residuals of the corrections. We also compare the quality of the telluric absorption correction achieved with moelcfit to the classical method based on a telluric standard star. (Abridged)

Context: The interaction of the light from astronomical objects with the constituents of the Earth's atmosphere leads to the formation of telluric absorption lines in ground-based collected spectra. Correcting for these lines, mostly affecting the red and infrared region of the spectrum, usually relies on observations of specific stars obtained close in time and airmass to the science targets, therefore using precious observing time. Aims: We present molecfit, a tool for correcting for telluric absorption lines based on synthetic modelling of the Earth's atmospheric transmission. Molecfit is versatile and can be used with data obtained with various ground-based telescopes and instruments. Methods: Molecfit combines a publicly available radiative transfer code, a molecular line database, atmospheric profiles, and various kernels to model the instrument line spread function. The atmospheric profiles are created by merging a standard atmospheric profile representative of a given observatory's climate, of local meteorological data, and of dynamically retrieved altitude profiles for temperature, pressure, and humidity. We discuss the various ingredients of the method, its applicability, and its limitations. We also show examples of telluric line correction on spectra obtained with a suite of ESO Very Large Telescope (VLT) instruments. Results: Compared to previous similar tools, molecfit takes the best results for temperature, pressure, and humidity in the atmosphere above the observatory into account. As a result, the standard deviation of the residuals after correction of unsaturated telluric lines is frequently better than 2% of the continuum. Conclusion: Molecfit is able to accurately model and correct for telluric lines over a broad range of wavelengths and spectral resolutions. (Abridged)

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.

We present a dense model grid with tailored input chemical composition appropriate for the Large Magellanic Cloud. We use a one-dimensional hydrodynamic stellar evolution code, which accounts for rotation, transport of angular momentum by magnetic fields, and stellar wind mass loss to compute our detailed models. We calculate stellar evolution models with initial masses of 70-500 Msun and with initial surface rotational velocities of 0-550 km/s, covering the core-hydrogen burning phase of evolution. We find our rapid rotators to be strongly influenced by rotationally induced mixing of helium, with quasi-chemically homogeneous evolution occurring for the fastest rotating models. Above 160 Msun, homogeneous evolution is also established through mass loss, producing pure helium stars at core hydrogen exhaustion independent of the initial rotation rate. Surface nitrogen enrichment is also found for slower rotators, even for stars that lose only a small fraction of their initial mass. For models above 150 MZAMS, and for models in the whole considered mass range later on, we find a considerable envelope inflation due to the proximity of these models to their Eddington limit. This leads to a maximum zero-age main sequence surface temperature of 56000 K, at 180 Msun, and to an evolution of stars in the mass range 50-100 Msun to the regime of luminous blue variables in the HR diagram with high internal Eddington factors. Inflation also leads to decreasing surface temperatures during the chemically homogeneous evolution of stars above 180 Msun. The cool surface temperatures due to the envelope inflation in our models lead to an enhanced mass loss, which prevents stars at LMC metallicity from evolving into pair-instability supernovae. The corresponding spin-down will also prevent very massive LMC stars to produce long-duration gamma-ray bursts, which might, however, originate from lower masses.

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

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.

Massive stars likely played an important role in the reionization of the Universe, and the formation of the first black holes. Massive stars in low-metallicity environments in the local Universe are reminiscent of their high redshift counterparts. In a previous paper, we reported on indications that the stellar winds of low-metallicity O stars may be stronger than predicted, which would challenge the current paradigm of massive star evolution. In this paper, we aim to extend our initial sample of six O stars in low-metallicity environments by four. We aim to derive their stellar and wind parameters, and compare these to radiation-driven wind theory and stellar evolution models. We have obtained intermediate-resolution VLT/X-Shooter spectra of our sample of stars. We derive the stellar parameters by fitting synthetic fastwind line profiles to the VLT/X-Shooter spectra using a genetic fitting algoritm. We compare our parameters to evolutionary tracks and obtain evolutionary masses and ages. We also investigate the effective temperature versus spectral type calibration for SMC and lower metallicities. Finally, we reassess the wind momentum versus luminosity diagram. The derived parameters of our target stars indicate stellar masses that reach values of up to 50 $M_{\odot}$. The wind strengths of our stars are, on average, stronger than predicted from radiation-driven wind theory and reminiscent of stars with an LMC metallicity. We discuss indications that the iron content of the host galaxies is higher than originally thought and is instead SMC-like. We find that the discrepancy with theory is lessened, but remains significant for this higher metallicity. This may imply that our current understanding of the wind properties of massive stars, both in the local universe as well as at cosmic distances, remains incomplete.

Multiplicity is one of the most fundamental observable properties of massive O-type stars and offers a promising way to discriminate between massive star formation theories. Nevertheless, companions at separations between 1 and 100 mas remain mostly unknown due to intrinsic observational limitations. [...] The Southern MAssive Stars at High angular resolution survey (SMASH+) was designed to fill this gap by providing the first systematic interferometric survey of Galactic massive stars. We observed 117 O-type stars with VLTI/PIONIER and 162 O-type stars with NACO/SAM, respectively probing the separation ranges 1-45 and 30-250mas and brightness contrasts of Delta H < 4 and Delta H < 5. Taking advantage of NACO's field-of-view, we further uniformly searched for visual companions in an 8''-radius down to Delta H = 8. This paper describes the observations and data analysis, reports the discovery of almost 200 new companions in the separation range from 1mas to 8'' and presents the catalog of detections, including the first resolved measurements of over a dozen known long-period spectroscopic binaries. Excluding known runaway stars for which no companions are detected, 96 objects in our main sample (DEC < 0 deg; H<7.5) were observed both with PIONIER and NACO/SAM. The fraction of these stars with at least one resolved companion within 200mas is 0.53. Accounting for known but unresolved spectroscopic or eclipsing companions, the multiplicity fraction at separation < 8'' increases to f_m = 0.91 +/- 0.03. The fraction of luminosity class V stars that have a bound companion reaches 100% at 30mas while their average number of physically connected companions within 8'' is f_c = 2.2 +/- 0.3. This demonstrates that massive stars form nearly exclusively in multiple systems. Additionally, the nine non-thermal (NT) radio emitters observed by SMASH+ are all resolved [...]

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.

It is now well established that a fraction of the massive (M>8 Msun) star population hosts strong, organised magnetic fields, most likely of fossil origin. The details of the generation and evolution of these fields are still poorly understood. The BinaMIcS project takes an important step towards the understanding of the interplay between binarity and magnetism during the stellar formation and evolution, and in particular the genesis of fossil fields, by studying the magnetic properties of close binary systems. The components of such systems are most likely formed together, at the same time and in the same environment, and can therefore help us to disentangle the role of initial conditions on the magnetic properties of the massive stars from other competing effects such as age or rotation. We present here the main scientific objectives of the BinaMIcS project, as well as preliminary results from the first year of observations from the associated ESPaDOnS and Narval spectropolarimetric surveys.

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.

Delta Cir is known as an O7.5 III eclipsing and spectroscopic binary with an eccentric orbit. Penny et al. discovered the presence of a third component in the IUE spectra. The eclipsing binary and the third body revolve around a common centre of gravity with a period of 1644 days in an eccentric orbit with the semimajor axis of 10 AU. We demonstrate the presence of the apsidal-line rotation with a period of about 141 years, which is considerably longer than its theoretically predicted value, based on the published radii of the binary components derived from the Hipparchos Hp light curve. However, our new solution of the same light curve resulted in smaller radii and a better agreement between the observed and predicted period of the apsidal- line advance. There are indications that the third body is a binary. The object was resolved by VLTI with the PIONIER combiner; in June 2012 the separation was 3.78 mas, magnitude difference in the H region 1.75 mag. This result means that (assuming the distance 770 pc) the inclination of the long orbit is 87.7 degrees.

Analyses of multi-epoch, high-resolution (R ~ 50.000) optical spectra of the O-type star HD 152246 (O9 IV according to the most recent classification), complemented by a limited number of earlier published radial velocities, led to the finding that the object is a hierarchical triple system, where a close inner pair (Ba-Bb) with a slightly eccentric orbit (e = 0.11) and a period of 6.0049 days revolves in a 470-day highly eccentric orbit (e = 0.865) with another massive and brighter component A. The mass ratio of the inner system must be low since we were unable to find any traces of the secondary spectrum. The mass ratio A/(Ba+Bb) is 0.89. The outer system has recently been resolved using long-baseline interferometry on three occasions. The interferometry confirms the spectroscopic results and specifies elements of the system. Our orbital solutions, including the combined radial-velocity and interferometric solution indicate an orbital inclination of the outer orbit of 112{\deg} and stellar masses of 20.4 and 22.8 solar masses. We also disentangled the spectra of components A and Ba and compare them to synthetic spectra from two independent programmes, TLUSTY and FASTWIND. In either case, the fit was not satisfactory and we postpone a better determination of the system properties for a future study, after obtaining observations during the periastron passage of the outer orbit (the nearest chance being March 2015). For the moment, we can only conclude that component A is an O9 IV star with v*sin(i) = 210 +\- 10 km/s and effective temperature of 33000 +\- 500 K, while component Ba is an O9 V object with v*sin(i) = 65 +/- 3 km/s and T_eff = 33600 +\- 600 K.

The enigmatic oxygen sequence Wolf-Rayet (WO) stars represent a very late stage in massive star evolution, although their exact nature is still under debate. The spectra of most of the WO stars have never been analysed through detailed modelling with a non-local thermodynamic equilibrium expanding atmosphere code. Here we present preliminary results of the first homogeneous analysis of the (apparently) single WOs.

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.

We report the discovery of the B[e] star VFTS 822 in the 30 Doradus star-forming region of the Large Magellanic Cloud, classified by optical spectroscopy from the VLT-FLAMES Tarantula Survey and complementary infrared photometry. VFTS 822 is a relatively low-luminosity (log $L$ = 4.04 $\pm$ 0.25 $L_{\odot}$) B8[e] star. In this Letter, we evaluate the evolutionary status of VFTS 822 and discuss its candidacy as a Herbig B[e] star. If the object is indeed in the pre-main sequence phase, it would present an exciting opportunity to measure mass accretion rates at low metallicity spectroscopically, to understand the effect of metallicity on accretion rates.

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.

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.

Correcting for the sky signature usually requires supplementary calibration data which are very expensive in terms of telescope time. In addition, the scheduling flexibility is restricted as these data have to be taken usually directly before/after the science observations due to the high variability of the telluric absorption which depends on the state and the chemical composition of the atmosphere at the time of observations. Therefore, a tool for sky correction, which does not require this supplementary calibration data, saves a significant amount of valuable telescope time and increases its efficiency. We developed a software package aimed at performing telluric feature corrections on the basis of synthetic absorption spectra.

Herbig Ae/Be stars are intermediate-mass pre-main sequence stars surrounded by circumstellar dust disks. Some are observed to produce jets, whose appearance as a sequence of shock fronts (knots) suggests a past episodic outflow variability. This "jet fossil record" can be used to reconstruct the outflow history. We present the first optical to near-infrared (NIR) VLT/X-shooter spectra of the jet from the Herbig Ae star HD 163296. We determine physical conditions in the knots, as well as their kinematic "launch epochs". Knots are formed simultaneously on either side of the disk, with a regular interval of ~16 yr. The velocity dispersion versus jet velocity and the energy input are comparable in both lobes. However, the mass loss rate, velocity, and shock conditions are asymmetric. We find Mjet/Macc ~ 0.01-0.1, consistent with magneto-centrifugal jet launching models. No evidence for dust is found in the high-velocity jet, suggesting it is launched within the sublimation radius (<0.5 au). The jet inclination measured from proper motions and radial velocities confirms it is perpendicular to the disk. A tentative relation is found between the structure of the jet and the photometric variability of the source. Episodes of NIR brightening were previously detected and attributed to a dusty disk wind. We report for the first time significant optical fadings lasting from a few days up to a year, coinciding with the NIR brightenings. These are likely caused by dust lifted high above the disk plane; this supports the disk wind scenario. The disk wind is launched at a larger radius than the high-velocity atomic jet, although their outflow variability may have a common origin. No significant relation between outflow and accretion variability could be established. Our findings confirm that this source undergoes periodic ejection events, which may be coupled with dust ejections above the disk plane.

Because the majority of massive stars are born as members of close binary systems, populations of massive main-sequence stars contain stellar mergers and products of binary mass transfer. We simulate populations of massive stars accounting for all major binary evolution effects based on the most recent binary parameter statistics and extensively evaluate the effect of model uncertainties. Assuming constant star formation, we find that $8^{+9}_{-4}\,\%$ of a sample of early type stars to be the product of a merger resulting from a close binary system. In total we find that $30^{+10}_{-15}\,\%$ of massive main-sequence stars are the product of binary interaction. We show that the commonly adapted approach to minimize the effects of binaries on an observed sample by excluding systems detected as binaries through radial velocity campaigns can be counterproductive. Systems with significant radial velocity variations are mostly pre-interaction systems. Excluding them substantially enhances the relative incidence of mergers and binary products in the non radial velocity variable sample. This poses a challenge for testing single stellar evolutionary models. It also raises the question of whether certain peculiar classes of stars, such as magnetic O-stars, are the result of binary interaction and it emphasizes the need to further study the effect of binarity on the diagnostics that are used to derive the fundamental properties (star-formation history, initial mass function, mass to light ratio) of stellar populations nearby and at high redshift.

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.

We present the results of a quantitative spectroscopic analysis of the oxygen-sequence Wolf- Rayet star DR1 in the low-metallicity galaxy IC 1613. Our models suggest that the strong oxygen emission lines are the result of the high temperature of this WO3 star and do not necessarily reflect a more advanced evolutionary stage than WC stars.

Massive stars rapidly change their masses through strong stellar winds and mass transfer in binary systems. We show that such mass changes leave characteristic signatures in stellar mass functions of young star clusters which can be used to infer their ages and to identify products of binary evolution. We model the observed present day mass functions of the young Galactic Arches and Quintuplet star clusters using our rapid binary evolution code. We find that shaping of the mass function by stellar wind mass loss allows us to determine the cluster ages to 3.5$\pm$0.7 Myr and 4.8$\pm$1.1 Myr, respectively. Exploiting the effects of binary mass exchange on the cluster mass function, we find that the most massive stars in both clusters are rejuvenated products of binary mass transfer, i.e. the massive counterpart of classical blue straggler stars. This resolves the problem of an apparent age spread among the most luminous stars exceeding the expected duration of star formation in these clusters. We perform Monte Carlo simulations to probe stochastic sampling, which support the idea of the most massive stars being rejuvenated binary products. We find that the most massive star is expected to be a binary product after 1.0$\pm$0.7 Myr in Arches and after 1.7$\pm$1.0 Myr in Quintuplet. Today, the most massive 9$\pm$3 stars in Arches and 8$\pm$3 in Quintuplet are expected to be such objects. Our findings have strong implications for the stellar upper mass limit and solve the discrepancy between the claimed 150 $\mathrm{M}_\odot$ limit and observations of fours stars with initial masses of 165-320 $\mathrm{M}_\odot$ in R136 and of SN 2007bi, which is thought to be a pair-instability supernova from an initial 250 $\mathrm{M}_\odot$ star. Using the stellar population of R136, we revise the upper mass limit to values in the range 200-500 $\mathrm{M}_\odot$.

The VLT-FLAMES Tarantula Survey (VFTS) has acquired multi-epoch spectroscopy of over 800 O, B and Wolf-Rayet stars in the 30 Doradus region with the aim to investigate a number of important questions related to the evolution of massive stars and of cluster dynamics. In this paper, I first provide an overview of the scientific results obtained by the VFTS consortium so far. I then review the constraints obtained on the multiplicity properties of massive stars in 30 Dor and compare them to our recent results from a Milky Way sample.

Context. Oxygen sequence Wolf-Rayet (WO) stars are thought to represent the final evolutionary stage of the most massive stars. The characteristic strong O vi emission possibly originates from an enhanced oxygen abundance in the stellar wind. Alternatively, the O vi emission can be caused by the high temperature of these stars, in which case the WO stars are the high-temperature extension of the more common carbon sequence Wolf-Rayet (WC) stars. Aims. By constraining the physical properties and evolutionary status of DR1, a WO star in the low-metallicity Local Group dwarf galaxy IC 1613 and one of only two objects of its class known in a SMC-like metallicity environment, we aim to investigate the nature of WO stars and their evolutionary connection with WC stars. Methods. We use the non-LTE atmosphere code cmfgen to model the observed spectrum of DR1 and to derive its stellar and wind parameters. We compare our values with other studies of WC and WO stars, as well as with the predictions of evolutionary models. We also model the surrounding nebula using the photo-ionization code cloudy. Results. The oxygen and carbon abundances that we obtain are comparable to values found for WC stars. The temperature and luminosity are, however, higher than those of WC stars. DR1 is embedded in the hottest known H ii region in the Local Group. The nebular properties can be consistently reproduced by cloudy models adopting DR1 as central ionizing source. Conclusions. Comparison of the abundances and temperature of DR1 with core helium-burning models show that DR1 is currently well into the second half of helium burning. If the properties of DR1 are representative for the WO class, it would imply that WO stars are the high-temperature and high-luminosity extension of the WC stars, and do not necessarily represent a later evolutionary stage.

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

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

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

Massive stars play a dominant role in the process of clustered star formation, with their feedback into the molecular cloud through ionizing radiation, stellar winds and outflows. The formation process of massive stars is poorly constrained because of their scarcity, the short formation timescale and obscuration. By obtaining a census of the newly formed stellar population, the star formation history of the young cluster and the role of the massive stars within it can be unraveled. We aim to reconstruct the formation history of the young stellar population of the massive star-forming region RCW 36. We study several dozens of individual objects, both photometrically and spectroscopically, look for signs of multiple generations of young stars and investigate the role of the massive stars in this process. We obtain a census of the physical parameters and evolutionary status of the young stellar population. Using a combination of near-infrared photometry and spectroscopy we estimate ages and masses of individual objects. We identify the population of embedded young stellar objects (YSO) by their infrared colors and emission line spectra. RCW 36 harbors a stellar population of massive and intermediate-mass stars located around the center of the cluster. Class 0/I and II sources are found throughout the cluster. The central population has a median age of 1.1 +/- 0.6 Myr. Of the stars which could be classified, the most massive ones are situated in the center of the cluster. The central cluster is surrounded by filamentary cloud structures; within these, some embedded and accreting YSOs are found. Our age determination is consistent with the filamentary structures having been shaped by the ionizing radiation and stellar winds of the central massive stars. The formation of a new generation of stars is ongoing, as demonstrated by the presence of embedded protostellar clumps, and two exposed jets.

Using XMM data, we study for the first time the X-ray emission of HM1 and IC2944/2948. Low-mass, pre-main-sequence objects with an age of a few Myr are detected, as well as a few background or foreground objects. Most massive stars in both clusters display the usual high-energy properties of that type of objects, though with log(Lx/Lbol) apparently lower in HM1 than in IC2944/2948. Compared with studies of other clusters, it seems that a low signal-to-noise ratio at soft energies, due to the high extinction, may be the main cause of this difference. In HM1, the two Wolf-Rayet stars show contrasting behaviors: WR89 is extremely bright, but much softer than WR87. It remains to be seen whether wind-wind collisions or magnetically confined winds can explain these emissions. In IC2944/2948, the X-ray sources concentrate around HD101205; a group of massive stars to the north of this object is isolated, suggesting that there exist two subclusters in the field-of-view.

R144 is a WN6h star in the 30 Doradus region. It is suspected to be a binary because of its high luminosity and its strong X-ray flux, but no periodicity could be established so far. Here, we present new Xshooter multi-epoch spectroscopy of R144 obtained at the ESO Very Large Telescope (VLT). We detect variability in position and/or shape of all the spectral lines. We measure radial velocity variations with an amplitude larger than 250 km/s in NIV and NV lines. Furthermore, the NIII and NV line Doppler shifts are anti-correlated and the NIV lines show a double-peaked profile on six of our seven epochs. We thus conclude that R144 is a double-lined spectroscopic binary. Possible orbital periods range from 2 to 6 months, although a period up to one year is allowed if the orbit is highly eccentric. We estimate the spectral types of the components to be WN5-6h and WN6-7h, respectively. The high luminosity of the system (log Lbol/Lsun ~ 6.8) suggests a present-day total mass content in the range of about 200 to 300 Msun, depending on the evolutionary stage of the components. This makes R144 the most massive binary identified so far, with a total mass content at birth possibly as large as 400 Msun. We briefly discuss the presence of such a massive object 60 pc away from the R136 cluster core in the context of star formation and stellar dynamics.

Context. HD 150136 is a triple hierarchical system and a non-thermal radio emitter. It is formed by an O3-3.5 V + O5.5-6 V close binary and a more distant O6.5-7 V tertiary. So far, only the inner orbital properties have been reliably constrained. Aims. To quantitatively understand the non-thermal emission process, accurate knowledge of the physical and orbital properties of the object is crucial. Here, we aim to investigate the orbital properties of the wide system and to constrain the inclinations of the inner and outer binaries, and with these the absolute masses of the system components. Methods. We used the PIONIER combiner at the Very Large Telescope Interferometer to obtain the very first interferometric measurements of HD 150136. We combined the interferometric observations with new and existing high resolution spectroscopic data to derive the orbital solution of the outer companion in the three-dimensional space. Results. The wide system is clearly resolved by PIONIER, with a projected separation on the plane of the sky of about 9 milli-arcsec. The best-fit orbital period, eccentricity, and inclination are 8.2 yr, 0.73 and 108 degr. We constrain the masses of the three stars of the system to 63 +/- 10, 40 +/- 6, and 33 +/- 12 Msun for the O3-3.5 V, O5.5-6 V and O6.5-7 V components. Conclusions. The dynamical masses agree within errors with the evolutionary masses of the components. Future interferometric and spectroscopic monitoring of HD 150136 should allow one to reduce the uncertainties to a few per cent only and to accurately constrain the distance to the system. This makes HD 150136 an ideal system to quantitatively test evolutionary models of high-mass stars as well as the physics of non-thermal processes occurring in O-type systems.

Globular clusters contain many stars with surface abundance patterns indicating contributions from hydrogen burning products, as seen in the anti-correlated elemental abundances of e.g. sodium and oxygen, and magnesium and aluminium. Multiple generations of stars can explain this phenomenon, with the second generation forming from a mixture of pristine gas and ejecta from the first generation. We show that massive binary stars may be a source of much of the material that makes this second generation of stars. Mass transfer in binaries is often non-conservative and the ejected matter moves slowly enough that it can remain inside a globular cluster and remain available for subsequent star formation. Recent studies show that there are more short-period massive binaries than previously thought, hence also more stars that interact and eject nuclear-processed material.

Rotation is thought to be a major factor in the evolution of massive stars, especially at low metallicity, with consequences for their chemical yields, ionizing flux and final fate. Determining the natal rotation-rate distribution of stars is of high priority given its importance as a constraint on theories of massive star formation and as input for models of stellar populations in the local Universe and at high redshift. Recently, it has become clear that the majority of massive stars interact with a binary companion before they die. We investigate how this affects the distribution of rotation rates. For this purpose, we simulate a massive binary-star population typical for our Galaxy assuming continuous star formation. We find that, because of binary interaction, 20^+5_-10% of all massive main-sequence stars have projected rotational velocities in excess of 200km/s. We evaluate the effect of uncertain input distributions and physical processes and conclude that the main uncertainties are the mass transfer efficiency and the possible effect of magnetic braking, especially if magnetic fields are generated or amplified during mass accretion and stellar mergers. The fraction of rapid rotators we derive is similar to that observed. If indeed mass transfer and mergers are the main cause for rapid rotation in massive stars, little room remains for rapidly rotating stars that are born single. This implies that spin down during star formation is even more efficient than previously thought. In addition, this raises questions about the interpretation of the surface abundances of rapidly rotating stars as evidence for rotational mixing. Furthermore, our results allow for the possibility that all early-type Be stars result from binary interactions and suggest that evidence for rotation in explosions, such as long gamma-ray bursts, points to a binary origin.

Jets around low- and intermediate-mass young stellar objects (YSOs) contain a fossil record of the recent accretion and outflow activity of their parent star-forming systems. We aim to understand whether the accretion/ejection process is similar across the entire stellar mass range of the parent YSOs. To this end we have obtained VLT/X-shooter spectra of HH 1042 and HH 1043, two newly discovered jets in the massive star-forming region RCW 36. HH 1042 is associated with the intermediate-mass YSO 08576nr292. Over 90 emission lines are detected in the spectra. High-velocity (up to 220 km/s) blue- and redshifted emission from a bipolar flow is observed in typical shock tracers. Low-velocity emission from the background cloud is detected in nebular tracers, including lines from high ionization species. We applied combined optical and infrared spectral diagnostic tools in order to derive the physical conditions (density, temperature, and ionization) in the jets. The measured mass outflow rates are Mjet ~ 10^-7 Msun/yr. We measure a high accretion rate for HH 1042 (Macc ~ 10^-6 Msun/yr) and Mjet/Macc ~ 0.1, comparable to low-mass sources and consistent with models for magneto-centrifugal jet launching. The knotted structure and velocity spread in both jets are interpreted as fossil signatures of a variable outflow rate. The mean velocities in both lobes of the jets are comparable, but the variations in Mjet and velocity in the two lobes are not symmetric, suggesting that the launching mechanism on either side of the accretion disk is not synchronized. For HH 1042, we have constructed an interpretative physical model with a stochastic or periodic outflow rate and a description of a ballistic flow as its constituents. The knotted structure and velocity spread can be reproduced qualitatively with the model, indicating that the outflow velocity varies on timescales on the order of 100 yr.

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.

Nearby companions alter the evolution of massive stars in binary systems. Using a sample of Galactic massive stars in nearby young clusters, we simultaneously measure all intrinsic binary characteristics relevant to quantify the frequency and nature of binary interactions. We find a large intrinsic binary fraction, a strong preference for short orbital periods and a flat distribution for the mass-ratios. Our results do not support the presence of a significant peak of equal-mass `twin' binaries. As a result of the measured distributions, we find that over seventy per cent of all massive stars exchange mass with a companion. Such a rate greatly exceeds previous estimates and implies that the majority of massive stars have their evolution strongly affected by interaction with a nearby companion.

We have performed new wide-field photometry of the young open cluster NGC 6231 to study the shape of the initial mass function (IMF) and mass segregation. We also investigated the reddening law toward NGC 6231 from optical to mid-infrared color excess ratios, and found that the total-to-selective extinction ratio is Rv = 3.2, which is very close to the normal value. But many early-type stars in the cluster center show large color excess ratios. We derived the surface density profiles of four member groups, and found that they reach the surface density of field stars at about 10', regardless of stellar mass. The IMF of NGC 6231 is derived for the mass range 0.8 -- 45 Msun. The slope of the IMF of NGC 6231 (Gamma = -1.1 +/- 0.1) is slightly shallower than the canonical value, but the difference is marginal. In addition, the mass function varies systematically, and is a strong function of radius - it is is very shallow at the center, and very steep at the outer ring suggesting the cluster is mass segregated. We confirm the mass segregation for the massive stars (m >~ 8 Msun) by a minimum spanning tree analysis. Using a Monte Carlo method, we estimate the total mass of NGC 6231 to be about 2.6 (+/- 0.6) x 10^3 Msun. We constrain the age of NGC 6231 by comparison with evolutionary isochrones. The age of the low-mass stars ranges from 1 to 7 Myr with a slight peak at 3 Myr. However the age of the high mass stars depends on the adopted models and is 3.5 +/- 0.5 Myr from the non- or moderately-rotating models of Brott et al. as well as the non-rotating models of Ekstr\"om et al. But the age is 4.0 -- 7.0 Myr if the rotating models of Ekstr\"om et al. are adopted. This latter age is in excellent agreement with the time scale of ejection of the high mass runaway star HD 153919 from NGC 6231, albeit the younger age cannot be entirely excluded.

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.

Detailed studies of resolved young massive star clusters are necessary to determine their dynamical state and evaluate the importance of gas expulsion and early cluster evolution. In an effort to gain insight into the dynamical state of the young massive cluster R136 and obtain the first measurement of its velocity dispersion, we analyse multi-epoch spectroscopic data of the inner regions of 30 Doradus in the Large Magellanic Cloud (LMC) obtained as part of the VLT-FLAMES Tarantula Survey. Following a quantitative assessment of the variability, we use the radial velocities of non-variable sources to place an upper limit of 6 km/s on the line-of-sight velocity dispersion of stars within a projected distance of 5 pc from the centre of the cluster. After accounting for the contributions of undetected binaries and measurement errors through Monte Carlo simulations, we conclude that the true velocity dispersion is likely between 4 and 5 km/s given a range of standard assumptions about the binary distribution. This result is consistent with what is expected if the cluster is in virial equilibrium, suggesting that gas expulsion has not altered its dynamics. We find that the velocity dispersion would be ~25 km/s if binaries were not identified and rejected, confirming the importance of the multi-epoch strategy and the risk of interpreting velocity dispersion measurements of unresolved extragalactic young massive clusters.

Although it has important ramifications for both the formation of star clusters and their subsequent dynamical evolution, rotation remains a largely unexplored characteristic of young star clusters (few Myr). Using multi-epoch spectroscopic data of the inner regions of 30 Doradus in the Large Magellanic Cloud (LMC) obtained as part of the VLT-FLAMES Tarantula Survey, we search for rotation of the young massive cluster R136. From the radial velocities of 36 apparently single O-type stars within a projected radius of 10 pc from the centre of the cluster, we find evidence, at the 95% confidence level, for rotation of the cluster as a whole. We use a maximum likelihood method to fit simple rotation curves to our data and find a typical rotational velocity of ~3 km/s. When compared to the low velocity dispersion of R136, our result suggests that star clusters may form with at least ~20% of the kinetic energy in rotation.

The presence of a nearby companion alters the evolution of massive stars in binary systems, leading to phenomena such as stellar mergers, X-ray binaries and gamma-ray bursts. Unambiguous constraints on the fraction of massive stars affected by binary interaction were lacking. We simultaneously measured all relevant binary characteristics in a sample of Galactic massive O stars and quantified the frequency and nature of binary interactions. Over seventy per cent of all massive stars will exchange mass with a companion, leading to a binary merger in one third of the cases. These numbers greatly exceed previous estimates and imply that binary interaction dominates the evolution of massive stars, with implications for populations of massive stars and their supernovae.

The O-type object 9 Sgr is a well-known synchrotron radio emitter. This feature is usually attributed to colliding-wind binary systems, but 9 Sgr was long considered a single star. We have conducted a long-term spectroscopic monitoring of this star to investigate its multiplicity and search for evidence for wind-wind interactions. Radial velocities are determined and analysed using various period search methods. Spectral disentangling is applied to separate the spectra of the components of the binary system. We derive the first ever orbital solution of 9 Sgr. The system is found to consist of an O3.5 V((f+)) primary and an O5- 5.5 V((f)) secondary moving around each other on a highly eccentric (e = 0.7), 8.6 year orbit. The spectra reveal no variable emission lines that could be formed in the wind interaction zone in agreement with the expected properties of the interaction in such a wide system. Our results provide further support to the paradigm of synchrotron radio emission from early-type stars being a manifestation of interacting winds in a binary system.

Whether massive stars can occasionally form in relative isolation or if they require a large cluster of lower-mass stars around them is a key test in the differentiation of star formation theories as well as how the initial mass function of stars is sampled. Previous attempts to find O-type stars that formed in isolation were hindered by the possibility that such stars are merely runaways from clusters, i.e., their current isolation does not reflect their birth conditions. We introduce a new method to find O-type stars that are not affected by such a degeneracy. Using the VLT-FLAMES Tarantula Survey and additional high resolution imaging we have identified stars that satisfy the following constraints: 1) they are O-type stars that are not detected to be part of a binary system based on RV time series analysis; 2) they are designated spectral type O7 or earlier ; 3) their velocities are within 1\sigma of the mean of OB-type stars in the 30 Doradus region, i.e. they are not runaways along our line-of-sight; 4) the projected surface density of stars does not increase within 3 pc towards the O-star (no evidence for clusters); 5) their sight lines are associated with gaseous and/or dusty filaments in the ISM, and 6) if a second candidate is found in the direction of the same filament with which the target is associated, both are required to have similar velocities. With these criteria, we have identified 15 stars in the 30 Doradus region, which are strong candidates for being high-mass stars that have formed in isolation. Additionally, we employed extensive MC stellar cluster simulations to confirm that our results rule out the presence of clusters around the candidates. Eleven of these are classified as Vz stars, possibly associated with the zero-age main sequence. We include a newly discovered W-R star as a candidate, although it does not meet all of the above criteria.

Using optical long baseline interferometry, we resolved for the first time the two wide components of HD167971, a candidate hierarchical triple system known to efficiently accelerate particles. Our multi-epoch VLTI observations provide direct evidence for a gravitational link between the O8 supergiant and the close eclipsing O + O binary. The separation varies from 8 to 15 mas over the three-year baseline of our observations, suggesting that the components evolve on a wide and very eccentric orbit (most probably e>0.5). These results provide evidence that the wide orbit revealed by our study is not coplanar with the orbit of the inner eclipsing binary. From our measurements of the near-infrared luminosity ratio, we constrain the spectral classification of the components in the close binary to be O6-O7, and confirm that these stars are likely main-sequence objects. Our results are discussed in the context of the bright non-thermal radio emission already reported for this system, and we provide arguments in favour of a maximum radio emission coincident with periastron passage. HD167971 turns out to be an efficient O-type particle accelerator that constitutes a valuable target for future high angular resolution radio imaging using VLBI facilities.

Context. HD150136 is one of the nearest systems harbouring an O3 star. Although this system was for a long time considered as binary, more recent investigations have suggested the possible existence of a third component. Aims. We present a detailed analysis of HD 150136 to confirm the triple nature of this system. In addition, we investigate the physical properties of the individual components of this system. Methods. We analysed high-resolution, high signal-to-noise data collected through multi-epoch runs spread over ten years. We applied a disentangling program to refine the radial velocities and to obtain the individual spectra of each star. With the radial velocities, we computed the orbital solution of the inner system, and we describe the main properties of the orbit of the outer star such as the preliminary mass ratio, the eccentricity, and the orbital-period range. With the individual spectra, we determined the stellar parameters of each star by means of the CMFGEN atmosphere code. Results. We offer clear evidence that HD 150136 is a triple system composed of an O3V((f\ast))-3.5V((f+)), an O5.5-6V((f)), and an O6.5-7V((f)) star. The three stars are between 0-3 Myr old. We derive dynamical masses of about 64, 40, and 35 Msun for the primary, the secondary and the third components by assuming an inclination of 49{\deg}. It currently corresponds to one of the most massive systems in our galaxy. The third star moves with a period in the range of 2950 to 5500 d on an outer orbit with an eccentricity of at least 0.3. This discovery makes HD 150136 the first confirmed triple system with an O3 primary star. However, because of the long orbital period, our dataset is not sufficient to constrain the orbital solution of the tertiary component with high accuracy.

As part of the VLT/X-shooter science verification, we obtained the first optical medium-resolution spectrum of a previously identified bright O-type object in NGC55, an LMC-like galaxy at a distance of \sim2.0 Mpc. Based on the stellar and nebular spectrum, we investigate the nature and evolutionary status of the central object(s) and its influence on the surrounding interstellar medium. We conclude that the source, NGC55_C1_31, is a composite object, likely a stellar cluster, which contains one or several hot (T_eff \simeq 50000 K) WN stars with a high mass-loss rate (\sim3 \times 10^{-5} M_\odot yr^{-1}) and a helium-rich composition (N_He/N_H = 0.8). The visual flux is dominated by OB-type (super)giant stars with T_eff \sim< 35000 K, solar helium abundance (N_He/N_H = 0.1), and mass-loss rate \sim2 \times 10^{-6} M_\odot yr^{-1}. The surrounding H II region has an electron density n_e < 10^2 cm^{-3} and an electron temperature T(OIII) \simeq 11500 \pm 600 K. The oxygen abundance of this region is [O/H] = 8.18 \pm 0.03 which corresponds to Z = 0.31 \pm 0.04 Z_\odot. We observed no significant gradients in T(OIII), n_e or [O/H] on a scale of 73 pc extending in four directions from the ionising source. The properties of the HII region can be reproduced by a CLOUDY model which uses the central cluster as ionising source, thus providing a self-consistent interpretation of the data. We also report on the serendipitous discovery of HeII nebular emission associated with the nearby source NGC55_C2_35, a feature usually associated with strong X-ray sources.

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.

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

As the brightest O-type X-ray source in the Carina nebula, HD 93250 (O4 III(fc)) is X-ray overluminous for its spectral type and has an unusually hard X-ray spectrum. Two different scenarios have been invoked to explain its X-ray properties: wind-wind interaction and magnetic wind confinement. Yet, HD 93250 shows absolutely constant radial velocities over time scales of years suggesting either a single star, a binary system seen pole-one or a very long period and/or highly eccentric system. Using the ESO Very Large Telescope Interferometer, we resolved HD 93250 as a close pair with similar components. We measured a near-infrared flux ratio of 0.8+/-0.1 and a separation of 1.5+/-0.2 x 10E-03 arcsec. At the distance of Carina, this corresponds to a projected physical distance of 3.5 A.U. While a quantitative investigation would require a full characterization of the orbit, the binary nature of HD 93250 allows us to qualitatively explain both its X-ray flux and hardness and its non-thermal radio emission in the framework of a colliding wind scenario. We also discuss various observational biases. We show that, due to line-blending of two similar spectral components, HD 93250 could have a period as short as one to several years despite the lack of measurable radial velocity variations.

We present the earliest results of our NACO/VLT sparse aperture masking (SAM) campaign to search for binarity in a sample of 60 O-type stars. We detect Delta Ks < 5 mag companions for 20-25% of our targets with separations in the range 30-100 mas (typically, 40 - 200 A.U.). Most of these companions were unknown, shedding thus new light on the multiplicity properties of massive stars in a separation and brightness regime that has been difficult to explore so far. Adding detections from other techniques (spectroscopy, interferometry, speckle, lucky imaging, AO), the fraction of O stars with at least one companion is 85% (51/60 targets). This is the largest multiplicity fraction ever found.

We present a spectroscopic analysis of VLT/X-Shooter observations of six O-type stars in the low-metallicity (Z ~ 1/7 Z\odot) galaxies IC 1613, WLM and NGC 3109. The stellar and wind parameters of these sources allow us, for the first time, to probe the mass-loss versus metallicity dependence of stellar winds below that of the Small Magellanic Cloud (at Z ~ 1/5Z\odot) by means of a modified wind momentum versus luminosity diagram. The wind strengths that we obtain for the objects in WLM and NGC 3109 are unexpectedly high and do not agree with theoretical predictions. The objects in IC 1613 tend towards a higher than expected mass-loss rate, but remain consistent with predictions within their error bars. We discuss potential systematic uncertainties in the mass-loss determinations to explain our results. However, if further scrutinization of these findings point towards an intrinsic cause for this unexpected sub-SMC mass-loss behavior, implications would include a higher than anticipated number of Wolf-Rayet stars and Ib/Ic supernovae in low-metallicity environments, but a reduced number of long-duration gamma-ray bursts produced through a single-star evolutionary channel.

Using an extended set of multi-epoch high-resolution high signal-to-noise ratio optical spectra, we readdress the multiplicity properties of the O-type stars in IC 2944 and in the Cen OB2 association. We present new evidence of binarity for five objects and we confirm the multiple nature of another two. We derive the first orbital solutions for HD 100099, HD 101436 and HD 101190 and we provide additional support for HD 101205 being a quadruple system. The minimal spectroscopic binary fraction in our sample is f_min = 0.57. Using numerical simulations, we show that the detection rate of our observational campaign is close to 90 per cent, leaving thus little room for undetected spectroscopic binary systems. The statistical properties of the O-star population in IC 2944 are similar, within the uncertainties, to the results obtained in the earlier papers in this series despite the fact that sample size effects limit the significance of the comparison. Using newly derived spectroscopic parallaxes, we reassess the distance to IC 2944 and obtained 2.3+/-0.3 kpc, in agreement with previous studies. We also confirm that, as far as the O stars are concerned, the IC 2944 cluster is most likely a single entity.

The distance of the very young open cluster Westerlund 2, which contains the very massive binary system WR 20a and is likely associated with a TeV source, has been the subject of much debate. We attempt a joint analysis of spectroscopic and photometric data of eclipsing binaries in the cluster to constrain its distance. A sample of 15 stars, including three eclipsing binaries (MSP 44, MSP 96, and MSP 223) was monitored with the FLAMES multi-object spectrograph. The spectroscopic data are analysed together with existing BV photometry. The analysis of the three eclipsing binaries clearly supports the larger values of the distance, around 8 kpc, and rules out values of about 2.4 - 2.8 kpc that have been suggested in the literature. Furthermore, our spectroscopic monitoring reveals no clear signature of binarity with periods shorter than 50 days in either the WN6ha star WR 20b, the early O-type stars MSP 18, MSP 171, MSP 182, MSP 183, MSP 199, and MSP 203, or three previously unknown mid O-type stars. The only newly identified candidate binary system is MSP 167. The absence of a binary signature is especially surprising for WR 20b and MSP 18, which were previously found to be bright X-ray sources. The distance of Westerlund 2 is confirmed to be around 8 kpc as previously suggested based on the spectrophotometry of its population of O-type stars and the analysis of the light curve of WR 20a. Our results suggest that short-period binary systems are not likely to be common, at least not among the population of O-type stars in the cluster.

VFTS 682 is located in an active star-forming region, at a projected distance of 29 pc from the young massive cluster R136 in the Tarantula Nebula of the Large Magellanic Cloud. It was previously reported as a candidate young stellar object, and more recently spectroscopically revealed as a hydrogen-rich Wolf-Rayet (WN5h) star. Our aim is to obtain the stellar properties, such as its intrinsic luminosity, and to investigate the origin of VFTS 682. To this purpose, we model optical spectra from the VLT-FLAMES Tarantula Survey with the non-LTE stellar atmosphere code CMFGEN, as well as the spectral energy distribution from complementary optical and infrared photometry. We find the extinction properties to be highly peculiar (RV ~4.7), and obtain a surprisingly high luminosity log(L/Lsun) = 6.5 \pm 0.2, corresponding to a present-day mass of ~150Msun. The high effective temperature of 52.2 \pm 2.5kK might be explained by chemically homogeneous evolution - suggested to be the key process in the path towards long gamma-ray bursts. Lightcurves of the object show variability at the 10% level on a timescale of years. Such changes are unprecedented for classical Wolf-Rayet stars, and are more reminiscent of Luminous Blue Variables. Finally, we discuss two possibilities for the origin of VFTS 682: (i) the star either formed in situ, which would have profound implications for the formation mechanism of massive stars, or (ii) VFTS 682 is a slow runaway star that originated from the dense cluster R136, which would make it the most massive runaway known to date.

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

We report the discovery that R139 in 30 Doradus is a massive spectroscopic binary system. Multi-epoch optical spectroscopy of R139 was obtained as part of the VLT-FLAMES Tarantula Survey, revealing a double-lined system. The two components are of similar spectral types; the primary exhibits strong C III 4650 emission and is classified as an O6.5 Iafc supergiant, while the secondary is an O6 Iaf supergiant. The radial-velocity variations indicate a highly eccentric orbit with a period of 153.9 days. Photometry obtained with the Faulkes Telescope South shows no evidence for significant variability within an 18 month period. The orbital solution yields lower mass limits for the components of M1sin^3 i = 78 \pm 8 Msun and M2sin^3 i = 66 \pm 7 Msun. As R139 appears to be the most massive binary system known to contain two evolved Of supergiants, it will provide an excellent test for atmospheric and evolutionary models.

We present observations of the embedded massive young stellar object (YSO) candidate 08576nr292, obtained with X-shooter and SINFONI on the ESO Very Large Telescope (VLT). The flux-calibrated, medium-resolution X-shooter spectrum (300 - 2500 nm) includes over 300 emission lines, but no (photospheric) absorption lines and is consistent with a reddened disk spectrum. Among the emission lines are three hydrogen series and helium lines, both permitted and forbidden metal lines, and CO first-overtone emission. A representative sample of lines with different morphologies is presented. The H alpha and Ca II triplet lines are very strong, with profiles indicative of outflow and - possibly - infall, usually observed in accreting stars. These lines include a blue-shifted absorption component at ~ -125 km/s. The He I and metal-line profiles are double-peaked, with a likely origin in a circumstellar disk. The forbidden lines, associated with outflow, have a single blue-shifted emission component centered at -125 km/s, coinciding with the absorption components in H alpha and Ca II. SINFONI H- and K-band integral-field spectroscopy of the cluster environment demonstrates that the [Fe II] emission is produced by a jet originating at the location of 08576nr292. Because the spectral type of the central object cannot be determined, its mass remains uncertain. We argue that 08576nr292 is an intermediate-mass YSO with a high accretion rate Mdot_acc ~ 10^-5 - 10^-6 Msun/yr. These observations demonstrate the potential of X-shooter and SINFONI to study in great detail an accretion disk-jet system, rarely seen around the more massive YSOs.

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.

Aims: The present study was conducted to determine the optical extinction curve for Cerro Paranal under typical clear-sky observing conditions, with an accuracy of 0.01 mag/airmass. Methods: The extinction curve of Paranal was obtained through low-resolution spectroscopy of 8 spectrophotometric standard stars observed with FORS1 mounted at the 8.2 m Very Large Telescope, covering a spectral range 3300-8000 A. A total of 600 spectra were collected on more than 40 nights distributed over six months, from October 2008 to March 2009. The average extinction curve was derived using a global fit algorithm, which allowed us to simultaneously combine all the available data. The main atmospheric parameters were retrieved using the LBLRTM radiative transfer code, which was also utilised to study the impact of variability of the main molecular bands of O2, O3, and H2O, and to estimate their column densities. Results: In general, the extinction curve of Paranal appears to conform to those derived for other astronomical sites in the Atacama desert, like La Silla and Cerro Tololo. However, a systematic deficit with respect to the extinction curve derived for Cerro Tololo before the El Chichon eruption is detected below 4000 A. We attribute this downturn to a non standard aerosol composition, probably revealing the presence of volcanic pollutants above the Atacama desert. An analysis of all spectroscopic extinction curves obtained since 1974 shows that the aerosol composition has been evolving during the last 35 years. The persistence of traces of non meteorologic haze suggests the effect of volcanic eruptions, like those of El Chichon and Pinatubo, lasts several decades. The usage of the standard CTIO and La Silla extinction curves implemented in IRAF and MIDAS produce systematic over/under-estimates of the absolute flux.

We present multi-instruments and multi-wavelengths observations of the famous LBV star Pistol Star. These observations are part of a larger program about early O stars at different metallicities. The Pistol star has been claimed as the most massive star known, with 250 solar masses. We present the preliminary results based on X-Shooter spectra, as well as the observations performed with the VLTI-AMBER and the VLT-NACO adaptive optics. The X-shooter spectrograph allows to obtain simultaneously a spectrum from the UV to the K-band with a resolving power of $\sim$15000. The preliminary results obtained indicate that Pistol Star has similar properties of Eta Car, including shells of matter, but also the binarity. Other objects of the program, here briefly presented, were selected for their particular nature: early O stars with mass discrepancies between stellar evolution models and observations, discrepancies with the wind momentum luminosity relation.

First results of near-IR adaptive optics (AO)-assisted imaging, interferometry, and spectroscopy of this Luminous Blue Variable (LBV) are presented. They suggest that the Pistol Star is at least double. If the association is physical, it would reinforce questions concerning the importance of multiplicity for the formation and evolution of extremely massive stars.

Binaries are excellent astrophysical laboratories that provide us with direct measurements of fundamental stellar parameters. Compared to single isolated star, multiplicity induces new processes, offering the opportunity to confront our understanding of a broad range of physics under the extreme conditions found in, and close to, astrophysical objects. In this contribution, we will discuss the parameter space occupied by massive binaries, and the observational means to investigate it. We will review the multiplicity fraction of OB stars within each regime, and in different astrophysical environments. In particular we will compare the O star spectroscopic binary fraction in nearby open clusters and we will show that the current data are adequately described by an homogeneous fraction of f~0.44. We will also summarize our current understanding of the observed parameter distributions of O+OB spectroscopic binaries. We will show that the period distribution is overabundant in short period binaries and that it can be described by a bi-modal Oepik law with a break point around P~10d. The distribution of the mass-ratios shows no indication for a twin population of equal mass binaries and seems rather uniform in the range 0.2< q=M_2/M_1<1.0.

We present the first ultraviolet (UV) and multi-epoch optical spectroscopy of 30 Dor 016, a massive O2-type star on the periphery of 30 Doradus in the Large Magellanic Cloud. The UV data were obtained with the Cosmic Origins Spectrograph on the Hubble Space Telescope as part of the Servicing Mission Observatory Verification program after Servicing Mission 4, and reveal #016 to have one of the fastest stellar winds known. From analysis of the CIV 1548-51 doublet we find a terminal velocity, v_infty=3450 +/- 50km/s. Optical spectroscopy is from the VLT-FLAMES Tarantula Survey, from which we rule out a massive companion (with 2d<P<1yr) to a confidence of 98%. The radial velocity of #016 is offset from the systemic value by -85km/s, suggesting that the star has traveled the 120pc from the core of 30 Doradus as a runaway, ejected via dynamical interactions.

We present adaptive optics (AO) near-infrared observations of the core of the Tr 14 cluster in the Carina region obtained with the ESO multi-conjugate AO demonstrator, MAD. Our campaign yields AO-corrected observations with an image quality of about 0.2 arcsec across the 2 arcmin field of view, which is the widest AO mosaic ever obtained. We detected almost 2000 sources spanning a dynamic range of 10 mag. The pre-main sequence (PMS) locus in the colour-magnitude diagram is well reproduced by Palla & Stahler isochrones with an age of 3 to 5 1E+05 yr, confirming the very young age of the cluster. We derive a very high (deprojected) central density n0~4.5(+/-0.5) \times 10^4 pc^-3 and estimate the total mass of the cluster to be about ~4.3^{+3.3}_{-1.5} \times 10^3 Msun, although contamination of the field of view might have a significant impact on the derived mass. We show that the pairing process is largely dominated by chance alignment so that physical pairs are difficult to disentangle from spurious ones based on our single epoch observation. Yet, we identify 150 likely bound pairs, 30% of these with a separation smaller than 0.5 arcsec (~1300AU). We further show that at the 2-sigma level massive stars have more companions than lower-mass stars and that those companions are respectively brighter on average, thus more massive. Finally, we find some hints of mass segregation for stars heavier than about 10 Msun. If confirmed, the observed degree of mass segregation could be explained by dynamical evolution, despite the young age of the cluster.

The Bulge is the least understood major stellar population of the Milky Way. Most of what we know about the formation and evolution of the Bulge comes from bright giant stars. The underlying assumption that giants represent all the stars, and accurately trace the chemical evolution of a stellar population, is under debate. In particular, recent observations of a few microlensed dwarf stars give a very different picture of the evolution of the Bulge from that given by the giant stars. [ABRIDGED] We perform a detailed elemental abundance analysis of dwarf stars in the Galactic bulge, based on high-resolution spectra that were obtained while the stars were optically magnified during gravitational microlensing events. [ABRIDGED] We present detailed elemental abundances and stellar ages for six new dwarf stars in the Galactic bulge. Combining these with previous events, here re-analysed with the same methods, we study a homogeneous sample of 15 stars, which constitute the largest sample to date of microlensed dwarf stars in the Galactic bulge. We find that the stars span the full range of metallicities from [Fe/H]=-0.72 to +0.54, and an average metallicity of <[Fe/H]>=-0.08+/-0.47, close to the average metallicity based on giant stars in the Bulge. Furthermore, the stars follow well-defined abundance trends, that for [Fe/H]<0 are very similar to those of the local Galactic thick disc. This suggests that the Bulge and the thick disc have had, at least partially, comparable chemical histories. At sub-solar metallicities we find the Bulge dwarf stars to have consistently old ages, while at super-solar metallicities we find a wide range of ages. Using the new age and abundance results from the microlensed dwarf stars we investigate possible formation scenarios for the Bulge.

With an orbital period of about 3.078 days, the double-lined spectroscopic binary HD 48099 is, still now, the only short-period O+O system known in the Mon OB2 association. Even though an orbital solution has already been derived for this system, few information are available about the individual stars. We present, in this paper, the results of a long-term spectroscopic campaign. We derive a new orbital solution and apply a disentangling method to recover the mean spectrum of each star. To improve our knowledge concerning both components, we determine their spectral classifications and their projected rotational velocities. We also constrain the main stellar parameters of both stars by using the CMFGEN atmosphere code and provide the wind properties for the primary star through the study of IUE spectra. This investigation reveals that HD 48099 is an O5.5 V((f))+O9 V binary with M_1 sin^3 i = 0.70 M_{\sun} and M_2 sin^3 i = 0.39 M_{\sun}, implying a rather low orbital inclination. This result, combined with both a large effective temperature and log g, suggests that the primary star (v sini ~ 91 km s^-1) is actually a fast rotator with a strongly clumped wind and a nitrogen abundance of about 8 times the solar value.

Based on a set of over 100 medium- to high-resolution optical spectra collected from 2003 to 2009, we investigate the properties of the O-type star population in NGC6611 in the core of the Eagle Nebula (M16). Using a much more extended data set than previously available, we revise the spectral classification and multiplicity status of the nine O-type stars in our sample. We confirm two suspected binaries and derive the first SB2 orbital solutions for two systems. We further report that two other objects are displaying a composite spectrum, suggesting possible long-period binaries. Our analysis is supported by a set of Monte-Carlo simulations, allowing us to estimate the detection biases of our campaign and showing that the latter do not affect our conclusions. The absolute minimal binary fraction in our sample is f_min=0.44 but could be as high as 0.67 if all the binary candidates are confirmed. As in NGC6231 (see Paper I), up to 75% of the O star population in NGC6611 are found in an O+OB system, thus implicitly excluding random pairing from a classical IMF as a process to describe the companion association in massive binaries. No statistical difference could be further identified in the binary fraction, mass-ratio and period distributions between NGC6231 and NGC6611, despite the difference in age and environment of the two clusters.

Many young extra-galactic clusters have a measured velocity dispersion that is too high for the mass derived from their age and total luminosity, which has led to the suggestion that they are not in virial equilibrium. Most of these clusters are confined to a narrow age range centred around 10 Myr because of observational constraints. At this age the cluster light is dominated by luminous evolved stars, such as red supergiants, with initial masses of ~13-22 Msun for which (primordial) binarity is high. In this study we investigate to what extent the observed excess velocity dispersion is the result of the orbital motions of binaries. We demonstrate that estimates for the dynamical mass of young star clusters, derived from the observed velocity dispersion, exceed the photometric mass by up-to a factor of 10 and are consistent with a constant offset in the square of the velocity dispersion. This can be reproduced by models of virialised star clusters hosting a massive star population of which ~25 is in binaries, with typical mass ratios of ~0.6 and periods of ~1000 days. We conclude that binaries play a pivotal role in deriving the dynamical masses of young (~10 Myr) moderately massive and compact (<1e5 Msun; > 1 pc) star clusters.

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 elemental abundances of 13 microlensed dwarf and subgiant stars in the Galactic bulge, which constitute the largest sample to date. We show that these stars span the full range of metallicity from Fe/H=-0.8 to +0.4, and that they follow well-defined abundance trends, coincident with those of the Galactic thick disc.

While it is well known that most massive stars are found to be part of binary or multiple systems, an accurate characterization of the statistical properties of these multiple objects is still lacking. In the present talk, we will review the current status of the field, emphasizing the need of using complementarity techniques to cover the large parameter space. We will also describe what we think is the place of interferometry in this context.

We present the results of a spectroscopic monitoring campaign devoted to three O-type stars in the Carina nebula. We derive the full SB2 orbital solution of the binary system Tr 16-112, an exceptional dissymmetrical system consisting of an O5.5-6V((f$^+$?p)) primary and a B2V-III secondary. We also report on low-amplitude brightness variations in Tr 16-112 that are likely due to the ellipsoidal shape of the O5.5-6 primary revolving in an eccentric orbit around the system's centre of mass. We detect for the first time a clear SB2 binary signature in the spectrum of HD 93343 (O8 + O8), although our data are not sufficient to establish an orbital solution. This system also displays low amplitude photometric modulations. On the other hand, no indication of multiplicity is found in the optical spectra of HD 93250. Finally, we discuss the general properties of multiple massive stars in the Carina OB1 association.

HD 152219 is a massive binary system with O9.5 III + B1-2 V/III components and a short orbital period of 4.2 d. Its primary component further displays clear line profile variability (LPV). The primary component being located within the pulsational instability domain predicted for high-luminosity stars, we previously suggested that the observed LPV could be associated with non-radial pulsations. The aim of the present work is to determine the nature of the observed LPV in the spectrum of the primary component of HD 152219.During a 4-night FEROS monitoring campaign, we collected a new set of 134 high signal-to-noise spectra. These new observations were then used to re-investigate the variability of different line profiles in the spectrum of HD 152219. Based on the present analysis, we discard the non-radial pulsations and point out the Rossiter-McLaughlin effect as the cause of the LPV in HD 152219. The upper limit on the amplitude of possible weak pulsations is set at a few parts per thousand of the continuum level.

Aims. We present the results obtained from a long-term spectroscopic campaign devoted to the multiplicity of O-type stars in the young open cluster NGC2244 and in the Mon OB2 association. Methods. Our spectroscopic monitoring was performed over several years, allowing us to probe different time-scales. For each star, several spectral diagnostic tools are applied, in order to search for line shifts and profile variations. We also measure the projected rotational velocity and revisit the spectral classification. Results. In our sample, several stars were previously considered as spectroscopic binaries, though only a few scattered observations were available. Our results now reveal a more complex situation. Our study identifies two new spectroscopic binaries (HD46149 in NGC2244 and HD46573 in MonOB2). The first object is a long-period double-lined spectroscopic binary, though the exact value of its period remains uncertain and the second object is classified as an SB1 system with a period of about 10.67 days but the time series of our observations do not enable us to derive a unique orbital solution for this system. We also classify another star as variable in radial velocity (HD46150) and we detect line profile variations in two rapid rotators (HD46056 and HD46485). Conclusions. This spectroscopic investigation places a firm lower limit (17%) on the binary fraction of O-stars in NGC2244 and reveals the lack of short-period O+OB systems in this cluster. In addition, a comparison of these new results with two other well-studied clusters (NGC6231 and IC1805) puts forward possible hints of a relation between stellar density and binarity, which could provide constraints on the theories about the formation and early evolution of hot stars.

Several colliding-wind massive binaries are known to be non-thermal emitters in the radio domain. This constitutes strong evidence for the fact that an efficient particle acceleration process is at work in these objects. The acceleration mechanism is most probably the Diffusive Shock Acceleration (DSA) process in the presence of strong hydrodynamic shocks due to the colliding-winds. In order to investigate the physics of this particle acceleration, we initiated a multiwavelength campaign covering a large part of the electromagnetic spectrum. In this context, the detailed study of the hard X-ray emission from these sources in the SIMBOL-X bandpass constitutes a crucial element in order to probe this still poorly known topic of astrophysics. It should be noted that colliding-wind massive binaries should be considered as very valuable targets for the investigation of particle acceleration in a similar way as supernova remnants, but in a different region of the parameter space.

We present here the analysis of an extensive set of high resolution optical spectra of Plaskett's star (HD 47129). We use a disentangling method to separate the individual spectra of each star. We derive a new orbital solution and discuss the spectral classification of both components. A Doppler tomography technique applied to the emission lines H alpha and He II 4686 yields a Doppler map that illustrates the wind interactions in the system. Finally, an atmosphere code is used to determine the different chemical abundances of the system components and the wind parameters. HD 47129 appears to be an O8 III/I + O7.5 III binary system in a post RLOF evolutionary stage, where matter has been transferred from the primary to the secondary star. The He overabundance of the secondary supports this scenario. In addition, the N overabundance and C underabundance of the primary component confirm previous results based on X-ray spectroscopy and indicate that the primary is an evolved massive star. Furthermore, the secondary star has a large rotational velocity that deforms its surface, leading to a non-uniform distribution in effective temperature. This could explain the variations in the equivalent widths of the secondary lines with phase. We suggest that the wind of the secondary star is confined near the equatorial plane because of its high rotational velocity, affecting the ram pressure equilibrium in the wind interaction zone.

The X-ray spectra of late type stars can generally be well fitted by a two temperature component model of the corona. We fnd that the temperature of both components are strong functions of stellar age, although the temperature of the hotter plasma in the corona shows a larger scatter and is probably affected by the activity of stars, such as flares. We confirm the power-law decay of the temperature of the hot plasma, but the temperature of the cool component decays linearly with log (age).

We present the results of a long-term high-resolution spectroscopy campaign on the O-type stars in NGC 6231. We revise the spectral classification and multiplicity of these objects and we constrain the fundamental properties of the O-star population. Almost three quarters of the O-type stars in the cluster are members of a binary system. The minimum binary fraction is 0.63, with half the O-type binaries having an orbital period of the order of a few days. The eccentricities of all the short-period binaries are revised downward, and henceforth match a normal period-eccentricity distribution. The mass-ratio distribution shows a large preference for O+OB binaries, ruling out the possibility that, in NGC 6231, the companion of an O-type star is randomly drawn from a standard IMF. Obtained from a complete and homogeneous population of O-type stars, our conclusions provide interesting observational constraints to be confronted with the formation and early-evolution theories of O stars.

We present the results of an optical and X-ray monitoring campaign on the short-period massive SB2 binary HD 152218. Combining our HiRes spectroscopic data with previous observations, we unveil the contradictions between the published orbital solutions. In particular, we solve the aliasing on the period and derive a value close to 5.604 d. Our eccentricity e = 0.259 +/- 0.006 is slightly lower than previously admitted. We show that HD 152218 is probably undergoing a relatively rapid apsidal motion of about 3deg/yr and we confirm the O9IV + O9.7V classification. We derive minimal masses of 15.82 +/- 0.26 Msol operator and 12.00 +/- 0.19 Msol operator and constrain the radius of the components to R1 = 10.3 +/- 1.3 Rsol and R2 = 7.8 +/- 1.7 Rsol. We also report the results of an XMM-Newton monitoring of the HD 152218 X-ray emission throughout its orbital motion. The averaged X-ray spectrum is relatively soft and it is well reproduced by a 2-T optically thin thermal plasma model with component temperatures about 0.3 and 0.7 keV. The system presents an increase of its X-ray flux by about 30% near apastron compared to periastron, which is interpreted as the signature of an ongoing wind-wind interaction process occurring within the wind acceleration region.

An increasing number of early-type (O and Wolf-Rayet) colliding wind binaries (CWBs) is known to accelerate particles up to relativistic energies. In this context, non-thermal emission processes such as inverse Compton (IC) scattering are expected to produce a high energy spectrum, in addition to the strong thermal emission from the shock-heated plasma. SIMBOL-X will be the ideal observatory to investigate the hard X-ray spectrum (above 10 keV) of these systems, i.e. where it is no longer dominated by the thermal emission. Such observations are strongly needed to constrain the models aimed at understanding the physics of particle acceleration in CWB. Such systems are important laboratories for investigating the underlying physics of particle acceleration at high Mach number shocks, and probe a different region of parameter space than studies of supernova remnants.

We present a spectroscopic analysis of four massive binary systems that are known or are good candidates to display the Struve-Sahade effect (defined as the apparent strengthening of the secondary spectrum of the binary when the star is approaching, and the corresponding weakening of the lines when it is receding). We use high resolution optical spectra to determine new orbital solutions and spectral types of HD 165052, HD 100213, HD 159176 and DH Cep. As good knowledge of the fundamental parameters of the considered systems is necessary to examine the Struve-Sahade effect. We then study equivalent width variations in the lines of both components of these binaries during their orbital cycle. In the case of these four systems, variations appear in the equivalent widths of some lines during the orbital cycle, but the definition given above can any longer be valid, since it is now clear that the effect modifies the primary spectrum as much as the secondary spectrum. Furthermore, the lines affected, and the way in which they are affected, depend on the considered system. For at least two of them (HD 100213 and HD 159176) these variations probably reflect the ellipsoidal variable nature of the system.

Aims: We analyze INTEGRAL-ISGRI data in order to probe the hard X-ray emission (above 20 keV) from point sources in the Cyg OB2 region and to investigate the putative non-thermal high-energy emission from early-type stars (Wolf-Rayet and O-type stars). Among the targets located in the field of view, we focus on the still unidentified EGRET source 3EG 2033+4118 that may be related to massive stars known to produce non-thermal emission in the radio domain, and on the wide colliding-wind binary WR 140. Methods: Using a large set of data obtained with the IBIS-ISGRI imager onboard INTEGRAL, we run the OSA software package in order to find point sources in the fully coded field of view of the instrument. Results: Our data do not allow the detection of a lower-energy counterpart of 3EG J2033+4118 nor of any other new point sources in the field of view, and we derive upper limits on the high-energy flux for a few targets: 3EG J2033+4118, TeV J2032+4130, WR140, WR146 and WR147. The results are discussed in the context of the multiwavelength investigation of these objects. Conclusions: The upper limits derived are valuable constraints for models aimed at understanding the acceleration of particles in non-thermal emitting massive stars, and of the still unidentified very-high gamma-ray source TeV J2032+4130.

We discuss the properties of the X-ray sources with faint optical counterparts in the very young open cluster NGC 6231. From their positions in the H-R diagram, we find that the bulk of these objects probably consists of low-mass pre-main sequence stars with masses in the range 0.3 to 3.0 M$_{\odot}$. The age distribution of these objects indicates that low-mass star formation in NGC 6231 started more than 10 Myr ago and culminated in a starburst-like event about 1 to 4 Myr ago when the bulk of the low-mass PMS stars as well as the massive cluster members formed. We find no evidence for a spatial age gradient that could point towards a sequential star formation process. Only a few X-ray sources have counterparts with a reddening exceeding the average value of the cluster or with infrared colours indicating the presence of a moderate near-IR excess. The X-ray spectra of the brightest PMS sources are best fitted by rather hard thermal plasma models and a significant fraction of these sources display flares in their light curve. The X-ray brightest flaring sources have decay times between 2 and 16 ks. The X-ray selected PMS stars in NGC 6231 have $\log{L_{\rm X}/L_{\rm bol}}$ values that increase strongly with decreasing bolometric luminosity and can reach a saturation level ($\log{L_{\rm X}/L_{\rm bol}} \sim -2.4$) for non-flaring sources and even more extreme values during flares.

Using a set of high-resolution spectra, we studied the physical and orbital properties of the O-type binary CPD-41 7733, located in the core of \ngc. We report the unambiguous detection of the secondary spectral signature and we derive the first SB2 orbital solution of the system. The period is 5.6815 +/- 0.0015 d and the orbit has no significant eccentricity. CPD-41 7733 probably consists of stars of spectral types O8.5 and B3. As for other objects in the cluster, we observe discrepant luminosity classifications while using spectroscopic or brightness criteria. Still, the present analysis suggests that both components display physical parameters close to those of typical O8.5 and B3 dwarfs. We also analyze the X-ray light curves and spectra obtained during six 30 ks XMM-Newton pointings spread over the 5.7 d period. We find no significant variability between the different pointings, nor within the individual observations. The CPD-41 7733 X-ray spectrum is well reproduced by a three-temperature thermal mekal model with temperatures of 0.3, 0.8 and 2.4 keV. No X-ray overluminosity, resulting e.g. from a possible wind interaction, is observed. The emission of CPD-41 7733 is thus very representative of typical O-type star X-ray emission.

Aims. The properties of the early-type stars in the core of the Westerlund2 cluster are examined in order to establish a link between the cluster and the very massive Wolf-Rayet binary WR20a as well as the HII complex RCW49. Methods. Photometric monitoring as well as spectroscopic observations of Westerlund2 are used to search for light variability and to establish the spectral types of the early-type stars in the cluster core. Results. The first light curves of the eclipsing binary WR20a in B and V filters are analysed and a distance of 8kpc is inferred. Three additional eclipsing binaries, which are probable late O or early B-type cluster members, are discovered, but none of the known early O-type stars in the cluster displays significant photometric variability above 1% at the 1-sigma level. The twelve brightest O-type stars are found to have spectral types between O3 and O6.5, significantly earlier than previously thought. Conclusions. The distance of the early-type stars in Westerlund2 is established to be in excellent agreement with the distance of WR20a, indicating that WR20a actually belongs to the cluster. Our best estimate of the cluster distance thus amounts to 8.0pm1.4kpc. Despite the earlier spectral types, the currently known population of early-type stars in Westerlund2 does not provide enough ionizing photons to account for the radio emission of the RCW49 complex. This suggests that there might still exist a number of embedded early O-stars in RCW49.

Since the discovery, with the EINSTEIN satellite, of strong X-ray emission associated with HD93162 (=WR25), this object has been predicted to be a colliding-wind binary system. However, radial-velocity variations that would prove the suspected binary nature have yet to be found. We spectroscopically monitored this object to investigate its possible variability to address this discordance. We compiled the largest available radial-velocity data set for this star to look for variations that might be due to binary motion. We derived radial velocities from spectroscopic data acquired mainly between 1994 and 2006, and searched these radial velocities for periodicities using different numerical methods. For the first time, periodic radial-velocity variations are detected. Our analysis definitively shows that the Wolf-Rayet star WR25 is an eccentric binary system with a probable period of about 208 days.

In this second paper, we pursue the analysis of the 180 ks XMM-Newton campaign towards the young open cluster NGC 6231 and we focus on its rich OB star population. We present a literature-based census of the OB stars in the field of view with more than one hundred objects, among which 30% can be associated with an X-ray source. All the O-type stars are detected in the X-ray domain as soft and reasonably strong emitters. In the 0.5-10.0 keV band, their X-ray luminosities scale with their bolometric luminosities as $\log L_\mathrm{X} - \log L_\mathrm{bol}=-6.912\pm0.153$. Such a scaling law holds in the soft (0.5-1.0 keV) and intermediate (1.0-2.5 keV) bands but breaks down in the hard band. While the two colliding wind binaries in our sample clearly deviate from this scheme, the remaining O-type objects show a very limited dispersion (40% or 20% according to whether `cool' dwarfs are included or not), much smaller than that obtained from previous studies. At our detection threshold and within our sample, the sole identified mechanism that produces significant modulations in the O star X-ray emission is related to wind interaction. The intrinsic X-ray emission of non-peculiar O-type stars seems thus constant for a given star and the level of its X-ray emission is accurately related to the its luminosity or, equivalently, to its wind properties. Among B-type stars, the detection rate is only about 25% in the sub-type range B0-B4 and remains mostly uniform throughout the different sub-populations while it drops significantly at later sub-types. The associated X-ray spectra are harder than those of O-type stars. Our analysis points towards the detected emission being associated with a physical PMS companion >... [see paper for the complete abstract]

This paper is the first of a series dedicated to the X-ray properties of the young open cluster NGC 6231. Our data set relies on an XMM-Newton campaign of a nominal duration of 180 ks and reveals that NGC 6231 is very rich in the X-ray domain too. Indeed, 610 X-ray sources are detected in the present field of view, centered on the cluster core. The limiting sensitivity of our survey is approximately 6E-15 erg/s/cm^2 but clearly depends on the location in the field of view and on the source spectrum. Using different existing catalogues, over 85% of the X-ray sources could be associated with at least one optical and/or infrared counterpart within a limited cross-correlation radius of 3 arcsec at maximum. The surface density distribution of the X-ray sources presents a slight N-S elongation. Once corrected for the spatial sensitivity variation of the EPIC instruments, the radial profile of the source surface density is well described by a King profile with a central density of about 8 sources per arcmin^2 and a core radius close to 3.1arcmin. The distribution of the X-ray sources seems closely related to the optical source distribution. The expected number of foreground and background sources should represent about 9% of the detected sources, thus strongly suggesting that most of the observed X-ray emitters are physically belonging to NGC 6231. Finally, beside a few bright but soft objects -- corresponding to the early-type stars of the cluster -- most of the sources are relatively faint (~5E-15 erg/s/cm^2) with an energy distribution peaked around 1.0-2.0 keV.

Considered as the core of the Sco OB1 association, the young open cluster NGC 6231 harbours a rich O-type star population. In 2001, the XMM-Newton satellite targeted the cluster for a nominal duration of about 180 ks. Thanks to the detector sensitivity, the EPIC cameras provided an unprecedented X-ray view of NGC 6231, revealing about 600 point-like sources. In this contribution, we review the main results that have been obtained thanks to this unprecedented data set. Concerning the O-type stars, we present the latest developments related to the so-called 'canonical' Lx-Lbol relation. The dispersion around this relation might actually be much smaller than previously thought. In our data set, the sole mechanism that yields a significant deviation from this scheme is wind interaction. It is also the sole mechanism that induces a significant variation of the early-type star X-ray flux. In a second part of this contribution, we probe the properties of the optically faint X-ray sources. Most of them are believed to be low mass pre-main sequence stars. Their analysis provides direct insight into the star formation history of the cluster.

We summarize the results of a multiwavelength observing campaign on the massive eclipsing binary CPD-41deg7742, another remarkable object in the young open cluster NGC 6231. Our campaign relies on high resolution echelle spectroscopy, narrow-band optical photometry, and XMM-Newton X-ray observations. Combined with the spectroscopic analysis, the light curve analysis provides a direct measurement of the masses and sizes of the system components. However, the most outstanding results come from the XMM-Newton observations. Our 180 ks campaign towards NGC 6231, and CPD-41deg7742, provides an unprecedented phase coverage of such a close early-type binary. The EPIC-MOS light curves almost fully cover the 2.4 day period of the system and the brightness of the object is sufficient to yield a time resolution as tight as 1 ks. The X-ray flux presents clear variations along the orbit, that we interpret as the signature of an unusual wind interaction. We indeed expect that, in this O+B system, the dominant primary wind crashes into the secondary surface, leading to a wind-photosphere interaction. As a strong support to our interpretation, we provide a geometrical model that associates an extra X-ray emission to the secondary inner surface. Though quite simple, the present model matches the main features of the X-ray light curve.

In the first paper of this series, we presented a detailed high-resolution spectroscopic study of CPD - 41degr7742, deriving for the first time an orbital solution for both components of the system. In this second paper, we focus on the analysis of the optical light curve and on recent XMM-Newton X-ray observations. In the optical, the system presents two eclipses, yielding an inclination i ~ 77degr. Combining the constraints from the photometry with the results of our previous work, we derive the absolute parameters of the system. We confirm that the two components of CPD - 41degr7742 are main sequence stars (O9 V + B1-1.5 V) with masses (M_1 ~ 18 Msol and M_2 ~ 10 Msol) and respective radii (R_1 ~ 7.5 Rsol and R_2 ~ 5.4 Rsol) close to the typical values expected for such stars. We also report an unprecedented set of X-ray observations that almost uniformly cover the 2.44-day orbital cycle. The X-ray emission from CPD - 41degr7742 is well described by a two-temperature thermal plasma model with energies close to 0.6 and 1.0 keV, thus slightly harder than typical early-type emission. The X-ray light curve shows clear signs of variability. The emission level is higher when the primary is in front of the secondary. During the high emission state, the system shows a drop of its X-ray emission that almost exactly matches the optical eclipse. We interpret the main features of the X-ray light curve as the signature of a wind-photosphere interaction, in which the overwhelming primary O9 star wind crashes into the secondary surface. Alternatively the light curve could result from a wind-wind interaction zone located near the secondary star surface. As a support to our interpretation, we provide a phenomenological geometric model that qualitatively reproduces the observed modulations of the X-ray emission.

We present the analysis of an extensive set of high resolution spectroscopic observations of HD93161, a visual binary with a separation of 2 arcsec. HD93161 A is a spectroscopic binary, with both components clearly detected throughout the orbit. The primary star is most probably of spectral type O8V, while the secondary is likely an O9V. We obtain the first orbital solution for this system, characterized by a period of 8.566\pm0.004 days. The minimum masses of the primary and secondary stars are 22.2\pm0.6 M_sol and 17.0\pm0.4 M_sol respectively. These values are quite large, suggesting a high inclination of the orbit. The second object, HD93161 B, displays an O6.5V((f)) spectral type and is thus slightly hotter than its neighbour. This star is at first sight single but presents radial velocity variations. We finally study HD93161 in the X-ray domain. No significant variability is detected. The X-ray spectrum is well described by a 2T model with kT_1~0.3keV and kT_2~0.7keV. The X-ray luminosity is rather moderate, without any large emission excess imputable to a wind interaction.

Based on a deep XMM-Newton observation of the young open cluster NGC 6231, we derive the main X-ray properties of its early-type star population. Among the 610 X-ray sources detected in the field, 42 are associated with early-type stars. We investigate their $L_\mathrm{X} / L_\mathrm{bol}$ relationship and we confirm the clear dichotomy between O- and B-type stars. The cut-off line between the two behaviours occurs at $L_\mathrm{bol} \sim 10^{38}$ erg s$^{-1}$ as previously proposed by \citet{BSC97}. The distinction between single and binary stars is not clear cut, except for the colliding wind system HD 152248. The X-ray detected B-stars in NGC 6231 appear to be more luminous than predicted from the \citeauthor{BSC97} relation. Though this suggests a bimodal distribution of the B-star X-ray emission, we caution however that these results might be biased by detection limits. Finally we investigate the X-ray variability of the detected sources and we find that about 40% of the X-ray emitters in the field of view present consistent signs of variability in the EPIC instruments. This fraction is much larger than previously thought. About one third of the early-type stars population, either singles or binaries, further shows variability. These variations do not seem to be exclusively related to binarity and could thus not be systematically interpreted in the framework of a wind-wind collision phenomenon.

We report the results of a multiwavelength study of the non-thermal radio emitter HD168112 (O5.5III(f^+)). The detailed analysis of two quasi-simultaneous XMM-Newton and VLA observations reveals strong variability of this star both in the X-ray and radio ranges. The X-ray observations separated by five months reveal a decrease of the X-ray flux of ~30%. The radio emission on the other hand increases by a factor 5-7 between the two observations obtained roughly simultaneously with the XMM-Newton pointings. The X-ray data reveal a hard emission that is most likely produced by a thermal plasma at kT ~2-3 keV while the VLA data confirm the non-thermal status of this star in the radio waveband. Comparison with archive X-ray and radio data confirms the variability of this source in both wavelength ranges over a yet ill defined time scale. The properties of HD168112 in the X-ray and radio domain point towards a binary system with a significant eccentricity and an orbital period of a few years. However, our optical spectra reveal no significant changes of the star's radial velocity suggesting that if HD168112 is indeed a binary, it must be seen under a fairly low inclination.

We analyse spectroscopic observations of WR20a revealing that this star is a massive early-type binary system with a most probable orbital period of \sim 3.675 days. Our spectra indicate that both components are most likely of WN6ha or O3If^*/WN6ha spectral type. The orbital solution for a period of 3.675 days yields extremely large minimum masses of 70.7 \pm 4.0 and 68.8 \pm 3.8 M_{\odot} for the two stars. These properties make WR20a a cornerstone system for the study of massive star evolution.

We report the first results of an XMM-Newton monitoring campaign of the NGC 6231 open cluster in the Sco OB1 association. This first paper focuses on the massive colliding wind binary HD 152248, which is the brightest X-ray source of the cluster. The campaign, with a total duration of 180 ksec, was split into six separate observations, following the orbital motion of HD 152248. The X-ray flux from this system presents a clear, asymmetric modulation with the phase and ranges from 0.73 to 1.18 10-12 erg s-1 cm-2 in the [0.5-10.0 keV] energy band. The maximum of the emission is reached slightly after apastron. The EPIC spectra are quite soft and peak around 0.8-0.9 keV. We characterize their shape using several combinations of mekal models and power-law spectra and we detect significant spectral variability in the [0.5-2.5 keV] energy band. We also perform 2-D hydrodynamical simulations using different sets of parameters that closely reproduce the physical and orbital configuration of the HD 152248 system at the time of the six XMM-Newton pointings. This allows a direct confrontation of the model predictions with the constraints deduced from the X-ray observations of the system. We show that the observed variation of the flux can be explained by a variation of the X-ray emission from the colliding wind zone, diluted by the softer X-ray contribution of the two O-type stars of the system. Our simulations also reveal that the interaction region of HD 152248 should be highly unstable, giving rise to shells of dense gas that are separated by low density regions. Finally, we perform a search for short-term variability in the light curves of the system and we show that trends are present within several of the 30 ksec exposures of our campaign. Further, most of these trends are in good .