Showing 1–50 of 97 results for author: Sander, A

New Wolf-Rayet wind yields and nucleosynthesis of Helium stars

Authors: Erin R. Higgins, Jorick S. Vink, Raphael Hirschi, Alison M. Laird, Andreas A. C. Sander

Abstract: Strong metallicity-dependent winds dominate the evolution of core He-burning, classical Wolf-Rayet (cWR) stars, which eject both H and He-fusion products such as 14N, 12C, 16O, 19F, 22Ne and 23Na during their evolution. The chemical enrichment from cWRs can be significant. cWR stars are also key sources for neutron production relevant for the weak s-process. We calculate stellar models of cWRs at… ▽ More Strong metallicity-dependent winds dominate the evolution of core He-burning, classical Wolf-Rayet (cWR) stars, which eject both H and He-fusion products such as 14N, 12C, 16O, 19F, 22Ne and 23Na during their evolution. The chemical enrichment from cWRs can be significant. cWR stars are also key sources for neutron production relevant for the weak s-process. We calculate stellar models of cWRs at solar metallicity for a range of initial Helium star masses (12-50M), adopting the recent hydrodynamical wind rates from Sander & Vink (2020). Stellar wind yields are provided for the entire post-main sequence evolution until core O-exhaustion. While literature has previously considered cWRs as a viable source of the radioisotope 26Al, we confirm that negligible 26Al is ejected by cWRs since it has decayed to 26Mg or proton-captured to 27Al. However, in Paper I, Higgins et al. (2023) we showed that very massive stars eject substantial quantities of 26Al, among other elements including N, Ne, and Na, already from the zero-age-main-sequence. Here, we examine the production of 19F and find that even with lower mass-loss rates than previous studies, our cWR models still eject substantial amounts of 19F. We provide central neutron densities (Nn) of a 30M cWR compared with a 32M post-VMS WR and confirm that during core He-burning, cWRs produce a significant number of neutrons for the weak s-process via the 22Ne(alpha,n)25Mg reaction. Finally, we compare our cWR models with observed [Ne/He], [C/He] and [O/He] ratios of Galactic WC and WO stars. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: Accepted for publication in MNRAS. 17 pages, 16 figures

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

Hydrodynamic simulation of Cygnus OB2: the absence of a cluster wind termination shock

Authors: Thibault Vieu, Cormac J. K. Larkin, Lucia Härer, Brian Reville, Andreas A. C. Sander, Varsha Ramachandran

Abstract: We perform a large-scale hydrodynamic simulation of a massive star cluster whose stellar population mimics that of the Cygnus OB2 association. The main-sequence stars are first simulated during 1.6 Myr, until a quasi-stationary state is reached. At this time the three Wolf-Rayet stars observed in Cygnus OB2 are added to the simulation, which continues to 2 Myr. Using a high-resolution grid in the… ▽ More We perform a large-scale hydrodynamic simulation of a massive star cluster whose stellar population mimics that of the Cygnus OB2 association. The main-sequence stars are first simulated during 1.6 Myr, until a quasi-stationary state is reached. At this time the three Wolf-Rayet stars observed in Cygnus OB2 are added to the simulation, which continues to 2 Myr. Using a high-resolution grid in the centre of the domain, we can resolve the most massive stars individually, which allows us to probe the kinetic structures at small (parsec) scales. We find that, although the cluster excavates a spherical "superbubble" cavity, the stellar population is too loosely distributed to blow a large-scale cluster wind termination shock, and that collective effects from wind-wind interactions are much less efficient than usually assumed. This challenges our understanding of the ultra-high energy emission observed from the region. △ Less

Submitted 19 June, 2024; originally announced June 2024.

Comments: 15 pages, 12 figures

Peering into the Wolf-Rayet phenomenon through [WO] and [WC] stars

Authors: Jesús A. Toalá, Helge Todt, Andreas A. C. Sander

Abstract: Spectroscopic observations have shown for decades that the Wolf-Rayet (WR) phenomenon is ubiquitous among stars with different initial masses. Although much effort to understand the winds from massive WR stars has been presented in the literature, not much has been done for such type of stars in the low-mass range. Here we present an attempt to understand the winds from [WR]-type stars using resul… ▽ More Spectroscopic observations have shown for decades that the Wolf-Rayet (WR) phenomenon is ubiquitous among stars with different initial masses. Although much effort to understand the winds from massive WR stars has been presented in the literature, not much has been done for such type of stars in the low-mass range. Here we present an attempt to understand the winds from [WR]-type stars using results from spectral analyses with the full non-LTE stellar atmosphere code PoWR. These results are put into context with the properties of massive WR stars. We found that WC+[WC] stars and WO+[WO] stars create independent sequences in the mass-loss rate ($\dot{M}$) and modified wind momentum ($D_\mathrm{mom}$) versus luminosity ($L$) diagrams. Our analysis indicates that even when the winds of WR and [WR] stars become optically thin, there is no breakdown of the general mass-loss trend, contrary to the observed ``weak wind phenomenon'' in OB stars. We report that all WR-type stars studied here broadly define single sequences in the wind efficiency ($η$) versus transformed mass-loss rate ($\dot{M}_\mathrm{t}$), the $\dot{M}_\mathrm{t}$-$T_\mathrm{eff}$ diagram, and the $(L, T_\mathrm{eff}, \dot{M})$ space, which suggest these to be fundamental properties of the WR phenomenon (regardless of the mass range), at least for WR-type stars of the O and C sequences. Our analytical estimations could drive computations of future stellar evolution models for WR-type stars. △ Less

Submitted 16 May, 2024; originally announced May 2024.

Comments: 11 pages; 5 Figures; 2 Tables; Accepted to MNRAS

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

Impact of main-sequence mass loss on the appearance, structure and evolution of Wolf-Rayet stars

Authors: Joris Josiek, Sylvia Ekström, Andreas A. C. Sander

Abstract: Stellar winds are one of the most important drivers of massive star evolution and a vital source of chemical, mechanical, and radiative feedback. Despite its significance, mass loss remains a major uncertainty in stellar evolution models. Particularly the interdependencies of different approaches with subsequent evolutionary stages and predicted observable phenomena are far from being systematical… ▽ More Stellar winds are one of the most important drivers of massive star evolution and a vital source of chemical, mechanical, and radiative feedback. Despite its significance, mass loss remains a major uncertainty in stellar evolution models. Particularly the interdependencies of different approaches with subsequent evolutionary stages and predicted observable phenomena are far from being systematically understood. In this study, we examine the impact of main sequence mass loss on the structure of massive stars throughout their evolution. A particular focus is placed on the consequences for entering the Wolf-Rayet (WR) regime and the subsequent evolution. Using the Geneva stellar evolution code, we compute grids of single, non-rotating stellar models at solar and Large Magellanic Cloud (LMC) metallicity of initial masses between 20 and 120 solar masses, with two representative prescriptions for high and low main sequence mass loss. We obtain detailed numerical predictions regarding the structure and evolution of massive stars, and infer the role of main sequence mass loss by comparison of the mass-loss rate prescriptions. We present implications for the overall evolutionary trajectory, including the evolution of WR stars, as well as the effect on stellar yields and stellar populations. Mass loss during the main sequence plays an important role due to its ability to affect the sequence and duration of all subsequent phases. We identify several distinct evolutionary paths for massive stars which are significantly influenced by the chosen main sequence mass-loss description. We also discuss the impact of uncertainties other than mass loss on the evolution, in particular those relating to convection. We further demonstrate that not just the total mass loss, but the specific mass-loss history throughout a star's life is a crucial determinant of many aspects, such as the resulting stellar yields. △ Less

Submitted 22 April, 2024; originally announced April 2024.

Comments: 15 pages, 12 figures, accepted for publication in A&A

The Wide-field Spectroscopic Telescope (WST) Science White Paper

Authors: Vincenzo Mainieri, Richard I. Anderson, Jarle Brinchmann, Andrea Cimatti, Richard S. Ellis, Vanessa Hill, Jean-Paul Kneib, Anna F. McLeod, Cyrielle Opitom, Martin M. Roth, Paula Sanchez-Saez, Rodolfo Smiljanic, Eline Tolstoy, Roland Bacon, Sofia Randich, Angela Adamo, Francesca Annibali, Patricia Arevalo, Marc Audard, Stefania Barsanti, Giuseppina Battaglia, Amelia M. Bayo Aran, Francesco Belfiore, Michele Bellazzini, Emilio Bellini , et al. (192 additional authors not shown)

Abstract: The Wide-field Spectroscopic Telescope (WST) is proposed as a new facility dedicated to the efficient delivery of spectroscopic surveys. This white paper summarises the initial concept as well as the corresponding science cases. WST will feature simultaneous operation of a large field-of-view (3 sq. degree), a high multiplex (20,000) multi-object spectrograph (MOS) and a giant 3x3 sq. arcmin integ… ▽ More The Wide-field Spectroscopic Telescope (WST) is proposed as a new facility dedicated to the efficient delivery of spectroscopic surveys. This white paper summarises the initial concept as well as the corresponding science cases. WST will feature simultaneous operation of a large field-of-view (3 sq. degree), a high multiplex (20,000) multi-object spectrograph (MOS) and a giant 3x3 sq. arcmin integral field spectrograph (IFS). In scientific capability these requirements place WST far ahead of existing and planned facilities. Given the current investment in deep imaging surveys and noting the diagnostic power of spectroscopy, WST will fill a crucial gap in astronomical capability and work synergistically with future ground and space-based facilities. This white paper shows that WST can address outstanding scientific questions in the areas of cosmology; galaxy assembly, evolution, and enrichment, including our own Milky Way; origin of stars and planets; time domain and multi-messenger astrophysics. WST's uniquely rich dataset will deliver unforeseen discoveries in many of these areas. The WST Science Team (already including more than 500 scientists worldwide) is open to the all astronomical community. To register in the WST Science Team please visit https://www.wstelescope.com/for-scientists/participate △ Less

Submitted 12 April, 2024; v1 submitted 8 March, 2024; originally announced March 2024.

Comments: 194 pages, 66 figures. Comments are welcome (wstelescope@gmail.com)

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

Chemical abundances and ionizing mechanisms in the star-forming double-ring of AM 0644-741 using MUSE data

Authors: V. M. A. Gómez-González, Y. D. Mayya, J. Zaragoza-Cardiel, G. Bruzual, S. Charlot, G. Ramos-Larios, L. M. Oskinova, A. A. C. Sander, S. Reyero Serantes

Abstract: We present the analysis of archival Very Large Telescope (VLT) Multi-Unit Spectroscopic Explorer (MUSE) observations of 179 HII regions in the star-forming double-ring collisional galaxy AM 0644-741 at 98.6 Mpc. We determined ionic abundances of He, N, O and Fe using the direct method for the brightest H II region (ID 39); we report $\log\rm{(\frac{N}{O})}=-1.3\pm0.2$ and… ▽ More We present the analysis of archival Very Large Telescope (VLT) Multi-Unit Spectroscopic Explorer (MUSE) observations of 179 HII regions in the star-forming double-ring collisional galaxy AM 0644-741 at 98.6 Mpc. We determined ionic abundances of He, N, O and Fe using the direct method for the brightest H II region (ID 39); we report $\log\rm{(\frac{N}{O})}=-1.3\pm0.2$ and $12+\log\rm{(\frac{O}{H})}=8.9\pm0.2$. We also find the so-called `blue-bump', broad He II $\lambda4686$, in the spectrum of this knot of massive star-formation; its luminosity being consistent with the presence of $\sim430$ Wolf-Rayet (WR) stars of the Nitrogen late-type. We determined the O abundances for 137 HII regions using the strong-line method; we report a median value of $12+\log\rm{(\frac{O}{H})}=8.5\pm0.8$. The location of three objects, including the WR complex, coincide with that of an Ultra Luminous X-ray source. Nebular He II is not detected in any H II region. We investigate the physical mechanisms responsible for the observed spectral lines using appropriate diagnostic diagrams and ionization models. We find that the H II regions are being photoionized by star clusters with ages $\sim2.5-20$ Myr and ionization potential $-3.5<$$\log\langle U\rangle$$<-3.0$. In these diagrams, a binary population is needed to reproduce the observables considered in this work. △ Less

Submitted 20 February, 2024; originally announced February 2024.

Comments: 20 pages. Accepted in MNRAS

Aperture and Resolution Effects on Ultraviolet Star-Forming Properties: Insights from Local Galaxies and Implications for High-Redshift Observations

Authors: Ilyse Clark, Danielle A. Berg, Claus Leitherer, Karla Z. Arellano-Cordova, Andreas A. C. Sander

Abstract: We present an analysis of the effects of spectral resolution and aperture scales on derived galaxy properties using far-ultraviolet (FUV) spectra of local star-forming galaxies from the International Ultraviolet Explorer (R~250, FOV~10"x20") and Cosmic Origins Spectrograph on the Hubble Space Telescope (R~15,000, FOV~2.5"). Using these spectra, we measured FUV luminosities, spectral slopes, dust a… ▽ More We present an analysis of the effects of spectral resolution and aperture scales on derived galaxy properties using far-ultraviolet (FUV) spectra of local star-forming galaxies from the International Ultraviolet Explorer (R~250, FOV~10"x20") and Cosmic Origins Spectrograph on the Hubble Space Telescope (R~15,000, FOV~2.5"). Using these spectra, we measured FUV luminosities, spectral slopes, dust attenuation, and equivalent widths. We find that galaxies with one dominant stellar cluster have FUV properties that are independent of aperture size, while galaxies with multiple bright clusters are sensitive to the total light fraction captured by the aperture. Additionally, we find significant correlations between the strength of stellar and interstellar absorption-lines and metallicity, indicating metallicity-dependent line-driven stellar winds and interstellar macroscopic gas flows shape the stellar and interstellar spectral lines, respectively. The observed line-strength versus metallicity relation of stellar-wind lines agrees with the prediction of population synthesis models for young starbursts. In particular, measurements of the strong stellar CIV 1548,1550 line provide an opportunity to determine stellar abundances as a complement to gas-phase abundances. We provide a relation between the equivalent width of the CIV line and the oxygen abundance of the galaxy. We discuss this relation in terms of the stellar-wind properties of massive stars. As the driving lines in stellar winds are mostly ionized iron species, the CIV line may eventually offer a method to probe alpha-element-to-iron ratios in star-forming galaxies once consistent models with non-solar abundance ratios are available. These results have important implications for the galaxy-scale, low-resolution observations of high-redshift galaxies from JWST (R~100-3,500). △ Less

Submitted 16 January, 2024; originally announced January 2024.

Comments: This paper has 31 pages total, 11 figures, and a figureset. Accepted for publication in ApJ

Optically-thick Structure in Early B Type Supergiant Stellar Winds at Low Metallicities

Authors: Timothy N. Parsons, Raman K. Prinja, Matheus Bernini-Peron, Alex W. Fullerton, Derck L. Massa, Lidia M. Oskinova, Daniel Pauli, Andreas A. C. Sander, Matthew J. Rickard

Abstract: Accurate determination of mass-loss rates from massive stars is important to understanding stellar and galactic evolution and enrichment of the interstellar medium. Large-scale structure and variability in stellar winds have significant effects on mass-loss rates. Time-series observations provide direct quantification of such variability. Observations of this nature are available for some Galactic… ▽ More Accurate determination of mass-loss rates from massive stars is important to understanding stellar and galactic evolution and enrichment of the interstellar medium. Large-scale structure and variability in stellar winds have significant effects on mass-loss rates. Time-series observations provide direct quantification of such variability. Observations of this nature are available for some Galactic early supergiant stars but not yet for stars in lower metallicity environments such as the Magellanic Clouds. We utilise ultraviolet spectra from the Hubble Space Telescope ULLYSES program to demonstrate that the presence of structure in stellar winds of supergiant stars at low metallicities may be discerned from single-epoch spectra. We find evidence that, for given stellar luminosities and mean stellar wind optical depths, structure is more prevalent at higher metallicities. We confirm, at Large Magellanic Cloud (0.5 Z_solar), Small Magellanic Cloud (0.2 Z_solar) and lower (0.14 -- 0.1 Z_solar) metallicities, earlier Galactic results that there does not appear to be correlation between the degree of structure in stellar winds of massive stars and stellar effective temperature. Similar lack of correlation is found with regard to terminal velocity of stellar winds. Additional and revised values for radial velocities of stars and terminal velocities of stellar winds are presented. Direct evidence of temporal variability, on timescales of several days, in stellar wind at low metallicity is found. We illustrate that narrow absorption components in wind-formed profiles of Galactic OB stellar spectra remain common in early B supergiant spectra at low metallicities, providing means for better constraining hot, massive star mass-loss rates. △ Less

Submitted 2 January, 2024; originally announced January 2024.

Comments: 18 pages plus Appendices, 9 figures plus 37 figures in Appendices. Accepted by MNRAS

A soft and transient ultraluminous X-ray source with 6-h modulation in the NGC 300 galaxy

Authors: A. Sacchi, P. Esposito, D. de Martino, R. Soria, G. L. Israel, A. A. C. Sander, L. Sidoli, D. A. H. Buckley, I. M. Monageng, A. Tiengo, M. Arca Sedda, C. Pinto, R. Di Stefano, M. Imbrogno, A. Carleo, G. Rivolta

Abstract: We investigate the nature of CXOU J005440.5-374320 (J0054), a peculiar bright ($\sim$$4\times10^{39}$ erg/s) and soft X-ray transient in the spiral galaxy NGC 300 with a 6-hour periodic flux modulation that was detected in a 2014 Chandra observation. Subsequent observations with Chandra and XMM-Newton, as well as a large observational campaign of NGC 300 and its sources performed with the Swift Ne… ▽ More We investigate the nature of CXOU J005440.5-374320 (J0054), a peculiar bright ($\sim$$4\times10^{39}$ erg/s) and soft X-ray transient in the spiral galaxy NGC 300 with a 6-hour periodic flux modulation that was detected in a 2014 Chandra observation. Subsequent observations with Chandra and XMM-Newton, as well as a large observational campaign of NGC 300 and its sources performed with the Swift Neil Gehrels Observatory, showed that this source exhibits recurrent flaring activity: four other outbursts were detected across $\sim$8 years of monitoring. Using data from the Swift/UVOT archive and from the XMM-Newton/OM and Gaia catalogues, we noted the source is likely associated with a bright blue optical/ultraviolet counterpart. This prompted us to perform follow-up observations with the Southern African Large Telescope in December 2019. With the multi-wavelength information at hand, we discuss several possibilities for the nature of J0054. Although none is able to account for the full range of the observed peculiar features, we found that the two most promising scenarios are a stellar-mass compact object in a binary system with a Wolf$-$Rayet star companion, or the recurrent tidal stripping of a stellar object trapped in a system with an intermediate-mass ($\sim1000$ $M_\odot$) black hole. △ Less

Submitted 24 November, 2023; originally announced November 2023.

Comments: 13 pages, 11 Figures, 3 Tables (the Table in appendix A will be available in the published version). Accepted for publication in A&A

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)

Clumping and X-Rays in cooler B supergiant stars

Authors: Matheus Bernini-Peron, Wagner L. F. Marcolino, Andreas A. C. Sander, Jean-Claude Bouret, Varsha Ramachandran, Julian Saling, Fabian R. N. Schneider, Lidia M. Oskinova, Francisco Najarro

Abstract: B supergiants (BSGs) are evolved stars with effective temperatures between 10 to 30 kK and are important to understand massive star evolution. Located on the edge of the line-driven wind regime, the study of their atmospheres is helpful to understand phenomena such as the bi-stability jump. Key UV features of their spectra have so far not been reproduced by models for types later than B1. Here, we… ▽ More B supergiants (BSGs) are evolved stars with effective temperatures between 10 to 30 kK and are important to understand massive star evolution. Located on the edge of the line-driven wind regime, the study of their atmospheres is helpful to understand phenomena such as the bi-stability jump. Key UV features of their spectra have so far not been reproduced by models for types later than B1. Here, we aim to remedy this situation via spectral analysis that accounts for wind clumping and X-rays. In addition, we investigate the evolutionary status of our sample stars based on the obtained stellar parameters. We determined parameters via quantitative spectroscopy using CMFGEN and PoWR codes. The models were compared to UV and optical data of four BSGs: HD206165, HD198478, HD53138, and HD164353. We also study the evolutionary status of our sample using GENEC and MESA tracks. When including clumping and X-rays, we find good agreements between synthetic and observed spectra for our sample stars. For the first time, we reproduced key lines in the UV. For that, we require a moderately clumped wind (f_infty > ~0.5). We also infer relative X-ray luminosities of ~10^-7.5 to 10^-8 -- lower than the typical ratio of 10^-7. Moreover, we find a possible mismatch between evolutionary and spectroscopic masses, which could be related to the mass-discrepancy problem present in other OB stars. Our results provide evidence that X-rays and clumping are needed to describe the winds of cool BSGs. However, their winds seem less structured than in earlier type stars. This aligns with observational X-rays and clumping constraints as well as recent hydrodynamical simulations. The BSGs' evolutionary status appears diverse: some objects are potentially post-red supergiants or merger products. The wind parameters provide evidence for a moderate mass-loss rate increase around the bi-stability jump. Abstract abridged △ Less

Submitted 21 July, 2023; v1 submitted 11 July, 2023; originally announced July 2023.

Comments: 27 pages, 22 figures, accepted for publication in A&A

Journal ref: A&A 677, A50 (2023)

Very Massive Stars and Pair-Instability Supernovae: Mass-loss Framework for low Metallicity

Authors: Gautham N. Sabhahit, Jorick S. Vink, Andreas A. C. Sander, Erin R. Higgins

Abstract: Very massive stars (VMS) up to 200-300 $M_\odot$ have been found in the Local Universe. If they would lose little mass they produce intermediate-mass black holes or pair-instability supernovae (PISNe). Until now, VMS modellers have extrapolated mass-loss vs. metallicity ($Z$) exponents from optically-thin winds, resulting in a range of PISN thresholds that might be unrealistically high in $Z$, as… ▽ More Very massive stars (VMS) up to 200-300 $M_\odot$ have been found in the Local Universe. If they would lose little mass they produce intermediate-mass black holes or pair-instability supernovae (PISNe). Until now, VMS modellers have extrapolated mass-loss vs. metallicity ($Z$) exponents from optically-thin winds, resulting in a range of PISN thresholds that might be unrealistically high in $Z$, as VMS develop optically-thick winds. We utilize the transition mass-loss rate of Vink and Gräfener (2012) that accurately predicts mass-loss rates of Of/WNh ("slash") stars that characterize the morphological transition from absorption-dominated O-type spectra to emission-dominated WNh spectra. We develop a wind efficiency framework, where optically thin winds transition to enhanced winds, enabling us to study VMS evolution at high redshift where individual stars cannot be resolved. We present a MESA grid covering $Z_\odot/2$ to $Z_\odot/100$. VMS above the transition evolve towards lower luminosity, skipping the cool supergiant phase but directly forming pure He stars at the end of hydrogen burning. Below the transition, VMS evolve as cooler luminous blue variables (LBVs) or yellow hypergiants (YHGs), naturally approaching the Eddington limit. Strong winds in this YHG/LBV regime -- combined with a degeneracy in luminosity -- result in a mass-loss runaway where a decrease in mass increases wind mass loss. Our models indicate an order-of-magnitude lower threshold than usually assumed, at $Z_\odot/20$ due to our mass-loss runaway. While future work on LBV mass loss could affect the PISN threshold, our framework will be critical for establishing definitive answers on the PISN threshold and galactic chemical evolution modelling. △ Less

Submitted 20 June, 2023; originally announced June 2023.

Comments: Accepted for publication in MNRAS, 18 pages, 8 figures, 2 appendix

From Dust to Nanodust: Resolving Circumstellar Dust from the Colliding-Wind Binary Wolf-Rayet (WR) 140

Authors: Ryan M. Lau, Jason Wang, Matthew J. Hankins, Thayne Currie, Vincent Deo, Izumi Endo, Olivier Guyon, Yinuo Han, Anthony P. Jones, Nemanja Jovanovic, Julien Lozi, Anthony F. J. Moffat, Takashi Onaka, Garreth Ruane, Andreas A. C. Sander, Samaporn Tinyanont, Peter G. Tuthill, Gerd Weigelt, Peredur M. Williams, Sebastien Vievard

Abstract: Wolf-Rayet (WR) 140 is the archetypal periodic dust-forming colliding-wind binary that hosts a carbon-rich WR (WC) star and an O-star companion with an orbital period of 7.93 years and an orbital eccentricity of 0.9. Throughout the past several decades, multiple dust-formation episodes from WR 140 have been observed that are linked to the binary orbit and occur near the time of periastron passage.… ▽ More Wolf-Rayet (WR) 140 is the archetypal periodic dust-forming colliding-wind binary that hosts a carbon-rich WR (WC) star and an O-star companion with an orbital period of 7.93 years and an orbital eccentricity of 0.9. Throughout the past several decades, multiple dust-formation episodes from WR 140 have been observed that are linked to the binary orbit and occur near the time of periastron passage. Given its predictable dust-formation episodes, WR 140 presents an ideal astrophysical laboratory for investigating the formation and evolution of dust in the hostile environment around a massive binary system. In this paper, we present near- and mid-infrared (IR) spectroscopic and imaging observations of WR 140 with Subaru/SCExAO+CHARIS, Keck/NIRC2+PyWFS, and Subaru/COMICS taken between 2020 June and Sept that resolve the circumstellar dust emission linked to its most recent dust-formation episode in 2016 Dec. Our spectral energy distribution (SED) analysis of WR 140's resolved circumstellar dust emission reveals the presence of a hot ($T_\mathrm{d}\sim1000$ K) near-IR dust component that is co-spatial with the previously known and cooler ($T_\mathrm{d}\sim500$ K) mid-IR dust component composed of $300-500$ Å-sized dust grains. We attribute the hot near-IR dust emission to the presence of nano-sized ("nanodust") grains and suggest they were formed from grain-grain collisions or the rotational disruption of the larger grain size population by radiative torques in the strong radiation field from the central binary. Lastly, we speculate on the astrophysical implications of nanodust formation around colliding-wind WC binaries, which may present an early source of carbonaceous nanodust in the interstellar medium. △ Less

Submitted 23 May, 2023; originally announced May 2023.

Comments: 21 pages, 8 figures, Accepted for publication in ApJ

A partially stripped massive star in a Be binary at low metallicity: A missing link towards Be X-ray binaries and double neutron star mergers

Authors: V. Ramachandran, J. Klencki, A. A. C. Sander, D. Pauli, T. Shenar, L. M. Oskinova, W. -R. Hamann

Abstract: Standard binary evolutionary models predict a significant population of core helium-burning stars that lost their hydrogen-rich envelope after mass transfer via Roche-lobe overflow. However, there is a scarcity of observations of such stripped stars in the intermediate mass regime (~1.5 - 8$ M_{\odot}$), which are thought to be prominent progenitors of SN Ib/c. Especially at low metallicity, a sig… ▽ More Standard binary evolutionary models predict a significant population of core helium-burning stars that lost their hydrogen-rich envelope after mass transfer via Roche-lobe overflow. However, there is a scarcity of observations of such stripped stars in the intermediate mass regime (~1.5 - 8$ M_{\odot}$), which are thought to be prominent progenitors of SN Ib/c. Especially at low metallicity, a significant fraction of these stars is expected to be only partially stripped, retaining a significant amount of hydrogen on their surfaces. For the first time, we discovered a partially stripped massive star in a binary with a Be-type companion located in the Small Magellanic Cloud (SMC) using a detailed spectroscopic analysis. The stripped-star nature of the primary is revealed by the extreme CNO abundance pattern and very high luminosity-to-mass ratio, which suggest that the primary is likely shell-hydrogen burning. Our target SMCSGS-FS 69 is the most luminous and most massive system among the known stripped star + Be binaries, with Mstripped ~3$ M_{\odot}$ and MBe ~17$ M_{\odot}$. Binary evolutionary tracks suggest an initial mass of Mini $\gtrsim 12 M_{\odot}$ for the stripped star and predict it to be in a transition phase towards a hot compact He star, which will eventually produce a stripped-envelope supernova. Our target marks the first representative of a so-far missing evolutionary stage in the formation pathway of Be X-ray binaries and double neutron star mergers. △ Less

Submitted 12 June, 2023; v1 submitted 23 May, 2023; originally announced May 2023.

Comments: Accepted for publication as A&A Letter to the Editor

Journal ref: A&A 674, L12 (2023)

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

Spectroscopic and evolutionary analyses of the binary system AzV 14 outline paths toward the WR stage at low metallicity

Authors: D. Pauli, L. M. Oskinova, W. -R. Hamann, D. M. Bowman, H. Todt, T. Shenar, A. A. C. Sander, C. Erba, V. M. A. Gómez-González, C. Kehrig, J. Klencki, R. Kuiper, A. Mehner, S. E. de Mink, M. S. Oey, V. Ramachandran, A. Schootemeijer, S. Reyero Serantes, A. Wofford

Abstract: The origin of the observed population of Wolf-Rayet (WR) stars in low-metallicity (low-Z) galaxies, such as the Small Magellanic Cloud (SMC), is not yet understood. Standard, single-star evolutionary models predict that WR stars should stem from very massive O-type star progenitors, but these are very rare. On the other hand, binary evolutionary models predict that WR stars could originate from pr… ▽ More The origin of the observed population of Wolf-Rayet (WR) stars in low-metallicity (low-Z) galaxies, such as the Small Magellanic Cloud (SMC), is not yet understood. Standard, single-star evolutionary models predict that WR stars should stem from very massive O-type star progenitors, but these are very rare. On the other hand, binary evolutionary models predict that WR stars could originate from primary stars in close binaries. We conduct an analysis of the massive O star, AzV 14, to spectroscopically determine its fundamental and stellar wind parameters, which are then used to investigate evolutionary paths from the O-type to the WR stage with stellar evolutionary models. Multi-epoch UV and optical spectra of AzV 14 are analyzed using the non-LTE stellar atmosphere code PoWR. An optical TESS light curve was extracted and analyzed using the PHOEBE code. The obtained parameters are put into an evolutionary context, using the MESA code. AzV 14 is a close binary system consisting of two similar main sequence stars with masses of 32 Msol. Both stars have weak stellar winds with mass-loss rates of log $\dot{M}$ = -7.7. Binary evolutionary models can explain the empirically derived stellar and orbital parameters. The model predicts that the primary will evolve into a WR star with T = 100 kK, while the secondary, which will accrete significant amounts of mass during the first mass transfer phase, will become a cooler WR star with T = 50 kK and are predicted to have compared to other WR stars increased oxygen abundances. This model prediction is supported by a spectroscopic analysis of a WR star in the SMC. We hypothesize that the populations of WR stars in low-Z galaxies may have bimodal temperature distributions. Hotter WR stars might originate from primary stars, while cooler WR stars are the evolutionary descendants of the secondary stars if they accreted a significant amount of mass. △ Less

Submitted 7 March, 2023; originally announced March 2023.

Comments: 21 pages (13 main body + 8 appendix), 16 figures, 9 tables

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

Exploring the influence of different velocity fields on Wolf-Rayet star spectra

Authors: Roel R. Lefever, Andreas A. C. Sander, Tomer Shenar, Luka G. Poniatowski, Karan Dsilva, Helge Todt

Abstract: Given their strong stellar winds, Wolf-Rayet (WR) stars exhibit emission line spectra that are predominantly formed in expanding atmospheric layers. The description of the wind velocity field $v(r)$ is therefore a crucial ingredient in the spectral analysis of WR stars, possibly influencing the determination of stellar parameters. In view of this, we perform a systematic study by simulating a sequ… ▽ More Given their strong stellar winds, Wolf-Rayet (WR) stars exhibit emission line spectra that are predominantly formed in expanding atmospheric layers. The description of the wind velocity field $v(r)$ is therefore a crucial ingredient in the spectral analysis of WR stars, possibly influencing the determination of stellar parameters. In view of this, we perform a systematic study by simulating a sequence of WR-star spectra for different temperatures and mass-loss rates using $β$-type laws with $0.5\leqβ\leq 20$. We quantify the impact of varying $v(r)$ by analysing diagnostic lines and spectral classifications of emergent model spectra computed with the Potsdam Wolf-Rayet (PoWR) code. We additionally cross-check these models with hydrodynamically consistent -- hydro -- model atmospheres. Our analysis confirms that the choice of the $β$-exponent has a strong impact on WR-star spectra, affecting line widths, line strengths and line profiles. In some parameter regimes, the entire range of WR subtypes could be covered. Comparison with observed WR stars and hydro models revealed that values of $β\gtrsim 8$ are unlikely to be realized in nature, but a range of $β$-values needs to be considered in spectral analysis. UV spectroscopy is crucial here to avoid an underestimation of the terminal velocity $v_\infty$. Neither single- nor double-$β$ descriptions yield an acceptable approximation of the inner wind when compared to hydro models. Instead, we find temperature shifts to lower $T_{2/3}$ when employing a hydro model. Additionally, there are further hints that round-lined profiles seen in several early WN stars are an effect from non-$β$ velocity laws. △ Less

Submitted 28 February, 2023; v1 submitted 27 February, 2023; originally announced February 2023.

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

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

The temperature dependency of Wolf-Rayet-type mass loss: An exploratory study for winds launched by the hot iron bump

Authors: A. A. C. Sander, R. R. Lefever, L. G. Poniatowski, V. Ramachandran, G. N. Sabhahit, J. S. Vink

Abstract: CONTEXT: The mass loss of He-burning stars, which are partially or completely stripped of their outer hydrogen envelope, is a catalyst of the cosmic matter cycle and decisive ingredient of massive star evolution. Yet, its theoretical fundament is only starting to emerge with major dependencies still to be uncovered. AIMS: A temperature or radius dependence is usually not included in descriptions… ▽ More CONTEXT: The mass loss of He-burning stars, which are partially or completely stripped of their outer hydrogen envelope, is a catalyst of the cosmic matter cycle and decisive ingredient of massive star evolution. Yet, its theoretical fundament is only starting to emerge with major dependencies still to be uncovered. AIMS: A temperature or radius dependence is usually not included in descriptions for the mass loss of classical Wolf-Rayet (cWR) stars, despite being crucial for other hot star wind domains. We thus aim to determine whether such a dependency will also be necessary for a comprehensive description of mass loss in the cWR regime. METHODS: Sequences of dynamically consistent atmosphere models were calculated with the hydrodynamic branch of the PoWR code along the temperature domain, using different choices for luminosity, mass, and surface abundances. For the first time, we allowed nonmonotonic velocity fields when solving the equation of motion. The resulting velocity structures were then interpolated for the comoving-frame radiative transfer, ensuring that the main wind characteristics were preserved. RESULTS: We find a strong dependence of the mass-loss rate with the temperature of the critical/sonic point which mainly reflects the different radii and resulting gravitational accelerations. Moreover, we obtain a relation between the observed effective temperature and the transformed mass-loss rate which seems to be largely independent of the underlying stellar parameters. The relation shifts for different clumping factors in the outer wind. Below a characteristic value of -4.5, the slope of this relation changes and the winds become transparent for He II ionizing photons. CONCLUSIONS: The mass loss of cWR stars is a high-dimensional problem but also shows inherent scalings which can be used to obtain an approximation of the observed effective temperature. (...) △ Less

Submitted 4 January, 2023; originally announced January 2023.

Comments: 16 pages + 5 page appendix, 17+9 figures, 3+2 tables. Accepted for publication in A&A

Journal ref: A&A 670, A83 (2023)

The Driving of Hot Star Winds

Authors: Andreas A. C. Sander

Abstract: In the regime of hot stars, winds were not seen as a common thing until the era of UV astronomy. Since we have access to the UV wavelength range, it has become clear that winds are not an exotic phenomenon limited to some special objects, but actually ubiquitous among hot and massive stars. The opacities due to spectral lines are the decisive ingredient that allows hot, massive stars to launch pow… ▽ More In the regime of hot stars, winds were not seen as a common thing until the era of UV astronomy. Since we have access to the UV wavelength range, it has become clear that winds are not an exotic phenomenon limited to some special objects, but actually ubiquitous among hot and massive stars. The opacities due to spectral lines are the decisive ingredient that allows hot, massive stars to launch powerful winds. While the fundamental principles of these so-called line-driven winds have been realized decades ago, their proper quantitative prediction is still a major challenge today. Established theoretical and empirical descriptions have allowed us to make major progress on all astrophysical scales. However, we are now reaching their limitations as we still lack various fundamental insights on the nature of hot star winds, thereby hampering us from drawing deeper conclusions, not least when dealing with stellar or sub-stellar companions. This has spawned a new generation of researchers searching for answers with a yet unprecedented level of detail in observational and new theoretical approaches. In these proceedings, the fundamental principles of driving hot star winds will be briefly reviewed. Starting from the classical CAK theory and its extensions, over Monte Carlo and recent comoving-frame-based simulations, the different methods to describe and model the acceleration of hot star winds will be introduced. The review continues with briefly discussing instabilities as well as qualitative and quantitative insights for OB- and Wolf-Rayet-star winds. Moreover, the challenges of companions and their impact on radiation-driven winds are outlined. △ Less

Submitted 21 November, 2022; originally announced November 2022.

Comments: 14 pages, 4 figures, to be published in the Proceedings of the International Astronomical Union for the IAU Symposium 370 "Winds of Stars and Exoplanets" (eds. A.A. Vidotto, L. Fossati, J.S. Vink)

Massive Stars in the Far and Extreme Ultraviolet

Authors: Andreas A. C. Sander

Abstract: From the main sequence to their late evolutionary stages, massive stars spend most of their life as hot stars. Due to their high effective temperatures, the maximum of their emitted flux falls into the ultraviolet (UV) regime. Consequently, these stars emit a significant number of photons with energies sufficiently high enough to ionize hydrogen and other elements. As simple as these fundamental c… ▽ More From the main sequence to their late evolutionary stages, massive stars spend most of their life as hot stars. Due to their high effective temperatures, the maximum of their emitted flux falls into the ultraviolet (UV) regime. Consequently, these stars emit a significant number of photons with energies sufficiently high enough to ionize hydrogen and other elements. As simple as these fundamental considerations are, as complex is a realistic estimate of the resulting ionizing fluxes, in particular for energies above 54 eV. Estimating the ionizing flux budget of hot stars requires accurate models of their spectral energy distributions (SEDs), covering in particular the far and extreme UV region. Modern atmosphere models that incorporate the so-called line-blanketing effect, i.e. taking into account the millions of lines from iron and other elements, yield a complex picture, illustrating that the SED of a hot, massive star usually deviates significantly from a blackbody. △ Less

Submitted 10 November, 2022; originally announced November 2022.

Comments: 3 pages, 4 figures, to be published in the International Astronomical Union Proceedings Series for the IAU GA 2022 (Astronomy in Focus, Focus Meeting 4, ed. José Espinosa)

Investigating the impact of different velocity fields on the spectral appearance of Wolf-Rayet stars

Authors: Roel R. Lefever, Tomer Shenar, Andreas A. C. Sander, Luka Poniatowski, Karan Dsilva, Helge Todt

Abstract: The emission line spectra of WR stars are often formed completely in the optically thick stellar wind. Hence, any assumption on the wind velocity law in a spectral analysis has a profound impact on the determination of the stellar parameters. By comparing Potsdam Wolf-Rayet (PoWR) model spectra calculated with different $β$ laws, we show that the velocity field heavily influences the spectra: by u… ▽ More The emission line spectra of WR stars are often formed completely in the optically thick stellar wind. Hence, any assumption on the wind velocity law in a spectral analysis has a profound impact on the determination of the stellar parameters. By comparing Potsdam Wolf-Rayet (PoWR) model spectra calculated with different $β$ laws, we show that the velocity field heavily influences the spectra: by using the appropriate $β$ laws, the entire range of late and early types can be covered with the same stellar model. △ Less

Submitted 13 September, 2022; originally announced September 2022.

Comments: 2 pages, 1 figure, to be published in the Proceedings of the International Astronomical Union for the IAU Symposium 361 "Massive Stars Near and Far"

The hydrogen clock to infer the upper stellar mass

Authors: Erin R. Higgins, Jorick S. Vink, Gautham N. Sabhahit, Andreas A. C. Sander

Abstract: The most massive stars dominate the chemical enrichment, mechanical and radiative feedback, and energy budget of their host environments. Yet how massive stars initially form and how they evolve throughout their lives is ambiguous. The mass loss of the most massive stars remains a key unknown in stellar physics, with consequences for stellar feedback and populations. In this work, we compare grids… ▽ More The most massive stars dominate the chemical enrichment, mechanical and radiative feedback, and energy budget of their host environments. Yet how massive stars initially form and how they evolve throughout their lives is ambiguous. The mass loss of the most massive stars remains a key unknown in stellar physics, with consequences for stellar feedback and populations. In this work, we compare grids of very massive star (VMS) models with masses ranging from 80-1000Msun, for a range of input physics. We include enhanced winds close to the Eddington limit as a comparison to standard O-star winds, with consequences for present-day observations of ~50-100Msun stars. We probe the relevant surface H abundances (Xs) to determine the key traits of VMS evolution compared to O stars. We find fundamental differences in the behaviour of our models with the enhanced-wind prescription, with a convergence on the stellar mass at 1.6 Myr, regardless of the initial mass. It turns out that Xs is an important tool in deciphering the initial mass due to the chemically homogeneous nature of VMS above a mass threshold. We use Xs to break the degeneracy of the initial masses of both components of a detached binary, and a sample of WNh stars in the Tarantula nebula. We find that for some objects, the initial masses are unrestricted and, as such, even initial masses of the order 1000Msun are not excluded. Coupled with the mass turnover at 1.6 Myr, Xs can be used as a 'clock' to determine the upper stellar mass. △ Less

Submitted 1 September, 2022; originally announced September 2022.

Comments: Accepted for publication in MNRAS, 14 figures

The enigmatic winds of Wolf-Rayet stars: Results from dynamically consistent atmosphere modelling

Authors: Andreas A. C. Sander

Abstract: Line-driven stellar winds are ubiquitous among hot massive stars. In some cases they can become so strong, that the whole star is cloaked by an optically thick wind. The strong outflow gives rise to large emission lines, defining the class of so-called Wolf-Rayet (WR) stars. While being major players in the evolution of massive stars, the formation of heavy black holes,and the distribution of elem… ▽ More Line-driven stellar winds are ubiquitous among hot massive stars. In some cases they can become so strong, that the whole star is cloaked by an optically thick wind. The strong outflow gives rise to large emission lines, defining the class of so-called Wolf-Rayet (WR) stars. While being major players in the evolution of massive stars, the formation of heavy black holes,and the distribution of elements, the occurrence and nature of WR winds is still quite enigmatic. A promising instrument towards a better theoretical understanding are stellar atmospheres allowing for a consistent inclusion of the hydrodynamics. By coupling stellar and wind parameters and the inclusion of a detailed non-LTE radiative transfer, they allow us to go beneath the observable layers and study the onset of WR-type winds. Establishing larger sets of models, we were able to make ground-breaking progress by identifying trends with mass and metallicity that deviate significantly from present empirical descriptions. Our modelling efforts reveal a complex picture for WR-type winds with strong, non-linear dependencies. Besides covering metallicity and mass, we further identify surface hydrogen as an important ingredient to retain WR-type mass loss at lower metallicity. Here, we present a summary of recent insights on the nature and onset of WR-type winds in massive stars including the consequences for stellar evolution, remaining open questions, and current efforts to overcome them. △ Less

Submitted 1 September, 2022; originally announced September 2022.

Comments: 6 pages, 3 figures, to be published in the Proceedings of the International Astronomical Union for the IAU Symposium 361 "Massive Stars Near and Far" (eds. N. St-Louis, J.S. Vink, J. Mackey)

Phase-resolved spectroscopic analysis of the eclipsing black hole X-ray binary M33 X-7: System properties, accretion, and evolution

Authors: V. Ramachandran, L. M. Oskinova, W. -R. Hamann, A. A. C. Sander, H. Todt, D. Pauli, T. Shenar, J. M. Torrejón, K. A. Postnov, J. M. Blondin, E. Bozzo, R. Hainich, D. Massa

Abstract: M33 X-7 is the only known eclipsing black hole high mass X-ray binary. The system is reported to contain a very massive O supergiant donor and a massive black hole in a short orbit. The high X-ray luminosity and its location in the metal-poor galaxy M33 make it a unique laboratory for studying the winds of metal-poor donor stars with black hole companions and it helps us to understand the potentia… ▽ More M33 X-7 is the only known eclipsing black hole high mass X-ray binary. The system is reported to contain a very massive O supergiant donor and a massive black hole in a short orbit. The high X-ray luminosity and its location in the metal-poor galaxy M33 make it a unique laboratory for studying the winds of metal-poor donor stars with black hole companions and it helps us to understand the potential progenitors of black hole mergers. Using phase-resolved simultaneous HST- and XMM-Newton-observations, we traced the interaction of the stellar wind with the black hole. Our comprehensive spectroscopic investigation of the donor star (X-ray+UV+optical) yields new stellar and wind parameters for the system that differs significantly from previous estimates. In particular, the masses of the components are considerably reduced to 38 for the O-star donor and 11.4 for the black hole. The O giant is overfilling its Roche lobe and shows surface He enrichment. The donor shows a densely clumped wind with a mass-loss rate that matches theoretical predictions. We investigated the wind-driving contributions from different ions and the changes in the ionization structure due to X-ray illumination. Toward the black hole, the wind is strongly quenched due to strong X-ray illumination. For this system, the standard wind-fed accretion scenario alone cannot explain the observed X-ray luminosity, pointing toward an additional mass overflow, which is in line with our acceleration calculations. The X-ray photoionization creates an He II emission region emitting $10^{47}$ ph/s. We computed binary evolutionary tracks for the system using MESA. Currently, the system is transitioning toward an unstable mass transfer phase, resulting in a common envelope of the black hole and donor. Since the mass ratio is q~3.3 and the period is short, the system is unlikely to survive the common envelope, but will rather merge. △ Less

Submitted 16 August, 2022; originally announced August 2022.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 667, A77 (2022)

Stellar wind properties of the nearly complete sample of O stars in the low metallicity young star cluster NGC346 in the SMC galaxy

Authors: M. J. Rickard, R. Hainich, W. -R. Hamann, L. M. Oskinova, R. K. Prinja, V. Ramachandran, D. Pauli, H. Todt, A. C. C. Sander, T. Shenar, Y. -H. Chu, J. S. Gallagher III

Abstract: Massive stars are among the main cosmic engines driving the evolution of star-forming galaxies. Their powerful ionising radiation and stellar winds inject a large amount of energy in the interstellar medium. Furthermore, mass-loss ($\dot{M}$) through radiatively driven winds plays a key role in the evolution of massive stars. Even so, the wind mass-loss prescriptions used in stellar evolution mode… ▽ More Massive stars are among the main cosmic engines driving the evolution of star-forming galaxies. Their powerful ionising radiation and stellar winds inject a large amount of energy in the interstellar medium. Furthermore, mass-loss ($\dot{M}$) through radiatively driven winds plays a key role in the evolution of massive stars. Even so, the wind mass-loss prescriptions used in stellar evolution models, population synthesis, and stellar feedback models often disagree with mass-loss rates empirically measured from the UV spectra of low metallicity massive stars. The most massive young star cluster in the low metallicity Small Magellanic Cloud galaxy is NGC346. This cluster contains more than half of all O stars discovered in this galaxy so far. A similar age, metallicity ($Z$), and extinction, the O stars in the NGC346 cluster are uniquely suited for a comparative study of stellar winds in O stars of different subtypes. We aim to use a sample of O stars within NGC346 to study stellar winds at low metallicity. We mapped the central 1" of NGC346 with the long-slit UV observations performed by the Space Telescope Imaging Spectrograph (STIS) on board of the {\em Hubble Space Telescope} and complemented these new datasets with archival observations. Multi-epoch observations allowed for the detection of wind variability. The UV dataset was supplemented by optical spectroscopy and photometry. The resulting spectra were analysed using a non-local thermal equilibrium model atmosphere code (PoWR) to determine wind parameters and ionising fluxes. The effective mapping technique allowed us to obtain a mosaic of almost the full extent of the cluster and resolve stars in its core. Among hundreds of extracted stellar spectra, 21 belong to O stars. Nine of them are classified as O stars for the first time. We analyse, in detail, the UV spectra of 19 O stars... (continued) △ Less

Submitted 19 July, 2022; originally announced July 2022.

Comments: 33 pages, 40 figures

Journal ref: A&A 666, A189 (2022)

Mass-loss implementation and temperature evolution of very massive stars

Authors: Gautham N. Sabhahit, Jorick S. Vink, Erin R. Higgins, Andreas A. C. Sander

Abstract: Very massive stars (VMS) dominate the physics of young clusters due to their ionising radiation and extreme stellar winds. It is these winds that determine their lifepaths until expiration. Observations in the Arches cluster show that VMS all have similar temperatures. The VLT-Flames Tarantula survey analysed VMS in the 30 Dor region of the LMC also finding a narrow range of temperatures, albeit a… ▽ More Very massive stars (VMS) dominate the physics of young clusters due to their ionising radiation and extreme stellar winds. It is these winds that determine their lifepaths until expiration. Observations in the Arches cluster show that VMS all have similar temperatures. The VLT-Flames Tarantula survey analysed VMS in the 30 Dor region of the LMC also finding a narrow range of temperatures, albeit at higher values - likely a metallicity effect. Using MESA, we study the main-sequence evolution of VMS with a new mass-loss recipe that switches from optically-thin O-star winds to optically-thick Wolf-Rayet type winds through the model-independent transition mass-loss rate of Vink & Gräfener. We examine the temperature evolution of VMS with mass loss that scales with the luminosity-over-mass (L/M) ratio and the Eddington parameter ($Γ_{\rm e}$), assessing the relevance of the surface hydrogen (H) abundance which sets the number of free electrons. We present grids of VMS models at Galactic and LMC metallicity and compare our temperature predictions with empirical results. Models with a steep $Γ_{\rm e}$-dependence evolve horizontally in the Hertzsprung-Russel (HR) diagram at nearly constant luminosities, requiring a delicate and unlikely balance between envelope inflation and enhanced mass loss over the entire VMS mass range. By contrast, models with a steep L/M-dependent mass loss are shown to evolve vertically in the HR-diagram at nearly constant Teff, naturally reproducing the narrow range of observed temperatures, as well as the correct trend with metallicity. This distinct behavior of a steeply dropping luminosity is a self-regulatory mechanism that keeps temperatures constant during evolution in the HR-diagram. △ Less

Submitted 26 July, 2022; v1 submitted 18 May, 2022; originally announced May 2022.

Comments: Accepted for publication in MNRAS. 18 pages, 16 figures, 1 appendix

The origin and impact of Wolf-Rayet-type mass loss

Authors: Andreas A. C. Sander, Jorick S. Vink, Erin R. Higgins, Tomer Shenar, Wolf-Rainer Hamann, Helge Todt

Abstract: Classical Wolf-Rayet (WR) stars mark an important stage in the late evolution of massive stars. As hydrogen-poor massive stars, these objects have lost their outer layers, while still losing further mass through strong winds indicated by their prominent emission line spectra. Wolf-Rayet stars have been detected in a variety of different galaxies. Their strong winds are a major ingredient of stella… ▽ More Classical Wolf-Rayet (WR) stars mark an important stage in the late evolution of massive stars. As hydrogen-poor massive stars, these objects have lost their outer layers, while still losing further mass through strong winds indicated by their prominent emission line spectra. Wolf-Rayet stars have been detected in a variety of different galaxies. Their strong winds are a major ingredient of stellar evolution and population synthesis models. Yet, a coherent theoretical picture of their strong mass-loss is only starting to emerge. In particular, the occurrence of WR stars as a function of metallicity (Z) is still far from being understood. To uncover the nature of the complex and dense winds of Wolf-Rayet stars, we employ a new generation of model atmospheres including a consistent solution of the wind hydrodynamics in an expanding non-LTE situation. With this technique, we can dissect the ingredients driving the wind and predict the resulting mass-loss for hydrogen-depleted massive stars. Our modelling efforts reveal a complex picture with strong, non-linear dependencies on the luminosity-to-mass ratio and Z with a steep, but not totally abrupt onset for WR-type winds in helium stars. With our findings, we provide a theoretical motivation for a population of helium stars at low Z, which cannot be detected via WR-type spectral features. Our study of massive He-star atmosphere models yields the very first mass-loss recipe derived from first principles in this regime. Implementing our first findings in stellar evolution models, we demonstrate how traditional approaches tend to overpredict WR-type mass loss in the young Universe. △ Less

Submitted 9 February, 2022; originally announced February 2022.

Comments: 6 pages, 3 figures, to appear in the proceedings of IAUS 366 "The Origin of Outflows in Evolved Stars"

The Link between Hot and Cool Outflows

Authors: Jorick S. Vink, A. A. C. Sander, E. R. Higgins, G. N. Sabhahit

Abstract: The link between hot and cool stellar outflows is shown to be critical for correctly predicting the masses of the most massive black holes (BHs) below the so-called pair-instability supernova (PISN) mass gap. Gravitational Wave (GW) event 190521 allegedly hosted an "impossibly" heavy BH of 85 Solar Masses. Here we show how our increased knowledge of both metallicity Z and temperature dependent mas… ▽ More The link between hot and cool stellar outflows is shown to be critical for correctly predicting the masses of the most massive black holes (BHs) below the so-called pair-instability supernova (PISN) mass gap. Gravitational Wave (GW) event 190521 allegedly hosted an "impossibly" heavy BH of 85 Solar Masses. Here we show how our increased knowledge of both metallicity Z and temperature dependent mass loss is critical for our evolutionary scenario of a low-Z blue supergiant (BSG) progenitor of an initially approx 100 Solar Masses star to work. We show using MESA stellar evolution modelling experiments that as long as we can keep such stars above 8000 K such low-Z BSGs can avoid strong winds, and keep a very large envelope mass intact before core collapse. This naturally leads to the Cosmic Time dependent maximum BH function below the PISN gap. △ Less

Submitted 11 February, 2022; v1 submitted 28 January, 2022; originally announced January 2022.

Comments: 6 pages, 4 figures. Talk in IAUS 366

The earliest O-type eclipsing binary in the Small Magellanic Cloud, AzV 476: A comprehensive analysis reveals surprisingly low stellar masses

Authors: D. Pauli, L. M. Oskinova, W. -R. Hamann, V. Ramachandran, H. Todt, A. A. C. Sander, T. Shenar, M. Rickard, J. Maíz Apellániz, R. Prinja

Abstract: Massive stars at low metallicity are among the main feedback agents in the early Universe and in present-day star forming galaxies. When in binaries, these stars are potential progenitors of gravitational-wave events. Knowledge of stellar masses is a prerequisite to understanding evolution and feedback of low-metallicity massive stars. Using abundant spectroscopic and photometric measurements of a… ▽ More Massive stars at low metallicity are among the main feedback agents in the early Universe and in present-day star forming galaxies. When in binaries, these stars are potential progenitors of gravitational-wave events. Knowledge of stellar masses is a prerequisite to understanding evolution and feedback of low-metallicity massive stars. Using abundant spectroscopic and photometric measurements of an outstandingly bright eclipsing binary, we compare its dynamic, spectroscopic, and evolutionary mass estimates and develop a binary evolution scenario. We comprehensively studied the eclipsing binary system, AzV 476, in the Small Magellanic Cloud. The light curve and radial velocities were analyzed to obtain the orbital parameters. The photometric and spectroscopic data in the UV and optical were analyzed using the Potsdam Wolf-Rayet model atmospheres. The obtained results are interpreted using binary-evolution tracks. AzV 476 consists of an O4IV-III((f))p primary and an O9.5:Vn secondary. Both components have similar current masses (~20 M$_{\odot}$) obtained from both the orbital and spectroscopic analysis. The wind mass-loss rate of log($\dot{M}$/(M$_{\odot}$/yr))=-6.2 of the primary is a factor of ten higher than a recent empirical prescription for single O stars in the SMC. Only close-binary evolution with mass transfer can reproduce the current stellar and orbital parameters. The binary evolutionary model reveals that the primary has lost about half of its initial mass and is already core helium burning. Our comprehensive analysis of AzV 476 yields a consistent set of parameters and suggests previous case B mass transfer. The derived stellar masses agree within their uncertainties. The moderate masses of AzV 476 underline the scarcity of bright massive stars in the SMC. The core helium burning nature of the primary indicates that stripped stars might be hidden among OB-type populations. △ Less

Submitted 22 January, 2022; originally announced January 2022.

Comments: 18 pages + 5 pages appendix, 17 figures, 8 tables. Accepted for publication in A&A

Journal ref: A&A 659, A9 (2022)

A meeting at z~3: Young massive galaxies and an AGN within 30kpc of the luminous QSO LBQS0302-0019

Authors: B. Husemann, G. Worseck, F. Arrigoni Battaia, A. A. C. Sander, T. Shanks

Abstract: Contrary to expectations from scenarios of black hole growth driven by galaxy interactions and mergers, dual active galactic nuclei (AGN) with kiloparsec separations are rarely observed and are very difficult to identify, in particular at high redshifts (i.e. z>2). Focussing on the recently discovered dual AGN system LBQS 0302-0019 at z=3.29, we seek to identify further group members in its enviro… ▽ More Contrary to expectations from scenarios of black hole growth driven by galaxy interactions and mergers, dual active galactic nuclei (AGN) with kiloparsec separations are rarely observed and are very difficult to identify, in particular at high redshifts (i.e. z>2). Focussing on the recently discovered dual AGN system LBQS 0302-0019 at z=3.29, we seek to identify further group members in its environment and to understand their formation history through deep high-angular-resolution imaging. We present deep Hubble Space Telescope (HST) Wide-field Camera 3 near-infrared imaging of LBQS 0302-0019. In combination with ground-based VLT/HAWK-I imaging, we infer accurate sizes, colours, ages, and stellar masses of companion galaxies. We clearly detect four companion objects close to LBQS 0302-0019 that also have faint signatures in the ground-based images. We constrain light-weighted ages and masses for the two most prominent companions, Jil1 and Jil2, to $t_\star=252_{-109}^{+222}$Myr with $\log(M_\star/[\mathrm{M}_\odot])= 11.2_{-0.1}^{+0.3}$ and $t_{\star}=19_{-14}^{+74}$Myr with $\log(M_\star/[\mathrm{M}_\odot])= 9.4_{-0.4}^{+0.9}$, respectively. The HST data also show that the obscured AGN, previously identified by strong nebular HeII emission, is associated with the young massive companion Jil2. Because very massive stars of the starburst cannot be solely responsible for the HeII emission, we strengthen our initial conclusion that Jil2 has been hosting an AGN. If the young starburst of Jil2 had been accompanied by sustained black hole growth, Jil2 may have contributed HeII-ionising flux to create the large HeII Ly$α$ proximity zone around LBQS 0302-0019. Hence, the duration of the current luminous AGN episode of LBQS 0302-0019 may have been overestimated. △ Less

Submitted 22 July, 2021; originally announced July 2021.

Comments: 7 pages, 3 figures and 1 table. Accepted for publication in A&A

Journal ref: A&A 653, A122 (2021)

Superadiabaticity and the metallicity independence of the Humphreys-Davidson limit

Authors: Gautham N. Sabhahit, Jorick S. Vink, Erin R. Higgins, Andreas A. C. Sander

Abstract: The Humphreys-Davidson (HD) limit sets the boundary between evolutionary channels of massive stars that either end their lives as red supergiants (RSGs) or as the hotter blue supergiants (BSGs) and Wolf-Rayet stars. Mixing in the envelopes of massive stars close to their Eddington limit is crucial for investigating the upper luminosity limit of the coolest supergiants. We study the effects of exce… ▽ More The Humphreys-Davidson (HD) limit sets the boundary between evolutionary channels of massive stars that either end their lives as red supergiants (RSGs) or as the hotter blue supergiants (BSGs) and Wolf-Rayet stars. Mixing in the envelopes of massive stars close to their Eddington limit is crucial for investigating the upper luminosity limit of the coolest supergiants. We study the effects of excess mixing in superadiabatic layers that are dominated by radiation pressure, and we critically investigate the effects of mixing and mass loss on the evolution of RSGs with log (Teff/K) < 3.68 - as a function of metallicity. Using MESA, we produce grids of massive star models at three metallicities: Galactic (Zsol), LMC (1/2 Zsol) and SMC (1/5 Zsol), with both high and low amounts of overshooting to study the upper luminosity limit of RSGs. We systematically study the effects of excess mixing in the superadiabatic layers of post-main sequence massive stars, overshooting above the hydrogen core and yellow supergiant (YSG) mass-loss rates on the fraction of core helium burning time spent as a RSG. We find that the excess mixing in the superadiabatic layers is stronger at lower metallicities, as it depends on the opacities in the hydrogen bump at log (Teff/K) ~ 4, which become more pronounced at lower metallicity. This shifts the cutoff luminosities to lower values at lower metallicities, thus balancing the first-order effect of mass loss. The opposing effects of mass loss and excess envelope mixing during post-main sequence evolution of stars with higher overshooting potentially results in a metallicity-independent upper luminosity limit. △ Less

Submitted 9 August, 2021; v1 submitted 5 July, 2021; originally announced July 2021.

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

Evolution of Wolf-Rayet stars as black hole progenitors

Authors: Erin R. Higgins, Andreas A. C. Sander, Jorick S. Vink, Raphael Hirschi

Abstract: Evolved Wolf-Rayet stars form a key aspect of massive star evolution, and their strong outflows determine their final fates. In this study, we calculate grids of stellar models for a wide range of initial masses at five metallicities (ranging from solar down to just 2% solar). We compare a recent hydrodynamically-consistent wind prescription with two earlier frequently-used wind recipes in stellar… ▽ More Evolved Wolf-Rayet stars form a key aspect of massive star evolution, and their strong outflows determine their final fates. In this study, we calculate grids of stellar models for a wide range of initial masses at five metallicities (ranging from solar down to just 2% solar). We compare a recent hydrodynamically-consistent wind prescription with two earlier frequently-used wind recipes in stellar evolution and population synthesis modelling, and we present the ranges of maximum final masses at core He-exhaustion for each wind prescription and metallicity Z. Our model grids reveal qualitative differences in mass-loss behaviour of the wind prescriptions in terms of "convergence". Using the prescription from Nugis & Lamers the maximum stellar black hole is found to converge to a value of 20-30Msun, independent of host metallicity, however when utilising the new physically-motivated prescription from Sander & Vink there is no convergence to a maximum black hole mass value. The final mass is simply larger for larger initial He-star mass, which implies that the upper black hole limit for He-stars below the pair-instability gap is set by prior evolution with mass loss, or the pair instability itself. Quantitatively, we find the critical Z for pair-instability (Z_PI) to be as high as 50% Zsolar, corresponding to the host metallicity of the LMC. Moreover, while the Nugis & Lamers prescription would not predict any black holes above the approx 130Msun pair-instability limit, with Sander & Vink winds included, we demonstrate a potential channel for very massive helium stars to form such massive black holes at ~2% Zsolar or below. △ Less

Submitted 25 May, 2021; originally announced May 2021.

Comments: Accepted for publication in MNRAS. 17 pages, 16 figures

Revisiting the archetypical wind accretor Vela X-1 in depth -- A case study of a well-known X-ray binary and the limits of our knowledge

Authors: Peter Kretschmar, Ileyk El Mellah, Silvia Martínez-Núñez, Felix Fürst, Victoria Grinberg, Andreas A. C. Sander, Jakob van den Eijnden, Nathalie Degenaar, Jesús Maíz-Apellániz, Francisco Jiménez Esteban, Mercedes Ramos-Lerate, Enrique Utrilla

Abstract: Context: Vela X-1 is one of the best studied X-ray binaries. Frequently though, specific values for its parameters have been used in subsequent studies without considering alternatives. Aims: We aim to provide a robust compilation and synthesis of the accumulated knowledge about Vela X-1 as a solid baseline for future studies and identify specific avenues of possible future research. Methods: We e… ▽ More Context: Vela X-1 is one of the best studied X-ray binaries. Frequently though, specific values for its parameters have been used in subsequent studies without considering alternatives. Aims: We aim to provide a robust compilation and synthesis of the accumulated knowledge about Vela X-1 as a solid baseline for future studies and identify specific avenues of possible future research. Methods: We explore the literature for Vela X-1 and on modelling efforts, describing the evolution of the system knowledge. We also add information derived from public data, especially the Gaia EDR3 release. Results: We update the distance to Vela X-1, the spectral classification for HD 77518 and find that the supergiant may be very close to filling its Roche lobe. Constraints on the clumpiness of the stellar wind have improved. The orbit is very well determined, but the uncertain inclination limits information on the neutron star mass. Estimates for the stellar wind have evolved towards lower velocities, supporting the idea of transient wind-captured disks around the neutron star. Hydrodynamic models and observations are consistent with an accretion wake trailing the neutron star. Conclusions: Vela X-1 is an excellent laboratory, but a lot of room remains to improve. Well-coordinated multi-wavelength observations and campaigns addressing the intrinsic variability are required. New opportunities will arise through new instrumentation. Models of the stellar wind should account for the orbital eccentricity and the non-spherical shape of HD 77581. Realistic multi-dimensional models of radiative transfer in the UV and X-rays are needed, but remain very challenging. Improved MHD models covering a wide range of scales would be required to improve understanding of the plasma-magnetosphere coupling. A full characterization of the accretion column remains another open challenge. (Abbreviated for arXiv) △ Less

Submitted 11 May, 2021; v1 submitted 27 April, 2021; originally announced April 2021.

Comments: 50 pages, 20 figures, 12 tables. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 652, A95 (2021)

Metallicity-dependent wind parameter predictions for OB stars

Authors: Jorick S. Vink, Andreas A. C. Sander

Abstract: Mass-loss rates and terminal wind velocities are key parameters that determine the kinetic wind energy and momenta of massive stars. Furthermore, accurate mass-loss rates determine the mass and rotational velocity evolution of mass stars, and their fates as neutron stars and black holes in function of metallicity (Z). Here we update our Monte Carlo mass-loss Recipe with new dynamically-consistent… ▽ More Mass-loss rates and terminal wind velocities are key parameters that determine the kinetic wind energy and momenta of massive stars. Furthermore, accurate mass-loss rates determine the mass and rotational velocity evolution of mass stars, and their fates as neutron stars and black holes in function of metallicity (Z). Here we update our Monte Carlo mass-loss Recipe with new dynamically-consistent computations of the terminal wind velocity -- as a function of Z. These predictions are particularly timely as the HST ULLYSES project will observe ultraviolet spectra with blue-shifted P Cygni lines of hundreds of massive stars in the low-Z Large and Small Magellanic Clouds, as well as sub-SMC metallicity hosts. Around 35 000 K, we uncover a weak-wind "dip" and we present diagnostics to investigate its physics with ULLYSES and X-Shooter data. We discuss how the dip may provide important information on wind-driving physics, and how this is of key relevance towards finding a new gold-standard for OB star mass-loss rates. For B supergiants below the Fe IV to III bi-stability jump, the terminal velocity is found to be independent of Z and M, while the mass-loss rate still varies as $\dot{M} \propto Z^{0.85}$. For O-type stars above the bi-stability jump we find a terminal-velocity dependence of $v_{\infty} \propto Z^{0.19}$ and the Z-dependence of the mass-loss rate is found to be as shallow as $\dot{M} \propto Z^{0.42}$, implying that to reproduce the `heavy' black holes from LIGO/VIRGO, the `low Z' requirement becomes even more stringent than was previously anticipated. △ Less

Submitted 9 April, 2021; v1 submitted 23 March, 2021; originally announced March 2021.

Comments: MNRAS Accepted. 12 pages, 13 figures. Notes on IDL/Python Mass-loss Recipe Usage added to the Appendix

Conditions in the WR 140 wind-collision region revealed by the 1.083-micron He I line profile

Authors: Peredur M. Williams, Watson P. Varricatt, André-Nicolas Chené, Michael F. Corcoran, Ted R. Gull, Kenji Hamaguchi, Anthony F. J. Moffat, Andrew M. T. Pollock, Noel D. Richardson, Christopher M. P. Russell, Andreas A. C. Sander, Ian R. Stevens, Gerd Weigelt

Abstract: We present spectroscopy of the P~Cygni profile of the 1.083-micron He I line in the WC7 + O5 colliding-wind binary (CWB) WR 140 (HD 193793), observed in 2008, before its periastron passage in 2009, and in 2016-17, spanning the subsequent periastron passage. Both absorption and emission components showed strong variations. The variation of the absorption component as the O5 star was occulted by the… ▽ More We present spectroscopy of the P~Cygni profile of the 1.083-micron He I line in the WC7 + O5 colliding-wind binary (CWB) WR 140 (HD 193793), observed in 2008, before its periastron passage in 2009, and in 2016-17, spanning the subsequent periastron passage. Both absorption and emission components showed strong variations. The variation of the absorption component as the O5 star was occulted by the wind-collision region (WCR) sets a tight constraint on its geometry. While the sightline to the O5 star traversed the WCR, the strength and breadth of the absorption component varied significantly on time-scales of days. An emission sub-peak was observed on all our profiles. The variation of its radial velocity with orbital phase was shown to be consistent with formation in the WCR as it swung round the stars in their orbit. Modelling the profile gave a measure of the extent of the sub-peak forming region. In the phase range 0.93-0.99, the flux in the sub-peak increased steadily, approximately inversely proportionally to the stellar separation, indicating that the shocked gas in the WCR where the line was formed was adiabatic. After periastron, the sub-peak flux was anomalously strong and varied rapidly, suggesting formation in clumps down-stream in the WCR. For most of the time, its flux exceeded the 2-10-keV X-ray emission, showing it to be a significant coolant of the shocked wind. △ Less

Submitted 18 February, 2021; originally announced February 2021.

Comments: 19 pages, 14 figures, accepted for publication in the MNRAS

The First Dynamical Mass Determination of a Nitrogen-rich Wolf-Rayet Star using a Combined Visual and Spectroscopic Orbit

Authors: Noel D. Richardson, Laura Lee, Gail Schaefer, Tomer Shenar, Andreas A. C. Sander, Grant M. Hill, Andrew G. Fullard, John D. Monnier, Narsireddy Anugu, Claire L Davies, Tyler Gardner, Cyprien Lanthermann, Stefan Kraus, Benjamin R. Setterholm

Abstract: We present the first visual orbit for the nitrogen-rich Wolf-Rayet binary, WR 133 (WN5o + O9I) based on observations made with the CHARA Array and the MIRC-X combiner. This orbit represents the first visual orbit for a WN star and only the third Wolf-Rayet star with a visual orbit. The orbit has a period of 112.8 d, a moderate eccentricity of 0.36, and a separation of $a$= 0.79 mas on the sky. We… ▽ More We present the first visual orbit for the nitrogen-rich Wolf-Rayet binary, WR 133 (WN5o + O9I) based on observations made with the CHARA Array and the MIRC-X combiner. This orbit represents the first visual orbit for a WN star and only the third Wolf-Rayet star with a visual orbit. The orbit has a period of 112.8 d, a moderate eccentricity of 0.36, and a separation of $a$= 0.79 mas on the sky. We combine the visual orbit with an SB2 orbit and Gaia parallax to find that the derived masses of the component stars are $M_{\rm WR}$ = $9.3\pm1.6 M_\odot$ and $M_{\rm O}$ = $22.6\pm 3.2 M_\odot$, with the large errors owing to the nearly face-on geometry of the system combined with errors in the spectroscopic parameters. We also derive an orbital parallax that is identical to the {\it Gaia}-determined distance. We present a preliminary spectral analysis and atmosphere models of the component stars, and find the mass-loss rate in agreement with polarization variability and our orbit. However, the derived masses are low compared to the spectral types and spectral model. Given the close binary nature, we suspect that WR 133 should have formed through binary interactions, and represents an ideal target for testing evolutionary models given its membership in the cluster NGC 6871. △ Less

Submitted 11 January, 2021; originally announced January 2021.

Comments: Accepted to ApJ Letters

Wind-envelope interaction as the origin of the slow cyclic brightness variations of luminous blue variables

Authors: Luca Grassitelli, Norbert Langer, Jonathan Mackey, Goetz Graefener, Nathan Grin, Andreas Sander, Jorick Vink

Abstract: Luminous blue variables (LBVs) are hot, very luminous massive stars displaying large quasi-periodic variations in brightness, radius,and photospheric temperature, on timescales of years to decades. The physical origin of this variability, called S Doradus cycle after its prototype, has remained elusive. Here, we study the feedback of stellar wind mass-loss on the envelope structure in stars near t… ▽ More Luminous blue variables (LBVs) are hot, very luminous massive stars displaying large quasi-periodic variations in brightness, radius,and photospheric temperature, on timescales of years to decades. The physical origin of this variability, called S Doradus cycle after its prototype, has remained elusive. Here, we study the feedback of stellar wind mass-loss on the envelope structure in stars near the Eddington limit. We perform a time-dependent hydrodynamic stellar evolutionary calculation, applying a stellar wind mass-loss prescription with a temperature-dependence inspired by the predicted systematic increase in mass-loss rates below 25 kK. We find that when the wind mass-loss rate crosses a well-defined threshold, a discontinuous change in the wind base conditions leads to a restructuring of the stellar envelope. The induced drastic radius and temperature changes, which occur on the thermal timescale of the inflated envelope, impose in turn mass-loss variations that reverse the initial changes, leading to a cycle that lacks a stationary equilibrium configuration. Our proof-of-concept model broadly reproduces the typical observational phenomenology of the S Doradus variability. We identify three key physical ingredients needed to trigger the instability: inflated envelopes in close proximity to the Eddington limit, a temperature range where decreasing opacities do not lead to an accelerating outflow, and a mass-loss rate that increases with decreasing temperature, crossing a critical threshold value within this temperature range. Our scenario and model provide testable predictions, and open the door for a consistent theoretical treatment of the LBV phase in stellar evolution, with consequences for their further evolution as single stars or in binary systems. △ Less

Submitted 3 December, 2020; v1 submitted 30 November, 2020; originally announced December 2020.

Comments: A&A Accepted, 14 pages, many colorful figures

Report number: AA/2020/38298

Journal ref: A&A 647, A99 (2021)

Maximum Black Hole mass across Cosmic Time

Authors: Jorick S. Vink, Erin R. Higgins, Andreas A. C. Sander, Gautham N. Sabhahit

Abstract: At the end of its life, a very massive star is expected to collapse into a black hole. The recent detection of an 85 Msun black hole from the gravitational wave event GW 190521 appears to present a fundamental problem as to how such heavy black holes exist above the approximately 50 Msun pair-instability limit where stars are expected to be blown to pieces with no remnant left. Using MESA, we show… ▽ More At the end of its life, a very massive star is expected to collapse into a black hole. The recent detection of an 85 Msun black hole from the gravitational wave event GW 190521 appears to present a fundamental problem as to how such heavy black holes exist above the approximately 50 Msun pair-instability limit where stars are expected to be blown to pieces with no remnant left. Using MESA, we show that for stellar models with non-extreme assumptions, 90..100 Msun stars at reduced metallicity (Z/Zsun < 0.1) can produce blue supergiant progenitors with core masses sufficiently small to remain below the fundamental pair-instability limit, yet at the same time lose an amount of mass via stellar winds that is small enough to end up in the range of an "impossible" 85 Msun black hole. The two key points are the proper consideration of core overshooting and stellar wind physics with an improved scaling of mass loss with iron (Fe) contents characteristic for the host galaxy metallicity. Our modelling provides a robust scenario that not only doubles the maximum black hole mass set by pair instability, but also allows us to probe the maximum stellar black hole mass as a function of metallicity and Cosmic time in a physically sound framework. △ Less

Submitted 16 March, 2021; v1 submitted 19 October, 2020; originally announced October 2020.

Comments: MNRAS Accepted. 9 pages, 3 figures. Prior Submission in September 2020 to Nature (see Research Square)

On the nature of massive helium star winds and Wolf-Rayet-type mass loss

Authors: Andreas A. C. Sander, Jorick S. Vink

Abstract: The mass-loss rates of massive helium stars are one of the major uncertainties in modern astrophysics. Regardless of whether they were stripped by a binary companion or managed to peel off their outer layers by themselves, the influence and final fate of helium stars -- in particular the resulting black hole mass -- highly depends on their wind mass loss as stripped-envelope objects. While empiric… ▽ More The mass-loss rates of massive helium stars are one of the major uncertainties in modern astrophysics. Regardless of whether they were stripped by a binary companion or managed to peel off their outer layers by themselves, the influence and final fate of helium stars -- in particular the resulting black hole mass -- highly depends on their wind mass loss as stripped-envelope objects. While empirical mass-loss constraints for massive helium stars have improved over the last decades, the resulting recipes are limited to metallicities with the observational ability to sufficiently resolve individual stars. Yet, theoretical efforts have been hampered by the complexity of Wolf-Rayet (WR) winds arising from the more massive helium stars. In an unprecedented effort, we calculate next-generation stellar atmosphere models resembling massive helium main sequence stars with Fe-bump driven winds up to $500\,M_\odot$ over a wide metallicity range between $2.0$ and $0.02\,Z_\odot$. We uncover a complex $Γ_\text{e}$-dependency of WR-type winds and their metallicity-dependent breakdown. The latter can be related to the onset of multiple scattering, requiring higher $L/M$-ratios at lower metallicity. Based on our findings, we derive the first ever theoretically-motivated mass-loss recipe for massive helium stars. We also provide estimates for LyC and He II ionizing fluxes, finding stripped helium stars to contribute considerably at low metallicity. In sharp contrast to OB-star winds, the mass loss for helium stars scales with the terminal velocity. While limited to the helium main sequence, our study marks a major step towards a better theoretical understanding of helium star evolution. △ Less

Submitted 3 September, 2020; originally announced September 2020.

Comments: 20 pages, 22+9 figures, accepted for publication in MNRAS; a python routine implementing the resulting mass-loss recipes can be found at https://armagh.space/asander

Optical and infrared study of the obscured B[e] supergiant high-mass X-ray binary IGR J16318-4848

Authors: Francis Fortin, Sylvain Chaty, Andreas Sander

Abstract: The supergiant High Mass X-ray Binary IGR J16318-4848 was the first source detected by the INTEGRAL satellite in 2003 and distinguishes itself by its high intrinsic absorption and B[e] phenomenon. It is the perfect candidate to study both binary interaction and the environment of supergiant B[e] stars. This study targets the local properties of IGR J16318-4848. We aim to clarify the geometry of th… ▽ More The supergiant High Mass X-ray Binary IGR J16318-4848 was the first source detected by the INTEGRAL satellite in 2003 and distinguishes itself by its high intrinsic absorption and B[e] phenomenon. It is the perfect candidate to study both binary interaction and the environment of supergiant B[e] stars. This study targets the local properties of IGR J16318-4848. We aim to clarify the geometry of this system, and distinguish different key emitting regions in the binary. We provide optical to near-infrared spectra from VLT/X-Shooter and analyse both fine structures of the lines and the broadband spectral energy distribution by adding archival mid-infrared Spitzer and Herschel data. We also performed a stellar atmosphere and wind modeling of the optical to near-infrared spectrum using the PoWR code. We determine the contribution of the irradiated inner edge of the dusty circumbinary disk, derive the velocity of an equatorial stellar wind, and suggest the compact object orbits within the cavity between the star and the disk. We report on flat-topped lines originating from a spherically symetric disk wind, along with the first detection of what is likely the polar component of the stellar wind. Stellar atmosphere and wind modeling shows that the central star may have a helium-enhanced atmosphere, likely because of its intense wind sheding part of its hydrogen envelope. Finally we compare the properties of IGR J16318-4848 with a similar source, CI Cameleopardis. △ Less

Submitted 10 April, 2020; originally announced April 2020.

Comments: ApJ in press, 20 pages, 16 figures

Why binary interaction does not necessarily dominate the formation of Wolf-Rayet stars at low metallicity

Authors: T. Shenar, A. Gilkis, J. S. Vink, H. Sana, A. A. C. Sander

Abstract: Classical Wolf-Rayet (WR) stars are massive, hydrogen depleted, post main-sequence stars that exhibit emission-line dominated spectra. For a given metallicity Z, stars exceeding a certain initial mass M_single(Z) can reach the WR phase through intrinsic mass-loss (single-star channel). Stars of lower masses can reach the WR phase via binary mass transfer (binary channel). It is commonly assumed th… ▽ More Classical Wolf-Rayet (WR) stars are massive, hydrogen depleted, post main-sequence stars that exhibit emission-line dominated spectra. For a given metallicity Z, stars exceeding a certain initial mass M_single(Z) can reach the WR phase through intrinsic mass-loss (single-star channel). Stars of lower masses can reach the WR phase via binary mass transfer (binary channel). It is commonly assumed that the binary channel dominates the formation of WR stars in environments with low Z such as the SMC and LMC. However, their reported WR binary fractions of 30-40% are comparable to that of the Galaxy. Here, we explain this apparent contradiction by considering the minimum initial mass M_spec(Z) needed for the stripped product to appear as a WR star. We calibrate M_spec(Z) using the lowest-luminosity WR stars in the Clouds and the Galaxy. A range of M_single(Z) values are explored using various evolution codes. We estimate the additional contribution of the binary channel by considering the interval [M_spec(Z), M_single(Z)], which characterises the initial-mass range in which binaries can form additional WR stars. Results: The WR-phenomenon ceases below luminosities of logL = 4.9, 5.25, and 5.6 [Lsun] in the Galaxy, LMC, and SMC, which translates to He-star masses of 7.5, 11, 17 Msun and initial masses of M_spec = 18, 23, 37 Msun. Stripped stars with lower initial masses in the respective galaxies would tend to not appear as WR stars. M_single lies in the range 20-30, 30-60, and > 40 Msun for the Galaxy, LMC, and SMC. We find that that the additional contribution of the binary channel is a non-trivial function of Z that cannot be conclusively claimed to be monotonically increasing with decreasing Z. Hence, one should not a-priori expect that binary interactions become increasingly important in forming WR stars at low Z, or that the WR binary fraction grows with decreasing Z. △ Less

Submitted 13 January, 2020; originally announced January 2020.

Comments: 10 pages, 5 figures, accepted to A&A on the 12.01.2020. Comments welcome!

Journal ref: A&A 634, A79 (2020)

The stellar and wind parameters of six prototypical HMXBs and their evolutionary status

Authors: R. Hainich, L. M. Oskinova, J. M. Torrejón, F. Fuerst, A. Bodaghee, T. Shenar, A. A. C. Sander, H. Todt, K. Spetzer, W. -R. Hamann

Abstract: High-mass X-ray binaries (HMXBs) are exceptional astrophysical laboratories that offer a rare glimpse into the physical processes that govern accretion on compact objects, massive-star winds, and stellar evolution. In a subset of the HMXBs, the compact objects accrete matter solely from winds of massive donor stars. These so-called wind-fed HMXBs are divided in persistent HMXBs and supergiant fast… ▽ More High-mass X-ray binaries (HMXBs) are exceptional astrophysical laboratories that offer a rare glimpse into the physical processes that govern accretion on compact objects, massive-star winds, and stellar evolution. In a subset of the HMXBs, the compact objects accrete matter solely from winds of massive donor stars. These so-called wind-fed HMXBs are divided in persistent HMXBs and supergiant fast X-ray transients (SFXTs) according to their X-ray properties. While it has been suggested that this dichotomy depends on the characteristics of stellar winds, they have been poorly studied. With this investigation, we aim to remedy this situation by systematically analyzing donor stars of wind-fed HMXBs that are observable in the UV, concentrating on those with neutron star (NS) companions. We obtained Swift X-ray data, HST UV spectra, and additional optical spectra for all our targets. Our multi-wavelength approach allows us to provide stellar and wind parameters for six donor stars (four wind-fed systems and two OBe X-ray binaries). The wind properties are in line with the predictions of the line-driven wind theory. Three of the donor stars are in an advanced evolutionary stage, while for some of the stars, the abundance pattern indicates that processed material might have been accreted. When passing by the NS in its tight orbit, the donor star wind has not yet reached its terminal velocity but it is still significantly slower; its speed is comparable with the orbital velocity of the NS companion. There are no systematic differences between the two types of wind-fed HMXBs (persistent versus transients) with respect to the donor stars. For the SFXTs in our sample, the orbital eccentricity is decisive for their transient X-ray nature. Based on the orbital parameters and the further evolution of the donor stars, the investigated HMXBs will presumably form Thorne-Żytkow objects in the future. △ Less

Submitted 8 January, 2020; originally announced January 2020.

Comments: 44 pages, 13 figures, 15 tables; accepted for publication in A&A

Journal ref: A&A 634, A49 (2020)

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)