Search

Stripped-envelope supernovae (SNe) are H-poor transients produced at the end of the life of massive stars that previously lost their H-rich envelope. Their progenitors are thought to be donor stars in mass-transferring binary systems, which were stripped of their H-rich envelopes some $10^6$yr before core collapse. A subset of the stripped-envelope SNe exhibit spectral and photometric features indicative of interaction between their ejecta and nearby circumstellar material (CSM). We examine whether mass transfer during, or shortly before, core collapse in massive binary systems can produce the CSM inferred from the observations of interacting H-poor SNe. We select 44 models from a comprehensive grid of detailed binary evolution models in which the mass donors are H-free and explode while transferring mass to a main-sequence companion. We find that in these models, mass transfer starts less than $\sim20$kyr before, and often continues until the core collapse of the donor star. Up to $0.8M_\odot$ of H-free material are removed from the donor star during this phase, which may produce a He-rich circumbinary material. We explore plausible assumptions for its spatial distribution at the time of explosion. When assuming that the CSM accumulates in a circumbinary disk, we find qualitative agreement with the supernova and CSM properties inferred from observed Type Ibn SNe, and to a lesser extent with constraints from Type Icn SNe. We find that our mass transferring stripped envelope SN progenitor models may produce up to $\sim$10% of all stripped envelope supernovae. The binary channel proposed in this work can qualitatively account for the observed key properties and rate of interacting H-poor SNe. Models for the evolution of the circumbinary material and the spectral evolution of exploding progenitors from this channel are needed to further test its significance., Comment: 28 pages, 23 figures. Submitted to Astronomy & Astrophysics. Abstract is abridged. Comments are welcome!

Context. The influence of magnetic fields on stellar evolution remains unresolved. It has been proposed that if there is a large-scale magnetic field in the stellar interior, torsional waves could arise, efficiently transporting angular momentum. In fact, the observed variations in the rotation periods of some magnetic stars may be attributed to these torsional waves' standing waves. Aims. To demonstrate the existence of torsional waves through modeling of the rotational period variations. Method. We conduct an eigenmode analysis of standing waves based on one-dimensional magnetohydrodynamic equations. The internal magnetic field structures are parametrically represented to treat poloidal fields with different degrees of central/surface concentration. The obtained frequencies are compared with the observed frequencies of the rotational period variations, thereby constraining the internal magnetic field structures. Results. The 67.6 years exhibited by CU Vir is reproduced for surface-concentrated magnetic field structures. The rotational period variations of all ten magnetic stars analyzed in this study are inconsistent with a centrally concentrated magnetic field. Conclusions. Torsional waves can reproduce the observations of rotational period variations. The large-scale magnetic fields within magnetic stars would be concentrated on the surface., Comment: 12 pages, 7 figures, 1 table. submitted to A&A

Among massive stars, binary interaction is the rule rather than the exception. The closest binaries, those with periods of less than about 10 days, undergo mass transfer during core-hydrogen burning, with many of them experiencing a nuclear-timescale contact phase. Current binary population synthesis models predict the mass-ratio distribution of contact binaries to be heavily skewed toward a mass ratio of unity, which is inconsistent with observations. It has been shown that effects of tidal deformation due to the Roche potential, as well as energy transfer in the common layers of a contact binary, alter the internal structure of close binary components. However, previous population studies neglected these effects. We model a population of massive binary stars that undergo mass transfer during core-hydrogen burning, while consistently considering the effects of tidal deformation and energy transfer in contact phases. We use the MESA binary-evolution code to compute large grids of models with primary star masses of $8\,M_\odot$ to $70\,M_\odot$ at Solar metallicity. We then perform a population synthesis study to predict distribution functions of the observational properties of close binary systems, focusing in particular on the mass and luminosity ratio distribution. We find that the effects of tidal deformation and energy transfer have a limited effect on the predicted mass-ratio distribution of massive contact binaries. Only a small fraction of the population has their mass ratio significantly shifted toward a more unequal configuration. However, we suggest that orbital hardening could affect the evolution of contact binaries and their progenitors, and we advocate for a homogeneous set of observed contact binary parameters., Comment: 12 pages, 10 Figures, (22 pages, 16 Figs. including appendices). Submitted to A&A

Massive contact binaries are fundamental test-beds to study properties of close binaries en route to stellar mergers. Progress in their understanding has been largely impeded due to a dearth in their observed counterparts. We accumulate 61 over-contact binary candidates from the OGLE survey of the Magellanic Clouds and form the largest known sample of such systems. We mine the OGLE-IV catalogue for B-type and earlier systems which match our synthetic light curves, and thereafter examine their properties in comparison with our binary populations. Simultaneously, we also investigate the degeneracy in the light curves between over-contact and near-contact binaries and attempt to distinguish these populations. We construct new synthetic populations of such systems from detailed $\texttt{MESA}$ models with initial parameters of M$_\textrm{total,i}=14-80\,\textrm{M}_{\odot}$ and P$_\textrm{i}=0.6-23\,$days, and compute analogous $\texttt{PHOEBE}$ light curves. We then extract OGLE systems that match our synthetic light curves and compare their numbers, and distributions in periods, absolute V-band magnitudes and colours (V-I) with our population predictions. We expect massive over-contact binaries to have smooth equal-depth light curves due to their equal-mass ($q\approx1$) and consequently, equal-temperature components. We mine 47 and 14 predominantly B-type over-contact OGLE binaries from the Large and Small Magellanic Clouds respectively, with periods of $\approx0.6-1\,$d, along with a few wider O-type systems. This sample is consistent with our theoretically expected numbers and colour-magnitude and period distributions. Additionally, we find that binaries nearing contact have periods as large as $4.5\,$d and mass ratios of $q\approx0.4-0.95$, some of which may masquerade as over-contact binaries and thereby lead to potential observational misclassifications., Comment: 14 pages, 11 figures, submitted to A\&A

Massive stars mainly form in close binaries, where their mutual interactions can profoundly alter their evolutionary paths. Evolved binaries consisting of a massive OB-type main-sequence star with a stripped helium star or a compact companion represent a crucial stage in the evolution towards double compact objects, whose mergers are (potentially) detectable via gravitational waves. The recent detection of X-ray quiet OB+black hole binaries and OB+stripped helium star binaries has set the stage for discovering more of these systems in the near future. In this work, based on 3670 detailed binary-evolution models and using empirical distributions of initial binary parameters, we compute the expected population of such evolved massive binaries in coeval stellar populations, including stars in star clusters and in galaxies with starburst activities, for ages up to 100 Myr. Our results are vividly illustrated in an animation that shows the evolution of these binaries in the color-magnitude diagram over time. We find that the number of OB+black hole binaries peaks around 10 Myr, and OB+neutron star binaries are most abundant at approximately 20 Myr. Both black holes and neutron stars can potentially be found in populations with ages up to 90 Myr. Additionally, we analyze the properties of such binaries at specific ages. We find that OB+helium stars and OB+black hole binaries are likely to be identifiable as single-lined spectroscopic binaries. Our research serves as a guide for future observational efforts to discover such binaries in young star clusters and starburst environments., Comment: 32 pages, 15 figures. Accepted for Publication in Astrophysical Journal Letters (ApJL)

Context. Recent astrometric and spectroscopic surveys of OB stars have revealed a few stellar-mass black holes (BHs) with orbital periods as low as 10 days. No X-ray counterpart has been detected, due to the absence of a radiatively efficient accretion disk around the BH. Yet, dissipative processes in the hot, dilute and strongly magnetized plasma around the BH (so-called BH corona) can still lead to non-thermal X-ray emission (e.g. synchrotron). Aims. We determine the X-ray luminosity distribution from BH+OB star binaries up to orbital periods of a few thousand days. Methods. We use detailed binary evolution models computed with MESA for initial primary masses of 10-90 $M_{\odot}$ and orbital periods from 1-3000 d. The X-ray luminosity is computed for a broad range of radiative efficiencies. Results. We show that particle acceleration through magnetic reconnection can heat the BH corona. A substantial fraction of the gravitational potential energy from the accreted plasma is converted into non-thermal X-ray emission. Our population synthesis analysis predicts at least 28 (up to 72) BH+OB star binaries in the Large Magellanic Cloud (LMC) to produce X-ray luminosity above 10$^{31}$ erg$\,$s$^{-1}$, observable through focused Chandra observations. We identify a population of SB1 systems in the LMC and HD96670 in the Milky Way comprising O stars with unseen companions of masses above 2.3 $M_{\odot}$ that aligns well with our predictions. The predicted luminosities of the OB companions to these X-ray-emitting BHs are 10$^{4.5-5.5}$ $L_{\odot}$. Conclusions. These results make the case for long-time exposure in X-rays of the stellar-mass BH candidates identified around OB stars. It will constrain the underlying population of X-ray-faint BHs, the evolution from single to double degenerate binaries, and the progenitors of gravitational wave mergers. (Abridged), Comment: A&A, 690, A256. DOI: https://doi.org/10.1051/0004-6361/202450940

The evolution and fate of massive stars are thought to be affected by rotationally induced internal mixing. The surface boron abundance is a sensitive tracer of this in early B-type main sequence stars. We test current stellar evolution models of massive main sequence stars which include rotational mixing through a systematic study of their predicted surface boron depletion. We construct a dense grid of rotating single star models using MESA, for which we employ a new nuclear network which follows all the stable isotopes up to silicon, including lithium, beryllium, boron, as well as the radioactive isotope aluminium-26. We also compile the measured physical parameters of the 90 Galactic early B-type stars with boron abundance information. We then compare each observed stars with our models through a Bayesian analysis, which yields the mixing efficiency parameter with which the star is reproduced the best, and the probability that it is represented by the stellar models. We find that about two-thirds of the sample stars are well represented by the stellar models, with the best agreement achieved for a rotational mixing efficiency of ~50% compared to the widely adopted value. The remaining one third of the stars, of which many are strongly boron depleted slow rotators, are largely incompatible with our models, for any rotational mixing efficiency. We investigate the observational incidence of binary companions and surface magnetic fields, and discuss their evolutionary implications. Our results confirm the concept of rotational mixing in radiative stellar envelopes. On the other hand, we find that a different boron depletion mechanism, and likely a different formation path, is required to explain about one-third of the sample stars. The large spread in the surface boron abundances of these stars may hold a clue to understanding their origin.

The majority of massive stars reside in binary systems, which are expected to experience mass transfer during their evolution. However, so far the conditions under which mass transfer leads to a common envelope, and thus possibly to a merging of both stars, are not well understood. Main uncertainties arise from the possible swelling of the mass gainer, and from angular momentum loss from the binary system, during non-conservative mass transfer. We have computed a dense grid of detailed models of stars accreting mass at constant rates, to determine their radius increase due to their thermal disequilibrium. While we find that models with faster than thermal timescale accretion generally expand, this expansion remains quite limited in the intermediate mass regime even for accretion rates which exceed the thermal timescale accretion rate by a factor of 100. Our models of massive accretion stars expand to extreme radii under those conditions. When the accretion rate exceed the Eddington accretion rate, our models expand dynamically. We have derived analytical fits to the radius evolution of our models and a prescription for the borderline between stable mass transfer and mergers for arbitrary accretion efficiencies. We then apply our results to grids of binary models adopting various constant mass transfer efficiencies and angular momentum budgets. We find that the former parameter has the stronger effect on the outcome of the Roche lobe overflow. Our results are consistent with detailed binary evolution models, and often lead to a smaller initial parameter space for stable mass transfer than other recipes in the literature. We use this method to investigate the origin of the Wolf-Rayet stars with O star companions in the Small Magellanic Cloud, and find that the efficiency of the mass transfer process which lead to the formation of the Wolf-Rayet star was likely below 50%., Comment: 19 pages

Massive binary evolution models are needed to predict massive star populations in star forming galaxies, the supernova diversity, and the number and properties of gravitational wave sources. Such models are often computed using so called rapid binary evolution codes, which approximate the evolution of the binary components based on detailed single star models. However, about one third of the interacting massive binary stars undergo mass transfer during core hydrogen burning (Case A mass transfer), whose outcome is difficult to derive from single star models. Here, we use a large grid of detailed binary evolution models for primaries in the initial mass range 10 to 40 Solar masses of LMC and SMC composition, to derive analytic fits for the key quantities needed in rapid binary evolution codes, i.e., the duration of core hydrogen burning, and the resulting donor star mass. Systems with shorter orbital periods produce up to 50% lighter stripped donors and have a up to 30% larger lifetime than wider systems. We find that both quantities depend strongly on the initial binary orbital period, but that the initial mass ratio and the mass transfer efficiency of the binary have little impact on the outcome. Our results are easily parameterisable and can be used to capture the effects of Case A mass transfer more accurately in rapid binary evolution codes., Comment: 10 pages

Over the last decade, evidence has accumulated that massive stars do not typically evolve in isolation but instead follow a tumultuous journey with a companion star on their way to core collapse. While Roche-lobe overflow appears instrumental for the production of a large fraction of supernovae (SNe) of Type Ib and Ic, variations in the initial orbital period Pinit of massive interacting binaries may also produce a wide diversity of case B, BC, or C systems, with preSN stars endowed from minute to massive H-rich envelopes. Focusing here on the explosion of the primary, donor star, originally of 12.6Msun, we use radiation-hydrodynamics and NLTE time-dependent radiative transfer to document the gas and radiation properties of such SNe, covering from Type Ib, IIb, II-L to II-P. Variations in Pinit are the root cause behind the wide diversity of our SN light curves, with single-peak, double-peak, fast-declining or plateau-like morphologies in the V band. The different ejecta structures, expansion rates, and relative abundances (e.g., H, He, 56Ni) are conducive to much diversity in spectral line shapes (absorption vs emission strength, width) and evolution. We emphasize that Halpha is a key tracer of these modulations, and that HeI7065 is an enduring optical diagnostic for the presence of He. Our grid of simulations fare well against representative SNe Ib, IIb, and IIP SNe, but interaction with circumstellar material, which is ignored in this work, is likely at the origin of the tension between our Type IIL SN models and observations (e.g., SN2006Y). Remaining discrepancies in our model rise time to bolometric maximum call for a proper account of both small-scale and large-scale structures in core-collapse SN ejecta. Discrepant Type IIP SN models, with a large plateau brightness but small line widths, may be cured by adopting more compact red-supergiant star progenitors., Comment: Submitted to A&A on Dec 22nd, 2023

New boron abundances or upper limits have been determined for 8 early-B stars in the young Galactic open cluster NGC 3293, using ultraviolet spectra obtained by the Hubble Space Telescope Cosmic Origins Spectrograph. With previous observations, there are now 18 early-B stars in this cluster with boron measurements. Six of the newly observed stars have projected rotational velocities greater than 200 km/s, allowing new constraints on rotationally driven mixing in main-sequence stars. When comparing to synthetic model populations, we find that the majority of our sample stars agree well with the predicted trends of stronger boron depletion for larger rotation and for larger mass or luminosity. Based on those, a smaller than the canonical rotational mixing efficiency,(fc = 0.0165 vs the more standard value of 0.033), appears to be required. However, our five most slowly rotating stars are not well explained by rotational mixing, and we speculate that they originate from binary mergers., Comment: 18 pages,7 figures, Submitted to AAS Journals

Blue supergiants are the brightest stars in their host galaxies and yet their evolutionary status has been a long-standing problem in stellar astrophysics. In this pioneering work, we present a large sample of 59 early B-type supergiants in the Large Magellanic Cloud with newly derived stellar parameters and identify the signatures of stars born from binary mergers among them. We simulate novel 1D merger models of binaries consisting of supergiants with hydrogen-free cores (primaries) and main-sequence companions (secondaries) and consider the effects of interaction of the secondary with the core of the primary. We follow the evolution of the new-born $16-40$ M$_{\odot}$ stars until core-carbon depletion, close to their final pre-explosion structure. Unlike stars which are born alone, stars born from such stellar mergers are blue throughout their core helium-burning phase and reproduce the surface gravities and Hertzsprung-Russel diagram positions of most of our sample. This indicates that the observed blue supergiants are structurally similar to merger-born stars. Moreover, the large nitrogen-to-carbon and oxygen ratios, and helium enhancements exhibited by at least half our data sample are uniquely consistent with our model predictions, leading us to conclude that a large fraction of blue supergiants are indeed products of binary mergers.

Many supernovae (SNe) imply an interaction of the SN ejecta with matter (CSM) surrounding the progenitor star. This suggests that many massive stars may undergo various degrees of envelope stripping shortly before exploding, and produce a considerable diversity in their pre-explosion CSM properties. We explore a generic set of ~100 detailed massive binary evolution models to characterize the amount of envelope stripping and the expected CSM configurations. Our binary models were computed with the MESA stellar evolution code, considering an initial primary star mass of 12.6 Msun, and focus on initial orbital periods above 500 d. We compute these models up to the time of the primary's iron core collapse. We find that Roche lobe overflow often leads to incomplete stripping of the mass donor, resulting in a large variety of pre-SN envelope masses. Many of our models' red supergiant (RSG) donors undergo core collapse during Roche lobe overflow, with mass transfer and thus system mass loss rates of up to 0.01 Msun/yr at that time. The corresponding CSM densities are similar to those inferred for Type IIn SNe like 1998S. In other cases, the mass transfer turns unstable, leading to a common envelope phase at such late time that the mass donor explodes before the common envelope is fully ejected or the system has merged. We argue that this may cause significant pre-SN variability, as for example in SN 2020tlf. Other models suggest a common envelope ejection just centuries before core collapse, which may lead to the strongest interactions, as in superluminous Type IIn SNe like 1994W, or 2006gy. Wide massive binaries offer a natural framework to understand a broad range of hydrogen-rich interacting SNe. On the other hand, the flash features observed in many Type IIP SNe, like in SN 2013fs, may indicate that RSG atmospheres are more extended than currently assumed., Comment: 29 pages, 16 figures. Accepted by A&A on 29.01.2024. Abstract is abridged. The main text and the appendix have changed following the comments from the anonymous referee

Context. It is common for massive stars to engage in binary interaction. In close binaries, the components can enter a contact phase, where both stars overflow their respective Roche lobes simultaneously. While there exist observational constraints on the stellar properties of such systems, the most detailed stellar evolution models that feature a contact phase are not fully reconcilable with those measurements. Aims. We aim to consistently model contact phases of binary stars in a 1D stellar evolution code. To this end, we develop the methodology to account for energy transfer in the common contact layers. Methods. We implement an approximative model for energy transfer between the components of a contact binary based on the von Zeipel theorem in the stellar evolution code MESA. We compare structure and evolution models with and without this transfer and analyze the implications for the observable properties of the contact phase. Results. Implementing energy transfer helps eliminating baroclinicity in the common envelope between the components of a contact binary, which, if present, would drive strong thermal flows. We find that accounting for energy transfer in massive contact binaries significantly alters the mass ratio evolution and can extend the lifetime of an unequal mass ratio contact system., Comment: 11 pages, 5 figures. Accepted to A&A

Context: The surface properties of rotating stars can vary from pole to equator, resulting in anisotropic stellar winds which are not included in the currently available evolutionary models. Aims: We develop a formalism to describe the mass and angular momentum loss of rotating stars which takes into account both the varying surface properties and distortion due to rotation. Methods: Adopting the mass-loss recipe for non-rotating stars, we assigned to each point on the surface of a rotating star an equivalent non-rotating star, for which the surface mass flux is given by the recipe. The global mass-loss and angular momentum loss rates are then given by integrating over the deformed stellar surface as appropriate. Evolutionary models were computed and our prescription is compared to the currently used simple mass-loss enhancement recipes for rotating stars. Results: We find that mass-loss rates are largely insensitive to rotation for models not affected by the bi-stability jump. For those affected by the bi-stability jump, the increase in mass-loss rates with respect to time is smoothed. As our prescription considers the variation of physical conditions over the stellar surface, the region affected by the bi-stability jump is able to grow gradually instead of the whole star suddenly being affected. Conclusion: We have provided an easy to implement and flexible, yet physically meaningful prescription for calculating mass and angular momentum loss rates of rotating stars in a one-dimensional stellar evolution code which compares favourably to more physically comprehensive models. The implementation of our scheme in the stellar evolution code MESA is available online: https://zenodo.org/record/7437006, Comment: Accepted for publication in A&A. 12 pages

We perform three-dimensional simulations of magnetorotational supernovae using a $39\,M_{\odot}$ progenitor star with two different initial magnetic field strengths of $10^{10}$ G and $10^{12}$ G in the core. Both models rapidly undergo shock revival and their explosion energies asymptote within a few hundred milliseconds to values of $\gtrsim 2\times10^{51}$ erg after conservatively correcting for the binding energy of the envelope. Magnetically collimated, non-relativistic jets form in both models, though the jets are subject to non-axisymmetric instabilities. The jets do not appear crucial for driving the explosion, as they only emerge once the shock has already expanded considerably. Our simulations predict moderate neutron star kicks of about $150\, \mathrm{km}\,\mathrm{s}^{-1}$, no spin-kick alignment, and rapid early spin-down that would result in birth periods of about $20\, \mathrm{ms}$, too slow to power an energetic gamma-ray burst jet. More than $0.2\,M_\odot$ of iron-group material are ejected, but we estimate that the mass of ejected $^{56}\mathrm{Ni}$ will be considerably smaller as the bulk of this material is neutron-rich. Explosive burning does not contribute appreciable amounts of $^{56}\mathrm{Ni}$ because the burned material originates from the slightly neutron-rich silicon shell. The iron-group ejecta also show no pronounced bipolar geometry by the end of the simulations. The models thus do not immediately fit the characteristics of observed hypernovae, but may be representative of other transients with moderately high explosion energies. The gravitational-wave emission reaches high frequencies of up to 2000 Hz and amplitudes of over 100 cm. The gravitational-wave emission is detectable out to distances of $\sim4$ Mpc in the planned Cosmic Explorer detector.

Spectroscopic observations of stars in young open clusters have revealed evidence for a dichotomous distribution of stellar rotational velocities, with 10-30% of stars rotating slowly and the remaining 70-90% rotating fairly rapidly. At the same time, high-precision multiband photometry of young star clusters shows a split main sequence band, which is again interpreted as due to a spin dichotomy. Recent papers suggest that extreme rotation is required to retrieve the photometric split. Our new grids of MESA models and the prevalent SYCLIST models show, however, that initial slow (0-35% of the linear Keplerian rotation velocities) and intermediate (50-65% of the Keplerian rotation velocities) rotation are adequate to explain the photometric split. These values are consistent with the recent spectroscopic measurements of cluster and field stars, and are likely to reflect the birth spin distributions of upper main-sequence stars. A fraction of the initially faster-rotating stars may be able to reach near-critical rotation at the end of their main-sequence evolution and produce Be stars in the turn-off region of young star clusters. However, we find that the presence of Be stars up to two magnitudes below the cluster turnoff advocates for a crucial role of binary interaction in creating Be stars. We argue that surface chemical composition measurements may help distinguish these two Be star formation channels. While only the most rapidly rotating, and therefore nitrogen-enriched, single stars can evolve into Be stars, slow pre-mass-transfer rotation and inefficient accretion allows for mild or no enrichment even in critically rotating accretion-induced Be stars. Our results shed new light on the origin of the spin distribution of young and evolved B-type main sequence stars., Comment: 18 pages, 14 figures. Accepted for publication in A&A

We probe how common extremely rapid rotation is among massive stars in the early universe by measuring the OBe star fraction in nearby metal-poor dwarf galaxies. We apply a new method that uses broad-band photometry to measure the galaxy-wide OBe star fractions in the Magellanic Clouds and three more distant, more metal-poor dwarf galaxies. We find OBe star fractions of ~20% in the Large Magellanic Cloud (0.5 Z_Solar), and ~30% in the Small Magellanic Cloud (0.2 Z_Solar) as well as in the so-far unexplored metallicity range from 0.1 Z_solar to 0.2 Z_solar occupied by the other three dwarf galaxies. Our results imply that extremely rapid rotation is common among massive stars in metal-poor environments such as the early universe., Comment: Conference proceedings for a talk in IAU Symposium 361: Massive Stars Near and Far, Ballyconnell, Ireland, 9-13 May 2022

Massive stars are predominantly found in binaries and higher order multiples. While the period and eccentricity distributions of OB stars are now well established across different metallicity regimes, the determination of mass-ratios has been mostly limited to double-lined spectroscopic binaries. As a consequence, the mass-ratio distribution remains subject to significant uncertainties. Open questions include the shape and extent of the companion mass-function towards its low-mass end and the nature of undetected companions in single-lined spectroscopic binaries. In this contribution, we present the results of a large and systematic analysis of a sample of over 80 single-lined O-type spectroscopic binaries (SB1s) in the Milky Way and in the Large Magellanic Cloud (LMC). We report on the developed methodology, the constraints obtained on the nature of SB1 companions, the distribution of O star mass-ratios at LMC metallicity and the occurrence of quiescent OB+black hole binaries., Comment: To appear in the proceedings of IAUS361: Massive stars near and far; 6 pages, 1 figure

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

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

Stripped-envelope stars can be observed as Wolf-Rayet (WR) stars, or as less luminous hydrogen-poor stars with low mass loss rates and transparent winds. Both types are potential progenitors of Type I core-collapse supernovae (SNe). We use grids of core-collapse models obtained from helium stars at different metallicities to study the effects of metallicity on the transients and remnants these stars produce. We characterise the surface and core properties of our core collapse models, and investigate their explodability employing three criteria. In cases where explosions are predicted, we estimate the ejecta mass, explosion energy, nickel mass and neutron star (NS) mass. Otherwise, we predict the mass of the resulting black hole (BH). We construct a simplified population model, and find that the properties SNe and compact objects depend strongly on metallicity. Ejecta masses and explosion energies for Type Ic SNe are best reproduced by models with Z=0.04 which exhibit strong winds during core helium burning. This implies that either their mass loss rates are underestimated, or that Type Ic SN progenitors experience mass loss through other mechanisms before exploding. The distributions of ejecta masses, explosion energies and nickel mass for Type Ib SNe are not well reproduced by progenitor models with WR mass loss, but are better reproduced if we assume no mass loss in progenitors with luminosities below the minimum WR star luminosity. We find that Type Ic SNe become more common as metallicity increases, and that the vast majority of progenitors of Type Ib SNe must be transparent-wind stripped-envelope stars. We find several models with pre-collapse CO-masses of up to $\sim 30 M_{\odot}$ may form $\sim 3 M_{\odot}$ BHs in fallback SNe. This may carry important consequences for our understanding of SNe, binary BH and NS systems, X-ray binary systems and gravitational wave transients., Comment: 30 pages, 20 figures, submitted to Astronomy & Astrophysics

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

Recent high-quality Hubble Space Telescope (HST) photometry shows that the main sequences (MS) stars of young star clusters form two discrete components in the color-magnitude diagram (CMD). Based on their distribution in the CMD, we show that stars of the blue MS component can be understood as slow rotators originating from stellar mergers. We derive the masses of the blue MS stars, and find that they follow a nearly flat mass function, which supports their unusual formation path. Our results imply that the cluster stars gain their mass in two different ways, by disk accretion leading to rapid rotation, contributing to the red MS, or by binary merger leading to slow rotation and populating the blue MS. We also derive the approximate merger time of the individual stars of the blue MS component, and find a strong early peak in the merger rate, with a lower level merger activity prevailing for tens of Myr. This supports recent binary formation models, and explains new velocity dispersion measurements for members of young star clusters. Our findings shed new light on the origin of the bi-modal mass, spin, and magnetic field distributions of main-sequence stars., Comment: Published in Nature Astronomy online. Link to the paper:https://rdcu.be/cGLpH

(abridged) When stripped from their hydrogen-rich envelopes, stars with initial masses between $\sim$7 and 11 M$_\odot$ develop massive degenerate cores and collapse. Depending on the final structure and composition, the outcome can range from a thermonuclear explosion, to the formation of a neutron star in an electron-capture supernova (ECSN). It has been recently demonstrated that stars in this mass range may initiate explosive oxygen burning when their central densities are still below $\rho_{\rm c} \lesssim 10^{9.6}$ g cm$^{-3}$. This makes them interesting candidates for type Ia supernovae -- which we call (C)ONe SNe Ia -- and might have broader implications for the formation of neutron stars via ECSNe. Here, we model the evolution of 252 helium stars with initial masses in the $0.8-3.5$ M$_\odot$ range, and metallicities between $Z=10^{-4}$ and $0.02$. We use these models to constrain the central densities, compositions and envelope masses at the time of explosive oxygen ignition. We further investigate the sensitivity of these properties to mass loss rate assumptions using additional models with varying wind efficiencies. We find that helium stars with masses between $\sim$1.8 and 2.7 M$_\odot$ evolve onto $1.35-1.37$ M$_\odot$ (C)ONe cores that initiate explosive burning at central densities between $\rm \log_{10}(\rho_c)\sim 9.3$ and 9.6. We constrain the amount of residual carbon retained after core carbon burning, and conclude that it plays a critical role in determining the final outcome: Chandrasekhar-mass degenerate cores that retain more than $\sim 0.005$ M$_\odot$ of carbon result in (C)ONe SNe Ia, while those with lower carbon mass become ECSNe. We find that (C)ONe SNe Ia are more likely to occur at high metallicities, whereas at low metallicities ECSNe dominate., Comment: Accepted in Astronomy & Astrophysics; v2 matches the accepted version. Comments are welcome

Massive stars that become stripped of their hydrogen envelope through binary interaction or winds can be observed either as Wolf-Rayet stars, if they have optically thick winds, or as transparent-wind stripped-envelope stars. We approximate their evolution through evolutionary models of single helium stars, and compute detailed model grids in the initial mass range 1.5 to 70 M$_{\odot}$ for metallicities between 0.01 and 0.04, from core helium ignition until core collapse. Throughout their lifetime, some stellar models expose the ashes of helium burning. We propose that models that have nitrogen-rich envelopes are candidate WN stars, while models with a carbon-rich surface are candidate WC stars during core helium burning, and WO stars afterwards. We measure metallicity dependance of the total lifetime of our models and the duration of their evolutionary phases. We propose an analytic estimate of the wind optical depth to distinguish models of Wolf-Rayet stars from transparent-wind stripped-envelope stars, and find that the luminosity ranges at which WN, WC and WO type stars can exist is a strong function of metallicity. We find that all carbon-rich models produced in our grids have optically thick winds and match the luminosity distribution of observed populations. We construct population models and predict the numbers of transparent-wind stripped-envelope stars and Wolf-Rayet stars, and derive their number ratios at different metallicities. We find that as metallicity increases, the number of transparent-wind stripped-envelope stars decreases and the number of Wolf-Rayet stars increases. At high metallicities WC and WO type stars become more common. We apply our population models to nearby galaxies, and find that populations are more sensitive to the transition luminosity between Wolf-Rayet stars and transparent-wind helium stars than to the metallicity dependent mass loss rates., Comment: 20 pages, 10 figures, accepted for publication Astronomy & Astrophysics

The first three observing runs with Advanced LIGO and Virgo have resulted in the detection of binary black hole mergers (BBH) with highly unequal mass components, which are difficult to reconcile with standard formation paradigms. The most representative of these is GW190814, a highly asymmetric merger between a 23 M$_{\odot}$ black hole and a 2.6 M$_{\odot}$ compact object. Here, we explore recent results suggesting that a sizeable fraction of stars with pre-collapse carbon-oxygen core masses above 10 M$_{\odot}$, and extending up to at least 30 M$_{\odot}$, may produce objects inside the so-called lower mass gap that bridges the division between massive pulsars and BHs in Galactic X-ray binaries. We demonstrate that such an explosion landscape would naturally cause a fraction of massive binaries to produce GW190814-like systems instead of symmetric-mass BBHs. We present examples of specific evolutionary channels leading to the formation of GW190814 and GW200210, a 24+2.8 M$_{\odot}$ merger discovered during the O3b observing run. We estimate the merger-rate density of these events in our scenario to be $\mathcal{O}$(5%) of the total BBH merger rate. Finally, we discuss the broader implications of this formation channel for compact object populations, and its possible relevance to less asymmetric merger events such as GW200105 and GW200115, Comment: A&A Letters, in press

Rotation contributes to internal mixing processes and observed variability in massive stars. A significant number of binary stars are not in strict synchronous rotation, including all eccentric systems. This leads to a tidally induced and time-variable differential rotation structure. We present a method for exploring the rotation structure of asynchronously rotating binaries. We solve the equations of motion of a 3D grid of volume elements located above the rigidly rotating core in the presence of gravitational, centrifugal, Coriolis, gas pressure and viscous forces to obtain the angular velocity as a function of the three spatial coordinates and time. We find that the induced rotation structure and its temporal variability depend on the degree of departure from synchronicity. In eccentric systems, the structure changes over the orbital cycle with maximum amplitudes occurring potentially at orbital phases other than periastron passage. We discuss the possible role of the time-dependent tidal flows in enhancing the mixing efficiency and speculate that, in this context, slowly rotating asynchronous binaries could have more efficient mixing than the analogous more rapidly rotating but tidally locked systems. We find that some observed nitrogen abundances depend on the orbital inclination, which, if real, would imply an inhomogeneous chemical distribution over the stellar surface or that tidally induced spectral line variability, which is strongest near the equator, affects the abundance determinations. Our models predict that, neglecting other angular momentum transfer mechanisms, a pronounced initial differential rotation structure converges toward average uniform rotation on the viscous timescale. We suggest that by taking into account the processes that are triggered by asynchronous rotation, a broader perspective of binary star structure, evolution and variability may be gleaned., Comment: 21 pages, 18 figures

The detection of the unusually heavy binary neutron star merger GW190425 marked a stark contrast to the mass distribution from known Galactic pulsars in double neutron star binaries and gravitational-wave source GW170817. We suggest here a formation channel for heavy binary neutron stars and light black hole - neutron star binaries in which massive helium stars, which had their hydrogen envelope removed during a common envelope phase, remain compact and avoid mass transfer onto the neutron star companion, possibly avoiding pulsar recycling. We present three-dimensional simulations of the supernova explosion of the massive stripped helium star and follow the mass fallback evolution and the subsequent accretion onto the neutron star companion. We find that fallback leads to significant mass growth in the newly formed neutron star. This can explain the formation of heavy binary neutron star systems such as GW190425, as well as predict the assembly of light black hole - neutron star systems such as GW200115. This formation avenue is consistent with the observed mass-eccentricity correlation of binary neutron stars in the Milky Way. Finally, avoiding mass transfer suggests an unusually long spin-period population of pulsar binaries in our Galaxy., Comment: 7 main pages, 3 main figures, plus appendices. Accepted for publication in ApJ Letters

We perform for the first time a 3D hydrodynamics simulation of the evolution of the last minutes pre-collapse of the oxygen shell of a fast-rotating massive star. This star has an initial mass of 38 M$_\odot$, a metallicity of $\sim$1/50 Z$_\odot$, an initial rotational velocity of 600 km s$^{-1}$, and experiences chemically homogeneous evolution. It has a silicon- and oxygen-rich (Si/O) convective layer at (4.7-17)$\times 10^{8}$ cm, where oxygen-shell burning takes place. The power spectrum analysis of the turbulent velocity indicates the dominance of the large-scale mode ($\ell \sim 3$), which has also been seen in non-rotating stars that have a wide Si/O layer. Spiral arm structures of density and silicon-enriched material produced by oxygen-shell burning appear in the equatorial plane of the Si/O shell. Non-axisymmetric, large-scale ($m \le 3$) modes are dominant in these structures. The spiral arm structures have not been identified in previous non-rotating 3D pre-supernova models. Governed by such a convection pattern, the angle-averaged specific angular momentum becomes constant in the Si/O convective layer, which is not considered in spherically symmetrical stellar evolution models. Such spiral arms and constant specific angular momentum might affect the ensuing explosion or implosion of the star., Comment: 6 pages, 6 figures, accepted for publication in MNRAS Letters

As potential progenitors of several exotic phenomena including gravitational wave sources, magnetic stars, and Be stars, close massive binary systems probe a crucial area of the parameter space in massive star evolution. Despite the importance of these systems, large uncertainties regarding the nature and efficiency of the internal mixing mechanisms still exist. In this work, we aim to provide robust observational constraints on the internal mixing processes by spectroscopically analyzing a sample of three massive overcontact binaries at different metallicities. Using optical phase-resolved spectroscopic data, we perform an atmosphere analysis using more traditional 1D techniques and using state-of-the-art 3D techniques. We compare and contrast the assumptions and results of each technique and investigate how the assumptions affect the final derived atmospheric parameters. We find that in all three cases, both components of system are highly overluminous indicating either efficient internal mixing of helium or previous non-conservative mass transfer. However, we do not find strong evidence of helium or CNO surface abundance changes usually associated with mixing. Additionally, we find that in unequal mass systems, the measured effective temperature and luminosity of the less massive component places it very close to the more massive component on the Hertzsprung-Russell diagram. These results were obtained independently using both of the techniques mentioned above, which suggests that these measurements are robust. The observed discrepancies between the temperature and the surface abundance measurements when compared to theoretical expectations indicate that unaccounted for additional physical mechanisms may be at play., Comment: Accepted for publication in A&A; 28 pages, 20 figures

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., Comment: A&A Accepted, 14 pages, many colorful figures

The majority of close massive binary stars with initial periods of a few days experience a contact phase, in which both stars overflow their Roche lobes simultaneously. We perform the first dedicated study of the evolution of massive contact binaries and provide a comprehensive prediction of their observed properties. We compute 2790 detailed binary models for the Large and Small Magellanic Clouds each, assuming a conservative mass transfer. The initial parameter space for both grids span total masses from 20 to 80$\,\textrm{M}_{\odot}$, orbital periods of 0.6 to 2 days and mass ratios of 0.6 to 1.0. We find that models that remain in contact over nuclear timescales evolve towards equal masses, echoing the mass ratios of their observed counterparts. Ultimately, the fate of our nuclear-timescale models is to merge on the main sequence. Our predicted period-mass ratio distributions of O-type contact binaries are similar for both galaxies, and we expect 10 such systems together in both Magellanic Clouds. While we can largely reproduce the observed distribution, we over-estimate the population of equal-mass contact binaries. This situation is somewhat remedied if we also account for binaries that are approaching contact. Our theoretical distributions work particularly well for contact binaries with periods $<$2 days and total masses $\lessapprox45\,\textrm{M}_{\odot}$. We expect stellar winds, non-conservative mass transfer and envelope inflation to have played a role in the formation of the more massive and longer-period contact binaries., Comment: 23 pages, 11 figures, accepted for publication in MNRAS

While magnetic fields have long been considered to be important for the evolution of magnetic non-degenerate stars and compact stars, it has become clear in recent years that actually all of the stars are deeply affected. This is particularly true regarding their internal angular momentum distribution, but magnetic fields may also influence internal mixing processes and even the fate of the star. We propose a new framework for stellar evolution simulations, in which the interplay between magnetic field, rotation, mass loss, and changes in the stellar density and temperature distributions are treated self-consistently. For average large-scale stellar magnetic fields which are symmetric to the axis of rotation of the star, we derive 1D evolution equations for the toroidal and poloidal components from the mean-field MHD equation by applying Alfven's theorem, and a conservative form of the angular momentum transfer due to the Lorentz force is formulated. We implement our formalism into a numerical stellar evolution code and simulate the magneto-rotational evolution of 1.5 M$_\odot$ stars. The Lorentz force aided by the $\Omega$ effect imposes torsional Alfven waves propagating through the magnetized medium, leading to near-rigid rotation within the Alfven timescale. Our models with different initial spins and B-fields can reproduce the main observed properties of Ap/Bp stars. Calculations continued to the red-giant regime show a pronounced core-envelope coupling, which reproduces the core and surface rotation periods determined by asteroseismic observations., Comment: Accepted for publication in A&A. 29 pages, 13+figures, 1 Table

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 merger of two white dwarfs (WD) is a natural outcome from the evolution of many binary stars. Recently, a WD merger product, IRAS 00500+6713, was identified. IRAS 00500+6713 consists of a central star embedded in a circular nebula. The analysis of the optical spectrum of the central star revealed that it is hot, hydrogen and helium free, and drives an extremely fast wind with a record breaking speed. The nebula is visible in infrared and in the [O III] line images. No nebula spectroscopy was obtained prior to our observations. Here we report the first deep X-ray imaging spectroscopic observations of IRAS 00500+6713. Both the central star and the nebula are detected in X-rays, heralding the WD merger products as a new distinct type of strong X-ray sources. Low-resolution X-ray spectra reveal large neon, magnesium, silicon, and sulfur enrichment of the central star and the nebula. We conclude that IRAS 00500+6713 resulted from a merger of an ONe and a CO WD, which supports earlier suggestion for a super-Chandrasekhar mass of this object. X-ray analysis indicates that the merger was associated with an episode of carbon burning and possibly accompanied by a SN Iax. In X-rays, we observe the point source associated with the merger product while the surrounding diffuse nebula is a supernova remnant. IRAS 00500+6713 will likely terminate its evolution with another peculiar Type I supernova, where the final core collapse to a neutron star might be induced by electron captures., Comment: accepted by A&A Letters, 9 pages including appendix

We analyze the properties of 42 rapidly rotating, low metallicity, quasi-chemically homogeneously evolving stellar models in the mass range between 4 and 45 $\,\mathrm{M}_\odot$ at the time of core collapse. Such models were proposed as progenitors for both superluminous supernovae (SLSNe) and long duration gamma-ray bursts (lGRBs), and the Type Ic-BL supernovae (SNe) that are associated with them. Our findings suggest that whether these models produce a magnetar driven SLSN explosion or a near-critically rotating black hole (BH) is not a monotonic function of the initial mass. Rather, their explodability varies non-monotonically depending on the late core evolution, once chemical homogeneity is broken. Using different explodability criteria we find that our models have a clear preference to produce SLSNe at lower masses, and lGRBs at higher masses; but find several exceptions, expecting lGRBs to form from stars as low as 10 $\,\mathrm{M}_\odot$, and SLSNe with progenitors as massive as 30 $\,\mathrm{M}_\odot$. In general, our models reproduce the predicted angular momenta, ejecta masses and magnetic field strengths at core collapse inferred for SLSNe and lGRBs, and suggest significant interaction with their circumstellar medium, particularly for explosions with low ejecta mass., Comment: Accepted for publication in ApJ. 17 pages, 11 figures, 2 tables

Much difficulty has so far prevented the emergence of a consistent scenario for the origin of Type Ib and Ic supernovae (SNe). Here, we follow a heuristic approach by examining the fate of helium stars in the mass range 4 to 12Msun, which presumably form in interacting binaries. The helium stars are evolved using stellar wind mass loss rates that agree with observations, and which reproduce the observed luminosity range of galactic WR stars, leading to stellar masses at core collapse in the range 3-5.5Msun. We then explode these models adopting an explosion energy proportional to the ejecta mass, roughly consistent with theoretical predictions. We impose a fixed 56Ni mass and strong mixing. The SN radiation from 3 to 100d is computed self-consistently starting from the input stellar models using the time-dependent non-local thermodynamic equilibrium radiative-transfer code CMFGEN. By design, our fiducial models yield similar light curves, with a rise time of ~20d and a peak luminosity of ~10^42.2erg/s, in line with representative SNe Ibc. The less massive progenitors retain a He-rich envelope and reproduce the color, line widths, and line strengths of a representative sample of SNe Ib, while stellar winds remove most of the helium in more massive progenitors, whose spectra match typical SNe Ic in detail. The transition between the predicted Ib-like and Ic-like spectra is continuous, but it is sharp, such that the resulting models essentially form a dichotomy. Further models computed with varying explosion energy, 56Ni mass, and long-term power injection from the remnant show that a moderate variation of these parameters can reproduce much of the diversity of SNe Ibc. We conclude that stars stripped by a binary companion can account for the vast majority of ordinary SNe Ib and Ic, and that stellar wind mass loss is the key to remove the helium envelope in SN Ic progenitors. [abridged], Comment: Accepted for publication in A&A

Steady-state circulation currents are predicted in tidally deformed binary stars, which are believed to be progenitors of double black-hole merger events. This work aims to quantitatively characterise the steady-state circulation currents in components of a tidally locked binary system and to explore the effects of such currents on numerical models. Previous results describing the circulation in a single rotating star and a binary star are used to deduce a new prescription for the internal circulation in tidally locked binaries. We explore the effect of this prescription numerically for binary systems with primary masses between 25 and 100 solar masses. When comparing circulation velocities in the radial direction for the single rotating star and binary star, it is found that the average circulation velocity in the binary star may be described as an enhancement to the circulation velocity in a single rotating star. This velocity enhancement is a simple function depending on the masses of the binary components and amounts to a factor of approximately two when the components have equal masses. It is found that the ehancement causes the formation of double helium stars through efficient mixing to occur for systems with higher initial orbital periods, lower primary masses and lower mass ratios, compared to the standard circulation scenario. Taking into account appropriate distributions for primary mass, initial period and mass ratio, models with enhanced mixing predict 2.4 times more double helium stars being produced in the parameter space than models without. We conclude that the effects of companion-induced circulation have strong implications for the formation of close binary black holes. Not only do the predicted detection rates increase but double black-hole systems with mass ratios as low as 0.8 may be formed when companion-induced circulation is taken into account., Comment: 17 pages, 10 figures, accepted for publication in A&A

Massive and very massive stars can play important roles in stellar populations by ejecting strong stellar winds and exploding in energetic phenomena. It is therefore imperative that their behavior be properly accounted for in synthetic model populations. We present nine grids of stellar evolutionary model sequences, together with finely resolved interpolated sequences and synthetic populations, of stars with 9-500 Msun and with metallicities ranging from Galactic metallicity down to 1/250 Zsun. The stellar models were computed with the Bonn evolutionary code with consistent physical ingredients, and covering core hydrogen- and core helium-burning phases. The interpolation and population synthesis were performed with our newly developed routine 'synStars'. Eight of the grids represent slowly rotating massive stars with a normal or classical evolutionary path, while one grid represents fast-rotating, chemically homogeneously evolving models. The grids contain data on stellar wind properties such as estimated wind velocity and kinetic energy of the wind, as well as common stellar parameters such as mass, radius, surface temperature, luminosity, mass-loss rate, and surface abundances of 34 isotopes. We also provide estimates of the helium and carbon-oxygen core mass for calculating the mass of stellar remnants. The Bonn Optimized Stellar Tracks (BoOST) project is published as simple tables that include stellar models, interpolated tracks, and synthetic populations. Covering the broadest mass and metallicity range of any published massive star evolutionary model sets to date, BoOST is ideal for further scientific applications such as star formation studies in both low- and high-redshift galaxies., Comment: Accepted by A&A

During the first three observing runs of the Advanced gravitational-wave detector network, the LIGO/Virgo collaboration detected several black hole binary (BHBH) mergers. As the population of detected BHBH mergers grows, it will become possible to constrain different channels for their formation. Here we consider the chemically homogeneous evolution (CHE) channel in close binaries, by performing population synthesis simulations that combine realistic binary models with detailed cosmological calculations of the chemical and star-formation history of the Universe. This allows us to constrain population properties, as well as cosmological and aLIGO/aVirgo detection rates of BHBH mergers formed through this pathway. We predict a BHBH merger rate at redshift zero of $5.8 \, \textrm{Gpc}^{-3} \textrm{yr}^{-1}$ through the CHE channel, to be compared with aLIGO/aVirgo's measured rate of ${53.2}_{-28.2}^{+55.8} \, \text{Gpc}^{-3} \text{yr}^{-1}$, and find that eventual merger systems have BH masses in the range $17 - 43 \, \textrm{M}_{\odot}$ below the pair-instability supernova (PISN) gap, and $>124 \, \textrm{M}_{\odot}$ above the PISN gap. We investigate effects of momentum kicks during black hole formation, and calculate cosmological and magnitude limited PISN rates. We also study the effects of high-redshift deviations in the star formation rate. We find that momentum kicks tend to increase delay times of BHBH systems, and our magnitude limited PISN rate estimates indicate that current deep surveys should be able to detect such events. Lastly, we find that our cosmological merger rate estimates change by at most $\sim 8\%$ for mild deviations of the star formation rate in the early Universe, and by up to $\sim 40\%$ for extreme deviations., Comment: 20 pages, 18 figures, 4 tables. References added, text clarifications added, results and figures unchanged

Type Ia supernovae (SNe Ia) are manifestations of stars deficient of hydrogen and helium disrupting in a thermonuclear runaway. While explosions of carbon-oxygen white dwarfs are thought to account for the majority of events, part of the observed diversity may be due to varied progenitor channels. We demonstrate that helium stars with masses between $\sim$1.8 and 2.5 M$_{\odot}$ may evolve into highly degenerate, near-Chandrasekhar mass cores with helium-free envelopes that subsequently ignite carbon and oxygen explosively at densities $\sim(1.8-5.9)\times 10^{9}$g cm$^{-3}$. This happens either due to core growth from shell burning (when the core has a hybrid CO/NeO composition), or following ignition of residual carbon triggered by exothermic electron captures on $^{24}$Mg (for a NeOMg-dominated composition). We argue that the resulting thermonuclear runaways is likely to prevent core collapse, leading to the complete disruption of the star. The available nuclear energy at the onset of explosive oxygen burning suffices to create ejecta with a kinetic energy of $\sim$10$^{51}$ erg, as in typical SNe Ia. Conversely, if these runaways result in partial disruptions, the corresponding transients would resemble SN Iax events similar to SN 2002cx. If helium stars in this mass range indeed explode as SNe Ia, then the frequency of events would be comparable to the observed SN Ib/c rates, thereby sufficing to account for the majority of SNe Ia in star-forming galaxies., Comment: accepted in Astronomy & Astrophysics; v2 matches the accepted version. Comments are welcome

Star clusters are the building blocks of galaxies. They are composed of stars of nearly equal age and chemical composition, allowing us to use them as chronometers and as testbeds for gauging stellar evolution. It has become clear recently that massive stars are formed preferentially in close binaries, in which mass transfer will drastically change the evolution of the stars. This is expected to leave a significant imprint in the distribution of cluster stars in the Hertzsprung-Russell diagram. Our results, based on a dense model grid of more than 50,000 detailed binary-evolution calculations, indeed show several distinct, coeval main-sequence (MS) components, most notably an extended MS turnoff region, and a group of near-critical rotating stars that is spread over a large luminosity range on the red side of the classical MS. We comprehensively demonstrate the time evolution of the features in an animation, and we derive analytic expressions to describe these features. We find quantitative agreement with results based on recent photometric and spectroscopic observations. We conclude that while other factors may also be at play, binary evolution has a major impact on the MS morphology of young star clusters., Comment: 12 pages, 5 figures, 1 animation, published in ApJL, 2020, 888, L12

Be stars are rapidly rotating B main sequence stars, which show line emission due to an outflowing disc. By studying the evolution of rotating single star models, we can assess their contribution to the observed Be star populations. We identify the main effects which are responsible for single stars to approach critical rotation as functions of initial mass and metallicity, and predict the properties of populations of rotating single stars. We perform population synthesis with single star models of initial masses ranging between 3 and 30 solar masses, initial equatorial rotation velocities between 0 and 600 kms$^{-1}$ at compositions representing the Milky Way, Large and Small Magellanic Clouds. These models include efficient core-envelope coupling mediated by internal magnetic fields and correspond to the maximum efficiency of Be star production. We predict Be star fractions and the positions of fast rotating stars in the colour-magnitude diagram. We identify stellar wind mass-loss and the convective core mass fraction as the key parameters which determine the time dependance of the stellar rotation rates. Using empirical distributions of initial rotational velocities,our single star models can reproduce the trends observed in Be star fractions with mass and metallicity. However,they fail to produce a significant number of stars rotating very close to critical. We also find that rapidly rotating Be stars in the Magellanic Clouds should have significant surface nitrogen enrichments, which may be in conflict with abundance determinations of Be stars. Single star evolution may explain the high number of Be stars if 70 to 80% of critical rotationwould be sufficient to produce the Be phenomenon. However even in this case, the unexplained presence of many Be stars far below the cluster turn-off indicates the importance of the binary channel for Be star production., Comment: in press

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

The massive O4.5 V + O5.5 V binary VFTS 352 in the Tarantula nebula is one of the shortest-period and most massive overcontact binaries known. Recent theoretical studies indicate that some of these systems could ultimately lead to the formation of gravitational waves via black hole binary mergers through the chemically homogeneous evolution pathway. By analyzing ultraviolet-optical phase-resolved spectroscopic data, we aim to constrain atmospheric and wind properties that could be later used to confront theoretical predictions from binary evolution. In particular, surface abundances are powerful diagnostics of the evolutionary status, mass transfer and the internal mixing processes. From a set of 32 VLT/FLAMES visual and 8 HST/COS ultraviolet spectra, we used spectral disentangling to separate the primary and secondary components. Using a genetic algorithm wrapped around the NLTE model atmosphere and spectral synthesis code FASTWIND, we perform an 11-parameter optimization to derive the atmospheric and wind parameters of both components, including the surface abundances of He, C, N, O and Si. We find that both components are hotter than expected compared to single-star evolutionary models indicating that additional mixing processes may be at play. However the derived chemical abundances do not show significant indications of mixing when adopting baseline values typical for the system environment., Comment: 24 pages, 4 tables, 12 figures, accepted for publication in ApJ

Context. The evolution of massive stars is strongly influenced by internal mixing processes such as semiconvection, convective core overshooting, and rotationally induced mixing. None of these is currently well constrained. Aims. We investigate models for massive stars in the Small Magellanic Cloud (SMC). We aim to constrain the various mixing efficiencies by comparing model results to observations. Methods. We use the stellar evolution code MESA to compute more than 60 grids of detailed evolutionary models for stars with initial masses of 9 to 100 Msol, in each grid assuming different combinations of mixing efficiencies. Results. We find that for most of the combinations of the mixing efficiencies, models in a wide mass range spend core-helium burning either only as blue supergiants, or only as red supergiants. The latter case corresponds to models that maintain a shallow slope of the hydrogen/helium (H/He) gradient separating the core and the envelope of the models. Only a small part of the mixing parameter space leads to models that produce a significant number of blue and red supergiants, which both exist abundantly in the SMC. Interestingly, these models contain steep H/He gradients, as is required to understand the hot, hydrogen-rich Wolf-Rayet stars in the SMC. We find that unless it is very fast, rotation has a limited effect on the H/He profiles in our models. Conclusions. While we use specific implementations of the considered mixing processes, they comprehensively probe the two first order structural parameters: the core mass and the H/He gradient in the core-envelope interface. Our results imply that, in massive stars, the H/He gradients above the helium cores become very steep. Our model grids can be tested with future observational surveys of the massive stars in the SMC, and thereby help to considerably reduce the uncertainties in models of massive star evolution., Comment: Accepted for publication in Astronomy and Astrophysics. 16 pages, 14 figures

We investigate how the proximity of supergiant donor stars to the Eddington-limit, and their advanced evolutionary stage, may influence the evolution of massive and ultra-luminous X-ray binaries with supergiant donor stars (SGXBs and ULXs). We construct models of massive stars with different internal hydrogen/helium gradients and different hydrogen-rich envelope masses, and expose them to slow mass loss to probe the response of the stellar radius. In addition, we compute the corresponding Roche-lobe overflow mass-transfer evolution with our detailed binary stellar evolution code, approximating the compact objects as point masses. We find that a hydrogen/helium gradient in the layers beneath the surface, as it is likely present in the well-studied donor stars of observed SGBXs, can enable nuclear timescale mass-transfer in SGXBs with a BH or a NS accretor, even for mass ratios in excess of 20. In our binary evolution models, the donor stars rapidly decrease their thermal equilibrium radius and can therefore cope with the inevitably strong orbital contraction imposed by the high mass ratio. Our results open a new perspective for understanding the large number of Galactic SGXBs, and their almost complete absence in the SMC. They may also offer a way to obtain more ULX systems, to find nuclear timescale mass-transfer in ULX systems even with neutron star accretors, and shed new light on the origin of the strong B-field in these neutron stars., Comment: 23 pages, 21 figures, we are thankful for any comments an this draft

Recent studies show that more than 70% of massive stars do not evolve as effectively single stars, but as members of interacting binary systems. The evolution of these stars is thus strongly altered compared to similar but isolated objects. We investigate the occurrence of parsec-scale mid-infrared nebulae around early-type stars. If they exist over a wide range of stellar properties, one possible overarching explanation is non-conservative mass transfer in binary interactions, or stellar mergers. For ~3850 stars (all OBA stars in the Bright Star Catalogue [BSC], Be stars, BeXRBs, and Be+sdO systems), we visually inspect WISE 22 $\mu$m images. Based on nebular shape and relative position, we distinguish five categories: offset bow shocks structurally aligned with the stellar space velocity, unaligned offset bow shocks, and centered, unresolved, and not classified nebulae. In the BSC, we find that 28%, 13%, and 0.4% of all O, B, and A stars, respectively, possess associated infrared (IR) nebulae. Additionally, 34/234 Be stars, 4/72 BeXRBs, and 3/17 Be+sdO systems are associated with IR nebulae. Aligned or unaligned bow shocks result from high relative velocities between star and interstellar medium (ISM) that are dominated by the star or the ISM, respectively. About 13% of the centered nebulae could be bow shocks seen head- or tail-on. For the rest, the data disfavor explanations as remains of parental disks, supernova remnants of a previous companion, and dust production in stellar winds. The existence of centered nebulae also at high Galactic latitudes strongly limits the global risk of coincidental alignments with condensations in the ISM. Mass loss during binary evolution seems a viable mechanism for the formation of at least some of these nebulae. In total, about 29% of the IR nebulae (2% of all OBA stars in the BSC) may find their explanation in the context of binary evolution., Comment: 21 pages, 5 tables, 6 figures, accepted for publication in A&A

We model the late evolution and mass loss history of rapidly rotating Wolf-Rayet stars in the mass range $5\,\rm{M}_{\odot}\dots 100\,\rm{M}_{\odot}$. We find that quasi-chemically homogeneously evolving single stars computed with enhanced mixing retain very little or no helium and are compatible with Type\,Ic supernovae. The more efficient removal of core angular momentum and the expected smaller compact object mass in our lower mass models lead to core spins in the range suggested for magnetar driven superluminous supernovae. Our more massive models retain larger specific core angular momenta, expected for long-duration gamma-ray bursts in the collapsar scenario. Due to the absence of a significant He envelope, the rapidly increasing neutrino emission after core helium exhaustion leads to an accelerated contraction of the whole star, inducing a strong spin-up, and centrifugally driven mass loss at rates of up to $10^{-2}\,\rm{M}_{\odot}~\rm{yr^{-1}}$ in the last years to decades before core collapse. Since the angular momentum transport in our lower mass models enhances the envelope spin-up, they show the largest relative amounts of centrifugally enforced mass loss, i.e., up to 25\% of the expected ejecta mass. Our most massive models evolve into the pulsational pair-instability regime. We would thus expect signatures of interaction with a C/O-rich circumstellar medium for Type~Ic superluminous supernovae with ejecta masses below $\sim 10\,\rm{M}_{\odot}$ and for the most massive engine-driven explosions with ejecta masses above $\sim 30\,\rm{M}_{\odot}$. Signs of such interaction should be observable at early epochs of the supernova explosion, and may be related to bumps observed in the light curves of superluminous supernovae, or to the massive circumstellar CO-shell proposed for Type~Ic superluminous supernova Gaia16apd., Comment: 20 pages, 15 figures, 2 tables Accepted for publication in ApJ