Search

Mass loss is one of the key parameters that determine stellar evolution. Despite the progress we have achieved over the last decades we still cannot match the observational derived values with theoretical predictions. Even worse, there are certain phases, such as the B[e] supergiants (B[e]SGs) and the Luminous Blue Variables (LBVs), where significant mass is lost through episodic or outburst activity. This leads to various structures forming around them that permit dust formation, making these objects bright IR sources. The ASSESS project aims to determine the role of episodic mass in the evolution of massive stars, by examining large numbers of cool and hot objects (such as B[e]SGs/LBVs). For this purpose, we initiated a large observation campaign to obtain spectroscopic data for ∼1000 IR-selected sources in 27 nearby galaxies. Within this project we successfully identified seven B[e] supergiants (one candidate) and four Luminous Blue Variables of which six and two, respectively, are new discoveries. We used spectroscopic, photometric, and light curve information to better constrain the nature of the reported objects. We particularly noted the presence of B[e]SGs at metallicity environments as low as 0.14 Z⊙.

The ERC «ASSESS» project aims to determine whether episodic mass loss is a dominant process in the evolution of the most massive stars by conducting the first extensive, multi-wavelength survey of evolved massive stars in the nearby Universe. The project hinges on the fact that mass-losing stars form dust and are bright in the mid-infrared. We investigate the properties of evolved targets, selected by mid-infrared archival photometry from Spitzer, in ~25 nearby galaxies and estimate the amount of ejected mass, which will constrain evolutionary models. In this work we present our first results from the northern hemisphere survey, which includes the galaxies NGC 6822 and IC 10.

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., Comment: A white paper submitted for the Voyage 2050 long-term plan in the ESA Science Programme. 21 pages, 1 table, 6 figures. arXiv admin note: substantial text overlap with arXiv:1903.05235

(abridged) Mass loss is a key parameter in the evolution of massive stars, with discrepancies between theory and observations and with unknown importance of the episodic mass loss. To address this we need increased numbers of classified sources stars spanning a range of metallicity environments. We aim to remedy the situation by applying machine learning techniques to recently available extensive photometric catalogs. We used IR/Spitzer and optical/Pan-STARRS, with Gaia astrometric information, to compile a large catalog of known massive stars in M31 and M33, which were grouped in Blue, Red, Yellow, B[e] supergiants, Luminous Blue Variables, Wolf-Rayet, and background galaxies. Due to the high imbalance, we implemented synthetic data generation to populate the underrepresented classes and improve separation by undersampling the majority class. We built an ensemble classifier using color indices. The probabilities from Support Vector Classification, Random Forests, and Multi-layer Perceptron were combined for the final classification. The overall weighted balanced accuracy is ~83%, recovering Red supergiants at ~94%, Blue/Yellow/B[e] supergiants and background galaxies at ~50-80%, Wolf-Rayets at ~45%, and Luminous Blue Variables at ~30%, mainly due to their small sample sizes. The mixing of spectral types (no strict boundaries in their color indices) complicates the classification. Independent application to IC 1613, WLM, and Sextans A galaxies resulted in an overall lower accuracy of ~70%, attributed to metallicity and extinction effects. The missing data imputation was explored using simple replacement with mean values and an iterative imputor, which proved more capable. We also found that r-i and y-[3.6] were the most important features. Our method, although limited by the sampling of the feature space, is efficient in classifying sources with missing data and at lower metallicitites., 27 pages, 18 figures, 6 tables; accepted for publication in A&A; typos and minor editing corrections

Observations of massive stars in young open clusters (< ~8 Myr) have shown that a majority of them are in binary systems, most of which will interact during their life. Populations of massive stars older than ~20 Myr allow us to probe the outcome of such interactions after many systems have experienced mass and angular momentum transfer. Using multi-epoch integral-field spectroscopy, we investigate the multiplicity properties of the massive-star population in NGC 330 (~40 Myr) in the Small Magellanic Cloud to search for imprints of stellar evolution on the multiplicity properties. From six epochs of VLT/MUSE observations supported by adaptive optics we extract spectra and measure radial velocities for stars brighter than F814W = 19. We identify single-lined spectroscopic binaries through significant RV variability as well as double-lined spectroscopic binaries, and quantify the observational biases for binary detection. The observed spectroscopic binary fraction is 13.2+/-2.0 %. Considering period and mass ratio ranges from log(P)=0.15-3.5, and q = 0.1-1.0, and a representative set of orbital parameter distributions, we find a bias-corrected close binary fraction of 34 +8 -7 %. This seems to decline for the fainter stars, which indicates either that the close binary fraction drops in the B-type domain, or that the period distribution becomes more heavily weighted towards longer orbital periods. Both fractions vary strongly in different regions of the color-magnitude diagram which probably reveals the imprint of the binary history of different groups of stars. We provide the first homogeneous RV study of a large sample of B-type stars at a low metallicity. The overall bias-corrected close binary fraction of B stars in NGC 330 is lower than the one reported for younger Galactic and LMC clusters. More data are needed to establish whether this result from an age or a metallicty effect., 19 page (incl. appendix), 10 figures, 4 tables, accepted for publication in A&A

We present a case study of using a novel method to identify red supergiant (RSG) candidates in NGC 6822, based on their 1.6 $\mu$m "H-bump". We collected 32 bands of photometric data for NGC 6822 ranging from optical to MIR. By using the theoretical spectra from MARCS, we demonstrate that there is a prominent difference around 1.6 $\mu$m ("H-bump") between low-surface-gravity (LSG) and high-surface-gravity (HSG) targets. Taking advantage of this feature, we identify efficient color-color diagrams (CCDs) of rzH and rzK to separate HSG and LSG targets from crossmatching of optical and NIR data. Moreover, synthetic photometry from ATLAS9 also give similar results. Further separating RSG candidates from the rest of the LSG candidates is done by using semi-empirical criteria on NIR CMDs and resulted in 323 RSG candidates. Meanwhile, the simulation of foreground stars from Besan\c{c}on models also indicates that our selection criteria is largely free from the contamination of Galactic giants. In addition to the "H-bump" method, we also use the traditional BVR method as a comparison and/or supplement, by applying a slightly aggressive cut to select as much as possible RSG candidates (358 targets). Furthermore, the Gaia astrometric solution is used to constrain the sample, where 181 and 193 targets were selected from the "H-bump" and BVR method, respectively. The percentages of selected targets in both methods are similar as $\sim$60\%, indicating the comparable accuracy of the two methods. In total, there are 234 RSG candidates after combining targets from both methods with 140 ($\sim$60\%) of them in common. The final RSG candidates are in the expected locations on the MIR CMDs, while the spatial distribution is also coincident with the FUV-selected star formation regions, suggesting the selection is reasonable and reliable., Comment: A&A accepted

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., Comment: Accepted to A&A The core of the paper 9 pages, 5 pages of Appendix

The extreme ultraviolet (EUV) spectra of distant star-forming 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. This makes EUV spectra poorly constrained. 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. Here 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 using mm recombination lines of HI, HeI and HeII. We find that in all cases the EUV spectra between 13.6 and 54.4 eV have similar frequency dependence: L_{\nu}~ \nu^{-4.5 +/- 0.4}. We compare the inferred values of the EUV spectral slopes with the values expected for a purely single stellar evolution model (Starburst99) and the Binary Population and Spectral Synthesis code (BPASS). We find that the observed spectral slope differs from the model predictions. This may imply that the fraction of interacting binaries in HII regions is substantially lower than assumed in BPASS. The technique demonstrated here allows one to deduce the EUV spectra of star forming regions providing critical insight into photon production rates at \lambda < 912 A and can serve as calibration to starburst synthesis models, improving our understanding of star formation in distant universe and the properties of ionizing flux during reionization., Comment: v1: 18 pages, 5 figures, 1 table, 1 appendix v2: BPASS version 2.1 -> 2.2.1, SB99 tracks explicitly added to figures, evolution with metallicity added, few points clarified. 20 pages, 6 figures, 1 table. Accepted to ApJ

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}$., 36 pages, 19 figures, 10 tables, published in A\&A

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., 19 pages, 17 figures, A&A accepted, catalogs are included in the ancillary files

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., 24 pages, 4 tables, 12 figures, accepted for publication in ApJ

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., Comment: 19 pages (incl. appendix), 21 figures, 3 tables, accepted for publication in A&A

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., Comment: 10 pages, 6 figures, accepted for publication in ApJ Letters

We present a clean, magnitude-limited (IRAC1 or WISE1$\leq$15.0 mag) multiwavelength source catalog for the LMC. The catalog was built upon crossmatching ($1''$) and deblending ($3''$) between the SEIP source list and Gaia DR2, with strict constraints on the Gaia astrometric solution to remove the foreground contamination. The catalog contains 197,004 targets in 52 different bands including 2 ultraviolet, 21 optical, and 29 infrared bands. Additional information about radial velocities and spectral/photometric classifications were collected from the literature. The bright end of our sample is mostly comprised of blue helium-burning stars (BHeBs) and red HeBs with inevitable contamination of main sequence stars at the blue end. After applying modified magnitude and color cuts based on previous studies, we identify and rank 2,974 RSG, 508 YSG, and 4,786 BSG candidates in the LMC in six CMDs. The comparison between the CMDs of the LMC and SMC indicates that the most distinct difference appears at the bright red end of the optical and near-infrared CMDs, where the cool evolved stars (e.g., RSGs, AGB, and RGs) are located, which is likely due to the effect of metallicity and SFH. Further quantitative comparison of colors of massive star candidates in equal absolute magnitude bins suggests that, there is basically no difference for the BSG candidates, but large discrepancy for the RSG candidates as LMC targets are redder than the SMC ones, which may be due to the combined effect of metallicity on both spectral type and mass-loss rate, and also the age effect. The $T_{\rm eff}$ of massive star populations are also derived from reddening-free color of $(J-K_{\rm S})_0$. The $T_{\rm eff}$ ranges are $35008000$ K for BSG population, with larger uncertainties towards the hotter stars., Comment: 18 pages, 14 figures, A&A accepted

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., Comment: 18 pages, 12 figures, A&A accepted

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., Comment: accepted to A&A Core of the paper 15 pages, Appendix: 41 pages

Metal-poor massive stars are assumed to be progenitors of certain supernovae, gamma-ray bursts, and compact object mergers that might contribute to the early epochs of the Universe with their strong ionizing radiation. However, this assumption remains mainly theoretical because individual spectroscopic observations of such objects have rarely been carried out below the metallicity of the Small Magellanic Cloud. Aims. Here we explore the predictions of the state-of-the-art theories of stellar evolution combined with those of stellar atmospheres about a certain type of metal-poor (0.02 Z·) hot massive stars, the chemically homogeneously evolving stars that we call Transparent Wind Ultraviolet INtense (TWUIN) stars. Methods. We computed synthetic spectra corresponding to a broad range in masses (20-130 M·) and covering several evolutionary phases from the zero-age main-sequence up to the core helium-burning stage. We investigated the influence of mass loss and wind clumping on spectral appearance and classified the spectra according to the Morgan-Keenan (MK) system. 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) are detected. 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 result of chemically homogeneous evolution and might therefore be 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 because they are hot without showing prominent emission lines during most of their lifetime. © 2019 ESO., This research was supported by grant 16-01116S (GACR). The Astronomical Institute Ondrejov is supported by the project RVO: 67985815. D. Sz. is grateful for the relevant discussion with G. Grafener and N. Langer, as well as to A. Szabo, as usual. A.A.C.S. is supported by the Deutsche Forschungsgemeinschaft (DFG) under grant HA 1455/26 and would like to thank STFC for funding under grant number ST/R000565/1. TS acknowledges funding from the European Research Council (ERC) under the European Union's DLV-772225-MULTIPLES Horizon 2020 research and innovation program.

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., Comment: Authors' version of a response to a technical comment published in Science; 8 pages, 1 figure

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], 20 pages (incl. appendix), 15 figures, 2 tables; accepted for publication in A&A

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., Authors' version of paper published in Science at http://science.sciencemag.org/content/359/6371/69 ; 15 pages main text (3 figures), 47 pages supplementary materials (10 figures, 3 tables)

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, 10 pages, 6 figures, accepted for publication in A&A

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., Comment: 7 pages, 4 figures; accepted for publication in MNRAS

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., Accepted for publication in A&A. Appendixes A and B have been truncated due to size limitations, the full version will be available on A&A

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\, Comment: 39 pages, 12 figures in the main text, 12 figures in the appendix. Accepted for publication in A&A. Tables C1-C5 will be available at the CDS. Appendix E has been truncated due to size limitations, the full version will be available on A&A

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., 5 pages; Accepted for publication in Astronomy and Astrophysics Letters

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., Comment: Accepted for publication in MNRAS Letters, 5 pages, 3 figures

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., Accepted for Publication in A&A, 16 pages, 17 figures and 4 tables

© ESO, 2017. Context. Theoretically, rotation-induced chemical mixing in massive stars has far reaching evolutionary consequences, affecting the sequence of morphological phases, lifetimes, nucleosynthesis, and supernova characteristics. Aims. 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 by analysing their surface nitrogen abundances. Methods. Using stellar and wind properties derived in a previous VFTS study we computed synthetic spectra for a set of up to 21 N ii-v lines in the optical spectral range, using the non-LTE atmosphere code FASTWIND. 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 ν e sini> 200 km s -1 . Results. We analysed the nitrogen abundance as a function of projected rotation rate ν e sini and confronted it with predictions of rotational mixing. 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. For the rapidly-spinning O-Type stars we can only provide upper limits on the nitrogen abundance, which are not in violation with theoretical expectations. Hence, the data cannot be used to test the physics of rotation induced mixing in the regime of high spin rates. Conclusions. While the surface abundances of 60-70 percent of presumed single O-Type giants to supergiants behave in conformity with expectations, at least 30-40 percent of our sample can not be understood in the current framework of rotational mixing for single stars. Even though we have excluded stars showing radial velocity variations, of our sample may have remained contaminated by post-interaction binary products. Hence, it is plausible that effects of binary interaction need to be considered to understand their surface properties. Alternatively, or in conjunction, the effects of magnetic fields or alternative mass-loss recipes may need to be invoked. ispartof: Astronomy & Astrophysics vol:600 status: published

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., 16 pages, 16 figures, paper accepted in Astronomy & Astrophysics

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., Accepted for publication in Astronomy & Astrophysics. 18 pages, 24 figures

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., Comment: 23 pages, 6 figures, accepted for publication in ApJ

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., Accepted for publication in Astronomy and Astrophysics. 10 pages, 8 figures

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., Comment: 19 pages, 12 images. Accepted for publication in Astronomy & Astrophysics

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., 6 pages, 3 figures, Proceedings IAU Symposium No. 307, 2014, 'New windows on massive stars: asteroseismology, interferometry, and spectropolarimetry'

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., 13 pages, 12 figures, accepted for publication in Astronomy & Astrophysics

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., 29 pages, 7 figures, Proceedings Joint Discussion 2. To be published in Hightlights of Astronomy, ed. T. Montmerle

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., Comment: 18 pages, 10 figures, accepted for publication in Astronomy & Astrophysics

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 $��= 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., Accepted for publication in ApJ

We present a new empirical prescription for the mass-loss rates of hydrogen-free Wolf-Rayet stars based on results of detailed spectral analyses of WC and WO stars. Compared to the prescription of Nugis & Lamers (2000), Ṁ is less sensitive to the surface helium abundance, implying a stronger mass loss at the late stages of Wolf-Rayet evolution. The winds of hydrogen-free WN stars have a strong metallicity dependence, while those of WC and WO stars have a very weak metallicity dependence.

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., 22 pages, 16 figures. Uses emulateapj style. Final three tables inserted as a figure to avoid issues with astro-ph's version of latex mangling the use of lscape. To be published in the Astrophysical Journal, January, 2011

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., Comment: 9 pages, 3 figures; accepted for publication in The Astrophysical Journal Letters

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., Comment: 2 pages, 1 figure, Proceedings IAU Symposium No. 307, 2014, 'New windows on massive stars: asteroseismology, interferometry, and spectropolarimetry'

The enigmatic oxygen-sequence Wolf-Rayet stars represent a rare stage in the evolution of massive stars. Their properties can provide unique constraints on the pre-supernova evolution of massive stars. This work presents the results of a quantitative spectroscopic analysis of the known single WO stars, with the aim to obtain the key stellar parameters and deduce their evolutionary state.X-Shooter spectra of the WO stars are modeled using the line-blanketed non-local thermal equilibrium atmosphere code cmfgen. The obtained stellar parameters show that the WO stars are very hot, with temperatures ranging from 150 kK to 210 kK. Their chemical composition is dominated by carbon (>50%), while the helium mass fraction is very low (down to 14%). Oxygen mass fractions reach as high as 25%. These properties can be reproduced with dedicated evolutionary models for helium stars, which show that the stars are post core-helium burning and very close to their eventual supernova explosion. The helium-star masses indicate initial masses or approximately 40 - 60M⊙.Thus, WO stars represent the final evolutionary stage of stars with estimated initial masses of 40 - 60M⊙. They are post core-helium burning and may explode as type Ic supernovae within a few thousand years.

Combined with the collecting power of the Very Large Telescope, X-shooter is the most sensitive medium-resolution spectrograph currently in operation, allowing us to perform quantitative spectroscopy of objects that, up to now, were deemed too faint. In addition, with its unique wavelength coverage, X-shooter provides access to an un-precedentedly large number of diagnostic lines. We review our recent work on massive stars in Local Group dwarf galaxies and in NGC 55, an irregular galaxy in the foreground of the Sculptor Group. The observations were obtained as part of the ESO Science Verification Programme and the NOVA-Dutch Guaranteed Time Observation Programme. The aim is to investigate the evolutionary status of high-mass stars in various environments and to test the theory of radiation-driven stellar winds at metallicities below that of the Small Magellanic Cloud.