Search

We provide line luminosities and spectroscopic templates of prominent optical emission lines of 133 Galactic Wolf-Rayet stars by exploiting Gaia DR3 parallaxes and optical spectrophotometry, and provide comparisons with 112 counterparts in the Magellanic Clouds. Average line luminosities of the broad blue (He II 4686, C III 4647,51, N III 4634,41, N V 4603,20) and yellow (C IV 5801,12) emission features for WN, WN/C, WC and WO stars have application in characterising the Wolf-Rayet populations of star-forming regions of distant, unresolved galaxies. Early-type WN stars reveal lower line luminosities in more metal poor environments, but the situation is less clear for late-type WN stars. LMC WC4-5 line luminosities are higher than their Milky Way counterparts, with line luminosities of Magellanic Cloud WO stars higher than Galactic stars. We highlight other prominent optical emission lines, N IV 3478,85 for WN and WN/C stars, O IV 3403,13 for WC and WO stars and O VI 3811,34 for WO stars. We apply our calibrations to representative metal-poor and metal-rich WR galaxies, IC 4870 and NGC 3049, respectively, with spectral templates also applied based on a realistic mix of subtypes. Finally, the global blue and C IV 5801,12 line luminosities of the Large (Small) Magellanic Clouds are 2.6e38 erg/s (9e36 erg/s) and 8.8e37 erg/s (4e36 erg/s), respectively, with the cumulative WR line luminosity of the Milky Way estimated to be an order of magnitude higher than the LMC., Comment: 15 pages, 10 figures, MNRAS in press. Appendix A (individual WR line luminosities), Appendix B (WR templates) available from https://zenodo.org/record/7573775

We report the discovery of a luminous late-time source at the position of the fast blue optical transient (FBOT) AT2018cow on images taken by the Hubble Space Telescope (HST) at 714 d and 1136 d after its explosion. This source is detected at both UV and optical wavelengths and has prominent H$\alpha$ emission. It has a very stable brightness between the two epochs and a very blue spectral energy distribution (SED) consistent with $f_\lambda$ $\propto$ $\lambda^{-4.1 \pm 0.1}$, i.e. the Rayleigh-Jeans tail of a hot blackbody with a very high temperature of log($T$/K) $>$ 4.6 and luminosity of log($L$/$L_\odot$) $>$ 7.0. This late-time source is unlikely to be an unrelated object in chance alignment, or due to a light echo of AT2018cow. Other possible scenarios also have some difficulties in explaining this late-time source, including companion star(s), star cluster, the survived progenitor star, interaction with circumstellar medium (CSM), magnetar, or delayed accretion in a tidal disruption event (TDE). Long-term and multi-wavelength monitoring will help to resolve its nature and finally reveal the origin of the "Cow"., Comment: 5 pages, 4 figures, accepted for publication in MNRAS Letters

The cluster R136 in the giant star-forming region 30 Doradus in the Large Magellanic Cloud (LMC) offers a unique opportunity to resolve a stellar population in a starburst-like environment. We obtain the near-infrared to ultraviolet extinction towards 50 stars in the core of R136, employing the `extinction without standards' method. To assure good fits over the full wavelength range, we combine and modify existing extinction laws. We detect a strong spatial gradient in the extinction properties across the core of R136, coinciding with a gradient in density of cold gas that is part of a molecular cloud lying northeast of the cluster. In line with previous measurements of R136 and the 30 Doradus region, we obtain a high total-to-relative extinction ($R_V = 4.38 \pm 0.87$). However, the high values of $R_V$ are accompanied by relatively strong extinction in the ultraviolet, contrary to what is observed for Galactic sightlines. The relatively strong ultraviolet extinction suggests that the properties of the dust towards R136 differ from those in the Milky Way. For $R_{V} \sim 4.4$, about three times fewer ultraviolet photons can escape from the ambient dust environment relative to the canonical Galactic value of $R_{V} \sim 3.1$ at the same $A_{V}$. Therefore, if dust in the R136 star-bursting environment is characteristic for cosmologically distant star-bursting regions, the escape fraction of ultraviolet photons from such regions is overestimated by a factor of three relative to the standard Milky Way assumption for the total-to-selective extinction. Furthermore, a comparison with average curves tailored to other regions of the LMC shows that large differences in ultraviolet extinction exist within this galaxy. Further investigation is required in order to decipher whether or not there is a relation between $R_V$ and ultraviolet extinction in the LMC., Accepted for publication in A&A (18 pages, 12 figures)

SN 2019yvr is the second Type Ib supernova (SN) with a possible direct detection of its progenitor (system); however, the spectral energy distribution (SED) of the pre-explosion source appears much cooler and overluminous than an expected helium-star progenitor. Using Hubble Space Telescope (HST) images and MUSE integral-field-unit (IFU) spectroscopy, we find the SN environment contains three episodes of star formation; the low ejecta mass suggests the SN progenitor is most likely from the oldest population, corresponding to an initial mass of 10.4$^{+1.5}_{-1.3}$ $M_\odot$. The pre-explosion SED can be reproduced by two components, one for the hot and compact SN progenitor and one for a cool and inflated yellow hypergiant (YHG) companion that dominates the brightness. Thus, SN 2019yvr could possibly be the first Type Ib/c SN for which the progenitor's binary companion is directly detected on pre-explosion images. Both the low progenitor mass and the YHG companion suggest significant binary interaction during their evolution. Similar to SN 2014C, SN 2019yvr exhibits a metamorphosis from Type Ib to Type IIn, showing signatures of interaction with hydrogen-rich circumstellar material (CSM) at >150 days; our result supports enhanced pre-SN mass loss as an important process for hydrogen-poor stars at the low mass end of core-collapse SN progenitors., 15 pages, 9 figures, accepted for publication in MNRAS

This paper reports an environmental analysis of 41 uniformly-selected stripped-envelope supernovae (SESNe) based on deep ultraviolet-optical images acquired by the Hubble Space Telescope. Young stellar populations are detected in most SN environments and their ages are derived with a hierarchical Bayesian approach. The age distributions are indistinguishable between Type IIb and Type Ib while that for Type Ic is systematically younger. This suggests that the Type Ic SN progenitors are more massive while the Type IIb and Type Ib SNe have very similar progenitor masses. Our result supports a hybrid envelope-stripping mechanism, in which the hydrogen envelopes of the SESN progenitors are stripped via a mass-insensitive process (e.g. binary interaction) while the helium envelopes are stripped via a mass-sensitive process (e.g. stellar wind of the post-binary interaction progenitor). We also provide progenitor constraints for three Type Ibn SNe and two broad-lined Type Ic SNe. All these results demonstrate the importance of the very diverse mass-loss processes in the origins of SESNe., 15 pages, 9 figures, accepted for publication in MNRAS

We aim to hunt for massive binaries hosting a black hole companion (OB+BH) and establish the natal mass-ratio distribution of massive stars at the subsolar metallicity environment of the Large Magellanic Cloud (LMC). We use the shift-and-add grid disentangling technique to characterize the hidden companions in 51 SB1 O-type and evolved B-type binaries in the LMC monitored in the framework of the Tarantula Massive Binary Monitoring (TMBM). Out of the 51 SB1 systems, 43 (84%) are found to have non-degenerate stellar companions, of which 28 are confident detections, and 15 are less certain (SB1: or SB2:). Of these 43 targets, one is found to be a triple (VFTS 64), and two are found to be quadruples (VFTS 120, 702). The remaining eight targets (16%) retain an SB1 classification. Aside from the unambiguous case of VFTS 243, analysed in detailed in a separate paper, we identify two additional OB+BH candidates: VFTS 514 and VFTS 779. Additional black holes may be present in the sample but at a lower probability. Our study firmly establishes a virtually flat natal mass-ratio distribution for O-type stars at LMC metallicity, covering the entire mass-ratio range (0.05 < q < 1) and periods in the range 0 < log P < 3 [d]. The nature of the OB+BH candidates should be verified through future monitoring, but the frequency of OB+BH candidates is generally in line with recent predictions at LMC metallicity., 41 pages (14 main article + 27 appendix), accepted to A&A

The presence of massive stars (MSs) in the region close to the Galactic Center (GC) poses several questions about their origin. The harsh environment of the GC favors specific formation scenarios, each of which should imprint characteristic kinematic features on the MSs. We present a 2D kinematic analysis of MSs in a GC region surrounding Sgr A* based on high-precision proper motions obtained with the Hubble Space Telescope. Thanks to a careful data reduction, well-measured bright stars in our proper-motion catalogs have errors better than 0.5 mas yr$^{-1}$. We discuss the absolute motion of the MSs in the field and their motion relative to Sgr A*, the Arches and the Quintuplet. For the majority of the MSs, we rule out any distance further than 3-4 kpc from Sgr A* using only kinematic arguments. If their membership to the GC is confirmed, most of the isolated MSs are likely not associated with either the Arches or Quintuplet clusters or Sgr A*. Only a few MSs have proper motions suggesting they are likely members of the Arches cluster, in agreement with previous spectroscopic results. Line-of-sight radial velocities and distances are required to shed further light on the origin of most of these massive objects. We also present an analysis of other fast-moving objects in the GC region, finding no clear excess of high-velocity escaping stars. We make our astro-photometric catalogs publicly available., Comment: 28 pages, 33 figures, 9 tables. Accepted for publication in MNRAS. [v3: Fixed bibliography]

We obtain distances to 383 Galactic Wolf-Rayet (WR) stars from $\textit{Gaia}$ DR2 parallaxes and Bayesian methods, with a prior based on H$\mathrm{\scriptsize{II}}$ regions and dust extinction. Distances agree with those from Bailer-Jones et al. for stars up to 2 kpc from the Sun, though deviate thereafter due to differing priors, leading to modest reductions in luminosities for recent WR spectroscopic results. We calculate visual and K-band absolute magnitudes, accounting for dust extinction contributions and binarity, and identify 187 stars with reliable absolute magnitudes. For WR and O stars within 2 kpc, we find a WR/O ratio of 0.09. The distances are used to generate absolute magnitude calibrations and obtain the $\mathrm{\textit{Gaia}}$ colour magnitude diagram for WR stars. Average $v^{\rm WR}$-band absolute magnitudes for WN stars range from -3.6 mag (WN3-4) to -7.0 mag (WN8-9ha), and -3.1 (WO2-4) to -4.6 mag (WC9), with standard deviations of $\sim$0.6 mag. Using H$\mathrm{\scriptsize{II}}$ region scale heights, we identify 31 WR stars at large (3$\sigma$, |z|$\geq$156 pc) distances from the mid-plane as potential runaways accounting for the Galactic warp, of which only 4 involve WN8-9 stars, contrary to previous claims., Comment: Amended section 3.3, table 4, table 5 and table 6 based on corrected Av to AK conversion. Minor updates to figures, WR/O ratio and WR stars with reliable results. 19 pages, 15 figures. Additional appendices updated (57 pages, 5 figures). To be published in MNRAS

Context: The star cluster R136 inside the LMC hosts a rich population of massive stars, including the most massive stars known. The strong stellar winds of these very luminous stars impact their evolution and the surrounding environment. We currently lack detailed knowledge of the wind structure that is needed to quantify this impact. Aims: To observationally constrain the stellar and wind properties of the massive stars in R136, in particular the parameters related to wind clumping. Methods: We simultaneously analyse optical and UV spectroscopy of 53 O-type and 3 WNh-stars using the FASTWIND model atmosphere code and a genetic algorithm. The models account for optically thick clumps and effects related to porosity and velocity-porosity, as well as a non-void interclump medium. Results: We obtain stellar parameters, surface abundances, mass-loss rates, terminal velocities and clumping characteristics and compare these to theoretical predictions and evolutionary models. The clumping properties include the density of the interclump medium and the velocity-porosity of the wind. For the first time, these characteristics are systematically measured for a wide range of effective temperatures and luminosities. Conclusions: We confirm a cluster age of 1.0-2.5 Myr and derive an initial stellar mass of $\geq 250 {\rm M}_\odot$ for the most massive star in our sample, R136a1. The winds of our sample stars are highly clumped, with an average clumping factor of $f_{\rm cl}=29\pm15$. We find tentative trends in the wind-structure parameters as a function of mass-loss rate, suggesting that the winds of stars with higher mass-loss rates are less clumped. We compare several theoretical predictions to the observed mass-loss rates and terminal velocities and find that none satisfactorily reproduces both quantities. The prescription of Krti\v{c}ka & Kub\'at (2018) matches best the observed mass-loss rates., Comment: Accepted for publication in A&A; Appendix I will not be included in the published version

We present a grid of stellar models at super-solar metallicity (Z = 0.020) extending the previous grids of Geneva models at solar and sub-solar metallicities. A metallicity of Z = 0.020 was chosen to match that of the inner Galactic disk. A modest increase of 43% (=0.02/0.014) in metallicity compared to solar models means that the models evolve similarly to solar models but with slightly larger mass loss. Mass loss limits the final total masses of the super-solar models to 35 M$_\odot$ even for stars with initial masses much larger than 100 M$_\odot$. Mass loss is strong enough in stars above 20 M$_\odot$ for rotating stars (25 M$_\odot$ for non-rotating stars) to remove the entire hydrogen-rich envelope. Our models thus predict SNII below 20 M$_\odot$ for rotating stars (25 M$_\odot$ for non-rotating stars) and SNIb (possibly SNIc) above that. We computed both isochrones and synthetic clusters to compare our super-solar models to the Westerlund 1 (Wd1) massive young cluster. A synthetic cluster combining rotating and non-rotating models with an age spread between log10 (age/yr) = 6.7 and 7.0 is able to reproduce qualitatively the observed populations of WR, RSG and YSG stars in Wd1, in particular their simultaneous presence at log10(L/L$_\odot$) = 5-5.5. The quantitative agreement is imperfect and we discuss the likely causes: synthetic cluster parameters, binary interactions, mass loss and their related uncertainties. In particular, mass loss in the cool part of the HRD plays a key role., Comment: 15 pages, 12 figures. Accepted for publication in MNRAS

The LIGO/Virgo gravitational--wave observatories have detected 50 BH-BH coalescences. This sample is large enough to have allowed several recent studies to draw conclusions about the branching ratios between isolated binaries versus dense stellar clusters as the origin of double BHs. It has also led to the exciting suggestion that the population is highly likely to contain primordial black holes. Here we demonstrate that such conclusions cannot yet be robust, because of the large current uncertainties in several key aspects of binary stellar evolution. These include the development and survival of a common envelope, the mass and angular momentum loss during binary interactions, mixing in stellar interiors, pair-instability mass loss and supernova outbursts. Using standard tools such as the population synthesis codes StarTrack and COMPAS and the detailed stellar evolution code MESA, we examine as a case study the possible future evolution of Melnick 34, the most massive known binary star system. We show that, despite its well-known orbital architecture, various assumptions regarding stellar and binary physics predict a wide variety of outcomes: from a close BH-BH binary (which would lead to a potentially detectable coalescence), through a wide BH-BH binary (which might be seen in microlensing observations), or a Thorne-Zytkow object, to a complete disruption of both objects by pair-instability supernovae. Thus since the future of massive binaries is inherently uncertain, sound predictions about the properties of BH-BH systems are highly challenging at this time. Consequently, drawing conclusions about the formation channels for the LIGO/Virgo BH-BH merger population is premature., Comment: comments welcomed

Stellar-mass black holes are the final remnants of stars born with more than 15 solar masses. Billions are expected to reside in the Local Group, yet only a few are known, mostly detected through X-rays emitted as they accrete material from a companion star. Here, we report on VFTS 243: a massive X-ray-faint binary in the Large Magellanic Cloud. With an orbital period of 10.4 d, it comprises an O-type star of 25 solar masses and an unseen companion of at least nine solar masses. Our spectral analysis excludes a non-degenerate companion at a 5σ confidence level. The minimum companion mass implies that it is a black hole. No other X-ray-quiet black hole is unambiguously known outside our Galaxy. The (near-)circular orbit and kinematics of VFTS 243 imply that the collapse of the progenitor into a black hole was associated with little or no ejected material or black-hole kick. Identifying such unique binaries substantially impacts the predicted rates of gravitational-wave detections and properties of core-collapse supernovae across the cosmos., This research has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (T.S. and H.S., grant agreement 772225: MULTIPLES). T.S. acknowledges support from the European Union’s Horizon 2020 programme under Marie Skłodowska-Curie grant agreement 101024605. This work is based on observations collected at the ESO under programme IDs 182.D-0222, 090.D-0323 and 092.D-0136. L.M. thanks the ESA and the Belgian Federal Science Policy Office (BELSPO) for their support in the framework of the PRODEX programme. P. Marchant acknowledges support from FWO junior postdoctoral fellowship 12ZY520N. We thank J. Hillier for making the CMFGEN code available. C.H. acknowledges support from the KU Leuven Research Council (grant C16/17/007: MAESTRO). P.A.C. acknowledges support from UK Science and Technology Facilities Council research grant ST/V000853/1. V.H.-B. acknowledges the support of the Natural Sciences and Engineering Research Council of Canada (NSERC) through grant RGPIN-2020-05990. M. Gieles acknowledges support from the Ministry of Science and Innovation through a Europa Excelencia grant (EUR2020-112157). The PoWR code, as well as the associated post-processing and visualization tool WRplot, has been developed under the guidance of W.-R. Hamann with substantial contributions from L. Koesterke, G. Gräfener, A. Sander, T.S. and other co-workers and students. This work has received funding from the ERC under the European Union’s Horizon 2020 research and innovation programme (grant agreement 945806). It is also supported by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy EXC 2181/1-390900948 (the Heidelberg STRUCTURES Excellence Cluster). S.T. acknowledges support from the Netherlands Research Council NWO (grants VENI 639.041.645, VIDI 203.061). A.H. and D.J.L. acknowledge support by the Spanish MCI through grant PGC-2018-0913741-B-C22 and the Severo Ochoa Programme through CEX2019-000920-S. J.M.A. acknowledges support from the Spanish Government Ministerio de Ciencia, Innovación y Universidades through grant PGC2018-095-049-B-C22.

Stars strongly impact their environment, and shape structures on all scales throughout the universe, in a process known as “feedback.” Due to the complexity of both stellar evolution and the physics of larger astrophysical structures, there remain many unanswered questions about how feedback operates and what we can learn about stars by studying their imprint on the wider universe. In this white paper, we summarize discussions from the Lorentz Center meeting “Bringing Stellar Evolution and Feedback Together” in 2022 April and identify key areas where further dialog can bring about radical changes in how we view the relationship between stars and the universe they live in.

The hot WN star WR2 (HD6327) has been claimed to have many singular characteristics. To explain its unusually rounded and relatively weak emission line profiles, it has been proposed that WR2 is rotating close to break-up with a magnetically confined wind. Alternatively, the line profiles could be explained by the dilution of WR2's spectrum by that of a companion. In this paper, we present a study of WR2 using near-infrared AO imaging and optical spectroscopy and polarimetry. Our spectra reveal the presence of weak photospheric absorption lines from a ~B2.5-4V companion, which however contributes only ~5-10% to the total light, suggesting that the companion is a background object. Therefore, its flux cannot be causing any significant dilution of the WR star's emission lines. The absence of intrinsic linear continuum polarization from WR2 does not support the proposed fast rotation. Our Stokes V spectrum was not of sufficient quality to test the presence of a moderately strong organized magnetic field but our new modelling indicates that to confine the wind the putative magnetic field must be significantly stronger than was previously suggested sufficiently strong as to make its presence implausible., 12 pages, 9 figures

Massive star evolution at low metallicity is closely connected to many fields in high-redshift astrophysics, but poorly understood. The Small Magellanic Cloud (SMC) is a unique laboratory to study this because of its metallicity of 0.2 Zsol, its proximity, and because it is currently forming stars. We used a spectral type catalog in combination with GAIA magnitudes to calculate temperatures and luminosities of bright SMC stars. By comparing these with literature studies, we tested the validity of our method, and using GAIA data, we estimated the completeness of stars in the catalog as a function of luminosity. This allowed us to obtain a nearly complete view of the most luminous stars in the SMC. When then compared with stellar evolution predictions. We also calculated the extinction distribution, the ionizing photon production rate, and the star formation rate. Our results imply that the SMS hosts only 30 very luminous main-sequence stars (M > 40 Msol; L > 10^5 Lsol), which are far fewer than expected from the number of stars in the luminosity range 3*10^4 < L/Lsol < 3*10^5 and from the typically quoted star formation rate in the SMC. Even more striking, we find that for masses above M > 20 Msol, stars in the first half of their hydrogen-burning phase are almost absent. This mirrors a qualitatively similar peculiarity that is known for the Milky Way and Large Magellanic Cloud. This amounts to a lack of hydrogen-burning counterparts of helium-burning stars, which is more pronounced for higher luminosities. We argue that a declining star formation rate or a steep initial mass function are unlikely to be the sole explanations for the dearth of young bright stars. Instead, many of these stars might be embedded in their birth clouds, although observational evidence for this is weak. We discuss implications for cosmic reionization and the top end of the initial mass function., 26 pages, 23 figures. Astronomy and Astrophysics, in press

We compared high-contrast near-infrared images of the core of R136 taken by VLT/SPHERE, in two epochs separated by 3.06 years. For the first time we monitored the dynamics of the detected sources in the core of R136 from a ground-based telescope with adaptive optics. The aim of these observations was to search for High prOper Motion cAndidates (HOMAs) in the central region of R136 (r, Accepted to be published in A&A/Letter. catalogue of the reliable-consistent sources are available online via CDS

We present the sharpest and deepest near infrared photometric analysis of the core of R136, a newly formed massive star cluster at the centre of the 30 Doradus star forming region in the Large Magellanic Cloud. We used the extreme adaptive optics of the SPHERE focal instrument implemented on the ESO Very Large Telescope and operated in its IRDIS imaging mode, for the second time with longer exposure time in the H- and K filters. Our aim was to (i) increase the number of resolved sources in the core of R136, and (ii) to compare with the first epoch to classify the properties of the detected common sources between the two epochs. Within the field of view (FOV) of 10.8"x12.1" (2.7pc x3.0pc), we detected 1499 sources in both H and K filters, for which 76% of these sources have visual companions closer than 0.2". The larger number of detected sources, enabled us to better sample the mass function (MF). The MF slopes are estimated at ages of 1, 1.5 and 2 Myr, at different radii, and for different mass ranges. The MF slopes for the mass range of 10-300 solar-mass are about 0.3 dex steeper than the mass range of 3-300 solar-mass, for the whole FOV and different radii. Comparing the JHK colours of 790 sources common in between the two epochs, 67% of detected sources in the outer region (r >3") are not consistent with evolutionary models at 1-2 Myr and with extinctions similar to the average cluster value, suggesting an origin from ongoing star formation within 30 Doradus, unrelated to R136., 24 pages, 20 figures, 6 tables. Accepted for publication in MNRAS

The recent discovery of a spectacular dust plume in the system 2XMM J160050.7-514245 (referred to as "Apep") suggested a physical origin in a colliding-wind binary by way of the "Pinwheel" mechanism. Observational data pointed to a hierarchical triple-star system, however several extreme and unexpected physical properties seem to defy the established physics of such objects. Most notably, a stark discrepancy was found in the observed outflow speed of the gas as measured spectroscopically in the line-of-sight direction compared to the proper motion expansion of the dust in the sky plane. This enigmatic behaviour arises at the wind base within the central Wolf-Rayet binary: a system that has so far remained spatially unresolved. Here we present an updated proper motion study deriving the expansion speed of Apep's dust plume over a two-year baseline that is four times slower than the spectroscopic wind speed, confirming and strengthening the previous finding. We also present the results from high-angular-resolution near-infrared imaging studies of the heart of the system, revealing a close binary with properties matching a Wolf-Rayet colliding-wind system. Based on these new observational constraints, an improved geometric model is presented yielding a close match to the data, constraining the orbital parameters of the Wolf-Rayet binary and lending further support to the anisotropic wind model., Comment: This article has been accepted for publication in the Monthly Notices of the Royal Astronomical Society. 17 pages, 9 figures, 6 tables

We present an optical analysis of 55 members of R136, the central cluster in the Tarantula Nebula of the Large Magellanic Cloud. Our sample was observed with STIS aboard the Hubble Space Telescope, is complete down to about 40\,$M_{\odot}$, and includes 7 very massive stars with masses over 100\,$M_{\odot}$. We performed a spectroscopic analysis to derive their physical properties. Using evolutionary models we find that the initial mass function (IMF) of massive stars in R136 is suggestive of being top-heavy with a power-law exponent $\gamma \approx 2 \pm 0.3$, but steeper exponents cannot be excluded. The age of R136 lies between 1 and 2\,Myr with a median age of around 1.6\,Myr. Stars more luminous than $\log L/L_{\odot} = 6.3$ are helium enriched and their evolution is dominated by mass loss, but rotational mixing or some other form of mixing could be still required to explain the helium composition at the surface. Stars more massive than 40\,$M_{\odot}$ have larger spectroscopic than evolutionary masses. The slope of the wind-luminosity relation assuming unclumped stellar winds is $2.41\pm0.13$ which is steeper than usually obtained ($\sim 1.8$). The ionising ($\log Q_0\,[{\rm ph/s}] = 51.4$) and mechanical ($\log L_{\rm SW}\,[{\rm erg/s}] = 39.1$) output of R136 is dominated by the most massive stars ($>100\,M_{\odot}$). R136 contributes around a quarter of the ionising flux and around a fifth of the mechanical feedback to the overall budget of the Tarantula Nebula. For a census of massive stars of the Tarantula Nebula region we combined our results with the VLT-FLAMES Tarantula Survey plus other spectroscopic studies. We observe a lack of evolved Wolf-Rayet stars and luminous blue and red supergiants., Comment: Accepted for publication in MNRAS, 21 pages, 14 figures, 2 tables plus 33 pages of supplementary material

Infrared imaging of the colliding-wind binary Apep has revealed a spectacular dust plume with complicated internal dynamics that challenges standard colliding-wind binary physics. Such challenges can be potentially resolved if a rapidly-rotating Wolf-Rayet star is located at the heart of the system, implicating Apep as a Galactic progenitor system to long-duration gamma-ray bursts. One of the difficulties in interpreting the dynamics of Apep is that the spectral composition of the stars in the system was unclear. Here we present visual to near-infrared spectra that demonstrate that the central component of Apep is composed of two classical Wolf-Rayet stars of carbon- (WC8) and nitrogen-sequence (WN4-6b) subtypes. We argue that such an assignment represents the strongest case of a classical WR+WR binary system in the Milky Way. The terminal line-of-sight wind velocities of the WC8 and WN4-6b stars are measured to be $2100 \pm 200$ and $3500 \pm 100$ km s$^{-1}$, respectively. If the mass-loss rate of the two stars are typical for their spectral class, the momentum ratio of the colliding winds is expected to be $\approx$ 0.4. Since the expansion velocity of the dust plume is significantly smaller than either of the measured terminal velocities, we explore the suggestion that one of the Wolf-Rayet winds is anisotropic. We can recover a shock-compressed wind velocity consistent with the observed dust expansion velocity if the WC8 star produces a significantly slow equatorial wind with a velocity of $\approx$530 km s$^{-1}$. Such slow wind speeds can be driven by near-critical rotation of a Wolf-Rayet star., Accepted for publication in MNRAS, 10 pages, 9 figures, 1 table

Galactic Wolf-Rayet (WR) star membership of star forming regions can be used to constrain the formation environments of massive stars. Here, we utilise $\textit{Gaia}$ DR2 parallaxes and proper motions to reconsider WR star membership of clusters and associations in the Galactic disk, supplemented by recent near-IR studies of young massive clusters. We find that only 18$-$36% of 553 WR stars external to the Galactic Centre region are located in clusters, OB associations or obscured star-forming regions, such that at least 64% of the known disk WR population are isolated, in contrast with only 13% of O stars from the Galactic O star Catalogue. The fraction located in clusters, OB associations or star-forming regions rises to 25$-$41% from a global census of 663 WR stars including the Galactic Centre region. We use simulations to explore the formation processes of isolated WR stars. Neither runaways, nor low mass clusters, are numerous enough to account for the low cluster membership fraction. Rapid cluster dissolution is excluded as mass segregation ensures WR stars remain in dense, well populated environments. Only low density environments consistently produce WR stars that appeared to be isolated during the WR phase. We therefore conclude that a significant fraction of WR progenitors originate in low density association-like surroundings which expand over time. We provide distance estimates to clusters and associations host to WR stars, and estimate cluster ages from isochrone fitting., Comment: 19 pages, 7 figures. Accepted for publication in MNRAS

Core-collapse supernovae (SNe) are the inevitable fate of most massive stars. Since most stars form in groups, SN progenitors can be constrained with information of their environments. It remains challenging to accurately analyse the various components in the environment and to correctly identify their relationships with the SN progenitors. Using a combined dataset of VLT/MUSE spatially-resolved integral-field-unit (IFU) spectroscopy and HST/ACS+WFC3 high-spatial resolution imaging, we present a detailed investigation of the environment of the Type II-P SN 2004dg and Type IIb SN 2012P. The two SNe occurred in a spiral arm of NGC 5806, where a star-forming complex is apparent with a giant H II region. By modelling the ionised gas, a compact star cluster and the resolved stars, we derive the ages and extinctions of stellar populations in the vicinity of the SNe. The various components are consistent with a sequence of triggered star formation as the spiral density wave swept through their positions. For SNe 2004dg and 2012P, we identify their host stellar populations and derive initial masses of $10.0^{+0.3}_{-0.2}~M_\odot$ and $15.2^{+2.0}_{-1.0}~M_\odot$ for their progenitors, respectively. Both results are consistent with those from pre-explosion images or nebular-phase spectroscopy. SN 2012P is spatially coincident but less likely to be coeval with the star-forming complex. As in this case, star formation bursts on small scales may appear correlated if they are controlled by any physical processes on larger scales; this may lead to a high probability of chance alignment between older SN progenitors and younger stellar populations., Comment: 21 pages, 17 figures, accepted for publication in MNRAS

NGC 300 X-1 and IC 10 X-1 are currently the only two robust extragalactic candidates for being Wolf-Rayet/black hole X-ray binaries, the Galactic analogue being Cyg X-3. These systems are believed to be a late product of high-mass X-ray binary evolution and direct progenitors of black hole mergers. From the analysis of Swift data, the orbital period of NGC 300 X-1 was found to be 32.8 h. We here merge the full set of existing data of NGC 300 X-1, using XMM-Newton, Chandra and Swift observations to derive a more precise value of the orbital period of 32.7932 ± 0.0029 h above a confidence level of 99.99%. This allows us to phase connect the X-ray light curve of the source with radial velocity measurements of He II lines performed in 2010. We show that, as for IC 10 X-1 and Cyg X-3, the X-ray eclipse corresponds to maximum of the blueshift of the He II lines, instead of the expected zero velocity. This indicates that for NGC 300 X-1 as well, the wind of the WR star is completely ionised by the black hole radiation and that the emission lines come from the region of the WR star that is in the shadow. We also present for the first time the light curve of two recent very long XMM-Newton observations of the source, performed on the 16th to 20th of December 2016.

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

Type Ibn supernovae (SNe Ibn) are intriguing stellar explosions whose spectra exhibit narrow helium lines with little hydrogen. They trace the presence of circumstellar material (CSM) formed via pre-SN eruptions of their stripped-envelope progenitors. Early work has generally assumed that SNe Ibn come from massive Wolf-Rayet (WR) stars via single star evolution. In this paper, we report ultraviolet (UV) and optical observations of two nearby Type Ibn SNe 2006jc and 2015G conducted with the Hubble Space Telescope (HST) at late times. A point source is detected at the position of SN 2006jc, and we confirm the conclusion of Maund et al. that it is the progenitor's binary companion. Its position on the Hertzsprung-Russell (HR) diagram corresponds to a star that has evolved off the main sequence (MS); further analysis implies a low initial mass for the companion star ($M_2$ $\le$ 11.9$^{+1.2}_{-0.8}$ $M_\odot$) and a secondary-to-primary initial mass ratio very close to unity ($q$ = $M_2/M_1$ $\sim$ 1); the SN progenitor's hydrogen envelope had been stripped through binary interaction. We do not detect the binary companion of SN 2015G. For both SNe, the surrounding stellar populations have relatively old ages and argue against any massive WR stars as their progenitors. These results suggest that SNe Ibn may have lower-mass origins in interacting binaries. As a result, they also provide evidence that the giant eruptions commonly seen in massive luminous blue variables (LBVs) can also occur in much lower-mass, stripped-envelope stars just before core collapse., Comment: 22 pages, 18 figures, accepted for publication in MNRAS

We observed seven epochs of spectropolarimetry in optical wavelengths for the Type IIb SN 2011hs, ranging from -3 to +40 days with respect to V -band maximum. A high degree of interstellar polarization was detected (up to ~3 percent), with a peak lying blueward of 4500A. Similar behaviours have been seen in some Type Ia SNe, but had never been observed in a Type IIb. We find that it is most likely the result of a relative enhancement of small silicate grains in the vicinity of the SN. Significant intrinsic continuum polarization was recovered at -3 and +2 days (p = 0.55 +\- 0.12 percent and p = 0.75 +\- 0.11 percent, respectively). We discuss the change of the polarization angle across spectral lines and in the continuum as diagnostics for the 3D structure of the ejecta. We see a gradual rotation by about -50 degree in the continuum polarization angle between -2 and +18 days after V - band maximum. A similar rotation in He I {\lambda}5876, H{\alpha} and the Ca II infrared triplet seems to indicate a strong influence of the global geometry on the line polarization features. The differences in the evolution of their respective loops on the Stokes q - u plane suggest that line specific geometries are also being probed. Possible interpretations are discussed and placed in the context of literature. We find that the spectropolarimetry of SN 2011hs is most similar to that of SN 2011dh, albeit with notable differences., Comment: Accepted for publication in MNRAS, 16 pages, 11 figures

Here we confirm Melnick 34, an X-ray bright star in the 30 Doradus region of the Large Magellanic Cloud, as an SB2 binary comprising WN5h+WN5h components. We present orbital solutions using 26 epochs of VLT/UVES spectra and 22 epochs of archival Gemini/GMOS spectra. Radial-velocity monitoring and automated template fitting methods both reveal a similar high eccentricity system with a mass ratio close to unity, and an orbital period in agreement with the 155.1 +/- 1 day X-ray light curve period previously derived by Pollock et al. Our favoured solution derived an eccentricity of 0.68 +/- 0.02 and mass ratio of 0.92 +/- 0.07, giving minimum masses of Ma_sin^{3}(i) = 65 +/- 7 Msun and Mb_sin^{3}(i) = 60 +/- 7 Msun. Spectral modelling using WN5h templates with CMFGEN reveals temperatures of T ~53 kK for each component and luminosities of log(La/Lsun) = 6.43 +/- 0.08 and log(Lb/Lsun) = 6.37 +/- 0.08, from which BONNSAI evolutionary modelling gives masses of Ma = 139 (+21,-18) Msun and Mb = 127 (+17,-17) Msun and ages of ~0.6 Myrs. Spectroscopic and dynamic masses would agree if Mk34 has an inclination of i ~50{\deg}, making Mk34 the most massive binary known and an excellent candidate for investigating the properties of colliding wind binaries. Within 2-3 Myrs, both components of Mk34 are expected to evolve to stellar mass black holes which, assuming the binary system survives, would make Mk34 a potential binary black hole merger progenitor and gravitational wave source., Comment: 21 pages, 18 figures, accepted for publication in MNRAS

We present observations and analysis of 18 stripped-envelope supernovae observed during 2013 -- 2018. This sample consists of 5 H/He-rich SNe, 6 H-poor/He-rich SNe, 3 narrow lined SNe Ic and 4 broad lined SNe Ic. The peak luminosity and characteristic time-scales of the bolometric light curves are calculated, and the light curves modelled to derive 56Ni and ejecta masses (MNi and Mej). Additionally, the temperature evolution and spectral line velocity-curves of each SN are examined. Analysis of the [O I] line in the nebular phase of eight SNe suggests their progenitors had initial masses $=2.8\pm{1.5}$ Msun which further strengthens the evidence that SE-SNe arise from low mass progenitors which are typically $$ and lack of clear bimodality in the distribution implies $, Comment: Accepted for publication in MNRAS, no changes to the previous submission

We present the spectroscopic analysis of 333 OB-type stars extracted from VLT-MUSE observations of the central 30 x 30 pc of NGC 2070 in the Tarantula Nebula on the Large Magellanic Cloud, the majority of which are analysed for the the first time. The distribution of stars in the spectroscopic Hertzsprung-Russell diagram (sHRD) shows 281 stars in the main sequence. We find two groups in the main sequence, with estimated ages of 2.1$\pm$0.8 and 6.2$\pm$2 Myr. A subgroup of 52 stars is apparently beyond the main sequence phase, which we consider to be due to emission-type objects and/or significant nebular contamination affecting the analysis. As in previous studies, stellar masses derived from the sHRD are systematically larger than those obtained from the conventional HRD, with the differences being largest for the most massive stars. Additionally, we do not find any trend between the estimated projected rotational velocity and evolution in the sHRD. The projected rotational velocity distribution presents a tail of fast rotators that resembles findings in the wider population of 30 Doradus. We use published spectral types to calibrate the HeI$\lambda$4921/HeII$\lambda$5411 equivalent-width ratio as a classification diagnostic for early-type main sequence stars when the classical blue-visible region is not observed. Our model-atmosphere analyses demonstrate that the resulting calibration is well correlated with effective temperature., Comment: 19 pages, 12 figures. Accepted in A&A

The blue compact dwarf galaxy NGC 5253 hosts a very young starburst containing twin nuclear star clusters, separated by a projected distance of 5 pc. One cluster (#5) coincides with the peak of the Hα emission and the other (#11) with a massive ultracompact H ii region. A recent analysis of these clusters shows that they have a photometric age of 1 ± 1 Myr, in apparent contradiction with the age of 3–5 Myr inferred from the presence of Wolf-Rayet features in the cluster #5 spectrum. We examine Hubble Space Telescope ultraviolet and Very Large Telescope optical spectroscopy of #5 and show that the stellar features arise from very massive stars (VMSs), with masses greater than 100 M ⊙, at an age of 1–2 Myr. We further show that the very high ionizing flux from the nuclear clusters can only be explained if VMSs are present. We investigate the origin of the observed nitrogen enrichment in the circumcluster ionized gas and find that the excess N can be produced by massive rotating stars within the first 1 Myr. We find similarities between the NGC 5253 cluster spectrum and those of metal-poor, high-redshift galaxies. We discuss the presence of VMSs in young, star-forming galaxies at high redshift; these should be detected in rest-frame UV spectra to be obtained with the James Webb Space Telescope. We emphasize that population synthesis models with upper mass cutoffs greater than 100 M ⊙ are crucial for future studies of young massive star clusters at all redshifts.

The upper mass limit of stars remains an open question in astrophysics. Here we discuss observations of the most massive stars (greater than 100 solar masses) in the local universe and how the observations fit in with theoretical predictions. In particular, the Large Magellanic Cloud plays host to numerous very massive stars, making it an ideal template to study the roles that environment, metallicity, and multiplicity play in the formation and evolution of the most massive stars. We will discuss the work that is instrumental in laying the groundwork for interpreting future observations by James Webb of starburst regions in the high redshift universe.

Despite their rarity, massive stars dominate the ecology of galaxies via their strong, radiatively-driven winds throughout their lives and as supernovae in their deaths. However, their evolution and subsequent impact on their environment can be significantly affected by the presence of a magnetic field. While recent studies indicate that about 7% of OB stars in the Milky Way host strong, stable, organised (fossil) magnetic fields at their surfaces, little is known about the fields of very massive stars, nor the magnetic properties of stars outside our Galaxy. We aim to continue searching for strong magnetic fields in a diverse set of massive and very massive stars (VMS) in the Large and Small Magellanic Clouds (LMC/SMC), and we evaluate the overall capability of FORS2 to usefully search for and detect stellar magnetic fields in extra-galactic environments. We have obtained FORS2 spectropolarimetry of a sample of 41 stars, which principally consist of spectral types B, O, Of/WN, WNh, and classical WR stars in the LMC and SMC. Four of our targets are Of?p stars; one of them was just recently discovered. Each spectrum was analysed to infer the longitudinal magnetic field. No magnetic fields were formally detected in our study, although Bayesian statistical considerations suggest that the Of?p star SMC159-2 is magnetic with a dipolar field of the order of 2.4 to 4.4kG. In addition, our first constraints of magnetic fields in VMS provide interesting insights into the formation of the most massive stars in the Universe., Comment: Accepted by A&A

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

The massive evolved Wolf-Rayet stars sometimes occur in colliding-wind binary systems in which dust plumes are formed as a result of the collision of stellar winds. These structures are known to encode the parameters of the binary orbit and winds. Here, we report observations of a previously undiscovered Wolf-Rayet system, 2XMM J160050.7-514245, with a spectroscopically determined wind speed of $\approx$3400 km s$^{-1}$. In the thermal infrared, the system is adorned with a prominent $\approx$12$''$ spiral dust plume, revealed by proper motion studies to be expanding at only $\approx$570 km s$^{-1}$. As the dust and gas appear coeval, these observations are inconsistent with existing models of the dynamics of such colliding wind systems. We propose that this contradiction can be resolved if the system is capable of launching extremely anisotropic winds. Near-critical stellar rotation is known to drive such winds, suggesting this Wolf-Rayet system as a potential Galactic progenitor system for long-duration gamma-ray bursts., Accepted to Nature Astronomy on 2018 Oct 1. 27 pages, 9 figures, and 4 tables

A previous spectroscopic study identified the very massive O2 III star VFTS 16 in the Tarantula Nebula as a runaway star based on its peculiar line-of-sight velocity. We use the Gaia DR2 catalog to measure the relative proper motion of VFTS 16 and nearby bright stars to test if this star might have been ejected from the central cluster, R136, via dynamical ejection. We find that the position angle and magnitude of the relative proper motion (0.338 +/- 0.046 mas/yr, or approximately 80 +\- 11 km/s) of VFTS 16 are consistent with ejection from R136 approximately 1.5 +/- 0.2 Myr ago, very soon after the cluster was formed. There is some tension with the presumed age of VFTS 16 that, from published stellar parameters, cannot be greater than 0.9 +0.3/-0.2 Myr. Older ages for this star would appear to be prohibited due to the absence of He I lines in its optical spectrum, since this sets a firm lower limit on its effective temperature. The dynamical constraints may imply an unusual evolutionary history for this object, perhaps indicating it is a merger product. Gaia DR2 also confirms that another very massive star in the Tarantula Nebula, VFTS 72 (alias BI253; O2 III-V(n)((f*)), is also a runaway on the basis of its proper motion as measured by Gaia. While its tangential proper motion (0.392 +/-0.062 mas/yr or 93 +/-15 km/s) would be consistent with dynamical ejection from R136 approximately 1 Myr ago, its position angle is discrepant with this direction at the 2$\sigma$ level. From their Gaia DR2 proper motions we conclude that the two ~100 solar mass O2 stars, VFTS 16 and VFTS72, are fast runaway stars, with space velocities of around 100 km/s relative to R136 and the local massive star population. The dynamics of VFTS16 are consistent with it having been ejected from R136, and this star therefore sets a robust lower limit on the age of the central cluster of ~1.3 Myr., Comment: Version 2 is now accepted by A&A with minor updates from V1

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

The empirical upper luminosity boundary $L_{\rm max}$ of cool supergiants, often referred to as the Humphreys-Davidson limit, is thought to encode information on the general mass-loss behaviour of massive stars. Further, it delineates the boundary at which single stars will end their lives stripped of their hydrogen-rich envelope, which in turn is a key factor in the relative rates of Type-II to Type-Ibc supernovae from single star channels. In this paper we have revisited the issue of $L_{\rm max}$ by studying the luminosity distributions of cool supergiants (SGs) in the Large and Small Magellanic Clouds (LMC/SMC). We assemble samples of cool SGs in each galaxy which are highly-complete above $\log L/L_{\odot}$=5.0, and determine their spectral energy distributions from the optical to the mid-infrared using modern multi-wavelength survey data. We show that in both cases $L_{\rm max}$ appears to be lower than previously quoted, and is in the region of $\log L/L_{\odot}$=5.5. There is no evidence for $L_{\rm max}$ being higher in the SMC than in the LMC, as would be expected if metallicity-dependent winds were the dominant factor in the stripping of stellar envelopes. We also show that $L_{\rm max}$ aligns with the lowest luminosity of single nitrogen-rich Wolf-Rayet stars, indicating of a change in evolutionary sequence for stars above a critical mass. From population synthesis analysis we show that the Geneva evolutionary models greatly over-predict the numbers of cool SGs in the SMC. We also argue that the trend of earlier average spectral types of cool SGs in lower metallicity environments represents a genuine shift to hotter temperatures. Finally, we use our new bolometric luminosity measurements to provide updated bolometric corrections for cool supergiants., Comment: 11 pages, 6 figs. Accepted for publication in MNRAS

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

Stars stripped of their envelope through interaction in a binary are generally not considered when accounting for ionizing radiation from stellar populations, despite the expectation that stripped stars emit hard ionizing radiation, form frequently and live 10-100 times longer than single massive stars. We compute the first grid of spectral models specially made for stars stripped in binaries for a range of progenitor masses (2-20$M_{\odot}$) and metallicities covering a wide range. For stripped stars with masses between 0.3-7$M_{\odot}$, we find high effective temperatures (20-100 kK, increasing with mass), small radii (0.2-1$R_{\odot}$) and high bolometric luminosities, comparable to that of their progenitor before stripping. The spectra show a continuous sequence that naturally bridges subdwarf-type stars and Wolf-Rayet like spectra. For intermediate masses we find hybrid spectral classes showing a mixture of absorption and emission lines. These appear for stars with mass loss rates of $10^{-8}-10^{-6}M_{\odot}$/yr, which have semi-transparent atmospheres. At low metallicity, substantial H-rich layers are left at the surface and we predict spectra that resemble O-type stars instead. We obtain spectra undistinguishable from subdwarfs for stripped stars with masses up to 1.7$M_{\odot}$, which questions whether the widely adopted canonical value of 0.47$M_{\odot}$ is uniformly valid. Increasing the observed sample of stripped stars will provide necessary tests for the physics of interaction, internal mixing and stellar winds. We investigate the feasibility to detect stripped stars next to an optically bright companion and recommend searches for their UV excess and possible emission lines, most notably HeII4686 in the optical and HeII1640 in the UV. Our models are publicly available for further investigations or inclusion in spectral synthesis simulations., Comment: Accepted for publication in A&A

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

We introduce VLT-MUSE observations of the central 2$'\times2'$ (30$\times$30 pc) of the Tarantula Nebula in the Large Magellanic Cloud. The observations provide an unprecedented spectroscopic census of the massive stars and ionised gas in the vicinity of R136, the young, dense star cluster located in NGC 2070, at the heart of the richest star-forming region in the Local Group. Spectrophotometry and radial-velocity estimates of the nebular gas (superimposed on the stellar spectra) are provided for 2255 point sources extracted from the MUSE datacubes, and we present estimates of stellar radial velocities for 270 early-type stars (finding an average systemic velocity of 271$\pm$41 km/s). We present an extinction map constructed from the nebular Balmer lines, with electron densities and temperatures estimated from intensity ratios of the [SII], [NII], and [SIII] lines. The interstellar medium, as traced by H$\alpha$ and [NII] $\lambda$6583, provides new insights in regions where stars are probably forming. The gas kinematics are complex, but with a clear bi-modal, blue- and red-shifted distribution compared to the systemic velocity of the gas centred on R136. Interesting point-like sources are also seen in the eastern cavity, western shell, and around R136; these might be related to phenomena such as runaway stars, jets, formation of new stars, or the interaction of the gas with the population of Wolf--Rayet stars. Closer inspection of the core reveals red-shifted material surrounding the strongest X-ray sources, although we are unable to investigate the kinematics in detail as the stars are spatially unresolved in the MUSE data. Further papers in this series will discuss the detailed stellar content of NGC 2070 and its integrated stellar and nebular properties., Comment: 50 pages, 13 figures. Accepted 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)

Previous analyses of the spectra of OB-type stars in the Magellanic Clouds have identified targets with low projected rotational velocities and relatively high nitrogen abundances. The VFTS obtained spectroscopy for B-type 434 stars. We have estimated atmospheric parameters and nitrogen abundances using TLUSTY model atmospheres for 54 B-type targets that appear to be single, have projected rotational velocities, Accepted in Astronomy and Astrophyics. Abstract significantly shortened - see paper for full abstract

We have carried out a spectroscopic variability survey of some of the most massive stars in the Arches cluster, using K-band observations obtained with SINFONI on the VLT. One target, F2, exhibits substantial changes in radial velocity; in combination with new KMOS and archival SINFONI spectra, its primary component is found to undergo radial velocity variation with a period of 10.483+/-0.002 d and an amplitude of ~350 km/s. A secondary radial velocity curve is also marginally detectable. We reanalyse archival NAOS-CONICA photometric survey data in combination with our radial velocity results to confirm this object as an eclipsing SB2 system, and the first binary identified in the Arches. We model it as consisting of an 82+/-12 M_sun WN8-9h primary and a 60+/-8 M_sun O5-6 Ia+ secondary, and as having a slightly eccentric orbit, implying an evolutionary stage prior to strong binary interaction. As one of four X-ray bright Arches sources previously proposed as colliding-wind massive binaries, it may be only the first of several binaries to be discovered in this cluster, presenting potential challenges to recent models for the Arches' age and composition. It also appears to be one of the most massive binaries detected to date; the primary's calculated initial mass of >~120 M_sun would arguably make this the most massive binary known in the Galaxy., 15 pages, 18 figures, accepted for publication in Astronomy & Astrophysics. Paper II in a series; for Paper I see arXiv:1803.09567

The Wolf-Rayet star Mk 34 was observed more than 50 times as part of the deep T-ReX Chandra ACIS-I X-ray imaging survey of the Tarantula Nebula in the Large Magellanic Cloud conducted between 2014 May and 2016 January. Its brightness showed one bright maximum and repeated faint minima which help define an X-ray recurrence time of $155.1\pm0.1$ days that is probably the orbital period of an eccentric binary system. The maximum immediately precedes the minimum in the folded X-ray light curve as confirmed by new Swift XRT observations. Notwithstanding its extreme median luminosity of $1.2\times10^{35}\mathrm{erg}~\mathrm{s}^{-1}$, which makes it over an order of magnitude brighter than comparable stars in the Milky Way, Mk 34 is almost certainly a colliding-wind binary system. Its spectrum shows phase-related changes of luminosity and absorption that are probably related to the orbital dynamics of two of the most massive stars known., Comment: 10 pages, 8 figures, 5 tables

IAU Commission 29 - Stellar Spectra has been one of the IAU commissions from the onset, until its dissolution at the most recent IAU General Assembly in Honolulu in 2015. This commission belonged to IAU Division G (“Stars and Stellar Physics”), the latter committed with fostering research in stellar astrophysics. Within the general field of stellar astrophysics, stellar spectroscopy plays a key role, as stellar spectra are a powerful tool providing a view into the detailed physical properties of stars and the physical processes occuring within them.