Showing 1–50 of 179 results for author: Crowther, P A

X-Shooting ULLYSES: Massive Stars at low metallicity -- IV. Spectral analysis methods and exemplary results for O stars

Authors: A. A. C. Sander, J. -C. Bouret, M. Bernini-Peron, J. Puls, F. Backs, S. R. Berlanas, J. M. Bestenlehner, S. A. Brands, A. Herrero, F. Martins, O. Maryeva, D. Pauli, V. Ramachandran, P. A. Crowther, V. M. A. Gómez-González, A. C. Gormaz-Matamala, W. -R. Hamann, D. J. Hillier, R. Kuiper, C. J. K. Larkin, R. R. Lefever, A. Mehner, F. Najarro, L. M. Oskinova, E. C. Schösser , et al. (4 additional authors not shown)

Abstract: CONTEXT: The spectral analysis of hot, massive stars is a fundamental astrophysical method to obtain their intrinsic properties and their feedback. Quantitative spectroscopy for hot, massive stars requires detailed numerical modeling of the atmosphere and an iterative treatment to obtain the best solution within a given framework. AIMS: We present an overview of different techniques for the quanti… ▽ More CONTEXT: The spectral analysis of hot, massive stars is a fundamental astrophysical method to obtain their intrinsic properties and their feedback. Quantitative spectroscopy for hot, massive stars requires detailed numerical modeling of the atmosphere and an iterative treatment to obtain the best solution within a given framework. AIMS: We present an overview of different techniques for the quantitative spectroscopy of hot stars employed within the X-Shooting ULLYSES collaboration, from grid-based approaches to tailored fits. By performing a blind test, we gain an overview about the similarities and differences of the resulting parameters. Our study aims to provide an overview of the parameter spread caused by different approaches. METHODS: For three different stars from the sample (SMC O5 star AzV 377, LMC O7 star Sk -69 50, and LMC O9 star Sk -66 171), we employ different atmosphere codes (CMFGEN, Fastwind, PoWR) and strategies to determine their best-fitting model. For our analyses, UV and optical spectra are used to derive the properties with some methods relying purely on optical data for comparison. To determine the overall spectral energy distribution, we further employ additional photometry from the literature. RESULTS: Effective temperatures for each of three sample stars agree within 3 kK while the differences in log g can be up to 0.2 dex. Luminosity differences of up to 0.1 dex result from different reddening assumptions, which seem to be larger for the methods employing a genetic algorithm. All sample stars are nitrogen-enriched. CONCLUSIONS: We find a reasonable agreement between the different methods. Tailored fitting tends to be able to minimize discrepancies obtained with more course or automatized treatments. UV spectral data is essential for the determination of realistic wind parameters. For one target (Sk -69 50), we find clear indications of an evolved status. △ Less

Submitted 3 July, 2024; originally announced July 2024.

Comments: 18+15 pages, 21+4 figures, under review at A&A, condensed abstract

X-Shooting ULLYSES: Massive stars at low metallicity -- V. Effect of metallicity on surface abundances of O stars

Authors: F. Martins, J. -C. Bouret, D. J. Hillier, S. A. Brands, P. A. Crowther, A. Herrero, F. Najarro, D. Pauli, J. Puls, V. Ramachandran, A. A. C. Sander, J. S. Vink, the XshootU collaboration

Abstract: Massive stars rotate faster, on average, than lower mass stars. Stellar rotation triggers hydrodynamical instabilities which transport angular momentum and chemical species from the core to the surface. Models of high-mass stars that include these processes predict that chemical mixing is stronger at lower metallicity. We aim to test this prediction by comparing the surface abundances of massive s… ▽ More Massive stars rotate faster, on average, than lower mass stars. Stellar rotation triggers hydrodynamical instabilities which transport angular momentum and chemical species from the core to the surface. Models of high-mass stars that include these processes predict that chemical mixing is stronger at lower metallicity. We aim to test this prediction by comparing the surface abundances of massive stars at different metallicities. We performed a spectroscopic analysis of single O stars in the Magellanic Clouds (MCs) based on the ULLYSES and XshootU surveys. We determined the fundamental parameters and helium, carbon, nitrogen, and oxygen surface abundances of 17 LMC and 17 SMC non-supergiant O6-9.5 stars. We complemented these determinations by literature results for additional MCs and also Galactic stars to increase the sample size and metallicity coverage. We investigated the differences in the surface chemical enrichment at different metallicities and compared them with predictions of three sets of evolutionary models. Surface abundances are consistent with CNO-cycle nucleosynthesis. The maximum surface nitrogen enrichment is stronger in MC stars than in Galactic stars. Nitrogen enrichment is also observed in stars with higher surface gravities in the SMC than in the Galaxy. This trend is predicted by models that incorporate chemical transport caused by stellar rotation. The distributions of projected rotational velocities in our samples are likely biased towards slow rotators. A metallicity dependence of surface abundances is demonstrated. The analysis of larger samples with an unbiased distribution of projected rotational velocities is required to better constrain the treatment of chemical mixing and angular momentum transport in massive single and binary stars. △ Less

Submitted 2 May, 2024; originally announced May 2024.

Comments: 15 pages + appendix. Accepted in Astronomy & Astrophysics

Mapping the core of the Tarantula Nebula with VLT-MUSE. III. A template for metal-poor starburst regions in the visual and far-ultraviolet

Authors: Paul A Crowther, N Castro

Abstract: We present the integrated VLT-MUSE spectrum of the central 2'x2' (30x30 pc$^{2}$) of NGC 2070, the dominant giant HII region of the Tarantula Nebula in the Large Magellanic Cloud, together withan empirical far-ultraviolet spectrum constructed via LMC template stars from the ULLYSES survey and Hubble Tarantula Treasury Project UV photometry. NGC 2070 provides a unique opportunity to compare results… ▽ More We present the integrated VLT-MUSE spectrum of the central 2'x2' (30x30 pc$^{2}$) of NGC 2070, the dominant giant HII region of the Tarantula Nebula in the Large Magellanic Cloud, together withan empirical far-ultraviolet spectrum constructed via LMC template stars from the ULLYSES survey and Hubble Tarantula Treasury Project UV photometry. NGC 2070 provides a unique opportunity to compare results from individual stellar populations (e.g. VLT FLAMES Tarantula Survey) in a metal-poor starburst region to the integrated results from the population synthesis tools Starburst99, Charlot & Bruzual and BPASS. The metallicity of NGC 2070 inferred from standard nebular strong line calibrations is 0.4$\pm$0.1 dex lower than obtained from direct methods. The Halpha inferred age of 4.2 Myr from Starburst99 is close to the median age of OB stars within the region, although individual stars span a broad range of 1-7 Myr. The inferred stellar mass is close to that obtained for the rich star cluster R136 within NGC 2070, although this contributes only 21% to the integrated far-UV continuum. HeII 1640 emission is dominated by classical WR stars and main sequence WNh+Of/WN stars. 18% of the NGC~2070 far UV continuum flux arises from very massive stars with >100 Msun, including several very luminous Of supergiants. None of the predicted population synthesis models at low metallicities are able to successfully reproduce the far-UV spectrum of NGC 2070. We attribute issues to the treatment of mass-loss in very massive stars, the lack of contemporary empirical metal-poor templates, plus WR stars produced via binary evolution. △ Less

Submitted 27 November, 2023; v1 submitted 13 November, 2023; originally announced November 2023.

Comments: 15 pages, 16 figures, plus Appendix. Accepted to MNRAS

Oxygen abundance of gamma Vel from [O III] 88um Herschel/PACS spectroscopy

Authors: Paul A Crowther, M J Barlow, P Royer, D J Hillier, J M Bestenlehner, P W Morris, R Wesson

Abstract: We present Herschel PACS spectroscopy of the [O III] 88.4um fine-structure line in the nearby WC8+O binary system gamma Vel to determine its oxygen abundance. The critical density of this line corresponds to several 10^5 R* such that it is spatially extended in PACS observations at the 336 pc distance to gamma Vel. Two approaches are used, the first involving a detailed stellar atmosphere analysis… ▽ More We present Herschel PACS spectroscopy of the [O III] 88.4um fine-structure line in the nearby WC8+O binary system gamma Vel to determine its oxygen abundance. The critical density of this line corresponds to several 10^5 R* such that it is spatially extended in PACS observations at the 336 pc distance to gamma Vel. Two approaches are used, the first involving a detailed stellar atmosphere analysis of gamma Vel using CMFGEN, extending to Ne ~ 10^0 cm^-3 in order to fully sample the line formation region of [O III] 88.4um. The second approach involves the analytical model introduced by Barlow et al. and revised by Dessart et al, additionally exploiting ISO LWS spectroscopy of [O III] 51.8um. We obtain higher luminosities for the WR and O components of gamma Vel with respect to De Marco et al, log L/L_sun = 5.31 and 5.56, respectively, due to the revised (higher) interferometric distance. We obtain an oxygen mass fraction of X_O = 1.0+/- 0.3% for an outer wind volume filling factor of f = 0.5+/-0.25, favouring either standard or slightly reduced Kunz et al. rates for the ^12C(alpha, gamma)^16O reaction from comparison with BPASS binary population synthesis models. We also revisit neon and sulphur abundances in the outer wind of gamma Vel from ISO SWS spectroscopy of [S IV] 10.5um [Ne II] 12.8um and [Ne III] 15.5um. The inferred neon abundance X_Ne = 2.0-0.6+0.4%, is in excellent agreement with BPASS predictions, while the.sulphur abundance of X_S = 0.04 +/- 0.01% agrees with the solar abundance, as expected for unprocessed elements. △ Less

Submitted 10 January, 2024; v1 submitted 23 October, 2023; originally announced October 2023.

Comments: 12 pages, 9 figures, plus Appendices, submitted to MNRAS, revised in response to referee report

Constraints on the multiplicity of the most massive stars known: R136 a1, a2, a3, and c

Authors: T. Shenar, H. Sana, P. A. Crowther, K. A. Bostroem, L. Mahy, F. Najarro, L. Oskinova, A. A. C. Sander

Abstract: The most massive stars known to date are R 136 a1, a2, a3, and c within the central cluster R 136a of the Tarantula nebula in the Large Magellanic Cloud (LMC), with reported masses in excess of 150-200$M_\odot$. However, the mass estimation of these stars relies on the assumption that they are single. We collected three epochs of spectroscopy for R 136 a1, a2, a3, and c with the Space Telescope Im… ▽ More The most massive stars known to date are R 136 a1, a2, a3, and c within the central cluster R 136a of the Tarantula nebula in the Large Magellanic Cloud (LMC), with reported masses in excess of 150-200$M_\odot$. However, the mass estimation of these stars relies on the assumption that they are single. We collected three epochs of spectroscopy for R 136 a1, a2, a3, and c with the Space Telescope Imaging Spectrograph (STIS) of the Hubble Space Telescope (HST) in the years 2020-2021 to probe potential radial-velocity (RV) variations. We combine these epochs with an additional HST/STIS observation taken in 2012. We use cross-correlation to quantify the RVs, and establish constraints on possible companions to these stars up to periods of ~10 yr. Objects are classified as binaries when the peak-to-peak RV shifts exceed 50 km/s, and when the RV shift is significant with respect to errors. R 136 a1, a2, and a3 do not satisfy the binary criteria and are thus classified as putatively single, although formal peak-to-peak RV variability on the level 40 km/s is noted for a3. Only R 136 c is classified as binary, in agreement with literature. We can generally rule out massive companions (M2 > ~50 Msun) to R 136 a1, a2, and a3 out to orbital periods of < 1 yr (separations < 5 au) at 95% confidence, or out to tens of years (separations < ~100 au) at 50% confidence. Highly eccentric binaries (e > ~0.9) or twin companions with similar spectra could evade detection down to shorter periods (> ~10 d), though their presence is not supported by the relative X-ray faintness of R 136 a1, a2, and a3. We derive a preliminary orbital solution with a 17.2 d period for the X-ray bright binary R 136 c, though more data are needed to conclusively derive its orbit. Our study supports a lower bound of 150-200 $M_\odot$ on the upper-mass limit at LMC metallicity △ Less

Submitted 22 September, 2023; originally announced September 2023.

Comments: Accepted to A&A

Journal ref: A&A 679, A36 (2023)

Spectroscopic analysis of hot, massive stars in large spectroscopic surveys with de-idealised models

Authors: J. M. Bestenlehner, T. Enßlin, M. Bergemann, P. A. Crowther, M. Greiner, M. Selig

Abstract: Upcoming large-scale spectroscopic surveys with e.g. WEAVE and 4MOST will provide thousands of spectra of massive stars, which need to be analysed in an efficient and homogeneous way. Usually, studies of massive stars are limited to samples of a few hundred objects which pushes current spectroscopic analysis tools to their limits because visual inspection is necessary to verify the spectroscopic f… ▽ More Upcoming large-scale spectroscopic surveys with e.g. WEAVE and 4MOST will provide thousands of spectra of massive stars, which need to be analysed in an efficient and homogeneous way. Usually, studies of massive stars are limited to samples of a few hundred objects which pushes current spectroscopic analysis tools to their limits because visual inspection is necessary to verify the spectroscopic fit. Often uncertainties are only estimated rather than derived and prior information cannot be incorporated without a Bayesian approach. In addition, uncertainties of stellar atmospheres and radiative transfer codes are not considered as a result of simplified, inaccurate or incomplete/missing physics or, in short, idealised physical models. Here, we address the question of "How to compare an idealised model of complex objects to real data?" with an empirical Bayesian approach and maximum a {\it posterior} approximations. We focus on application to large scale optical spectroscopic studies of complex astrophysical objects like stars. More specifically, we test and verify our methodology on samples of OB stars in 30 Doradus region of the Large Magellanic Clouds using a grid of FASTWIND model atmospheres. Our spectroscopic model de-idealisation analysis pipeline takes advantage of the statistics that large samples provide by determining the model error to account for the idealised stellar atmosphere models, which are included into the error budget. The pipeline performs well over a wide parameter space and derives robust stellar parameters with representative uncertainties. △ Less

Submitted 12 September, 2023; originally announced September 2023.

Comments: Submitted to MNRAS, 21 pages, 9 figures

Stellar properties of observed stars stripped in binaries in the Magellanic Clouds

Authors: Y. Gotberg, M. R. Drout, A. P. Ji, J. H. Groh, B. A. Ludwig, P. A. Crowther, N. Smith, A. de Koter, S. E. de Mink

Abstract: Massive stars (~8-25Msun) stripped of their hydrogen-rich envelopes via binary interaction are thought to be the main progenitors for merging neutron stars and stripped-envelope supernovae. We recently presented the discovery of the first set of such stripped stars in a companion paper. Here, we fit the spectra of ten stars with new atmosphere models in order to constrain their stellar properties… ▽ More Massive stars (~8-25Msun) stripped of their hydrogen-rich envelopes via binary interaction are thought to be the main progenitors for merging neutron stars and stripped-envelope supernovae. We recently presented the discovery of the first set of such stripped stars in a companion paper. Here, we fit the spectra of ten stars with new atmosphere models in order to constrain their stellar properties precisely. We find that the stellar properties align well with the theoretical expectations from binary evolution models for helium-core burning envelope-stripped stars. The fits confirm that the stars have high effective temperatures (Teff~50-100kK), high surface gravities (log g ~5), and hydrogen-poor/helium-rich surfaces (X(H, surf)~0-0.4) while showing for the first time a range of bolometric luminosities (10^3-10^5 Lsun), small radii (~0.5-1Rsun), and low Eddington factors (Gamma_e~0.006-0.4). Using these properties, we derive intermediate current masses (~1-8Msun), which suggest that their progenitors were massive stars (~5-25Msun) and that a subset will reach core-collapse, leaving behind neutron stars or black holes. Using the model fits, we also estimate the emission rates of ionizing photons for these stars, which agree well with previous model expectations. Further, by computing models for a range of mass-loss rates, we find that the stellar winds are weaker than predicted by any existing scheme (Mdot(wind)<~ 1e-9 Msun/yr). The properties of this first sample of intermediate mass helium stars suggest they both contain progenitors of type Ib and IIb supernovae, and provide important benchmarks for binary evolution and population synthesis models. △ Less

Submitted 30 June, 2023; originally announced July 2023.

Comments: Accepted for publication in ApJ. Spectral models available via zenodo as described in text. Please also see a companion manuscript detailing the discovery of these systems in today's arXiv posting

X-Shooting ULLYSES: massive stars at low metallicity. I. Project Description

Authors: Jorick S. Vink, A. Mehner, P. A. Crowther, A. Fullerton, M. Garcia, F. Martins, N. Morrell, L. M. Oskinova, N. St-Louis, A. ud-Doula, A. A. C. Sander, H. Sana, J. -C. Bouret, B. Kubatova, P. Marchant, L. P. Martins, A. Wofford, J. Th. van Loon, O. Grace Telford, Y. Gotberg, D. M. Bowman, C. Erba, V. M. Kalari, M. Abdul-Masih, T. Alkousa , et al. (56 additional authors not shown)

Abstract: Observations of individual massive stars, super-luminous supernovae, gamma-ray bursts, and gravitational-wave events involving spectacular black-hole mergers, indicate that the low-metallicity Universe is fundamentally different from our own Galaxy. Many transient phenomena will remain enigmatic until we achieve a firm understanding of the physics and evolution of massive stars at low metallicity… ▽ More Observations of individual massive stars, super-luminous supernovae, gamma-ray bursts, and gravitational-wave events involving spectacular black-hole mergers, indicate that the low-metallicity Universe is fundamentally different from our own Galaxy. Many transient phenomena will remain enigmatic until we achieve a firm understanding of the physics and evolution of massive stars at low metallicity (Z). The Hubble Space Telescope has devoted 500 orbits to observe 250 massive stars at low Z in the ultraviolet (UV) with the COS and STIS spectrographs under the ULLYSES program. The complementary ``X-Shooting ULLYSES'' (XShootU) project provides enhanced legacy value with high-quality optical and near-infrared spectra obtained with the wide-wavelength coverage X-shooter spectrograph at ESO's Very Large Telescope. We present an overview of the XShootU project, showing that combining ULLYSES UV and XShootU optical spectra is critical for the uniform determination of stellar parameters such as effective temperature, surface gravity, luminosity, and abundances, as well as wind properties such as mass-loss rates in function of Z. As uncertainties in stellar and wind parameters percolate into many adjacent areas of Astrophysics, the data and modelling of the XShootU project is expected to be a game-changer for our physical understanding of massive stars at low Z. To be able to confidently interpret James Webb Space Telescope (JWST) spectra of the first stellar generations, the individual spectra of low Z stars need to be understood, which is exactly where XShootU can deliver. △ Less

Submitted 1 June, 2023; v1 submitted 10 May, 2023; originally announced May 2023.

Comments: Accepted in A&A - 35 Pages, 12 Figures, 4 Tables, 2 Large Tables

Journal ref: A&A 675, A154 (2023)

X-Shooting ULLYSES: Massive stars at low metallicity. III. Terminal wind speeds of ULLYSES massive stars

Authors: C. Hawcroft, H. Sana, L. Mahy, J. O. Sundqvist, A. de Koter, P. A. Crowther, J. M. Bestenlehner, S. A. Brands, A. David-Uraz, L. Decin, C. Erba, M. Garcia, W. -R. Hamann, A. Herrero, R. Ignace, N. D. Kee, B. Kubátová, R. Lefever, A. Moffat, F. Najarro, L. Oskinova, D. Pauli, R. Prinja, J. Puls, A. A. C. Sander , et al. (4 additional authors not shown)

Abstract: The winds of massive stars have an impact on stellar evolution and on the surrounding medium. The maximum speed reached by these outflows, the terminal wind speed, is a global wind parameter and an essential input for models of stellar atmospheres and feedback. With the arrival of the ULLYSES programme, a legacy UV spectroscopic survey with HST, we have the opportunity to quantify the wind speeds… ▽ More The winds of massive stars have an impact on stellar evolution and on the surrounding medium. The maximum speed reached by these outflows, the terminal wind speed, is a global wind parameter and an essential input for models of stellar atmospheres and feedback. With the arrival of the ULLYSES programme, a legacy UV spectroscopic survey with HST, we have the opportunity to quantify the wind speeds of massive stars at sub-solar metallicity (in the Large and Small Magellanic Clouds, 0.5Z and 0.2Z) at an unprecedented scale. We empirically quantify the wind speeds of a large sample of OB stars, including supergiants, giants, and dwarfs at sub-solar metallicity. Using these measurements, we investigate trends of terminal wind speed with a number of fundamental stellar parameters, namely effective temperature, metallicity, and surface escape velocity. We empirically determined the terminal wind speed for a sample of 149 OB stars in the Magellanic Clouds either by directly measuring the maximum velocity shift of the absorption component of the Civ 1548-1550 line profile, or by fitting synthetic spectra produced using the Sobolev with exact integration method. Stellar parameters were either collected from the literature, obtained using spectral-type calibrations, or predicted from evolutionary models. We find strong trends of terminal wind speed with effective temperature and surface escape speed when the wind is strong enough to cause a saturated P Cygni profile in Civ 1548-1550. We find evidence for a metallicity dependence on the terminal wind speed proportional to Z^0.22+-0.03 when we compared our results to previous Galactic studies. Our results suggest that effective temperature rather than surface escape speed should be used as a straightforward empirical prediction of terminal wind speed and that the observed metallicity dependence is steeper than suggested by earlier works. △ Less

Submitted 24 March, 2023; v1 submitted 21 March, 2023; originally announced March 2023.

Comments: 21 pages, 16 figures, 8 tables. Accepted in A&A

Extinction towards the cluster R136 in the Large Magellanic Cloud: An extinction law from the near-infrared to the ultraviolet

Authors: Sarah A. Brands, Alex de Koter, Joachim M. Bestenlehner, Paul A. Crowther, Lex Kaper, Saida M. Caballero-Nieves, Götz Gräfener

Abstract: 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 an… ▽ More 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. △ Less

Submitted 16 March, 2023; originally announced March 2023.

Comments: Accepted for publication in A&A (18 pages, 12 figures)

Journal ref: A&A 673, A132 (2023)

Bringing Stellar Evolution & Feedback Together: Summary of proposals from the Lorentz Center Workshop, 2022

Authors: Sam Geen, Poojan Agrawal, Paul A. Crowther, B. W. Keller, Alex de Koter, Zsolt Keszthelyi, Freeke van de Voort, Ahmad A. Ali, Frank Backs, Lars Bonne, Vittoria Brugaletta, Annelotte Derkink, Sylvia Ekström, Yvonne A. Fichtner, Luca Grassitelli, Ylva Götberg, Erin R. Higgins, Eva Laplace, Kong You Liow, Marta Lorenzo, Anna F. McLeod, Georges Meynet, Megan Newsome, G. André Oliva, Varsha Ramachandran , et al. (12 additional authors not shown)

Abstract: 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 whit… ▽ More 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 April 2022, and identify key areas where further dialogue can bring about radical changes in how we view the relationship between stars and the universe they live in. △ Less

Submitted 31 January, 2023; originally announced January 2023.

Comments: Accepted to the Publications of the Astronomical Society of the Pacific

Line Luminosities of Galactic and Magellanic Cloud Wolf-Rayet stars

Authors: Paul A. Crowther, G. Rate, Joachim M. Bestenlehner

Abstract: 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 fe… ▽ More 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. △ Less

Submitted 3 February, 2023; v1 submitted 26 January, 2023; originally announced January 2023.

Comments: 15 pages, 10 figures, MNRAS in press. Appendix A (individual WR line luminosities), Appendix B (WR templates) available from https://zenodo.org/record/7573775

The VLT-FLAMES Tarantula Survey: Observational evidence for two distinct populations of massive runaway stars in 30 Doradus

Authors: H. Sana, O. H. Ramírez-Agudelo, V. Hénault-Brunet, L. Mahy, L. A. Almeida, A. de Koter, J. M. Bestenlehner, C. J. Evans, N. Langer, F. R. N. Schneider, P. A. Crowther, S. E. de Mink, A. Herrero, D. J. Lennon, M. Gieles, J. Maíz Apellániz, M. Renzo, E. Sabbi, J. Th. van Loon, J. S. Vink

Abstract: Two main scenarios have been proposed for origin of massive runaway stars -- dynamical ejection or release from a binary at the first core collapse -- but their relative contribution remains debated. Using two large spectroscopic campaigns towards massive stars in 30 Doradus, we aim to provide observational constraints on the properties of the O-type runaway population in the most massive active… ▽ More Two main scenarios have been proposed for origin of massive runaway stars -- dynamical ejection or release from a binary at the first core collapse -- but their relative contribution remains debated. Using two large spectroscopic campaigns towards massive stars in 30 Doradus, we aim to provide observational constraints on the properties of the O-type runaway population in the most massive active star-forming region in the Local group. We use RV measurements of the O-type star populations in 30 Doradus obtained by the VLT-FLAMES Tarantula Survey and the Tarantula Massive Binary Monitoring to identify single and binary O-type runaways. We discuss their rotational properties and qualitatively compare observations with expectations of ejection scenarios. We identify 23 single and one binary O-type runaway objects, most of them outside the main star-forming regions in 30 Doradus. We find an overabundance of rapid rotators (vsini > 200km/s) among the runaway population, providing an explanation of the overabundance of rapidly rotating stars in the 30 Doradus field. Considerations of the projected rotation rates and runaway line-of-sight (los) velocities reveal a conspicuous absence of rapidly rotating (vsini > 210k/ms), fast moving (v_{los} > 60km/s) runaways, and suggest the presence of two different populations of runaway stars: a population of rapidly-spinning but slowly moving runaways and a population of fast moving but slowly rotating ones. These are detected with a ratio close to 2:1 in our sample. We argue that slowly moving but rapidly spinning runaways result from binary ejections, while rapidly moving but slowly spinning runaways could result from dynamical ejections. Given that detection biases will more strongly impact the slow-moving population, our results suggest that the binary evolution scenario dominates the current massive runaway population in 30 Doradus. △ Less

Submitted 24 November, 2022; originally announced November 2022.

Comments: Accepted for publication in A&A Letters; 9 pages, 5 figures

Journal ref: A&A 668, L5 (2022)

A UV census of the environments of stripped-envelope supernovae

Authors: Ning-Chen Sun, Justyn R. Maund, Paul A. Crowther

Abstract: 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… ▽ More 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. △ Less

Submitted 2 March, 2023; v1 submitted 12 September, 2022; originally announced September 2022.

Comments: 15 pages, 9 figures, accepted for publication in MNRAS

X-ray properties of early-type stars in the Tarantula Nebula from T-ReX

Authors: Paul A Crowther, Patrick S Broos, Leisa K Townsley, Andy M T Pollock, Katie A Tehrani, Marc Gagne

Abstract: We reassess the historical $L_{X}/L_{Bol}$ relation for early-type stars from a comparison between T-ReX, the Chandra ACIS X-ray survey of the Tarantula Nebula in the LMC, and contemporary spectroscopic analysis of massive stars obtained primarily from VLT/FLAMES, VLT/MUSE and HST/STIS surveys. For 107 sources in common (some host to multiple stars), the majority of which are bolometrically lumino… ▽ More We reassess the historical $L_{X}/L_{Bol}$ relation for early-type stars from a comparison between T-ReX, the Chandra ACIS X-ray survey of the Tarantula Nebula in the LMC, and contemporary spectroscopic analysis of massive stars obtained primarily from VLT/FLAMES, VLT/MUSE and HST/STIS surveys. For 107 sources in common (some host to multiple stars), the majority of which are bolometrically luminous (40% exceed $10^6 L_{\odot}$), we find an average $\log L_{X} /L_{Bol} = -6.90 \pm 0.65$. Excluding extreme systems Mk 34 (WN5h+WN5h), R140a (WC4+WN6+) and VFTS 399 (O9 IIIn+?), plus four WR sources with anomalously hard X-ray components (R130, R134, R135, Mk 53) and 10 multiple sources within the spatially crowded core of R136a, $\log L_{X}/L_{Bol} = -7.00 \pm 0.49$, in good agreement with Galactic OB stars. No difference is found between single and binary systems, nor between O, Of/WN and WR stars, although there does appear to be a trend towards harder X-ray emission from O dwarfs, through O (super)giants, Of/WN stars and WR stars. The majority of known OB stars in the Tarantula are not detected in the T-ReX point source catalogue, so we have derived upper limits for all undetected OB stars for which log $L_{Bol}/L_{\odot} \geq 5.0$. A survival analysis using detected and upper-limit log $L_{X}/L_{Bol}$ values indicates no significant difference between luminous O stars in the LMC and the Carina Nebula. This analysis suggests that metallicity does not strongly influence $L_{X}/L_{Bol}$. Plasma temperatures for single, luminous O stars in the Tarantula ($\overline{kT_{m}}=1.0$ keV) are higher than counterparts in Carina ($\overline{kT_{m}}=0.5$ keV). △ Less

Submitted 21 July, 2022; v1 submitted 18 July, 2022; originally announced July 2022.

Comments: 14 pages, 8 figures, accepted for MNRAS plus Appendices providing X-ray plasma fits (A1), properties of early-type stars (A2), LMC baseline abundances (B1), plus X-ray upper limits for other luminous early-type stars (S1)

ULLYSES and Complementary Surveys of Massive Stars in the Magellanic Clouds

Authors: Paul A Crowther

Abstract: An overview is provided of the scientific goals of the Magellanic Cloud component of the STScI Directors Discretionary UV initiative ULLYSES, together with the complementary spectroscopic survey XShootU (VLT/Xshooter) and other ancillary datasets. Together, ULLYSES and XShootU permit the first comprehensive, homogeneous study of wind densities and velocities in metal-poor massive stars, plus UV/op… ▽ More An overview is provided of the scientific goals of the Magellanic Cloud component of the STScI Directors Discretionary UV initiative ULLYSES, together with the complementary spectroscopic survey XShootU (VLT/Xshooter) and other ancillary datasets. Together, ULLYSES and XShootU permit the first comprehensive, homogeneous study of wind densities and velocities in metal-poor massive stars, plus UV/optical spectroscopic libraries for population synthesis models and a large number of interstellar sight-lines towards the Magellanic Clouds. △ Less

Submitted 18 July, 2022; originally announced July 2022.

Comments: 10 pages, 9 figures, review paper for IAU Symp 361 "Massive Stars Near and Far" (eds. Nicole St-Louis, Jorick Vink, Jonathan Mackey)

An X-ray quiet black hole born with a negligible kick in a massive binary within the Large Magellanic Cloud

Authors: Tomer Shenar, Hugues Sana, Laurent Mahy, Kareem El-Badry, Pablo Marchant, Norbert Langer, Calum Hawcroft, Matthias Fabry, Koushik Sen, Leonardo A. Almeida, Michael Abdul-Masih, Julia Bodensteiner, Paul A. Crowther, Mark Gieles, Mariusz Gromadzki, Vincent Henault-Brunet, Artemio Herrero, Alex de Koter, Patryk Iwanek, Szymon Kozłowski, Daniel J. Lennon, Jesus Maız Apellaniz, Przemysław Mroz, Anthony F. J. Moffat, Annachiara Picco , et al. (13 additional authors not shown)

Abstract: 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 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-ty… ▽ More 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 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-sigma 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. △ Less

Submitted 15 July, 2022; originally announced July 2022.

Comments: Accepted to Nature Astronomy, 64 pages, 15 figures, 2 tables; ESO press release: https://www.eso.org/public/news/eso2210/; Nat Asr paper URL: https://www.nature.com/articles/s41550-022-01730-y

The Tarantula Massive Binary Monitoring VI: Characterisation of hidden companions in 51 single-lined O-type binaries, a flat mass-ratio distribution, and black-hole binary candidates

Authors: Tomer Shenar, Hugues Sana, Laurent Mahy, Jesus Maiz Apellaniz, Paul A. Crowther, Mariusz Gromadzki, Artemio Herrero, Norbert Langer, Pablo Marchant, Fabian R. N. Schneider, Koushik Sen, Igor Soszynski, S. Toonen

Abstract: 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 T… ▽ More 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. △ Less

Submitted 5 October, 2022; v1 submitted 15 July, 2022; originally announced July 2022.

Comments: 41 pages (14 main article + 27 appendix), accepted to A&A

Journal ref: A&A 665, A148 (2022)

A hot and luminous source at the site of the fast transient AT2018cow at 2-3 years after its explosion

Authors: Ning-Chen Sun, Justyn R. Maund, Paul A. Crowther, Liang-Duan Liu

Abstract: 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$α$ emission. It has a very stable brightness between the two epochs and a very blue spectral energy distr… ▽ More 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$α$ emission. It has a very stable brightness between the two epochs and a very blue spectral energy distribution (SED) consistent with $f_λ$ $\propto$ $λ^{-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". △ Less

Submitted 3 March, 2022; originally announced March 2022.

Comments: 5 pages, 4 figures, accepted for publication in MNRAS Letters

The R136 star cluster dissected with Hubble Space Telescope/STIS. III. The most massive stars and their clumped winds

Authors: Sarah A. Brands, Alex de Koter, Joachim M. Bestenlehner, Paul A. Crowther, Jon O. Sundqvist, Joachim Puls, Saida M. Caballero-Nieves, Michael Abdul-Masih, Florian A. Driessen, Miriam García, Sam Geen, Götz Gräfener, Calum Hawcroft, Lex Kaper, Zsolt Keszthelyi, Norbert Langer, Hugues Sana, Fabian R. N. Schneider, Tomer Shenar, Jorick S. Vink

Abstract: 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 propert… ▽ More 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čka & Kubát (2018) matches best the observed mass-loss rates. △ Less

Submitted 22 February, 2022; originally announced February 2022.

Comments: Accepted for publication in A&A; Appendix I will not be included in the published version

Journal ref: A&A 663, A36 (2022)

Grids of stellar models with rotation VII: Models from 0.8 to 300 M$_\odot$ at super-solar metallicity (Z = 0.020)

Authors: Norhasliza Yusof, Raphael Hirschi, Patrick Eggenberger, Sylvia Ekström, Cyril Georgy, Yves Sibony, Paul A. Crowther, Georges Meynet, Hasan Abu Kassim, Wan Aishah Wan Harun, André Maeder, Jose H. Groh, Eoin Farrell, Laura Murphy

Abstract: 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 ma… ▽ More 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. △ Less

Submitted 21 January, 2022; originally announced January 2022.

Comments: 15 pages, 12 figures. Accepted for publication in MNRAS

Melnick 33Na: a very massive colliding wind binary system in 30 Doradus

Authors: Joachim M. Bestenlehner, Paul A. Crowther, Patrick S. Broos, Andrew M. T. Pollock, Leisa K. Townsley

Abstract: We present spectroscopic analysis of the luminous X-ray source Melnick 33Na (Mk 33Na, HSH95 16) in the LMC 30 Doradus region (Tarantula Nebula), utilising new time-series VLT/UVES spectroscopy. We confirm Mk 33Na as a double-lined O-type spectroscopic binary with a mass ratio $q = 0.63 \pm 0.02$, $e = 0.33 \pm 0.01$ and orbital period of $18.3 \pm 0.1$ days, supporting the favoured period from X-r… ▽ More We present spectroscopic analysis of the luminous X-ray source Melnick 33Na (Mk 33Na, HSH95 16) in the LMC 30 Doradus region (Tarantula Nebula), utilising new time-series VLT/UVES spectroscopy. We confirm Mk 33Na as a double-lined O-type spectroscopic binary with a mass ratio $q = 0.63 \pm 0.02$, $e = 0.33 \pm 0.01$ and orbital period of $18.3 \pm 0.1$ days, supporting the favoured period from X-ray observations obtained via the Tarantula -- Revealed by X-rays (T-ReX) survey. Disentangled spectra of each component provide spectral types of OC2.5 If* and O4 V for the primary and secondary respectively - unusually for an O supergiant the primary exhibits strong CIV 4658 emission and weak NV 4603-20, justifying the OC classification. Spectroscopic analysis favours extreme physical properties for the primary ($T_{\rm eff} = 50$ kK, $\log L/L_{\odot} = 6.15$) with system components of $M_{1} = 83 \pm 19 M_{\odot}$ and $M_{2} = 48 \pm 11 M_{\odot}$ obtained from evolutionary models, which can be reconciled with results from our orbital analysis (e.g. $M_{1} \sin^3 i = 20.0 \pm 1.2 M_{\odot}$) if the system inclination is $\sim 38^{\circ}$ and it has an age of 0.9 to 1.6 Myr. This establishes Mk 33Na as one of the highest mass binary systems in the LMC, alongside other X-ray luminous early-type binaries Mk34 (WN5h+WN5h), R144 (WN5/6h+WN6/7h) and especially R139 (O6.5\,Iafc+O6\,Iaf). △ Less

Submitted 30 November, 2021; originally announced December 2021.

Comments: Accepted for publication in MNRAS, 16 pages and 14 figures plus supplementary material

An environmental analysis of the Type Ib SN 2019yvr and the possible presence of an inflated binary companion

Authors: Ning-Chen Sun, Justyn R. Maund, Paul A. Crowther, Ryosuke Hirai, Amir Kashapov, Ji-Feng Liu, Liang-Duan Liu, Emmanouil Zapartas

Abstract: 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 epis… ▽ More 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. △ Less

Submitted 23 December, 2021; v1 submitted 11 November, 2021; originally announced November 2021.

Comments: 15 pages, 9 figures, accepted for publication in MNRAS

The Uncertain Future of Massive Binaries Obscures the Origin of LIGO/Virgo Sources

Authors: K. Belczynski, A. Romagnolo, A. Olejak, J. Klencki, D. Chattopadhyay, S. Stevenson, M. Coleman Miller, J. -P. Lasota, P. A. Crowther

Abstract: 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.… ▽ More 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. △ Less

Submitted 24 August, 2021; originally announced August 2021.

Comments: comments welcomed

The Tarantula Massive Binary Monitoring V. R 144: a wind-eclipsing binary with a total mass > 140 Msun

Authors: T. Shenar, H. Sana, P. Marchant, B. Pablo, N. Richardson, A. F. J. Moffat, T. Van Reeth, R. H. Barba, D. M. Bowman, P. Broos, P. A. Crowther, J. S. Clark, A. de Koter, S. E. de Mink, K. Dsilva, G. Graefener, I. D. Howarth, N. Langer, L. Mahy, J. Maiz Apellaniz, A. M. T. Pollock, F. R. N. Schneider, L. Townsley, J. S. Vink

Abstract: R 144 is the visually brightest WR star in the Large Magellanic Cloud (LMC). R 144 was reported to be a binary, making it potentially the most massive binary thus observed. We perform a comprehensive spectral, photometric, orbital, and polarimetric analysis of R 144. R 144 is an eccentric (e=0.51) 74.2-d binary comprising two relatively evolved (age~2 Myr), H-rich WR stars. The hotter primary (W… ▽ More R 144 is the visually brightest WR star in the Large Magellanic Cloud (LMC). R 144 was reported to be a binary, making it potentially the most massive binary thus observed. We perform a comprehensive spectral, photometric, orbital, and polarimetric analysis of R 144. R 144 is an eccentric (e=0.51) 74.2-d binary comprising two relatively evolved (age~2 Myr), H-rich WR stars. The hotter primary (WN5/6h, T=50 kK) and the cooler secondary (WN6/7h,T=45kK) have nearly equal masses. The combination of low rotation and H-depletion observed in the system is well reproduced by contemporary evolution models that include boosted mass-loss at the upper-mass end. The systemic velocity of R 144 and its relative isolation suggest that it was ejected as a runaway from the neighbouring R 136 cluster. The optical light-curve shows a clear orbital modulation that can be well explained as a combination of two processes: excess emission stemming from wind-wind collisions and double wind eclipses. Our light-curve model implies an orbital inclination of i=60.4+-1.5deg, resulting in accurately constrained dynamical masses of 74+-4 and 69+-4 Msun. Assuming that both binary components are core H-burning, these masses are difficult to reconcile with the derived luminosities (logL1,2 = 6.44, 6.39 [Lsun]), which correspond to evolutionary masses of the order of 110 and 100Msun, respectively. Taken at face value, our results imply that both stars have high classical Eddington factors of Gamma_e = 0.78+-0.1. If the stars are on the main sequence, their derived radii (~25Rsun) suggest that they are only slightly inflated, even at this high Eddington factor. Alternatively, the stars could be core-He burning, strongly inflated from the regular size of classical Wolf-Rayet stars (~1Rsun), a scenario that could help resolve the observed mass discrepancy. △ Less

Submitted 28 October, 2021; v1 submitted 7 April, 2021; originally announced April 2021.

Comments: Accepted to A&A on the April 06 2021. 17 pages + 5 page appendix. Version 28.10.2021: corrected typo in Eq. (2); numerator and denominator were inversed on RHS (only typo; no impact on results)

Journal ref: A&A 650, A147 (2021)

Mapping the core of the Tarantula Nebula with VLT-MUSE II. The spectroscopic Hertzsprung-Russell diagram of OB stars in NGC 2070

Authors: N. Castro, P. A. Crowther, C. J. Evans, J. S. Vink, J. Puls, A. Herrero, M. Garcia, F. J. Selman, M. M. Roth, S. Simón-Díaz

Abstract: 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 s… ▽ More 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$λ$4921/HeII$λ$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. △ Less

Submitted 5 February, 2021; originally announced February 2021.

Comments: 19 pages, 12 figures. Accepted in A&A

Journal ref: A&A 648, A65 (2021)

A dearth of young and bright massive stars in the Small Magellanic Cloud

Authors: A. Schootemeijer, N. Langer, D. Lennon, C. J. Evans, P. A. Crowther, S. Geen, I. Howarth, A. de Koter, K. M. Menten, J. S. Vink

Abstract: 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 lu… ▽ More 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. △ Less

Submitted 10 December, 2020; originally announced December 2020.

Comments: 26 pages, 23 figures. Astronomy and Astrophysics, in press

Journal ref: A&A 646, A106 (2021)

Toward a better understanding of supernova environments: a study of SNe 2004dg and 2012P in NGC 5806 with HST and MUSE

Authors: Ning-Chen Sun, Justyn R. Maund, Paul A. Crowther, Xuan Fang, Emmanouil Zapartas

Abstract: 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 i… ▽ More 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. △ Less

Submitted 7 April, 2021; v1 submitted 27 November, 2020; originally announced November 2020.

Comments: 21 pages, 17 figures, accepted for publication in MNRAS

2D kinematics of massive stars near the Galactic Center

Authors: Mattia Libralato, Daniel J. Lennon, Andrea Bellini, Roeland van der Marel, Simon J. Clark, Francisco Najarro, Lee R. Patrick, Jay Anderson, Luigi R. Bedin, Paul A. Crowther, Selma E. de Mink, Christopher J. Evans, Imants Platais, Elena Sabbi, Sangmo Tony Sohn

Abstract: 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… ▽ More 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. △ Less

Submitted 13 November, 2020; v1 submitted 21 October, 2020; originally announced October 2020.

Comments: 28 pages, 33 figures, 9 tables. Accepted for publication in MNRAS. [v3: Fixed bibliography]

The R136 star cluster dissected with Hubble Space Telescope/STIS. II. Physical properties of the most massive stars in R136

Authors: Joachim M. Bestenlehner, Paul A. Crowther, Saida M. Caballero-Nieves, Fabian R. N. Schneider, Sergio Simon-Diaz, Sarah A. Brands, Alex de Koter, Goetz Graefener, Artemio Herrero, Norbert Langer, Daniel J. Lennon, Jesus Maiz Apellaniz, Joachim Puls, Jorick S. Vink

Abstract: 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 evolution… ▽ More 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 $γ\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. △ Less

Submitted 10 September, 2020; originally announced September 2020.

Comments: Accepted for publication in MNRAS, 21 pages, 14 figures, 2 tables plus 33 pages of supplementary material

The extreme colliding-wind system Apep: resolved imagery of the central binary and dust plume in the infrared

Authors: Y. Han, P. G. Tuthill, R. M. Lau, A. Soulain, J. R. Callingham, P. M. Williams, P. A. Crowther, B. J. S. Pope, B. Marcote

Abstract: 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… ▽ More 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. △ Less

Submitted 13 August, 2020; originally announced August 2020.

Comments: This article has been accepted for publication in the Monthly Notices of the Royal Astronomical Society. 17 pages, 9 figures, 6 tables

Unlocking Galactic Wolf-Rayet stars with $\textit{Gaia}$ DR2 II: Cluster and Association membership

Authors: Gemma Rate, Paul A. Crowther, Richard J. Parker

Abstract: 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 exter… ▽ More 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. △ Less

Submitted 5 May, 2020; originally announced May 2020.

Comments: 19 pages, 7 figures. Accepted for publication in MNRAS

Two Wolf-Rayet stars at the heart of colliding-wind binary Apep

Authors: J. R. Callingham, P. A. Crowther, P. M. Williams, P. G. Tuthill, Y. Han, B. J. S. Pope, B. Marcote

Abstract: 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 diff… ▽ More 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. △ Less

Submitted 1 May, 2020; originally announced May 2020.

Comments: Accepted for publication in MNRAS, 10 pages, 9 figures, 1 table

The changing-type SN 2014C may come from an 11-$M_\odot$ star stripped by binary interaction and violent eruption

Authors: Ning-Chen Sun, Justyn R. Maund, Paul A. Crowther

Abstract: SN 2014C was an unprecedented supernova (SN) that displayed a metamorphosis from Type Ib to Type IIn over $\sim$200 days. This transformation is consistent with a helium star having exploded in a cavity surrounded by a dense shell of the progenitor's stripped hydrogen envelope. For at least 5 years post-explosion, the ejecta continued to interact with an outer, extended component of circumstellar… ▽ More SN 2014C was an unprecedented supernova (SN) that displayed a metamorphosis from Type Ib to Type IIn over $\sim$200 days. This transformation is consistent with a helium star having exploded in a cavity surrounded by a dense shell of the progenitor's stripped hydrogen envelope. For at least 5 years post-explosion, the ejecta continued to interact with an outer, extended component of circumstellar medium (CSM) that was ejected even before the dense shell. It is still unclear, however, what kind of progenitor could have undergone such a complicated mass-loss history before it produced this peculiar SN. In this paper, we report a new analysis of SN 2014C's host star cluster based on data from the Hubble Space Telescope (HST). By carefully fitting its spectral energy distribution (SED), we derive a precise cluster age of 20.0$^{+3.5}_{-2.6}$ Myr, which corresponds to the progenitor's lifetime assuming coevolution. Combined with binary stellar evolution models, we find that SN 2014C's progenitor may have been an $\sim$11-$M_\odot$ star in a relatively wide binary system. The progenitor's envelope was partially stripped by Case C or Case BC mass transfer via binary interaction, followed by a violent eruption that ejected the last hydrogen layer before terminal explosion. Thus, SN 2014C, in common with SNe 2006jc and 2015G, may be a third example that violent eruptions, with mass-loss rates matching luminous blue variable (LBV) giant eruptions, can also occur in much lower-mass massive stars if their envelopes are partially or completely stripped in interacting binaries. △ Less

Submitted 29 July, 2020; v1 submitted 20 March, 2020; originally announced March 2020.

Comments: 19 pages, 13 figures, minor revision, accepted for publication in MNRAS

A search for strong magnetic fields in massive and very massive stars in the Magellanic Clouds

Authors: S. Bagnulo, G. A. Wade, Y. Naze, J. H. Grunhut, M. E. Shultz, D. J. Asher, P. A. Crowther, C. J. Evans, A. David-Uraz, I. D. Howarth, N. Morrell, M. S. Munoz, C. Neiner, J. Puls, M. K. Szymanski, J. S. Vink

Abstract: 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, o… ▽ More 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. △ Less

Submitted 27 February, 2020; originally announced February 2020.

Comments: Accepted by A&A

Journal ref: A&A, 635, A163 (15p, 2020)

Unlocking Galactic Wolf-Rayet stars with $\textit{Gaia}$ DR2 I: Distances and absolute magnitudes

Authors: Gemma Rate, Paul A. Crowther

Abstract: 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 spectro… ▽ More 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$σ$, |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. △ Less

Submitted 11 February, 2020; v1 submitted 20 December, 2019; originally announced December 2019.

Comments: 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

Properties of OB star-black hole systems derived from detailed binary evolution models

Authors: N. Langer, C. Schürmann, K. Stoll, P. Marchant, D. J. Lennon, L. Mahy, S. E. de Mink, M. Quast, W. Riedel, H. Sana, P. Schneider, A. Schootemeijer, Chen Wang, L. A. Almeida, J. M. Bestenlehner, J. Bodensteiner, N. Castro, S. Clark, P. A. Crowther, P. Dufton, C. J. Evans, L. Fossati, G. Gräfener, L. Grassitelli, N. Grin , et al. (16 additional authors not shown)

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

Submitted 10 April, 2020; v1 submitted 20 December, 2019; originally announced December 2019.

Comments: 19 pages, 14 figures, Astronomy and Astrophysics, in press

Journal ref: A&A 638, A39 (2020)

Origins of Type Ibn SNe 2006jc/2015G in interacting binaries and implications for pre-SN eruptions

Authors: Ning-Chen Sun, Jusytn R. Maund, Ryosuke Hirai, Paul A. Crowther, Philipp Podsiadlowski

Abstract: 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 ultraviole… ▽ More 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. △ Less

Submitted 4 December, 2019; v1 submitted 17 September, 2019; originally announced September 2019.

Comments: 22 pages, 18 figures, accepted for publication in MNRAS

Investigating the origin of the spectral line profiles of the Hot Wolf-Rayet Star WR2

Authors: A. -N. Chené, N. St-Louis, A. F. J. Moffat, O. Schnurr, P. A. Crowther, G. A. Wade, N. D. Richardson, C. Baranec, C. A. Ziegler, N. M. Law, R. Riddle, G. A. Rate, É. Artigau, E. Alecian, BinaMIcS collaboration

Abstract: 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… ▽ More 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. △ Less

Submitted 14 May, 2019; originally announced May 2019.

Comments: 12 pages, 9 figures

Journal ref: MNRAS, 484, 5834 (2019)

Weighing Melnick 34: the most massive binary system known

Authors: Katie A. Tehrani, Paul A. Crowther, Joachim M. Bestenlehner, Stuart P. Littlefair, A. M. T. Pollock, Richard J. Parker, Olivier Schnurr

Abstract: 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 cl… ▽ More 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°, 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. △ Less

Submitted 15 January, 2019; originally announced January 2019.

Comments: 21 pages, 18 figures, accepted for publication in MNRAS

Investigating the properties of stripped-envelope supernovae, what are the implications for their progenitors?

Authors: S. J. Prentice, C. Ashall, P. A. James, L. Short, P. A. Mazzali, D. Bersier, P. A. Crowther, C. Barbarino, T. -W. Chen, C. M. Copperwheat, M. J. Darnley, L. Denneau, N. Elias-Rosa, M. Fraser, L. Galbany, A. Gal-Yam, J. Harmanen, D. A. Howell, G. Hosseinzadeh, C. Inserra, E. Kankare, E. Karamehmetoglu, G. P. Lamb, M. Limongi, K. Maguire , et al. (19 additional authors not shown)

Abstract: 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). Additionall… ▽ More 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 $<20$ Msun. The bolometric light curve properties are examined in combination with those of other SE events from the literature. The resulting dataset gives the Mej distribution for 80 SE-SNe, the largest such sample in the literature to date, and shows that SNe Ib have the lowest median Mej, followed by narrow lined SNe Ic, H/He-rich SNe, broad lined SNe Ic, and finally gamma-ray burst SNe. SNe Ic-6/7 show the largest spread of Mej, ranging from $\sim 1.2 - 11$ Msun, considerably greater than any other subtype. For all SE-SNe $<$Mej$>=2.8\pm{1.5}$ Msun which further strengthens the evidence that SE-SNe arise from low mass progenitors which are typically $<5$ Msun at the time of explosion, again suggesting Mzams $<25$ Msun. The low $<$Mej$>$ and lack of clear bimodality in the distribution implies $<30$ Msun progenitors and that envelope stripping via binary interaction is the dominant evolutionary pathway of these SNe. △ Less

Submitted 13 December, 2018; v1 submitted 10 December, 2018; originally announced December 2018.

Comments: Accepted for publication in MNRAS, no changes to the previous submission

Anisotropic winds in Wolf-Rayet binary identify potential gamma-ray burst progenitor

Authors: J. R. Callingham, P. G. Tuthill, B. J. S. Pope, P. M. Williams, P. A. Crowther, M. Edwards, B. Norris, L. Kedziora-Chudczer

Abstract: 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… ▽ More 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. △ Less

Submitted 16 November, 2018; originally announced November 2018.

Comments: Accepted to Nature Astronomy on 2018 Oct 1. 27 pages, 9 figures, and 4 tables

Response to comment on "An excess of massive stars in the local 30 Doradus starburst"

Authors: F. R. N. Schneider, H. Sana, C. J. Evans, J. M. Bestenlehner, N. Castro, L. Fossati, G. Gräfener, N. Langer, O. H. Ramírez-Agudelo, C. Sabín-Sanjulián, S. Simón-Díaz, F. Tramper, P. A. Crowther, A. de Koter, S. E. de Mink, P. L. Dufton, M. Garcia, M. Gieles, V. Hénault-Brunet, A. Herrero, R. G. Izzard, V. Kalari, D. J. Lennon, J. Maíz Apellániz, N. Markova , et al. (7 additional authors not shown)

Abstract: 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 Doradu… ▽ More 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. △ Less

Submitted 25 July, 2018; originally announced July 2018.

Comments: Authors' version of a response to a technical comment published in Science; 8 pages, 1 figure

The VLT-FLAMES Tarantula Survey. XXIX. Massive star formation in the local 30 Doradus starburst

Authors: F. R. N. Schneider, O. H. Ramírez-Agudelo, F. Tramper, J. M. Bestenlehner, N. Castro, H. Sana, C. J. Evans, C. Sabín-Sanjulián, S. Simón-Díaz, N. Langer, L. Fossati, G. Gräfener, P. A. Crowther, S. E. de Mink, A. de Koter, M. Gieles, A. Herrero, R. G. Izzard, V. Kalari, R. S. Klessen, D. J. Lennon, L. Mahy, J. Maíz Apellániz, N. Markova, J. Th. van Loon , et al. (2 additional authors not shown)

Abstract: 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 a… ▽ More 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] △ Less

Submitted 10 July, 2018; originally announced July 2018.

Comments: 20 pages (incl. appendix), 15 figures, 2 tables; accepted for publication in A&A

Journal ref: A&A 618, A73 (2018)

Probing the rotational velocity of Galactic WO stars with spectropolarimetry

Authors: H. F. Stevance, R. Ignace, P. A. Crowther, J. R. Maund, B. Davies, G. Rate

Abstract: Oxygen sequence Wolf-Rayet stars (WO) are thought to be the final evolution phase of some high mass stars, as such they may be the progenitors of type Ic SNe as well as potential progenitors of broad-lined Ic and long gamma-ray bursts. We present the first spectropolarimetric observations of the Galactic WO stars WR93b and WR102 obtained with FORS1 on the VLT. We find no sign of a line effect, whi… ▽ More Oxygen sequence Wolf-Rayet stars (WO) are thought to be the final evolution phase of some high mass stars, as such they may be the progenitors of type Ic SNe as well as potential progenitors of broad-lined Ic and long gamma-ray bursts. We present the first spectropolarimetric observations of the Galactic WO stars WR93b and WR102 obtained with FORS1 on the VLT. We find no sign of a line effect, which could be expected if these stars were rapid rotators. We also place constraints on the amplitude of a potentially undetected line effect. This allows us to derive upper limits on the possible intrinsic continuum polarisation, and find P$_{\rm cont}$ < 0.077 percent and P$_{\rm cont}$ < 0.057 percent for WR93b and WR102, respectively. Furthermore, we derive upper limits on the rotation of our WO stars by considering our results in the context of the wind compression effect. We estimate that for an edge-on case the rotational velocity of WR93b is v$_{\rm rot}$ < 324 km/s while for WR102 v$_{\rm rot}$ < 234 km/s. These correspond to values of v$_{\rm rot}$/v$_{\rm crit}$ <19 percent and <10 percent, respectively, and values of log(j)<18.0 cm$^2$/s for WR93b and <17.6 cm^2 /s for WR102. The upper limits found on v$_{\rm rot}$/v$_{\rm crit}$ and log(j) for our WO stars are therefore similar to the estimates calculated for Galactic WR stars that do show a line effect. Therefore, although the presence of a line effect in single WR stars is indicative of fast rotation, the absence of a line effect does not rule out significant rotation, even when considering the edge-on scenario. △ Less

Submitted 10 July, 2018; v1 submitted 5 July, 2018; originally announced July 2018.

Comments: Accepted for Publication in MNRAS (9 pages, 3 figures)

Gaia DR2 reveals a very massive runaway star ejected from R136

Authors: D. J. Lennon, C. J. Evans, R. P. van der Marel, J. Anderson, I. Platais, A. Herrero, S. E. de Mink, H. Sana, E. Sabbi, L. R. Bedin, P. A. Crowther, N. Langer, M. Ramos Lerate, A. del Pino, M. Renzo, S. Simón-Díaz, F. R. N. Schneider

Abstract: 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… ▽ More 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$σ$ 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. △ Less

Submitted 24 July, 2018; v1 submitted 21 May, 2018; originally announced May 2018.

Comments: Version 2 is now accepted by A&A with minor updates from V1

Journal ref: A&A 619, A78 (2018)

The luminosities of cool supergiants in the Magellanic Clouds, and the Humphreys-Davidson limit revisited

Authors: Ben Davies, Paul A. Crowther, Emma R. Beasor

Abstract: 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 T… ▽ More 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. △ Less

Submitted 17 April, 2018; originally announced April 2018.

Comments: 11 pages, 6 figs. Accepted for publication in MNRAS

The Arches cluster revisited: II. A massive eclipsing spectroscopic binary in the Arches cluster

Authors: M. E. Lohr, J. S. Clark, F. Najarro, L. R. Patrick, P. A. Crowther, C. J. Evans

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 an… ▽ More 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. △ Less

Submitted 16 April, 2018; originally announced April 2018.

Comments: 15 pages, 18 figures, accepted for publication in Astronomy & Astrophysics. Paper II in a series; for Paper I see arXiv:1803.09567

Journal ref: A&A 617, A66 (2018)

The VLT-FLAMES Tarantula Survey: XXVIII. Nitrogen abundances for apparently single dwarf and giant B-type stars with small projected rotational velocities

Authors: P. L. Dufton, A. Thompson, P. A. Crowther, C. J. Evans, F. R. N. Schneider, A. de Koter, S. E. de Mink, R. Garland, N. Langer, D. J. Lennon, C. M. McEvoy, O. H. Ramírez-Agudelo, H. Sana, S. Símon Díaz, W. D. Taylor, J. S. Vink

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

Submitted 5 April, 2018; originally announced April 2018.

Comments: Accepted in Astronomy and Astrophyics. Abstract significantly shortened - see paper for full abstract

Journal ref: A&A 615, A101 (2018)

The 155-day X-ray cycle of the very massive Wolf-Rayet star Melnick 34 in the Large Magellanic Cloud

Authors: A. M. T. Pollock, P. A. Crowther, K. Tehrani, Patrick S. Broos, Leisa K. Townsley

Abstract: 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 eccentr… ▽ More 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. △ Less

Submitted 2 March, 2018; originally announced March 2018.

Comments: 10 pages, 8 figures, 5 tables