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.

We present the most comprehensive RSG sample for the SMC up to now, including 1,239 RSG candidates. The initial sample is derived based on a source catalog for the SMC with conservative ranking. Additional spectroscopic RSGs are retrieved from the literature, as well as RSG candidates selected from the inspection of CMDs. We estimate that there are in total $\sim$ 1,800 or more RSGs in the SMC. We purify the sample by studying the infrared CMDs and the variability of the objects, though there is still an ambiguity between AGBs and RSGs. There are much less RSGs candidates ($\sim$4\%) showing PAH emission features compared to the Milky Way and LMC ($\sim$15\%). The MIR variability of RSG sample increases with luminosity. We separate the RSG sample into two subsamples ("risky" and "safe") and identify one M5e AGB star in the "risky" subsample. Most of the targets with large variability are also the bright ones with large MLR. Some targets show excessive dust emission, which may be related to previous episodic mass loss events. We also roughly estimate the total gas and dust budget produced by entire RSG population as $\rm \sim1.9^{+2.4}_{-1.1}\times10^{-6}~M_{\odot}/yr$ in the most conservative case. Based on the MIST models, we derive a linear relation between $T_{\rm eff}$ and observed $\rm J-K_S$ color with reddening correction for the RSG sample. By using a constant bolometric correction and this relation, the Geneva evolutionary model is compared with our RSG sample, showing a good agreement and a lower initial mass limit of $\sim$7 $\rm M_\odot $ for the RSG population. Finally, we compare the RSG sample in the SMC and the LMC. Despite the incompleteness of LMC sample in the faint end, the result indicates that the LMC sample always shows redder color (except for the $\rm IRAC1-IRAC2$ and $\rm WISE1-WISE2$ colors due to CO absorption) and larger variability than the SMC sample.

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.

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.

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.

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.

We present the concept of BlueMUSE, a blue-optimised, medium spectral resolution, panoramic integral field spectrograph based on the MUSE concept and proposed for the Very Large Telescope. With an optimised transmission down to 350 nm, a larger FoV (1.4 x 1.4 arcmin$^2$) and a higher spectral resolution compared to MUSE, BlueMUSE will open up a new range of galactic and extragalactic science cases allowed by its specific capabilities, beyond those possible with MUSE. For example a survey of massive stars in our galaxy and the Local Group will increase the known population of massive stars by a factor $>$100, to answer key questions about their evolution. Deep field observations with BlueMUSE will also significantly increase samples of Lyman-alpha emitters, spanning the era of Cosmic Noon. This will revolutionise the study of the distant Universe: allowing the intergalactic medium to be detected unambiguously in emission, enabling the study of the exchange of baryons between galaxies and their surroundings. By 2030, at a time when the focus of most of the new large facilities (ELT, JWST) will be on the infra-red, BlueMUSE will be a unique facility, outperforming any ELT instrument in the Blue/UV. It will have a strong synergy with ELT, JWST as well as ALMA, SKA, Euclid and Athena.

SN2010da/NGC 300 ULX-1 was first detected as a supernova impostor in May 2010 and was recently discovered to be a pulsating ultraluminous X-ray source. In this letter, we present VLT/X-shooter spectra of this source obtained in October 2018, covering the wavelength range 350-2300 nm. The $J$- and $H$-bands clearly show the presence of a red supergiant donor star that is best matched by a MARCS stellar atmosphere with $T_{\rm eff} = 3650 - 3900$ K and $\log(L_{\rm bol}/L_{\odot}) = 4.25\pm0.10$, which yields a stellar radius $R = 310 \pm 70 R_{\odot}$. To fit the full spectrum, two additional components are required: a blue excess that can be fitted either by a hot blackbody (T $\gtrsim 20,000$ K) or a power law (spectral index $\alpha \approx 4$) and is likely due to X-ray emission reprocessed in the outer accretion disk or the donor star; and a red excess that is well fitted by a blackbody with a temperature of $\sim 1100$ K, and is likely due to warm dust in the vicinity of SN2010da. The presence of a red supergiant in this system implies an orbital period of at least 0.8-2.1 years, assuming Roche lobe overflow. Given the large donor-to-compact object mass ratio, orbital modulations of the radial velocity of the red supergiant are likely undetectable. However, the radial velocity amplitude of the neutron star is large enough (up to 40-60 km s$^{-1}$) to potentially be measured in the future, unless the system is viewed at a very unfavorable inclination.

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.

We present a clean, magnitude-limited (IRAC1 or WISE1 $\leq$ 15.0 mag) multiwavelength source catalog for the SMC with 45,466 targets in total, with the purpose of building an anchor for future studies, especially for the massive star populations at low-metallicity. The catalog contains data in 50 different bands including 21 optical and 29 infrared bands, ranging from the ultraviolet to the far-infrared. Additionally, radial velocities and spectral classifications were collected from the literature, as well as infrared and optical variability statistics were retrieved from different projects. The catalog was essentially built upon a $1''$ crossmatching and a $3''$ deblending between the SEIP source list and Gaia DR2 photometric data. Further constraints on the proper motions and parallaxes from Gaia DR2 allowed us to remove the foreground contamination. We estimated that about 99.5\% of the targets in our catalog were most likely genuine members of the SMC. By using the evolutionary tracks and synthetic photometry from MIST and the theoretical $\rm J-K_S$ color cuts, we identified 1,405 RSG, 217 YSG and 1,369 BSG candidates in the SMC in five different CMDs, where attention should also be paid to the incompleteness of our sample. We ranked the candidates based on the intersection of different CMDs. A comparison between the models and observational data shows that the lower limit of initial mass for the RSGs population may be as low as 7 or even 6 $M_{\odot}$ and the RSG is well separated from the AGB population even at faint magnitude, making RSGs a unique population connecting the evolved massive and intermediate stars, since stars with initial mass around 6 to 8 $M_{\odot}$ are thought to go through a second dredge-up to become AGBs. We encourage the interested reader to further exploit the potential of our catalog.

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.

Context. Metal-poor massive stars are supposed to be progenitors of certain supernovae, gamma-ray bursts and compact object mergers, potentially contributing to the early epochs of the Universe with their strong ionizing radiation. However, they remain mainly theoretical as individual spectroscopic observations of such objects have rarely been carried out below the metallicity of the SMC. Aims. This work aims at exploring what our state-of-the-art theories of stellar evolution combined with those of stellar atmospheres predict about a certain type of metal-poor (0.02 Z$_{\odot}$) hot massive stars, the chemically homogeneously evolving ones, called TWUIN stars. Methods. Synthetic spectra corresponding to a broad range in masses (20-130 M$_{\odot}$) and covering several evolutionary phases from the zero-age main-sequence up to the core helium-burning stage were computed. Results. We find that TWUIN stars show almost no emission lines during most of their {core hydrogen-burning} lifetimes. Most metal lines are completely absent, including nitrogen. During their core helium-burning stage, lines switch to emission and even some metal lines (oxygen and carbon, but still almost no nitrogen) show up. Mass loss and clumping play a significant role in line-formation in later evolutionary phases, particularly during core helium-burning. Most of our spectra are classified as an early O type giant or supergiant, and we find Wolf-Rayet stars of type WO in the core helium-burning phase. Conclusions. An extremely hot, early O type star observed in a low-metallicity galaxy could be the outcome of chemically homogeneous evolution $-$ and therefore the progenitor of a long-duration gamma-ray burst or a type Ic supernova. TWUIN stars may play an important role in reionizing the Universe due to their being hot without showing prominent emission lines during the majority of their lifetimes.

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.

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]

We present a timing analysis of multiple XMM-Newton and NuSTAR observations of the ultra-luminous pulsar NGC 7793 P13 spread over its 65d variability period. We use the measured pulse periods to determine the orbital ephemeris, confirm a long orbital period with P_orb = 63.9 (+0.5,-0.6) d, and find an eccentricity of e <= 0.15. The orbital signature is imprinted on top of a secular spin-up, which seems to get faster as the source becomes brighter. We also analyse data from dense monitoring of the source with Swift and find an optical photometric period of 63.9 +/- 0.5 d and an X-ray flux period of 66.8 +/- 0.4 d. The optical period is consistent with the orbital period, while the X-ray flux period is significantly longer. We discuss possible reasons for this discrepancy, which could be due to a super-orbital period caused by a precessing accretion disk or an orbital resonance. We put the orbital period of P13 into context with the orbital periods implied for two other ultra-luminous pulsars, M82 X-2 and NGC 5907 ULX and discuss possible implications for the system parameters.

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.

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.

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

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

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.

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.

HD 50138 is a B[e] star surrounded by a large amount of circumstellar gas and dust. Its spectrum shows characteristics which may indicate either a pre- or a post-main-sequence system. Mapping the kinematics of the gas in the inner few au of the system contributes to a better understanding of its physical nature. We present the first high spatial and spectral resolution interferometric observations of the Br-gamma line of HD~50138, obtained with VLTI/AMBER. The line emission originates from a region more compact (up to 3 au) than the continuum-emitting region. Blue- and red-shifted emission originates from the two different hemispheres of an elongated structure perpendicular to the polarization angle. The velocity of the emitting medium decreases radially. An overall offset along the NW direction between the line- and continuum-emitting regions is observed. We compare the data with a geometric model of a thin Keplerian disk and a spherical halo on top of a Gaussian continuum. Most of the data are well reproduced by this model, except for the variability, the global offset and the visibility at the systemic velocity. The evolutionary state of the system is discussed; most diagnostics are ambiguous and may point either to a post-main-sequence or a pre-main-sequence nature.

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

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

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.

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.

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.

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.

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

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 Suzaku and the Rossi X-ray Timing Explorer, we have conducted a series of four simultaneous observations of the galactic black hole candidate Cyg X-1 in what were historically faint and spectrally hard low states. Additionally, all of these observations occurred near superior conjunction with our line of sight to the X-ray source passing through the dense phases of the focused wind from the mass donating secondary. One of our observations was also simultaneous with observations by the Chandra-High Energy Transmission Grating. These latter spectra are crucial for revealing the ionized absorption due to the secondary's focused wind. Such absorption is present and must be accounted for in all four spectra. These simultaneous data give an unprecedented view of the 0.8-300 keV spectrum of Cyg X-1, and hence bear upon both corona and X-ray emitting jet models of black hole hard states. Three models fit the spectra well: coronae with thermal or mixed thermal/non-thermal electron populations, and jets. All three models require a soft component that we fit with a low temperature disk spectrum with an inner radius of only a few tens of GM/c^2. All three models also agree that the known spectral break at 10\,keV is not solely due to the presence of reflection, but each gives a different underlying explanation for the augmentation of this break. Thus whereas all three models require that there is a relativistically broadened Fe line, the strength and inner radius of such a line is dependent upon the specific model, {thus making premature line-based estimates of the black hole spin in the Cyg X-1 system. We look at the relativistic line in detail, accounting for the narrow Fe emission and ionized absorption detected by HETG. Although the specific relativistic parameters of the line are continuum-dependent, none of the broad line fits allow for an inner disk radius that is >40 GM/c^2.