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

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 VIII. Stellar and wind parameters of newly revealed stripped stars in Be binaries

Authors: V. Ramachandran, A. A. C. Sander, D. Pauli, J. Klencki, F. Backs, F. Tramper, M. Bernini-Peron, P. Crowther, W. -R. Hamann, R. Ignace, R. Kuiper, S. Oey, L. M. Oskinova, T. Shenar, H. Todt, J. S. Vink, L. Wang, A. Wofford, the XShootU collaboration

Abstract: On the route towards merging neutron stars and stripped-envelope supernovae, binary population synthesis predicts a large number of post-interaction systems with massive stars that have stripped off their outer layers. Yet, observations of such stars in the intermediate-mass regime below the Wolf-Rayet masses are rare. Using X-Shooting ULLYSES (XShootU) data, we discovered three partially stripped… ▽ More On the route towards merging neutron stars and stripped-envelope supernovae, binary population synthesis predicts a large number of post-interaction systems with massive stars that have stripped off their outer layers. Yet, observations of such stars in the intermediate-mass regime below the Wolf-Rayet masses are rare. Using X-Shooting ULLYSES (XShootU) data, we discovered three partially stripped star + Be/Oe binaries in the Magellanic Clouds. We analyzed the UV and optical spectra using the PoWR model atmosphere code by superimposing model spectra corresponding to each component. The estimated current masses of the partially stripped stars fall within the intermediate mass range of 4-8 $M_{\odot}$. These objects are overluminous for their stellar masses, matching core He-burning luminosities. Their Be/Oe secondaries have much higher masses than their stripped primaries (mass ratio > 2). All three partially stripped stars show significant nitrogen enrichment and carbon and oxygen depletion on their surfaces. Additionally, one of our sample stars exhibits significant helium enrichment. Our study provides the first comprehensive determination of the wind parameters of partially stripped stars in the intermediate mass range. The wind mass-loss rates of these stars are found to be on the order of $10^{-7} M_\odot$ yr$^{-1}$, which is over ten times higher than that of OB stars of the same luminosity. Current evolutionary models characterizing this phase typically employ OB or WR mass-loss rates, which underestimate or overestimate stripped stars' mass-loss rates by an order of magnitude. Binary evolution models indicate that the observed primaries had initial masses of 12-17 $M_{\odot}$, making them potential candidates for stripped-envelope supernovae that form neutron stars. If they survive the explosion, these systems may become Be X-ray binaries and later double neutron stars. △ Less

Submitted 25 June, 2024; originally announced June 2024.

Comments: Accepted for publication in Astronomy & Astrophysics

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

X-Shooting ULLYSES: Massive Stars at Low Metallicity

Authors: Jorick S. Vink, Paul Crowther, Alex Fullerton, Miriam Garcia, Fabrice Martins, Nidia Morrell, Lida Oskinova, Nicole St. Louis, Asif ud-Doula, Andreas Sander, Hugues Sana, Jean-Claude Bouret, Brankica Kubatova, Pablo Marchant, Lucimara P. Martins, Aida Wofford, Jacco van Loon, O. Grace Telford, Ylva Götberg, Dominic Bowman, Christi Erba, Venu Kalari, The XShootU Collaboration

Abstract: The Hubble Space Telescope has devoted 500 orbits to observing 250 massive stars with low metallicity in the ultraviolet (UV) range within the framework of the ULLYSES program. The X-Shooting ULLYSES (XShootU) project enhances the legacy value of this UV dataset by providing high-quality optical and near-infrared spectra, which are acquired using the wide-wavelength-coverage X-shooter spectrograph… ▽ More The Hubble Space Telescope has devoted 500 orbits to observing 250 massive stars with low metallicity in the ultraviolet (UV) range within the framework of the ULLYSES program. The X-Shooting ULLYSES (XShootU) project enhances the legacy value of this UV dataset by providing high-quality optical and near-infrared spectra, which are acquired using the wide-wavelength-coverage X-shooter spectrograph at ESO's Very Large Telescope. XShootU emphasises the importance of combining UV with optical spectra for the consistent determination of key stellar parameters such as effective temperature, surface gravity, luminosity, abundances, and wind characteristics including mass-loss rates as a function of metallicity. Since uncertainties in these parameters have implications across various branches of astrophysics, the data and modelling generated by the XShootU project are poised to significantly advance our understanding of massive stars at low metallicity. This is particularly crucial for confidently interpreting JWST data of the earliest stellar generations, making XShootU a unique resource for comprehending individual spectra of low-metallicity stars. △ Less

Submitted 30 April, 2024; originally announced May 2024.

Comments: 6 pages, 6 figures. ESO Large Programme Overview

Journal ref: ESO Messenger, 2024

X-Shooting ULLYSES: Massive Stars at low metallicity II. DR1: Advanced optical data products for the Magellanic Clouds

Authors: H. Sana, F. Tramper, M. Abdul-Masih, R. Blomme, K. Dsilva, G. Maravelias, L. Martins, A. Mehner, A. Wofford, G. Banyard, C. L. Barbosa, J. Bestenlehner, C. Hawcroft, D. John Hillier, H. Todt, C. J. K. Larkin, L. Mahy, F. Najarro, V. Ramachandran, M. C. Ramirez-Tannus, M. M. Rubio-Diez, A. A. C. Sander, T. Shenar, J. S. Vink, F. Backs , et al. (12 additional authors not shown)

Abstract: Using the medium resolution spectrograph X-shooter, spectra of 235 OB and Wolf-Rayet (WR) stars in sub-solar metallicity environments have been secured. [...]This second paper focuses on the optical observations of 232 Magellanic Clouds targets. It describes the uniform reduction of the UVB (300 - 560 nm) and VIS (550 - 1020 nm) XShootU data as well as the preparation of advanced data products [..… ▽ More Using the medium resolution spectrograph X-shooter, spectra of 235 OB and Wolf-Rayet (WR) stars in sub-solar metallicity environments have been secured. [...]This second paper focuses on the optical observations of 232 Magellanic Clouds targets. It describes the uniform reduction of the UVB (300 - 560 nm) and VIS (550 - 1020 nm) XShootU data as well as the preparation of advanced data products [...] . The data reduction of the raw data is based on the ESO CPL X-shooter pipeline. We paid particular attention to the determination of the response curves [...] We implemented slit-loss correction, absolute flux calibration, (semi-)automatic rectification to the continuum, and a correction for telluric lines. The spectra of individual epochs were corrected for the barycentric motion, re-sampled and co-added, and the spectra from the two arms were merged into a single flux calibrated spectrum covering the entire optical range with maximum signal-to-noise ratio. [...] We provide three types of data products: (i) two-dimensional spectra for each UVB and VIS exposure; (ii) one-dimensional UVB and VIS spectra before and after response-correction, as well as after applying various processing, including absolute flux calibration, telluric removal, normalisation and barycentric correction; and (iii) co-added flux-calibrated and rectified spectra over the full optical range, for which all available XShootU exposures were combined. For many of the targets, the final signal-to-noise ratio per resolution element is above 200 in both the UVB and the VIS co-added spectra. The reduced data and advanced scientific data products will be made available to the community upon publication of this paper. [...] △ Less

Submitted 26 February, 2024; originally announced February 2024.

Comments: 22 pages, 19 figures; accepted for publication in Astronomy & Astrophysics. Links to online tables and databases will be included upon publication

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

A resolved study of the inner regions of nearby galaxies with an excess of young massive stars: missing link in the AGN-starburst connection?

Authors: Guinevere Kauffmann, Iker Millán-Irigoyen, Paul Crowther, Claudia Maraston

Abstract: We have selected galaxies with very high levels of H$α$ emission (EQW(H$α$)$>$700 Å.) in their central regions from the final data release of the MaNGA survey . Our study focuses on 14 very well-resolved nearby galaxies with stellar masses in the range $9.5 < \log M_*/(M_{\odot}) < 11.5$. We investigate a variety of procedures for selecting galaxy regions that are likely to harbour excess populati… ▽ More We have selected galaxies with very high levels of H$α$ emission (EQW(H$α$)$>$700 Å.) in their central regions from the final data release of the MaNGA survey . Our study focuses on 14 very well-resolved nearby galaxies with stellar masses in the range $9.5 < \log M_*/(M_{\odot}) < 11.5$. We investigate a variety of procedures for selecting galaxy regions that are likely to harbour excess populations of young massive stars, finding that selection in the 2-dimensional space of extinction-corrected H$α$ EQW and [SIII]/[SII] line ratio produces the best results. By comparing stacked spectra covering these regions with stacked spectra covering normal starburst regions with 100Å$<$EQW(H$α$)$<$200Å, we obtain the following main results: 1) Clear signatures of excess Wolf Rayet stars are found in half of the H$α$ excess regions, 2) Galaxy regions containing excess Wolf Rayet stars are more often associated with the presence of high-ionization emission lines characteristic of accreting black holes. Excess [NeIII] is detected in 4 out of 8 of the WR regions and there are tentative [FeX] detections in 2 galaxies. 3) Regions of the galaxy with excess Wolf Rayet stars are located where the interstellar medium has larger ionized gas turbulent velocities and higher neutral gas overdensities. We make a first attempt to constrain changes in the high mass end of the stellar initial mass function (IMF) using the HR-pyPopStar evolutionary population synthesis models that include high wavelength-resolution theoretical atmosphere libraries for Wolf Rayet stars. △ Less

Submitted 9 October, 2023; originally announced October 2023.

Comments: accepted for publication in MNRAS

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 study of 1000 galaxies with unusually young and massive stars in the SDSS: a search for hidden black holes

Authors: Guinevere Kauffmann, Claudia Maraston, Johan Comparat, Paul Crowther

Abstract: We select 1076 galaxies with extinction-corrected Halpha equivalent widths too large to be explained with a Kroupa (2001) IMF, and compare these with a control sample of galaxies that is matched in stellar mass, redshift and 4000 AA break strength, but with normal Halpha equivalent widths. Our goal is to study how processes such as black hole growth and energetic feedback processes from massive st… ▽ More We select 1076 galaxies with extinction-corrected Halpha equivalent widths too large to be explained with a Kroupa (2001) IMF, and compare these with a control sample of galaxies that is matched in stellar mass, redshift and 4000 AA break strength, but with normal Halpha equivalent widths. Our goal is to study how processes such as black hole growth and energetic feedback processes from massive stars differ between galaxies with extreme central Halpha emission and galaxies with normal young central stellar populations. The stellar mass distribution of Halpha excess galaxies is peaked at 3 \times 10^10 Msun and almost all fall well within the star-forming locus in the [OIII]/Hbeta versus [NII]/Halpha BPT disgram. Halpha excess galaxies are twice as likely to exhibit Halpha line asymmetries and 1.55 times more likely to be detected at 1 GHz in the VLA FIRST survey compared to control sample galaxies. The radio luminosity per unit stellar mass decreases with the stellar age of the system. Using stacked spectra, we demonstrate that [NeV] emission is not present in the very youngest of the radio-quiet Halpha excess galaxies with detectable Wolf-Rayet features, suggesting that black hole growth has not yet commenced in such systems. [NeV] emission is detected in Halpha excess galaxies with radio detections and the strength of the line correlates with the radio luminosity. This is the clearest indication for a possible population of black holes that may be forming in a subset of the Halpha excess population. △ Less

Submitted 8 April, 2022; originally announced April 2022.

Comments: 17 pages, 22 figures, accepted in MNRAS

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

Extreme adaptive optics astrometry of R136. Searching for high proper motion stars

Authors: Zeinab Khorrami, M. Langlois, F. Vakili, P. C. Clark, A. S. M. Buckner, M. Gonzalez, P. Crowther, R. Wunsch, J. Palous, A. Boccaletti, S. Lumsden, E. Moraux

Abstract: We compared high-contrast near-infrared images of the core of R136 taken by VLT/SPHERE, in two epochs separated by 3.06 years. For the first time we monitored the dynamics of the detected sources in the core of R136 from a ground-based telescope with adaptive optics. The aim of these observations was to search for High prOper Motion cAndidates (HOMAs) in the central region of R136 (r<6") where it… ▽ More We compared high-contrast near-infrared images of the core of R136 taken by VLT/SPHERE, in two epochs separated by 3.06 years. For the first time we monitored the dynamics of the detected sources in the core of R136 from a ground-based telescope with adaptive optics. The aim of these observations was to search for High prOper Motion cAndidates (HOMAs) in the central region of R136 (r<6") where it has been challenging for other instruments. Two bright sources (K<15mag and V<16mag) are located near R136a1 and R136c (massive WR stars) and have been identified as potential HOMAs. These sources have significantly shifted in the images with respect to the mean shift of all reliable detected sources and their neighbours, and six times their own astrometric errors. We calculate their proper motions to be 1.36\pm0.22 mas/yr (321\pm52 km/s) and 1.15\pm0.11 mas/yr (273\pm26 km/s). We discuss different possible scenarios to explain the magnitude of such extreme proper motions, and argue for the necessity to conduct future observations to conclude on the nature of HOMAs in the core of R136. △ Less

Submitted 26 April, 2021; originally announced April 2021.

Comments: Accepted to be published in A&A/Letter. catalogue of the reliable-consistent sources are available online via CDS

Journal ref: A&A 649, L8 (2021)

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)

High contrast and resolution near infrared photometry of the core of R136

Authors: Zeinab Khorrami, Maud Langlois, Paul C. Clark, Farrokh Vakili, Anne S. M. Buckner, Marta Gonzalez, Paul Crowther, Richard Wunsch, Jan Palous, Stuart Lumsden, Estelle Moraux

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

Submitted 11 February, 2021; originally announced February 2021.

Comments: 24 pages, 20 figures, 6 tables. Accepted for publication in MNRAS

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)

The Tarantula Massive Binary Monitoring: III. Atmosphere analysis of double-lined spectroscopic systems

Authors: L. Mahy, H. Sana, M. Abdul-Masih, L. A. Almeida, N. Langer, T. Shenar, A. de Koter, S. E. de Mink, S. de Wit, N. J. Grin, C. J. Evans, A. F. J. Moffat, F. R. N. Schneider, R. Barbá, J. S. Clark, P. Crowther, G. Gräfener, D. J. Lennon, F. Tramper, J. S. Vink

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

Submitted 14 December, 2019; originally announced December 2019.

Comments: accepted to A&A Core of the paper 15 pages, Appendix: 41 pages

Massive stars in the Tarantula Nebula: A Rosetta Stone for Extragalactic Supergiant HII Regions

Authors: Paul Crowther

Abstract: A review of the properties of the Tarantula Nebula (30 Doradus) in the Large Magellanic Cloud is presented, primarily from the perspective of its massive star content. The proximity of the Tarantula and its accessibility to X-ray through radio observations permit it to serve as a Rosetta Stone amongst extragalactic supergiant HII regions since one can consider both its integrated characteristics a… ▽ More A review of the properties of the Tarantula Nebula (30 Doradus) in the Large Magellanic Cloud is presented, primarily from the perspective of its massive star content. The proximity of the Tarantula and its accessibility to X-ray through radio observations permit it to serve as a Rosetta Stone amongst extragalactic supergiant HII regions since one can consider both its integrated characteristics and the individual properties of individual massive stars. Recent surveys of its high mass stellar content, notably the VLT FLAMES Tarantula Survey (VFTS), are reviewed, together with VLT/MUSE observations of the central ionizing region NGC 2070 and HST/STIS spectroscopy of the young dense cluster R136, provide a near complete Hertzsprung-Russell diagram of the region, and cumulative ionizing output. Several high mass binaries are highlighted, some of which have been identified from a recent X-ray survey. Brief comparisons with the stellar content of giant HII regions in the Milky Way (NGC 3372) and Small Magellanic Cloud (NGC 346) are also made, together with Green Pea galaxies and star forming knots in high-z galaxies. Finally, the prospect of studying massive stars in metal poor galaxies is evaluated. △ Less

Submitted 8 November, 2019; v1 submitted 5 November, 2019; originally announced November 2019.

Comments: 23 pages, 12 figures, review paper accepted for "Luminous Stars in Nearby Galaxies" special issue of Galaxies (guest editor Roberta Humphreys)

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

BlueMUSE: Project Overview and Science Cases

Authors: Johan Richard, Roland Bacon, Jérémy Blaizot, Samuel Boissier, Alessandro Boselli, NicolasBouché, Jarle Brinchmann, Norberto Castro, Laure Ciesla, Paul Crowther, Emanuele Daddi, Stefan Dreizler, Pierre-Alain Duc, David Elbaz, Benoit Epinat, Chris Evans, Matteo Fossati, Michele Fumagalli, Miriam Garcia, Thibault Garel, Matthew Hayes, Angela Adamo, Artemio Herrero, Emmanuel Hugot, Andrew Humphrey , et al. (37 additional authors not shown)

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

Submitted 9 November, 2021; v1 submitted 4 June, 2019; originally announced June 2019.

Comments: 60 pages, 22 figures, minor updates

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)

The 3D shape of Type IIb SN 2011hs

Authors: H. F. Stevance, J. R. Maund, D. Baade, J. Bruten, A. Cikota, P. Höflich, L. Wang, J. C. Wheeler, A. Clocchiatti, J. Spyromilio, F. Patat, Y. Yang, P. Crowther

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

Submitted 22 January, 2019; originally announced January 2019.

Comments: Accepted for publication in MNRAS, 16 pages, 11 figures