Search

Very massive stars with masses over 100 Msun are key objects in the Universe for our understanding of chemical and energetic feedback in the Universe, but their evolution and fate are almost entirely determined by their wind mass loss. We aim to determine the mass of the most massive star known in the Local Group R136a1. For this we compute the first hydrodynamically consistent non-local thermodynamical equilibrium atmosphere models for both R136a1 (WN5h) as well as the binary system R144 (WN5/6h+WN6/7h) in the Tarantula nebula. Using the Potsdam Wolf-Rayet code, we simultaneously empirically derive and theoretically predict mass-loss rates and wind velocities. By fitting synthetic spectra derived from these models to multi-wavelength observations, we constrain the stellar and wind properties of R144 and R136a1. We first determine the clumping stratification required by our hydro-models to fit the spectra of R144 by using the available dynamical mass estimates for the two components. We then utilise this clumping stratification in hydrodynamic models of R136a1 and estimate a mass of $M_\mathrm{Hydro}$ of 233 Msun. Remarkably, the estimated mass is close to and entirely consistent with chemical homogeneous mass relations. This present-day mass of 233 Msun provides a lower limit to the initial stellar mass, that could be far higher due to previous wind mass loss., Comment: Accepted for publication in A&A, 20 pages, 18 figures, 4 tables

Using Gaia DR3 we derive new distances and luminosities for a sample of Galactic B supergiants which were thought to be post main-sequence (MS) objects from their HR diagram location beyond the terminal-age MS (TAMS). When applying the newer Gaia distances in addition to enhanced amounts of core-boundary mixing, aka convective overshooting, we show that these Galactic B supergiants are likely enclosed within the MS band, indicating an evolutionary stage of steady core hydrogen burning. We discuss the importance of considering enhanced overshooting and how vectors in the mass-luminosity plane (ML-plane) can be used to disentangle the effects of wind mass loss from interior mixing. We finish with the key message that any proposed solution to the BSG problem should consider not only an explanation for the sheer number of B supergiants inside the Hertzsprung gap, but should at the same time also account for the steep drop in rotation rates identified at spectral type B1 -- corresponding to an effective temperature of $\sim$21 kK, and for which two distinct families of solutions have been proposed., Comment: 10 pages, 3 figures, Comments welcome, Review in Circumstellar Matter in Hot Star Systems (for Galaxies)

We present the discovery of twelve metal-poor and distant pre-main sequence (PMS) candidates in the dwarf irregular galaxy Wolf-Lundmark-Melotte (WLM)~968 kpc away, at a present-day metallicity of [Fe/H] around $-$0.9. These candidates have masses between 1.25-5 M$_{\odot}$, with ages <10 Myr, and exhibit significant near-infrared excesses at 2.5 and 4.3 $\mu$m. They are concentrated within a cluster roughly 10 pc (2'') across situated in the HII region [HM95]-9. These are the most distant and metal-poor PMS stars known, and can offer new quantitative insights into star formation at low-metallicities., Comment: Accepted for publication in ApJ Letters. Comments welcome. 6 pages, 4 figures, and 1 table

We use the inertial-inflow model of massive star formation to describe the formation of globular clusters (GCs) in turbulent molecular clouds. A key aspect of this model is that the maximum stellar mass scales linearly with cloud mass, such that extremely massive stars (EMSs, $10^{3-4}\,{\rm M}_\odot$) form in massive GCs ($\gtrsim10^5\,M_\odot$). The total wind mass loss is dominated by accreting EMSs (aEMSs), whose wind mass-loss rates become comparable to their accretion rates ($\gtrsim10^{-2}\,{\rm M}_\odot\,{\rm yr}^{-1}$). These winds pollute the intra-cluster medium with hot-hydrogen burning yields during GC formation. We propose a parameterised model for the evolution of the stellar mass function during GC formation ($\sim 1-2\,{\rm Myr}$), accounting for gas inflow, wind mass loss and mixing of aEMS yields with pristine gas that has initial proto-GC abundances. Low-mass stars ($\lesssim1\,{\rm M}_\odot$) form continuously from this mixed gas and their abundances resemble observed abundance trends with GC mass and metallicity, specifically: (i) the helium spread in a typical GC is small ($\Delta Y \simeq 0.01$) and increases with GC mass; (ii) the fraction of polluted stars increases with GC mass and metallicity; (iii) the extent of the Mg-Al anticorrelations is more pronounced in metal-poor and massive GCs. We conclude that GCs formed with a population of EMSs, that may leave behind black holes with masses above the pair-instability gap ($\gtrsim120\,{\rm M}_\odot$) and that the nitrogen-rich galaxies found by the James Webb Space Telescope ({\it JWST}) are in the first phase of galaxy evolution, when a large fraction of star formation occurs in the form of GCs., Comment: 28 pages, 19 figures, submitted to MNRAS. New idea for addressing the "globular cluster multiple population problem". Comments welcome

While stellar expansion after core-hydrogen exhaustion related to thermal imbalance has been documented for decades, the physical phenomenon of stellar inflation that occurs close to the Eddington limit has only come to the fore in recent years. We aim to elucidate the differences between these physical mechanisms for stellar radius enlargement, especially as additional terms such as `bloated' and `puffed-up' stars have been introduced in the recent massive star literature. We employ single and binary star MESA structure and evolution models for both constant mass, as well as models allowing for the mass to change, due to winds or binary interaction. We find cases that were previously attributed to stellar inflation in fact to be due to stellar expansion. We also highlight that while the opposite effect of expansion is contraction, the removal of an inflated zone should not be referred to as contraction but {\it deflation} as the star is still in thermal balance., Comment: Accepted for publication in A&A, 19 pages, 14 figures

Recent studies of high-redshift galaxies using JWST, such as GN-z11 revealed highly elevated levels of nitrogen (N). This phenomenon extends to gravitationally-lensed galaxies like the Sunburst Arc at z = 2.37, as well as to globular clusters (GCs). We propose that this originates from the presence of very massive stars (VMSs) with masses ranging from 100 to 1000\,\Msun. The He {\sc ii} observed in the Sunburst Arc could also stem from the disproportionately large contribution of VMSs. We build an entirely new Framework for massive star evolution which is no longer set by Dutch or other mass-loss "recipes" but which take the physics of $\Gamma$ or $L/M$-dependent winds into account. We discuss the mass-loss kink and the transition mass-loss rate between optically thin and thick winds, before we study the evaporative mass-loss history of VMSs. Our novel evolution models exhibit vertical evolution in the HR-diagram from the zero-age main sequence due to a self-regulatory effect driven by their wind-dominated nature, and we discuss what wind physics sets the stellar upper-mass limit. Our estimate for the Sunburst Arc in Vink (2023) suggests that the significant amounts of N found in star-forming galaxies likely arise from VMSs. We evaluate the strengths and weaknesses of previous hypotheses, including fast rotating massive stars and supermassive stars (SMSs), and we conclude that only our VMS model satisfies the relevant criteria. Finally, we advocate for the inclusion of VMSs in population synthesis and chemical evolution models, emphasizing the need for a self-consistent wind approach, which currently does not exist. Even minor inaccuracies in mass-loss rates dramatically impact the stellar evolution of VMS, as well as their ionizing and chemical feedback., Comment: 6 pages, 4 figures, proceedings IAUS 391: "The first chapters of our cosmic history with JWST"

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., Comment: Accepted for publication in MNRAS. 17 pages, 16 figures

We analyse the current knowledge and uncertainties in detailed stellar evolution and wind modelling to evaluate the mass of the most massive stellar black hole (BH) at solar metallicity. Contrary to common expectations that it is the most massive stars that produce the most massive BHs, we find that the maximum $M_{\rm BH}^{\rm Max} \simeq 30 \pm 10\,M_{\odot}$ is found in the canonical intermediate range between $M_{\rm ZAMS} \simeq 30$ and $50\,M_{\odot}$ instead. The prime reason for this seemingly counter-intuitive finding is that very massive stars (VMS) have increasingly high mass-loss rates, leading to substantial mass evaporation before they expire as stars, ending as lighter BHs than their canonical O-star counterparts., Comment: A&A Letters, accepted. 2 Figures

Stellar winds form an integral part of astronomy. The solar wind affects Earth's magnetosphere, while the winds of hot massive stars are highly relevant for galactic feedback through their mechanical wind energy. In different parts of the stellar HR diagram different forces dominate. On the hot side of the HRD radiative forces on ionised gas particles are active, while on the cool side molecular and dust opacities take over. Moreover, due to the convective envelopes, alternative physical ingredients may start to dominate. I will describe the basic equation of motion and give a few examples, mostly focusing on the winds from massive stars. Here mass-loss rates significantly affect the stellar evolution all the way to core collapse as a supernova and/or black hole formation event., Comment: This is a pre-print of a chapter for the Encyclopedia of Astrophysics (edited by I. Mandel, section editor F.R.N. Schneider) to be published by Elsevier as a Reference Module

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., Comment: 6 pages, 6 figures. ESO Large Programme Overview

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, Comment: 194 pages, 66 figures. Comments are welcome (wstelescope@gmail.com)

Traditionally, the pair instability (PI) mass gap is located between 50\,and 130\,$M_{\odot}$, with stellar mass black holes (BHs) expected to "pile up" towards the lower PI edge. However, this lower PI boundary is based on the assumption that the star has already lost its hydrogen (H) envelope. With the announcement of an "impossibly" heavy BH of 85\,$M_{\odot}$ as part of GW\,190521 located inside the traditional PI gap, we realised that blue supergiant (BSG) progenitors with small cores but large Hydrogen envelopes at low metallicity ($Z$) could directly collapse to heavier BHs than had hitherto been assumed. The question of whether a single star can produce such a heavy BH is important, independent of gravitational wave events. Here, we systematically investigate the masses of stars inside the traditional PI gap by way of a grid of 336 detailed MESA stellar evolution models calculated across a wide parameter space, varying stellar mass, overshooting, rotation, semi-convection, and $Z$. We evolve low $Z$ stars in the range $10^{-3} < Z / Z_{\odot} < Z_{\rm SMC}$, making no prior assumption regarding the mass of an envelope, but instead employing a wind mass loss recipe to calculate it. We compute critical Carbon-Oxygen and Helium core masses to determine our lower limit to PI physics, and we provide two equations for $M_{\text{core}}$ and $M_{\text{final}}$ that can also be of use for binary population synthesis. Assuming the H envelope falls into the BH, we confirm the maximum BH mass below PI is $M_{\text{BH}} \simeq 93.3$ $M_{\odot}$. Our grid allows us to populate the traditional PI gap, and we conclude that the distribution of BHs above the traditional boundary is not solely due to the shape of the initial mass function (IMF), but also to the same stellar interior physics (i.e. mixing) that which sets the BH maximum., Comment: 24 pages, 14 figures. Accepted in MNRAS

Recent studies of high-redshift galaxies with James Webb Space Telescope (JWST), such as GN-z11 at $z=10.6,$ show unexpectedly significant amounts of nitrogen (N) in their spectra. As this phenomenology appears to extend to gravitionally lensed galaxies at Cosmic noon such as the Sunburst Arc at $z=2.37$, as well as globular clusters overall, we suggest that the common ingredient among them are very massive stars (VMSs) with zero-age main sequence (ZAMS) masses in the range of 100-1000 $M_{\odot}$. The He II in the Sunburst Arc might also be the result of the disproportionally large contribution of VMS to the total stellar contribution. We analyse the pros and cons of the previous suggestions, including classical Wolf-Rayet (cWR) stars and supermassive stars (SMSs), to conclude that only our VMS alternative ticks all the relevant boxes. We discuss the VMS mass-loss history via their peculiar vertical evolution in the HR diagram resulting from a self-regulatory effect of these wind-dominated VMSs and we estimate that the large amounts of N present in star-forming galaxies may indeed result from VMSs. We conclude that VMSs should be included in population synthesis and chemical evolution models. Moreover, that it is critical for this to be done self-consistently, as a small error in their mass-loss rates would have dramatic consequences for their stellar evolution, as well as their ionising and chemical feedback., Comment: Accepted as an A&A Letter. 5 pages, 2 figures

Red supergiants (RSG) are key objects for the evolution of massive stars and their endpoints, but uncertainties in their underlying mass-loss mechanism have thus far prevented an appropriate framework for massive star evolution. We analyse an empirical mass loss"kink" feature uncovered by Yang et al., and we highlight its similarity to hot star radiation-driven wind models and observations at the optically thin/thick transition point. We motivate a new RSG mass-loss prescription that depends on the Eddington factor Gamma (including both a steep L dependence and an inverse steep M dependence). We subsequently implement this new RSG mass-loss prescription in the stellar evolution code MESA. We find that our physically motivated mass-loss behaviour naturally reproduces the Humphreys-Davidson limit without a need for any ad-hoc tweaks. It also resolves the RSG supernova "problem". We argue that a universal behaviour of radiation-driven winds across the HR diagram, independent of the exact source of opacity, is a key feature of the evolution of the most massive stars., Comment: Accepted Letter in Astronomy & Astrophysics (A&A). 4 pages. 3 figures

In this introductory chapter of the Special Issue entitled `The Structure and Evolution of Stars', we highlight the recent major progress made in our understanding in the physics that governs stellar interiors. In so doing, we combine insight from observations, 1D evolutionary modelling and 2+3D rotating (magneto)hydrodynamical simulations. Therefore, a complete and compelling picture of the necessary ingredients in state-of-the-art stellar structure theory and areas in which improvements still need to be made are contextualised. Additionally, the over-arching perspective that links all the themes of subsequent chapters is presented., Comment: Accepted for publication to Galaxies. Introductory chapter to Special Issue: https://www.mdpi.com/journal/galaxies/special_issues/star_structure_evolution

The most massive stars provide an essential source of recycled material for young clusters and galaxies. While very massive stars (VMS, M>100M) are relatively rare compared to O stars, they lose disproportionately large amounts of mass already from the onset of core H-burning. VMS have optically thick winds with elevated mass-loss rates in comparison to optically thin standard O-star winds. We compute wind yields and ejected masses on the main sequence, and we compare enhanced mass-loss rates to standard ones. We calculate solar metallicity wind yields from MESA stellar evolution models in the range 50 - 500M, including a large nuclear network of 92 isotopes, investigating not only the CNO-cycle, but also the Ne-Na and Mg-Al cycles. VMS with enhanced winds eject 5-10 times more H-processed elements (N, Ne, Na, Al) on the main sequence in comparison to standard winds, with possible consequences for observed anti-correlations, such as C-N and Na-O, in globular clusters. We find that for VMS 95% of the total wind yields is produced on the main sequence, while only ~5% is supplied by the post-main sequence. This implies that VMS with enhanced winds are the primary source of 26Al, contrasting previous works where classical Wolf-Rayet winds had been suggested to be responsible for Galactic 26Al enrichment. Finally, 200M stars eject 100 times more of each heavy element in their winds than 50M stars, and even when weighted by an IMF their wind contribution is still an order of magnitude higher than that of 50M stars., Comment: Accepted for publication in MNRAS. 14 pages, 10 figures

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\"afener (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., Comment: Accepted for publication in MNRAS, 18 pages, 8 figures, 2 appendix

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., Comment: Accepted in A&A - 35 Pages, 12 Figures, 4 Tables, 2 Large Tables

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., Comment: Accepted to the Publications of the Astronomical Society of the Pacific

WEAVE, the new wide-field, massively multiplexed spectroscopic survey facility for the William Herschel Telescope, will see first light in late 2022. WEAVE comprises a new 2-degree field-of-view prime-focus corrector system, a nearly 1000-multiplex fibre positioner, 20 individually deployable 'mini' integral field units (IFUs), and a single large IFU. These fibre systems feed a dual-beam spectrograph covering the wavelength range 366$-$959\,nm at $R\sim5000$, or two shorter ranges at $R\sim20\,000$. After summarising the design and implementation of WEAVE and its data systems, we present the organisation, science drivers and design of a five- to seven-year programme of eight individual surveys to: (i) study our Galaxy's origins by completing Gaia's phase-space information, providing metallicities to its limiting magnitude for $\sim$3 million stars and detailed abundances for $\sim1.5$ million brighter field and open-cluster stars; (ii) survey $\sim0.4$ million Galactic-plane OBA stars, young stellar objects and nearby gas to understand the evolution of young stars and their environments; (iii) perform an extensive spectral survey of white dwarfs; (iv) survey $\sim400$ neutral-hydrogen-selected galaxies with the IFUs; (v) study properties and kinematics of stellar populations and ionised gas in $z<0.5$ cluster galaxies; (vi) survey stellar populations and kinematics in $\sim25\,000$ field galaxies at $0.3\lesssim z \lesssim 0.7$; (vii) study the cosmic evolution of accretion and star formation using $>1$ million spectra of LOFAR-selected radio sources; (viii) trace structures using intergalactic/circumgalactic gas at $z>2$. Finally, we describe the WEAVE Operational Rehearsals using the WEAVE Simulator., Comment: 41 pages, 27 figures, accepted for publication by MNRAS; updated version including information on individual grants in a revised Acknowledgements section, corrections to the affiliation list, and an updated references list

The ages of stars have historically relied on isochrone fitting of standardised grids of models. While these stellar models have provided key constraints on observational samples of massive stars, they inherit many systematic uncertainties, mainly in the internal mixing mechanisms applied throughout the grid, fundamentally undermining the isochrone method. In this work, we utilise the M-L plane of Higgins & Vink as a method of determining stellar age, with mixing-corrected models applying a calibrated core overshooting alpha_ov and rotation rate to fit the observational data. We provide multiple test-beds to showcase our new method, while also providing comparisons to the commonly-used isochrone method, highlighting the dominant systematic errors. We reproduce the evolution of individual O stars, and analyse the wider sample of O and B supergiants from the VLT-FLAMES Tarantula Survey, providing dedicated models with estimates for alpha_ov, Omega/Omega_crit, and ultimately stellar ages. The M-L plane highlights a large discrepancy in the spectroscopic masses of the O supergiant sample. Furthermore the M-L plane also demonstrates that the evolutionary masses of the B supergiant sample are inappropriate. Finally, we utilise detached eclipsing binaries, VFTS 642 and VFTS 500, and present their ages resulting from their precise dynamical masses, offering an opportunity to constrain their interior mixing. For the near-TAMS system, VFTS 500, we find that both components require a large amount of core overshooting (alpha_ov ~ 0.5), implying an extended main-sequence width. We hence infer that the vast majority of B supergiants are still burning hydrogen in their cores., Comment: 13 pages, 13 figures. Accepted for publication in MNRAS. Comments welcome

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., Comment: Accepted for publication in MNRAS, 14 figures

The origins of coalescing binary black holes (BBHs) detected by the advanced LIGO/Virgo are still under debate, and clues may be present in the joint mass-spin distribution of these merger events. Here we construct phenomenological models containing two sub-populations to investigate the BBH population detected in gravitational wave observations. We find that our models can explain the GWTC-3 data rather well, and several constraints to our model are required by the data: first, the maximum mass for the component with a stellar-origin, $m_{\rm max}$, is $39.1^{+2.4}_{-2.7}M_{\odot}$ at 90\% credibility; second, about $15\%$ of the mergers happen in dynamical environments, in which $7-16\%$ of events are hierarchical mergers, and these BHs have an average spin magnitude significantly larger than the first-generation mergers, with ${\rm d}\mu_{\rm a} > 0.4 $ at $99\%$ credibility; third, the dynamical component BHs tend to pair with each other with larger total mass and higher mass ratio. An independent analysis focusing on spins is also carried out, and we find that the spin amplitude of component BHs can be divided into two groups according to a division mass $m_{\rm d} = 46.1^{+5.6}_{-5.1}M_{\odot}$. These constraints can be naturally explained by current formation channels, and our results suggest that some of the observed events were likely from AGN disks., Comment: 21 pages, 8 figures. Published in The Astrophysical Journal Letters

The Vera C. Rubin Legacy Survey of Space and Time holds the potential to revolutionize time domain astrophysics, reaching completely unexplored areas of the Universe and mapping variability time scales from minutes to a decade. To prepare to maximize the potential of the Rubin LSST data for the exploration of the transient and variable Universe, one of the four pillars of Rubin LSST science, the Transient and Variable Stars Science Collaboration, one of the eight Rubin LSST Science Collaborations, has identified research areas of interest and requirements, and paths to enable them. While our roadmap is ever-evolving, this document represents a snapshot of our plans and preparatory work in the final years and months leading up to the survey's first light., Comment: 202 pages (in book format) 34 figures plus chapter heading figures (13)

The sharpest optical images of the R136 cluster in the Large Magellanic Cloud are presented, allowing for the first time to resolve members of the central core, including R136a1, the most massive star known. These data were taken using the Gemini speckle imager Zorro in medium-band filters with effective wavelengths similar to BVRI achieving angular resolutions between 30-40 mas. All stars previously known in the literature, having $V<$16 mag within the central $2''\times2''$ were recovered. Visual companions ($\geq$40 mas; 2000 au) were detected for the WN5h stars R136 a1, and a3. Photometry of the visual companion of a1 suggests it is of mid O spectral type. Based on new photometric luminosities using the resolved Zorro imaging, the masses of the individual WN5h stars are estimated to be between 150-200 $M_{\odot}$, lowering significantly the present-day masses of some of the most massive stars known. These mass estimates are critical anchor points for establishing the stellar upper-mass function., Comment: Accepted for publication in ApJ; 14 pages, 8 figures, 3 tables

Stellar-mass black holes are the final remnants of stars born with more than 15 solar masses. Billions are expected to reside in the Local Group, yet only few are known, mostly detected through X-rays emitted as they accrete material from a companion star. Here, we report on VFTS 243: a massive X-ray faint binary in the Large Magellanic Cloud. With an orbital period of 10.4-d, it comprises an O-type star of 25 solar masses and an unseen companion of at least nine solar masses. Our spectral analysis excludes a non-degenerate companion at a 5-sigma confidence level. The minimum companion mass implies that it is a black hole. No other X-ray quiet black hole is unambiguously known outside our Galaxy. The (near-)circular orbit and kinematics of VFTS 243 imply that the collapse of the progenitor into a black hole was associated with little or no ejected material or black-hole kick. Identifying such unique binaries substantially impacts the predicted rates of gravitational-wave detections and properties of core-collapse supernovae across the Cosmos., Comment: Accepted to Nature Astronomy, 64 pages, 15 figures, 2 tables; ESO press release: https://www.eso.org/public/news/eso2210/; Nat Asr paper URL: https://www.nature.com/articles/s41550-022-01730-y

The winds of massive stars are important for their direct impact on the interstellar medium, and for their influence on the final state of a star prior to it exploding as a supernova. However, the dynamics of these winds is understood primarily via their illumination from a single central source. The Doppler shift seen in resonance lines is a useful tool for inferring these dynamics, but the mapping from that Doppler shift to the radial distance from the source is ambiguous. Binary systems can reduce this ambiguity by providing a second light source at a known radius in the wind, seen from orbitally modulated directions. From the nature of the collision between the winds, a massive companion also provides unique additional information about wind momentum fluxes. Since massive stars are strong ultraviolet (UV) sources, and UV resonance line opacity in the wind is strong, UV instruments with a high resolution spectroscopic capability are essential for extracting this dynamical information. Polarimetric capability also helps to further resolve ambiguities in aspects of the wind geometry that are not axisymmetric about the line of sight, because of its unique access to scattering direction information. We review how the proposed MIDEX-scale mission Polstar can use UV spectropolarimetric observations to critically constrain the physics of colliding winds, and hence radiatively-driven winds in general. We propose a sample of 20 binary targets, capitalizing on this unique combination of illumination by companion starlight, and collision with a companion wind, to probe wind attributes over a range in wind strengths. Of particular interest is the hypothesis that the radial distribution of the wind acceleration is altered significantly, when the radiative transfer within the winds becomes optically thick to resonance scattering in multiple overlapping UV lines., Comment: 26 pages, 12 figures, Review in a topical collection series of Astrophysics and Space Sciences on the proposed Polstar satellite. arXiv admin note: substantial text overlap with arXiv:2111.11552

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\"afener. We examine the temperature evolution of VMS with mass loss that scales with the luminosity-over-mass (L/M) ratio and the Eddington parameter ($\Gamma_{\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 $\Gamma_{\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., Comment: Accepted for publication in MNRAS. 18 pages, 16 figures, 1 appendix

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

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., Comment: 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 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., Comment: 6 pages, 4 figures. Talk in IAUS 366

One objective of the Polstar spectropolarimetry mission is to characterize the degree of nonconservative mass transfer that occurs at various stages of binary evolution, from the initial mass reversal to the late Algol phase. The proposed instrument combines spectroscopic and polarimetric capabilities, where the spectroscopy can resolve Doppler shifts in UV resonance lines with 10 km/s precision, and polarimetry can resolve linear polarization with 1e-3 precision or better. The spectroscopy will identify absorption by mass streams seen in projection against the stellar disk as a function of orbital phase, hot accretion spots, as well as scattering from extended splash structures, circumbinary disks, and other flows in and above/below the orbital plane (e.g. jets) that fail to be transferred conservatively. The polarimetry affects more the light coming from material not seen against the stellar disk, allowing the geometry of the scattering to be tracked, resolving ambiguities left by the spectroscopy and light-curve information. For example, nonconservative mass streams ejected in the polar direction will produce polarization of the opposite sign from conservative transfer accreting in the orbital plane. Also, time domain coverage over a range of phases of the binary orbit are well supported by the Polstar observing strategy. Combining these elements will significantly improve our understanding of the mass transfer process and the amount of mass that can escape from the system, an important channel for changing the final mass, and ultimate supernova, of the large number of massive stars found in binaries at close enough separation to undergo interaction., Comment: Refereed paper in Ap&SS 367:Topical Collection, 23 pages, 10 figures

The most massive stars are thought to lose a significant fraction of their mass in a steady wind during the main-sequence and blue supergiant phases. This in turn sets the stage for their further evolution and eventual supernova, with consequences for ISM energization and chemical enrichment. Understanding these processes requires accurate observational constraints on the mass-loss rates of the most luminous stars, which can also be used to test theories of stellar wind generation. In the past, mass-loss rates have been characterized via collisional emission processes such as H$\alpha$ and free-free radio emission, but these so-called "density squared" diagnostics require correction in the presence of widespread clumping. Recent observational and theoretical evidence points to the likelihood of a ubiquitously high level of such clumping in hot-star winds, but quantifying its effects requires a deeper understanding of the complex dynamics of radiatively driven winds. Furthermore, large-scale structures arising from surface anisotropies and propagating throughout the wind can further complicate the picture by introducing further density enhancements, affecting mass-loss diagnostics. Time series spectroscopy of UV resonance lines with high resolution and high signal-to-noise are required to better understand this complex dynamics, and help correct "density squared" diagnostics of mass-loss rates. The proposed Polstar mission easily provides the necessary resolution at the sound-speed scale of 20 km s$^{-1}$, on three dozen bright targets with signal-to noise an order of magnitude above that of the celebrated IUE MEGA campaign, via continuous observations that track structures advecting through the wind in real time., Comment: white paper

As sources of chemical enrichment, ionizing radiation and energetic feedback, massive stars drive the ecology of their host galaxies despite their relative rarity, additionally to yielding compact remnants, which can generate gravitational waves. The evolution of massive stars is crucially informed by their detailed mass-loss history; however, wind structures on a variety of scales cause important uncertainties on their mass-loss rates. Binary systems can place further constraints on the mass-loss properties of massive stars, especially colliding-wind binaries. In this paper, we review how the proposed MIDEX-scale mission Polstar can critically constrain the physics of colliding winds (and hence radiatively-driven winds in general) with ultraviolet spectropolarimetric observations, providing an unprecedented improvement on the accuracy of the determination of both mass-loss rates and the velocity structure of the winds of massive stars. We propose a sample of 17 targets that will allow us to study a variety of wind-colliding systems spanning a large parameter space using the spatial information yielded by both spectroscopic and polarimetric data obtained with Polstar., Comment: 21 pages, 12 figures 1 table, Whitepaper

Polstar is a proposed NASA MIDEX mission that would feature a high resolution UV spectropolarimeter capable of measure all four Stokes parameters onboard a 60cm telescope. The mission would pioneer the field of time-domain UV spectropolarimetry. Time domain UV spectropolarimetry offers the best resource to determine the geometry and physical conditions of protoplanetary disks from the stellar surface to <5 AU. We detail two key objectives that a dedicated time domain UV spectropolarimetry survey, such as that enabled by Polstar, could achieve: 1) Test the hypothesis that magneto-accretion operating in young planet-forming disks around lower-mass stars transitions to boundary layer accretion in planet-forming disks around higher mass stars; and 2) Discriminate whether transient events in the innermost regions of planet-forming disks of intermediate mass stars are caused by inner disk mis-alignments or from stellar or disk emissions., Comment: 16 pages, 7 figures

Massive stars have strong stellar winds that direct their evolution through the upper Hertzsprung-Russell diagram and determine the black hole mass function. Secondly, wind strength dictates the atmospheric structure that sets the ionising flux. Thirdly, the wind directly intervenes with the stellar envelope structure, which is decisive for both single star and binary evolution, affecting predictions for gravitational wave events. Key findings of current hot-star research include: * The traditional line-driven wind theory is being updated with Monte Carlo and co-moving frame computations, revealing a rich multi-variate behaviour of the mass-loss rate dM/dt in terms of M, L, Eddington Gamma, Teff , and chemical composition Z. Concerning the latter, dM/dt is shown to depend on the iron (Fe) opacity, making Wolf-Rayet populations, and gravitational wave events dependent on host galaxy Z. * On top of smooth mass-loss behaviour, there are several transitions in the Hertzsprung-Russell diagram, involving bi-stability jumps around Fe recombination temperatures, leading to quasi-stationary episodic, and not necessarily eruptive, Luminous Blue Variable and pre-SN mass loss. * Moreover, there are kinks. At 100 Solar Masses, a high Gamma mass-loss transition implies that hydrogen-rich very massive stars have higher mass-loss rates than commonly considered. At the other end of the mass spectrum, low-mass stripped helium stars no longer appear as Wolf-Rayet stars, but as optically-thin stars. These stripped stars, in addition to very massive stars, are two newly identified sources of ionising radiation that could play a key role in local star formation as well as at high-redshift., Comment: Invited Annual Review of Astronomy and Astrophysics (ARA&A 2022) - 49 pages, 13 figures (now Link to Review in Advance first posted online on May 13, 2022)

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., Comment: Accepted for publication in MNRAS. 15 pages, 13 figures

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., Comment: Accepted for publication in MNRAS. 17 pages, 16 figures

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., Comment: MNRAS Accepted. 12 pages, 13 figures. Notes on IDL/Python Mass-loss Recipe Usage added to the Appendix

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., Comment: MNRAS Accepted. 9 pages, 3 figures. Prior Submission in September 2020 to Nature (see Research Square)

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

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 $\Gamma_\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., Comment: 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

With the launch of JWST and the upcoming installation of extremely large telescopes, the first galaxies in our Universe will finally be revealed. Their light will be dominated by massive stars, which peak in in the ultra-violet (UV) part of the electromagnetic spectrum. Star formation is the key driver of the evolution of our Universe. At young ages, within 10 Million years, both high and low mass stars generate complex UV emission processes which are poorly understood yet are vital for interpreting high red-shift line emission. For these reasons, the Hubble Space Telescope (HST) will devote 1000 orbits to obtaining a UV Legacy Library of Young Stars as Essential Standards (ULLYSES). The purpose of this Overview is to outline the basic physical principles driving UV emission processes from local (within 100 parsecs of) star formation, ranging from huge star-forming complexes containing hundreds of massive and very-massive stars (VMS), such as 30 Doradus (the Tarantula Nebula) in the neighbouring Magellanic Clouds (only 50 kpc away), to galaxies near and far, out to the epoch of Cosmic Reionization., Comment: 12 pages, 5 figures. Comments welcome. Accepted Editorial Overview for Special Issue in Galaxies: "Star Formation in the Ultraviolet", Ed. Jorick Vink

Current massive star evolution grids are not able to simultaneously reproduce the empirical upper luminosity limit of red supergiants, the Humphrey-Davidson (HD) limit at high and low metallicity. In this study, we provide a better understanding of what drives massive star evolution to blue and red supergiant phases, with the ultimate aim of reproducing the HD limit at varied metallicities (Z). For solar, LMC, and SMC Z, we develop eight grids of MESA models for the mass range 20-60M to probe the effect of semiconvection and overshooting. We compare rotating and non-rotating models with efficient (alpha_semi = 100) and inefficient semi-convection (alpha_semi = 0.1), with high and low core overshooting (alpha_ov of 0.1 or 0.5). The red and blue supergiant evolutionary phases are investigated by comparing the fraction of core He-burning lifetimes spent in each phase. We find that the extension of the convective core by overshooting alpha_ov = 0.5 has an effect on the post-MS evolution which can disable semiconvection leading to more RSGs, but a lack of BSGs. We therefore implement alpha_ov = 0.1 which switches on semiconvective mixing, though for standard alpha_semi = 1, would result in an HD limit which is higher than observed at low Z. Therefore, we need to implement very efficient semiconvection of alpha_semi = 100 which reproduces the HD limit at log L ~ 5.5 for the Magellanic Clouds while simultaneously reproducing the Galactic HD limit of log L ~ 5.8 naturally. The effect of semiconvection is not active at high Z due to the depletion of the envelope structure by strong mass loss such that semiconvective regions could not form. Z-dependent mass loss plays an indirect, yet decisive role in setting the HD limit as a function of Z. For a combination of efficient semiconvection and low overshooting with standard Z-dependent mass loss, we find a natural HD limit at all metallicities., Comment: 11pages, 7 figures, accepted for publication in A&A

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., Comment: 19 pages, 14 figures, Astronomy and Astrophysics, in press

Current burning issues in stellar physics, for both hot and cool stars, concern their magnetism. In hot stars, stable magnetic fields of fossil origin impact their stellar structure and circumstellar environment, with a likely major role in stellar evolution. However, this role is complex and thus poorly understood as of today. It needs to be quantified with high-resolution UV spectropolarimetric measurements. In cool stars, UV spectropolarimetry would provide access to the structure and magnetic field of the very dynamic upper stellar atmosphere, providing key data for new progress to be made on the role of magnetic fields in heating the upper atmospheres, launching stellar winds, and more generally in the interaction of cool stars with their environment (circumstellar disk, planets) along their whole evolution. UV spectropolarimetry is proposed on missions of various sizes and scopes, from POLLUX on the 15-m telescope LUVOIR to the Arago M-size mission dedicated to UV spectropolarimetry., Comment: 25 pages, 11 figures, white paper submitted to ESA Voyage 2050 call

We find that applying a theoretical wind mass-loss rate from Monte Carlo radiative transfer models for hydrogen-deficient stars results in significantly more leftover hydrogen following stable mass transfer through Roche-lobe overflow than when we use an extrapolation of an empirical fit for Galactic Wolf-Rayet stars, for which a negligible amount of hydrogen remains in a large set of binary stellar evolution computations. These findings have implications for modelling progenitors of Type Ib and Type IIb supernovae. Most importantly, our study stresses the sensitivity of the stellar evolution models to the assumed mass-loss rates and the need to develop a better theoretical understanding of stellar winds., Comment: Published in MNRAS

We aim to constrain massive star evolution models using the unique testbed eclipsing binary HD166734 with new grids of MESA stellar evolution models, adopting calibrated prescriptions of overshooting, mass loss, and rotation. We introduce a novel tool: the "mass-luminosity plane" or "M-L plane", as an equivalent to the traditional HR diagram, utilising it to reproduce the testbed binary HD166734 with newly calibrated MESA stellar evolution models for single stars. We can only reproduce the Galactic binary system with an enhanced amount of core overshooting (alpha = 0.5), mass loss, and rotational mixing. We can utilise the gradient in the M-L plane to constrain the amount of mass loss to 0.5 - 1.5 times the standard Vink et al. 2001 prescriptions, and we can exclude extreme reduction or multiplication factors. The extent of the vectors in the M-L plane leads us to conclude that the amount of core overshooting is larger than is normally adopted in contemporary massive star evolution models. We furthermore conclude that rotational mixing is mandatory to get the nitrogen abundance ratios between the primary and secondary components to be correct (3:1) in our testbed binary system. Our calibrated grid of models, alongside our new M-L plane approach, present the possibility of a widened main sequence due to an increased demand for core overshooting. The increased amount of core overshooting is not only needed to explain the extended main sequence, but the enhanced overshooting is also needed to explain the location of the upper-luminosity limit of the red supergiants. Finally, the increased amount of core overshooting has -- via the compactness parameter -- implications for supernova explodibility., Comment: Revised version, 14 pages, 19 figures, accepted for publication by Astronomy & Astrophysics

One of the largest surprises from the LIGO results regarding the first gravitational wave detection (GW 150914) was the fact the black holes (BHs) were "heavy", of order 30 - 40 solar masses. The most promising explanation for this obesity is that the BH-BH merger occurred at low metallicity (Z): when the iron (Fe) contents is lower this is expected to result in weaker mass loss during the Wolf-Rayet (WR) phase. We therefore critically evaluate the claims for the reasons of heavy BHs as a function of Z in the literature. Furthermore, weaker stellar winds might lead to more rapid stellar rotation, allowing WR and BH progenitor evolution in a chemically homogeneous manner. However, there is as yet no empirical evidence for more rapid rotation amongst WR stars in the low Z environment of the Magellanic Clouds. Due to the intrinsic challenge of determining WR rotation rates from emission lines, the most promising avenue to constrain rotation-rate distributions amongst various WR subgroups is through the utilisation of their emission lines in polarised light. We thus provide an overview of linear spectro-polarimetry observations of both single and binary WRs in the Galaxy, as well as the Large and Small Magellanic Clouds, at 50% and 20% of solar Z, respectively. Initial results suggest that the route of chemically homogeneous evolution (CHE) through stellar rotation is challenging, whilst the alternative of a post-LBV or common envelope evolution is more likely., Comment: 6 pages, 3 figures, oral contribution: IAU Symposium 346 "High Mass X-ray Binaries: illuminating the passage from massive binaries to merging compact objects", Vienna, Austria, 27-31 August 2018

Massive star evolution is dominated by key physical processes such as mass loss, convection and rotation, yet these effects are poorly constrained, even on the main sequence. We utilise a detached, eclipsing binary HD166734 as a testbed for single star evolution to calibrate new MESA stellar evolution grids. We introduce a novel method of comparing theoretical models with observations in the 'Mass-Luminosity Plane', as an equivalent to the HRD (see Higgins & Vink, 2018). We reproduce stellar parameters and abundances of HD166734 with enhanced overshooting (alpha=0.5), mass loss and rotational mixing. When comparing the constraints of our testbed to the systematic grid of models we find that a higher value of alpha = 0.5 (rather than 0.1) results in a solution which is more likely to evolve to a neutron star than a black hole, due to a lower value of the compactness parameter.