Showing 1–50 of 363 results for author: Langer, N

Whispering in the dark: faint X-ray emission from black holes with OB star companions

Authors: Koushik Sen, Ileyk El Mellah, Norbert Langer, Xiao-Tian Xu, Martin Quast, Daniel Pauli

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

Submitted 12 June, 2024; originally announced June 2024.

Comments: Submitted to A&A on May 31, 2024. Comments and feedback are very welcome

An absence of binary companions to Wolf-Rayet stars in the Small Magellanic Cloud: implications for mass loss and black hole masses at low metallicity

Authors: A. Schootemeijer, T. Shenar, N. Langer, N. Grin, H. Sana, G. Gräfener C. Schürmann, C. Wang, X. -T. Xu

Abstract: In order to predict the black hole mass distributions at high redshift, we need to understand whether very massive single stars ($M>40$ M$_\odot$) at low metallicity $Z$ lose their hydrogen-rich envelopes, like their metal-rich counterparts, or whether a binary companion is required to achieve this. To test this, we undertake a deep spectroscopic search for binary companions of the seven apparentl… ▽ More In order to predict the black hole mass distributions at high redshift, we need to understand whether very massive single stars ($M>40$ M$_\odot$) at low metallicity $Z$ lose their hydrogen-rich envelopes, like their metal-rich counterparts, or whether a binary companion is required to achieve this. To test this, we undertake a deep spectroscopic search for binary companions of the seven apparently single Wolf-Rayet (WR) stars in the Small Magellanic Cloud (SMC; $Z \simeq 1/5 Z_\odot$). For each of them, we acquired six high-quality VLT-UVES spectra spread over 1.5 years. By using the narrow N V lines in these spectra, we monitor radial velocity (RV) variations to search for binary motion. We find low RV variations between 6 and 23 km/s for the seven WR stars, with a median standard deviation of $5$ km/s. Our Monte Carlo simulations imply probabilities below ~5% for any of our target WR stars to have a binary companion more massive than ~5 M$_\odot$ at orbital periods of less than a year. We estimate that the probability that all our target WR stars have companions with orbital periods shorter than 10 yr is below ~10$^{-5}$, and argue that the observed modest RV variations may originate from intrinsic atmosphere or wind variability. Our findings imply that metal-poor massive stars born with $M \gtrsim 40$ M$_\odot$ can lose most of their hydrogen-rich envelopes via stellar winds or eruptive mass loss, which strongly constrains their initial mass - black hole mass relation. We also identify two of our seven target stars (SMC AB1 and SMC AB11) as runaway stars with a peculiar radial velocity of ~80 km/s. Moreover, with all five previously detected WR binaries in the SMC exhibiting orbital periods of below 20 d, a puzzling absence of intermediate-to-long-period WR binaries has emerged, with strong implications for the outcome of massive binary interaction at low metallicity. △ Less

Submitted 3 June, 2024; originally announced June 2024.

Comments: Accepted for publication in Astronomy & Astrophysics

Boron depletion in Galactic early B-type stars reveals two different main sequence star populations

Authors: Harim Jin, Norbert Langer, Daniel J. Lennon, Charles R. Proffitt

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

Submitted 28 May, 2024; originally announced May 2024.

Rotational synchronisation of B-type binaries in 30 Doradus

Authors: D. J. Lennon, P. L. Dufton, J. I. Villaseñor, N. Langer, C. J. Evans, H. Sana, W. D. Taylor

Abstract: The spin evolution of stars in close binary systems can be strongly affected by tides. We investigate the rotational synchronisation of the stellar components for 69 SB1 systems and 14 SB2 B-type systems in the 30 Doradus region of the Large Magellanic Cloud using observations from the VFTS and BBC surveys. Their orbital periods range from a few to a few hundred days, while estimated primary masse… ▽ More The spin evolution of stars in close binary systems can be strongly affected by tides. We investigate the rotational synchronisation of the stellar components for 69 SB1 systems and 14 SB2 B-type systems in the 30 Doradus region of the Large Magellanic Cloud using observations from the VFTS and BBC surveys. Their orbital periods range from a few to a few hundred days, while estimated primary masses for these systems are in the range 5-20 Msun with mass ratio ranges of approximately q=0.03-0.5 and q=0.6-1.0 for the SB1 and SB2 systems, respectively. Projected rotational velocities of the stellar components have been compared with their synchronous velocities derived from the orbital periods. We find that effectively all systems with orbital period of more than 10 days must be asynchronous, whilst all the systems with periods of less than 3 days are likely synchronised. In terms of the stellar fractional radius (r), our results imply that all systems with r<0.1 are asynchronous, with those having r>0.2 probably being synchronised. For the apparently synchronised systems our results are more consistent with synchronisation at the mean orbital angular velocity rather than with that at periastron. △ Less

Submitted 27 May, 2024; originally announced May 2024.

Comments: 13 pages, 5 figures. submitted to A&A

Exploring the borderline between stable mass transfer and mergers in close binary evolution

Authors: Christoph Schürmann, Norbert Langer

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

Submitted 12 April, 2024; originally announced April 2024.

Comments: 19 pages

Analytic approximations for massive close post-mass transfer binary systems

Authors: Christoph Schürmann, Norbert Langer, Joana A. Kramer, Pablo Marchant, Chen Wang, Koushik Sen

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

Submitted 12 April, 2024; originally announced April 2024.

Comments: 10 pages

A sequence of Type Ib, IIb, II-L, and II-P supernovae from binary-star progenitors of varying initial separation

Authors: Luc Dessart, Claudia P. Gutierrez, Andrea Ercolino, Harim Jin, Norbert Langer

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

Submitted 20 February, 2024; originally announced February 2024.

Comments: Submitted to A&A on Dec 22nd, 2023

Journal ref: A&A 685, A169 (2024)

Boron Abundances in Early B Dwarfs of the Galactic Open Cluster NGC 3293

Authors: Charles R. Proffitt, Harim Jin, Simone Daflon, Daniel J. Lennon, Norbert Langer, Katia Cunha, Talawanda Monroe

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

Submitted 2 May, 2024; v1 submitted 23 January, 2024; originally announced January 2024.

Comments: 18 pages,7 figures, Submitted to AAS Journals

Evidence for stellar mergers of evolved massive binaries: blue supergiants in the Large Magellanic Cloud

Authors: Athira Menon, Andrea Ercolino, Miguel A. Urbaneja, Daniel J. Lennon, Artemio Herrero, Ryosuke Hirai, Norbert Langer, Abel Schootemeijer, Emmanouil Chatzopoulos, Juhan Frank, Sagiv Shiber

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

Submitted 9 November, 2023; originally announced November 2023.

Sulfur isotope ratios in the Large Magellanic Cloud

Authors: Y. Gong, C. Henkel, K. M. Menten, C. -H. R. Chen, Z. Y. Zhang, Y. T. Yan, A. Weiss, N. Langer, J. Z. Wang, R. Q. Mao, X. D. Tang, W. Yang, Y. P. Ao, M. Wang

Abstract: Sulfur isotope ratios have emerged as a promising tool for tracing stellar nucleosynthesis, quantifying stellar populations, and investigating the chemical evolution of galaxies. While extensively studied in the Milky Way, in extragalactic environments they remain largely unexplored. We focus on investigating the sulfur isotope ratios in the Large Magellanic Cloud (LMC) to gain insights into sulfu… ▽ More Sulfur isotope ratios have emerged as a promising tool for tracing stellar nucleosynthesis, quantifying stellar populations, and investigating the chemical evolution of galaxies. While extensively studied in the Milky Way, in extragalactic environments they remain largely unexplored. We focus on investigating the sulfur isotope ratios in the Large Magellanic Cloud (LMC) to gain insights into sulfur enrichment in this nearby system and to establish benchmarks for such ratios in metal-poor galaxies. We conducted pointed observations of CS and its isotopologues toward N113, one of the most prominent star-formation regions in the LMC, utilizing the Atacama Pathfinder EXperiment 12~m telescope. We present the first robust detection of C$^{33}$S in the LMC by successfully identifying two C$^{33}$S transitions on a large scale of $\sim$5 pc. Our measurements result in an accurate determination of the $^{34}$S/$^{33}$S isotope ratio, which is 2.0$\pm$0.2. Our comparative analysis indicates that the $^{32}$S/$^{33}$S and $^{34}$S/$^{33}$S isotope ratios are about a factor of 2 lower in the LMC than in the Milky Way. Our findings suggest that the low $^{34}$S/$^{33}$S isotope ratio in the LMC can be attributed to a combination of the age effect, low metallicity, and star formation history. △ Less

Submitted 18 October, 2023; v1 submitted 26 September, 2023; originally announced September 2023.

Comments: 6 pages, 1 figures, 2 tables, accepted for publication in A&A, adjusted to the final version

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

The young massive SMC cluster NGC 330 seen by MUSE III. Stellar parameters and rotational velocities

Authors: J. Bodensteiner, H. Sana, P. L. Dufton, C. Wang, N. Langer, G. Banyard, L. Mahy, A. de Koter, S. E. de Mink, C. J. Evans, Y. Götberg, V. Hénault-Brunet, L. R. Patrick, F. R. N. Schneider

Abstract: The origin of initial rotation rates of stars, and how a star's surface rotational velocity changes during the evolution, either by internal angular momentum transport or due to interactions with a binary companion, remain open questions in stellar astrophysics. Here, we aim to derive the physical parameters and study the distribution of (projected) rotational velocities of B-type stars in the 35… ▽ More The origin of initial rotation rates of stars, and how a star's surface rotational velocity changes during the evolution, either by internal angular momentum transport or due to interactions with a binary companion, remain open questions in stellar astrophysics. Here, we aim to derive the physical parameters and study the distribution of (projected) rotational velocities of B-type stars in the 35 Myr-old, massive cluster NGC 330 in the Small Magellanic Cloud. NGC 330 is in an age range where the number of post-interaction binaries is predicted to be high near the cluster turnoff (TO). We develop a simultaneous photometric and spectroscopic grid-fitting method adjusting atmosphere models on multi-band Hubble Space Telescope photometry and Multi Unit Spectroscopic Explorer spectroscopy. This allows us to homogeneously constrain the physical parameters of over 250 B and Be stars, brighter than mF814W = 18.8 mag. The rotational velocities of Be stars in NGC 330 are significantly higher than the ones of B stars. The rotational velocities vary as a function of the star's position in the color-magnitude diagram, qualitatively following predictions of binary population synthesis. A comparison to younger clusters shows that stars in NGC 330 rotate more rapidly on average. The rotational velocities of the 35 Myr old population in NGC 330 quantitatively agree with predictions for a stellar population that underwent significant binary interactions: the bulk of the B stars could be single stars or primaries in pre-interaction binaries. The rapidly spinning Be stars could be mass and angular momentum gainers in previous interactions, while those Be stars close to the TO may be spun-up single stars. The slowly rotating, apparently single stars above the TO could be merger products. The different vsini-characteristics of NGC 330 compared to younger populations can be understood in this framework. △ Less

Submitted 28 August, 2023; originally announced August 2023.

Comments: 18 pages (incl. appendix), 15 figures, 3 tables, accepted for publication in A&A

Journal ref: A&A 680, A32 (2023)

A massive helium star with a sufficiently strong magnetic field to form a magnetar

Authors: T. Shenar, G. A. Wade, P. Marchant, S. Bagnulo, J. Bodensteiner, D. M. Bowman, A. Gilkis, N. Langer, A. -N. Chene, L. Oskinova, T. Van Reeth, H. Sana, N. St-Louis, A. Soares de Oliveira, H. Todt, S. Toonen

Abstract: Magnetars are highly magnetized neutron stars; their formation mechanism is unknown. Hot helium-rich stars with spectra dominated by emission lines are known as Wolf-Rayet stars. We observe the binary system HD 45166 using spectropolarimetry, finding that it contains a Wolf-Rayet star with a mass of 2 solar masses and a magnetic field of 43 kilogauss. Stellar evolution calculations indicate that t… ▽ More Magnetars are highly magnetized neutron stars; their formation mechanism is unknown. Hot helium-rich stars with spectra dominated by emission lines are known as Wolf-Rayet stars. We observe the binary system HD 45166 using spectropolarimetry, finding that it contains a Wolf-Rayet star with a mass of 2 solar masses and a magnetic field of 43 kilogauss. Stellar evolution calculations indicate that this component will explode as a type Ib or IIb supernova, and the strong magnetic field favors a magnetar remnant. We propose that the magnatized Wolf-Rayet star formed by the merger of two lower mass helium stars. △ Less

Submitted 16 August, 2023; originally announced August 2023.

Comments: Published in Science on the 18 August 2023. Radial velocities, spectra, and software available in: https://zenodo.org/record/8042656 ESO press release: www.eso.org/public/news/eso2313

Interacting supernovae from wide massive binary systems

Authors: Andrea Ercolino, Harim Jin, Norbert Langer, Luc Dessart

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

Submitted 18 March, 2024; v1 submitted 3 August, 2023; originally announced August 2023.

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

The B-type Binaries Characterisation Programme II. VFTS 291: A stripped star from a recent mass transfer phase

Authors: J. I. Villaseñor, D. J. Lennon, A. Picco, T. Shenar, P. Marchant, N. Langer, P. L. Dufton, F. Nardini, C. J. Evans, J. Bodensteiner, S. E. de Mink, Y. Götberg, I. Soszyński, W. D. Taylor, H. Sana

Abstract: Recent studies of massive binaries with putative black hole companions have uncovered a phase of binary evolution that has not been observed before, featuring a bloated stripped star that very recently ceased transferring mass to a main-sequence companion. In this study, we focus on the candidate system VFTS 291, a binary with an orbital period of 108 d and a high semi-amplitude velocity (… ▽ More Recent studies of massive binaries with putative black hole companions have uncovered a phase of binary evolution that has not been observed before, featuring a bloated stripped star that very recently ceased transferring mass to a main-sequence companion. In this study, we focus on the candidate system VFTS 291, a binary with an orbital period of 108 d and a high semi-amplitude velocity ($K_1=93.7\pm0.2$ km s$^{-1}$). Through our analysis of the disentangled spectra of the two components, together with dynamical and evolutionary arguments, we identify a narrow-lined star of ~1.5-2.5 $M_\odot$ dominating the spectrum, and an early B-type main-sequence companion of $13.2\pm1.5$ $M_\odot$. The low mass of the narrow-lined star, and the high mass ratio, suggest that VFTS 291 is a post-mass-transfer system, with the narrow-lined star being bloated and stripped of its hydrogen-rich envelope, sharing many similarities with other recently discovered stripped stars. Our finding is supported by our detailed binary evolution models, which indicate that the system can be well explained by an initial configuration consisting of an 8.1 $M_\odot$ primary with an 8 $M_\odot$ companion in a 7 d orbital period. While some open questions remain, particularly concerning the surface helium enrichment of the stripped star and the rotational velocity of the companion, we expect that high-resolution spectroscopy may help reconcile our estimates with theory. Our study highlights the importance of multi-epoch spectroscopic surveys to identify and characterize binary interaction products, and provides important insights into the evolution of massive binary stars. △ Less

Submitted 11 September, 2023; v1 submitted 15 July, 2023; originally announced July 2023.

Comments: 25 pages, 17 figures. Published by MNRAS, Advance Access

A helium-burning white dwarf binary as a supersoft X-ray source

Authors: J. Greiner, C. Maitra, F. Haberl, R. Willer, J. M. Burgess, N. Langer, J. Bodensteiner, D. A. H. Buckley, I. M. Monageng, A. Udalski, H. Ritter, K. Werner, P. Maggi, R. Jayaraman, R. Vanderspek

Abstract: Type Ia supernovae are cosmic distance indicators, and the main source of iron in the Universe, but their formation paths are still debated. Several dozen supersoft X-ray sources, in which a white dwarf accretes hydrogen-rich matter from a non-degenerate donor star, have been observed and suggested as Type Ia supernovae progenitors. However, observational evidence for hydrogen, which is expected t… ▽ More Type Ia supernovae are cosmic distance indicators, and the main source of iron in the Universe, but their formation paths are still debated. Several dozen supersoft X-ray sources, in which a white dwarf accretes hydrogen-rich matter from a non-degenerate donor star, have been observed and suggested as Type Ia supernovae progenitors. However, observational evidence for hydrogen, which is expected to be stripped off the donor star during the supernova explosion, is lacking. Helium-accreting white dwarfs, which would circumvent this problem, have been predicted for more than 30 years, also including their appearance as supersoft X-ray sources, but have so far escaped detection. Here we report a supersoft X-ray source with an accretion disk whose optical spectrum is completely dominated by helium, suggesting that the donor star is hydrogen-free. We interpret the luminous and supersoft X-rays as due to helium burning near the surface of the accreting white dwarf. The properties of our system provides evidence for extended pathways towards Chandrasekhar mass explosions based on helium accretion, in particular for stable burning in white dwarfs at lower accretion rates than expected so far. This may allow to recover the population of the sub-energetic so-called Type Iax supernovae, up to 30% of all Type Ia supernovae, within this scenario. △ Less

Submitted 23 March, 2023; originally announced March 2023.

Journal ref: Nature 615 (2023), 605

Modeling contact binaries, II. The effect of energy transfer

Authors: Matthias Fabry, Pablo Marchant, Norbert Langer, Hugues Sana

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

Submitted 16 March, 2023; originally announced March 2023.

Comments: 11 pages, 5 figures. Accepted to A&A

Journal ref: A&A 672, A175 (2023)

A model of anisotropic winds from rotating stars for evolutionary calculations

Authors: Ben Hastings, Norbert Langer, Joachim Puls

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

Submitted 21 February, 2023; originally announced February 2023.

Comments: Accepted for publication in A&A. 12 pages

Journal ref: A&A 672, A60 (2023)

Reverse Algols and hydrogen-rich Wolf-Rayet stars from very massive binaries

Authors: K. Sen, N. Langer, D. Pauli, G. Gräfener, A. Schootemeijer, H. Sana, T. Shenar, L. Mahy, C. Wang

Abstract: Massive star feedback affects the evolution of galaxies, where the most massive stars may have the largest impact. The majority of massive stars are born as members of close binary systems. Here, we investigate detailed evolutionary models of very massive binaries (30$\dots$90$M_{\odot}$) with Large Magellanic Cloud (LMC) metallicity. We identify four effects defying the conventional knowledge of… ▽ More Massive star feedback affects the evolution of galaxies, where the most massive stars may have the largest impact. The majority of massive stars are born as members of close binary systems. Here, we investigate detailed evolutionary models of very massive binaries (30$\dots$90$M_{\odot}$) with Large Magellanic Cloud (LMC) metallicity. We identify four effects defying the conventional knowledge of binary evolution, which are all related to the proximity of the models to the Eddington limit. We find that the majority of systems undergo mass transfer during core hydrogen burning. During the ensuing nuclear timescale evolution, many mass donors remain more massive than their companions (``reverse Algols''), and nuclear timescale mass transfer may be interrupted or absent altogether. Furthermore, due to the elevated luminosity-to-mass ratio, many of the core-hydrogen burning donors may develop Wolf-Rayet-type winds, at luminosities where single stars would not. We identify observational counterparts of very massive reverse Algol binaries in the LMC, and discuss their contribution to the observed hydrogen-rich Wolf-Rayet stars. We argue that an understanding of very massive Algol systems is key to predicting the advanced evolution of very massive binaries, including their ability to evolve into observable gravitational wave sources. △ Less

Submitted 9 February, 2023; originally announced February 2023.

Comments: 15 pages, 5 figures, proposed for acceptance in A&A. Animation at https://zenodo.org/record/4068428#.Y-SrOtLMKV4

Journal ref: A&A 672, A198 (2023)

Hubble Space Telescope survey of Magellanic Cloud star clusters. Photometry and astrometry of 113 clusters and early results

Authors: A. P. Milone, G. Cordoni, A. F. Marino, F. D'Antona, A. Bellini, M. Di Criscienzo, E. Dondoglio, E. P. Lagioia, N. Langer, M. V. Legnardi, M. Libralato, H. Baumgardt, M. Bettinelli, Y. Cavecchi, R. de Grijs, L. Deng, B. Hastings, C. Li, A. Mohandasan, A. Renzini, E. Vesperini, C. Wang, T. Ziliotto, M. Carlos, G. Costa , et al. (7 additional authors not shown)

Abstract: In the past years, we have undertaken an extensive investigation of LMC and SMC star clusters based on HST data. We present photometry and astrometry of stars in 101 fields observed with the WFC/ACS, UVIS/WFC3 and NIR/WFC3 cameras. These fields comprise 113 star clusters. We provide differential-reddening maps and illustrate various scientific outcomes that arise from the early inspection of the p… ▽ More In the past years, we have undertaken an extensive investigation of LMC and SMC star clusters based on HST data. We present photometry and astrometry of stars in 101 fields observed with the WFC/ACS, UVIS/WFC3 and NIR/WFC3 cameras. These fields comprise 113 star clusters. We provide differential-reddening maps and illustrate various scientific outcomes that arise from the early inspection of the photometric catalogs. In particular, we provide new insights on the extended main-sequence turn-off (eMSTO) phenomenon: i) We detected eMSTOs in two clusters, KMHK361 and NGC265, which had no previous evidence of multiple populations. This finding corroborates the conclusion that the eMSTO is a widespread phenomenon among clusters younger than ~2 Gyr. ii) The homogeneous color-magnitude diagrams (CMDs) of 19 LMC clusters reveal that the distribution of stars along the eMSTO depends on cluster age. iii) We discovered a new feature along the eMSTO of NGC1783, which consists of a distinct group of stars going on the red side of the eMSTO in CMDs composed of ultraviolet filters. Furthermore, we derived the proper motions of stars in the fields of view of clusters with multi-epoch images. Proper motions allowed us to separate the bulk of bright field stars from cluster members and investigate the internal kinematics of stellar populations in various LMC and SMC fields. As an example, we analyze the field around NGC346 to disentangle the motions of its stellar populations, including NGC364 and BS90, young and pre-MS stars in the star-forming region associated with NGC346, and young and old field stellar populations of the SMC. Based on these results and the fields around five additional clusters, we find that young SMC stars exhibit elongated proper-motion distributions that point toward the LMC, thus bringing new evidence for a kinematic connection between the LMC and SMC. △ Less

Submitted 15 December, 2022; originally announced December 2022.

Comments: 37 pages, 27 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 672, A161 (2023)

Three dimensional magnetorotational core-collapse supernova explosions of a 39 solar mass progenitor star

Authors: Jade Powell, Bernhard Mueller, David R. Aguilera-Dena, Norbert Langer

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

Submitted 4 May, 2023; v1 submitted 30 November, 2022; originally announced December 2022.

The initial spin distribution of B-type stars revealed by the split main sequences of young star clusters

Authors: Chen Wang, Ben Hastings, Abel Schootemeijer, Norbert Langer, Selma E. de Mink, Julia Bodensteiner, Antonino Milone, Stephen Justham, Pablo Marchant

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

Submitted 7 December, 2022; v1 submitted 28 November, 2022; originally announced November 2022.

Comments: 18 pages, 14 figures. Accepted for publication in A&A

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

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

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

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

Submitted 24 November, 2022; originally announced November 2022.

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

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

Searching for compact objects in the single-lined spectroscopic binaries of the young Galactic cluster NGC 6231

Authors: G. Banyard, L. Mahy, H. Sana, J. Bodensteiner, J. I. Villaseńor, K. Sen, N. Langer, S. de Mink, A. Picco, T. Shenar

Abstract: Recent evolutionary computations predict that a few percent of massive OB stars in binary systems should have a dormant BH companion. Despite several reported X-ray quiet OB+BH systems over the last couple of years, finding them with certainty remains challenging. These have great importance as they can be gravitational wave (GW) source progenitors, and are landmark systems in constraining superno… ▽ More Recent evolutionary computations predict that a few percent of massive OB stars in binary systems should have a dormant BH companion. Despite several reported X-ray quiet OB+BH systems over the last couple of years, finding them with certainty remains challenging. These have great importance as they can be gravitational wave (GW) source progenitors, and are landmark systems in constraining supernova kick physics. This work aims to characterise the hidden companions to the single-lined spectroscopic binaries (SB1s) in the B star population of the young open Galactic cluster NGC 6231 to find candidate systems for harbouring compact object companions. With the orbital solutions for each SB1 previously constrained, we applied Fourier spectral disentangling to multi-epoch optical VLT/FLAMES spectra of each target to extract a potential signature of a faint companion, and to identify newly disentangled double-lined spectroscopic binaries (SB2s). For targets where the disentangling does not reveal any signature of a stellar companion, we performed atmospheric and evolutionary modelling on the primary to obtain constraints on the unseen companion. Seven newly classified SB2 systems with mass ratios down to near 0.1 were identified. From the remaining targets, for which no faint companion could be extracted from the spectra, four are found to have companion masses in the predicted mass ranges of neutron stars (NSes) and BHes. Two of these have companion masses between 1 and 3.5 $M_{\odot}$, making them potential hosts of NSes (or lower mass main sequence stars). The other two are between 2.5 to 8 $M_{\odot}$ and 1.6 and 26 $M_{\odot}$, respectively, and so are identified as candidates for harbouring BH companions. However, unambiguous identification of these systems as X-ray quiet compact object harbouring binaries requires follow up observations. △ Less

Submitted 13 October, 2022; originally announced October 2022.

Comments: Submitted to A&A

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

Massive stars in metal-poor dwarf galaxies are often extreme rotators

Authors: Abel Schootemeijer, Danny J. Lennon, Miriam Garcia, Norbert Langer, Ben Hastings, Christoph Schürmann

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

Submitted 4 October, 2022; originally announced October 2022.

Comments: Conference proceedings for a talk in IAU Symposium 361: Massive Stars Near and Far, Ballyconnell, Ireland, 9-13 May 2022

The Nature of Unseen Companions in Massive Single-Line Spectroscopic Binaries

Authors: Hugues Sana, Michael Abdul-Masih, Gareth Banyard, Julia Bodensteiner, Dominic M. Bowman, Karan Dsilva, C. Eldridge, Matthias Fabry, Abigail J. Frost, Calum Hawcroft, Soetkin Janssens, Laurent Mahy, Pablo Marchant, Norbert Langer, Timothy Van Reeth, Koushik Sen, Tomer Shenar

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

Submitted 13 September, 2022; originally announced September 2022.

Comments: To appear in the proceedings of IAUS361: Massive stars near and far; 6 pages, 1 figure

A census of OBe stars in nearby metal-poor dwarf galaxies reveals a high fraction of extreme rotators

Authors: A. Schootemeijer, D. J. Lennon, M. Garcia, N. Langer, B. Hastings, C. Schuermann

Abstract: The Early Universe, together with many nearby dwarf galaxies, is deficient in heavy elements. The evolution of massive stars in such environments is thought to be affected by rotation. Extreme rotators amongst them tend to form decretion disks and manifest themselves as OBe stars. We use a combination of U B, GAIA, Spitzer, and Hubble Space Telescope photometry to identify the complete populations… ▽ More The Early Universe, together with many nearby dwarf galaxies, is deficient in heavy elements. The evolution of massive stars in such environments is thought to be affected by rotation. Extreme rotators amongst them tend to form decretion disks and manifest themselves as OBe stars. We use a combination of U B, GAIA, Spitzer, and Hubble Space Telescope photometry to identify the complete populations of massive OBe stars - one hundred to thousands in number - in five nearby dwarf galaxies. This allows us to derive the galaxy-wide fractions of main sequence stars that are OBe stars (f_OBe), and how it depends on absolute magnitude, mass, and metallicity (Z). We find f_OBe = 0.22 in the Large Magellanic Cloud (0.5 Z_Sun), increasing to f_OBe = 0.31 in the Small Magellanic Cloud (0.2 Z_Sun). In the so far unexplored metallicity regime below 0.2 Z_Sun, in Holmberg I, Holmberg II, and Sextans A, we also obtain high OBe star fractions of 0.27, 0.27, and 0.27, respectively. These high OBe star fractions, and the strong contribution in the stellar mass range which dominates the production of supernovae, shed new light on the formation channel of OBe stars, as well as on the preference of long-duration gamma-ray bursts and superluminous supernovae to occur in metal-poor galaxies. △ Less

Submitted 11 September, 2022; originally announced September 2022.

Comments: 19 pages, 23 figures - accepted for publication in Astronomy & Astrophysics

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

Hunting for red supergiant binaries: UVIT photometry of the SMC

Authors: L. R. Patrick, D. Thilker, D. J. Lennon, L. Bianchi, A. Schootemeijer, R. Dorda, N. Langer, I. Negueruela

Abstract: We present UVIT/Astrosat UV photometry of the RSG population of the Small Cloud galaxy (SMC). As RSGs are extremely faint in the far-UV, these observations directly probe potential companion stars. From a sample of 861 SMC RSGs, we find 88 have detections at far-UV wavelengths: a clear signature of binarity. Stellar parameters are determined for both components, which allows us to study - for the… ▽ More We present UVIT/Astrosat UV photometry of the RSG population of the Small Cloud galaxy (SMC). As RSGs are extremely faint in the far-UV, these observations directly probe potential companion stars. From a sample of 861 SMC RSGs, we find 88 have detections at far-UV wavelengths: a clear signature of binarity. Stellar parameters are determined for both components, which allows us to study - for the first time - the mass-ratio (q) distribution of RSG binary systems. We find a flat mass-ratio distribution best describes the observations up to M{RSG}~15 Msun. We account for our main observing bias (i.e. the limiting magnitude of the UVIT survey) to determine the intrinsic RSG binary fraction of 18.8+/-1.5 %, for mass-ratios in the range 0.3 < q < 1.0 and orbital periods approximately in the range 3 < log P[ days] < 8. △ Less

Submitted 9 September, 2022; originally announced September 2022.

Comments: IAU Symposium 361 contribution

Open questions in massive star research

Authors: N. Langer

Abstract: In discussing open question in the field of massive stars, I consider their evolution from birth to death. After touching upon massive star formation, which may be bi-modal and not lead to a zero-age main sequence at the highest masses, I consider the consequences of massive stars being close to their Eddington limit. Then, when discussing the effects of a binary companion, I highlight the importa… ▽ More In discussing open question in the field of massive stars, I consider their evolution from birth to death. After touching upon massive star formation, which may be bi-modal and not lead to a zero-age main sequence at the highest masses, I consider the consequences of massive stars being close to their Eddington limit. Then, when discussing the effects of a binary companion, I highlight the importance of massive Algols and contact binaries for understanding the consequences of mass transfer, and the role of binaries in forming Wolf-Rayet stars. Finally, a discussion on pair instability supernovae and of superluminous supernovae is provided. △ Less

Submitted 9 September, 2022; originally announced September 2022.

Comments: to appear in Proceedings for IAU Symposium 361: Massive Stars Near and Far, held in Ballyconnell, Ireland, 9-13 May 2022, N. St-Louis, J. S. Vink & J. Mackey, eds

A synthetic population of Wolf-Rayet stars in the LMC based on detailed single and binary star evolution models

Authors: D. Pauli, N. Langer, D. R. Aguilera-Dena, C. Wang, P. Marchant

Abstract: Without doubt, mass transfer in close binary systems contributes to the populations of Wolf-Rayet (WR) stars in the Milky Way and the Magellanic Clouds. However, the binary formation channel is so far not well explored. We want to remedy this by exploring large grids of detailed binary and single star evolution models computed with the publicly available MESA code, for a metallicity appropriate fo… ▽ More Without doubt, mass transfer in close binary systems contributes to the populations of Wolf-Rayet (WR) stars in the Milky Way and the Magellanic Clouds. However, the binary formation channel is so far not well explored. We want to remedy this by exploring large grids of detailed binary and single star evolution models computed with the publicly available MESA code, for a metallicity appropriate for the Large Magellanic Cloud (LMC). The binary models are calculated through Roche-lobe overflow and mass transfer, until the initially more massive star exhausts helium in its core. We distinguish models of WR and helium stars based on the estimated stellar wind optical depth. We use these models to build a synthetic WR population, assuming constant star formation. Our models can reproduce the WR population of the LMC to significant detail, including the number and luminosity functions of the main WR subtypes. We find that for binary fractions of 100% (50%), all LMC WR stars below $10^6\,L_{\odot}$ ($10^{5.7}\,L_{\odot}$) are stripped binary mass donors. We also identify several insightful mismatches. With a single star fraction of 50\%, our models produce too many yellow supergiants, calling either for a larger initial binary fraction, or for enhanced mass-loss near the Humphreys-Davidson limit. Our models predict more long-period WR binaries than observed, arguably due to an observational bias towards short periods. Our models also underpredict the shortest-period WR binaries, which may have implications for understanding the progenitors of double black hole mergers. The fraction of binary produced WR stars may be larger than often assumed, and outline the risk to mis-calibrate stellar physics when only single star models are used to reproduce the observed WR stars. △ Less

Submitted 22 August, 2022; originally announced August 2022.

Comments: 15 pages + 13 pages appendix, 14 figures, 2 tables

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

The spins of stripped B stars support magnetic internal angular momentum transport

Authors: C. Schürmann, N. Langer, X. Xu, C. Wang

Abstract: In order to predict the spins of stellar remnants we need to understand the evolution of the internal rotation of stars, and to identify at which stage the rotation of the contracting cores of evolved stars decouples from their expanding envelopes. The donor stars of mass transferring binaries lose almost their entire envelope and may thus offer a direct view on their core rotation. After the mass… ▽ More In order to predict the spins of stellar remnants we need to understand the evolution of the internal rotation of stars, and to identify at which stage the rotation of the contracting cores of evolved stars decouples from their expanding envelopes. The donor stars of mass transferring binaries lose almost their entire envelope and may thus offer a direct view on their core rotation. After the mass transfer event they contract and fade rapidly, although they are well observable when caught in the short-lived B-star phase. The B-type primary of the galactic binary system LB-1, which was originally suggested to contain a massive black hole, is nicely explained as a stripped star accompanied by a fainter Be star. The narrow absorption lines in the primary's spectrum signify extremely slow rotation, atypical of B-type main-sequence stars. Here we investigate the evolution of mass donors in generic grids of detailed binary evolution models, where both stars include differential rotation, internal angular momentum transport, and spin-orbit coupling. Whereas the mass gainers are typically spun-up during the mass transfer, we find that the spins of the stripped donor models depend sensitively on the employed mechanism for internal angular momentum transport. Purely hydrodynamic transport cannot explain the observed slow rotation, while models including magnetic angular momentum transport are able to reproduce the observed rotation of LB-1 and similar stars, independent of the initial rotation rate. In such models the spin of the white dwarfs that emerge at the end of the evolution is independent of the mass stripping. We find evidence that the mass transfer in LB-1 was moderately non-conservative. △ Less

Submitted 25 November, 2022; v1 submitted 5 August, 2022; originally announced August 2022.

Comments: 18 pages, 13 figures. Accepted in Astronomy and Astrophysics

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

Identifying quiescent compact objects in massive Galactic single-lined spectroscopic binaries

Authors: L. Mahy, H. Sana, T. Shenar, K. Sen, N. Langer, P. Marchant, M. Abdul-Masih, G. Banyard, J. Bodensteiner, D. M. Bowman, K. Dsilva, M. Fabry, C. Hawcroft, S. Janssens, T. Van Reeth, C. Eldridge

Abstract: Aims. To search for these rare objects, we study 32 Galactic O-type stars that were reported as SB1s in the literature. In our sample we include Cyg X-1, which is known to host an accreting stellar-mass BH, and HD 74194, a supergiant fast X-ray transient, in order to validate our methodology. The final goal is to characterise the nature of the unseen companions to determine if they are MS stars, s… ▽ More Aims. To search for these rare objects, we study 32 Galactic O-type stars that were reported as SB1s in the literature. In our sample we include Cyg X-1, which is known to host an accreting stellar-mass BH, and HD 74194, a supergiant fast X-ray transient, in order to validate our methodology. The final goal is to characterise the nature of the unseen companions to determine if they are MS stars, stripped helium stars, triples, or compact objects such as neutron stars or stellar-mass BHs. Methods. After measuring radial velocities and deriving orbital solutions for all the systems in our sample, we performed spectral disentangling to extract putative signatures of faint secondary companions from the composite spectra. We derived stellar parameters for the visible stars and estimated the mass ranges of the secondary stars using the binary mass function. Variability observed in the photometric TESS light curves was also searched for indications of the presence of putative companions, degenerate or not. Results. In 17 of the 32 systems reported as SB1s, we extract secondary signatures, down to mass ratios of ~0.15. For the 17 newly detected double-lined spectroscopic binaries (SB2s), we derive physical properties of the individual components and discuss why they have not been detected as such before. Among the remaining systems, we identify nine systems with possible NS or low-mass MS companions. For Cyg X-1 and HD 130298, we are not able to extract any signatures for the companions, and the minimum masses of their companions are estimated to be about 7Msun. Our simulations show that secondaries with such a mass should be detectable from our dataset, no matter their nature: MS stars, stripped helium stars or even triples. While this is expected for Cyg X-1, confirming our methodology, our simulations also strongly suggest that HD 130298 could be another candidate to host a stellar-mass BH. △ Less

Submitted 15 July, 2022; originally announced July 2022.

Comments: 20 pages (31 pages of Appendix), 19 figures, A&A in press. The abstract is shorter than in the paper, and some figures, and tables have been cut in the paper

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

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

Abstract: Stellar-mass black holes are the final remnants of stars born with more than 15 solar masses. Billions are expected to reside in the Local Group, yet only few are known, mostly detected through X-rays emitted as they accrete material from a companion star. Here, we report on VFTS 243: a massive X-ray faint binary in the Large Magellanic Cloud. With an orbital period of 10.4-d, it comprises an O-ty… ▽ More Stellar-mass black holes are the final remnants of stars born with more than 15 solar masses. Billions are expected to reside in the Local Group, yet only few are known, mostly detected through X-rays emitted as they accrete material from a companion star. Here, we report on VFTS 243: a massive X-ray faint binary in the Large Magellanic Cloud. With an orbital period of 10.4-d, it comprises an O-type star of 25 solar masses and an unseen companion of at least nine solar masses. Our spectral analysis excludes a non-degenerate companion at a 5-sigma confidence level. The minimum companion mass implies that it is a black hole. No other X-ray quiet black hole is unambiguously known outside our Galaxy. The (near-)circular orbit and kinematics of VFTS 243 imply that the collapse of the progenitor into a black hole was associated with little or no ejected material or black-hole kick. Identifying such unique binaries substantially impacts the predicted rates of gravitational-wave detections and properties of core-collapse supernovae across the Cosmos. △ Less

Submitted 15 July, 2022; originally announced July 2022.

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

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

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

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

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

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

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

Red supergiant stars in binary systems. I. Identification and characterisation in the Small Magellanic Cloud from the UVIT ultraviolet imaging survey

Authors: L. R. Patrick, D. Thilker, D. J. Lennon, L. Bianchi, A. Schootemeijer, R. Dorda, N. Langer, I. Negueruela

Abstract: We aim to identify and characterise binary systems containing red supergiant (RSG) stars in the Small Magellanic Cloud (SMC) using a newly available ultraviolet (UV) point source catalogue obtained using the Ultraviolet Imaging Telescope (UVIT) on board AstroSat. We select a sample of 560 SMC RSGs based on photometric and spectroscopic observations at optical wavelengths and cross-match this with… ▽ More We aim to identify and characterise binary systems containing red supergiant (RSG) stars in the Small Magellanic Cloud (SMC) using a newly available ultraviolet (UV) point source catalogue obtained using the Ultraviolet Imaging Telescope (UVIT) on board AstroSat. We select a sample of 560 SMC RSGs based on photometric and spectroscopic observations at optical wavelengths and cross-match this with the far-UV point source catalogue using the UVIT F172M filter, finding 88 matches down to m$_{F172M}$=20.3 ABmag, which we interpret as hot companions to the RSGs. Stellar parameters (luminosities, effective temperatures and masses) for both components in all 88 binary systems are determined and we find mass distributions in the ranges 6.1 to 22.3 Solar masses for RSGs and 3.7 to 15.6 Solar masses for their companions. The most massive RSG binary system in the SMC has a combined mass of 32 $\pm$4 M$_\odot$, with a mass ratio (q) of 0.92. By simulating observing biases, we find an intrinsic multipliciy fraction of 18.8 $\pm$ 1.5% for mass ratios in the range 0.3 < q < 1.0 and orbital periods approximately in the range 3 < log P [days] < 8. By comparing our results with those of a similar mass on the main-sequence, we determine the fraction of single stars to be ~20% and argue that the orbital period distribution declines rapidly beyond log P ~ 3.5. We study the mass-ratio distribution of RSG binary systems and find that a uniform distribution best describes the data below 14 M$_\odot$. Above 15 M$_\odot$, we find a lack of high mass-ratio systems. △ Less

Submitted 25 April, 2022; originally announced April 2022.

Comments: Accepted for publication in MNRAS; 17 pages

Stripped-envelope stars in different metallicity environments. II. Type I supernovae and compact remnants

Authors: David R. Aguilera-Dena, Bernhard Müller, John Antoniadis, Norbert Langer, Luc Dessart, Alejandro Vigna-Gómez, Sung-Chul Yoon

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

Submitted 31 March, 2022; originally announced April 2022.

Comments: 30 pages, 20 figures, submitted to Astronomy & Astrophysics

Journal ref: A&A 671, A134 (2023)

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

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

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

Submitted 22 February, 2022; originally announced February 2022.

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

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

Stellar mergers as the origin of the blue main-sequence band in young star clusters

Authors: Chen Wang, Norbert Langer, Abel Schootemeijer, Antonino Milone, Ben Hastings, Xiao-Tian Xu, Julia Bodensteiner, Hugues Sana, Norberto Castro, D. J. Lennon, Pablo Marchant, A. de Koter, Selma E. de Mink

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

Submitted 11 February, 2022; originally announced February 2022.

Comments: Published in Nature Astronomy online. Link to the paper:https://rdcu.be/cGLpH

Thermonuclear and Electron-Capture Supernovae from Stripped-Envelope Stars

Authors: Savvas Chanlaridis, John Antoniadis, David R. Aguilera-Dena, Götz Gräfener, Norbert Langer, Nikolaos Stergioulas

Abstract: (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… ▽ More (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 $ρ_{\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}(ρ_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. △ Less

Submitted 7 October, 2022; v1 submitted 3 January, 2022; originally announced January 2022.

Comments: Accepted in Astronomy & Astrophysics; v2 matches the accepted version. Comments are welcome

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

Stripped-Envelope Stars in Different Metallicity Environments I. Evolutionary Phases, Classification and Populations

Authors: David R. Aguilera-Dena, Norbert Langer, John Antoniadis, Daniel Pauli, Luc Dessart, Alejandro Vigna-Gómez, Götz Gräfener, Sung-Chul Yoon

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

Submitted 7 February, 2022; v1 submitted 13 December, 2021; originally announced December 2021.

Comments: 20 pages, 10 figures, accepted for publication Astronomy & Astrophysics

Journal ref: A&A 661, A60 (2022)

The VLT-FLAMES survey of massive stars. NGC2004#115: A triple system hosting a possible short period B+BH binary

Authors: D. J. Lennon, P. L. Dufton, J. I. Villaseñor, C. J. Evans, N. Langer, R. Saxton, I. M. Monageng, S. Toonen

Abstract: The star NGC2004#115 in the LMC, originally classified as an (SB1) Be spectroscopic binary, bears some morphological resemblance to the Galactic systems LB-1 and HR 6819, both of which are proposed as either Be+black hole (BH) or Be+stripped He-star systems. Two data-sets (ESO/VLT and SALT) of multi-epoch optical spectra of NGC 2004#115, separated by baseline of $\sim$20 years, lead us conclude it… ▽ More The star NGC2004#115 in the LMC, originally classified as an (SB1) Be spectroscopic binary, bears some morphological resemblance to the Galactic systems LB-1 and HR 6819, both of which are proposed as either Be+black hole (BH) or Be+stripped He-star systems. Two data-sets (ESO/VLT and SALT) of multi-epoch optical spectra of NGC 2004#115, separated by baseline of $\sim$20 years, lead us conclude it is a triple system hosting an inner binary with a period of 2.92 d, eccentricity $\sim$0.0 and mass function $\sim$0.07 $M_\odot$. The inner binary harbours a B-type star (the primary) with projected rotational velocity of 10km/s, and luminosity $\log L/L_\odot$=3.87, contributing $\sim$60% of the V-band light to the system. The secondary is not detected, while the tertiary, which contributes 40% of the light, is tentatively identified as a less luminous B-type star with high projected rotational velocity. No ellipsoidal light variability is detected, with stringent limits being set by MACHO and Gaia data. Assuming the primary to be a main sequence star yields a mass of 8.6$ M_\odot$, while the additional assumption of synchronous rotation constrains the inclination to be almost pole-on with i~9 degrees, implying the secondary is a BH with a mass of $\sim$25 $M_\odot$. A low mass stripped star with similar luminosity is ruled out as a potential solution as its mass implies a Roche radius that is substantially smaller than the stellar radius. The outer period likely exceeds 120 days and, while the disk-like emission is variable (it is almost absent in the SALT dataset), it may be associated with the inner binary rather than the rapidly rotating tertiary. XMM-Newton provides an upper limit of 5x$10^{33}$ ergs/s on the X-ray flux, consistent with, though not constraining of, the system hosting a quiescent B+BH binary. A number of caveats to this scenario are discussed in the paper. △ Less

Submitted 30 November, 2021; v1 submitted 23 November, 2021; originally announced November 2021.

Comments: 19 pages, 21 figures. Accepted for publication in Astronomy and Astrophysics Version 2

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

Uncovering astrometric black hole binaries with massive main-sequence companions with Gaia

Authors: S. Janssens, T. Shenar, H. Sana, S. Faigler, N. Langer, P. Marchant, T. Mazeh, C. Schürmann, S. Shahaf

Abstract: The hunt for compact objects is on. Rarely seen massive binaries with a compact object are a crucial phase in the evolution towards compact object mergers. In Gaia data release 3 (DR3), the first Gaia astrometric orbital solutions for binary sources will become available. We investigate how many black holes (BH) with massive main-sequence dwarf companions (OB+BH binaries) are expected to be detect… ▽ More The hunt for compact objects is on. Rarely seen massive binaries with a compact object are a crucial phase in the evolution towards compact object mergers. In Gaia data release 3 (DR3), the first Gaia astrometric orbital solutions for binary sources will become available. We investigate how many black holes (BH) with massive main-sequence dwarf companions (OB+BH binaries) are expected to be detected as binaries in DR3 and at the end of the nominal 5-yr mission (DR4). We estimate the fraction of identifiable OB+BH binaries and discuss the distributions of the masses of both components and the orbital periods. We study the impact of different BH-formation scenarios. Using tailored models for the massive star population, which assume a direct collapse and no kick upon BH formation (the fiducial case), we estimate the fraction of OB+BH systems that Gaia will detect as binaries. A distance distribution according to that of the second Alma Luminous Star catalogue (ALSII) is assumed. We investigate how many of the systems detected as binaries are identifiable as OB+BH binaries, using a method based on astrometric data. In the fiducial case we conservatively estimate that 77% of the OB+BH binaries in ALSII will be detected as binaries in DR3, of which 89% are identifiable as OB+BH binaries. This leads to a total of around 190 OB+BH binaries, a 20-fold increase in the known sample of OB+BH binaries, covering an uncharted parameter space of long-period binaries. The size and properties of the identifiable OB+BH population will contain crucial observational constraints to improve our understanding of BH formation. In DR4, the detected fraction will increase to 85%, of which 82% will be identifiable. Hence, an additional ~5 systems could be identified, which are expected to have either very short or long periods. The fractions become smaller for different BH-formation scenarios. (truncated) △ Less

Submitted 11 November, 2021; originally announced November 2021.

Comments: 22 pages, 18 figures, 2 tables, accepted for publication in A&A

Journal ref: A&A 658, A129 (2022)

Detailed models of interacting short-period massive binary stars

Authors: K. Sen, N. Langer, P. Marchant, A. Menon, S. E. de Mink, A. Schootemeijer, C. Schürmann, L. Mahy, B. Hastings, K. Nathaniel, H. Sana, C. Wang, X. T. Xu

Abstract: About a quarter of massive binary stars undergo mass transfer while both stars burn hydrogen at their cores, first on the thermal and then on the nuclear timescale. The nuclear timescale mass transfer leads to observable counterparts: the semi-detached so-called massive Algol binaries. However, comprehensive model predictions for these systems are sparse. We study them using a large grid of ~10,00… ▽ More About a quarter of massive binary stars undergo mass transfer while both stars burn hydrogen at their cores, first on the thermal and then on the nuclear timescale. The nuclear timescale mass transfer leads to observable counterparts: the semi-detached so-called massive Algol binaries. However, comprehensive model predictions for these systems are sparse. We study them using a large grid of ~10,000 detailed binary evolution models calculated with the stellar evolution code MESA, covering initial donor masses between 10-40 M$_{\odot}$ and initial orbital periods above 1.4 d, at a metallicity suitable for the Large Magellanic Cloud (LMC). Our models imply ~30, or ~3% of the ~1,000 core hydrogen burning O-star binaries in the LMC to be currently in the semi-detached phase. Our donor models are up to 25-times more luminous than single stars of identical mass and effective temperature, which agrees with the observed Algols. A comparison of our models with the observed orbital periods and mass ratios implies rather conservative mass transfer in some systems, while very inefficient one in others. This is generally well reproduced by our spin-dependent mass transfer algorithm, except for the lowest considered masses. The observations reflect the slow increase of the surface nitrogen enrichment of the donors during the semi-detached phase all the way to CNO equilibrium. We also investigate the properties of our models after core hydrogen depletion of the donor star, when these models correspond to Wolf-Rayet/helium+OB star binaries. A dedicated spectroscopic survey of massive Algol systems may allow to derive the dependence of the efficiency of thermal timescale mass transfer on the binary parameters, as well as the efficiency of semiconvective mixing in the stellar interior. This would be a crucial step towards reliable binary models up to the formation of supernovae and compact objects. △ Less

Submitted 9 December, 2021; v1 submitted 5 November, 2021; originally announced November 2021.

Comments: Accepted for publication in A&A (December 6, 2021)

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

Explodability fluctuations of massive stellar cores enable asymmetric compact object mergers such as GW190814

Authors: John Antoniadis, David R. Aguilera-Dena, Alejandro Vigna-Gómez, Michael Kramer, Norbert Langer, Bernhard Müller, Thomas M. Tauris, Chen Wang, Xiao-Tian Xu

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

Submitted 21 December, 2021; v1 submitted 4 October, 2021; originally announced October 2021.

Comments: A&A Letters, in press

Journal ref: A&A 657, L6 (2022)

The B-type Binaries Characterisation Programme I. Orbital solutions for the 30 Doradus population

Authors: J. I. Villaseñor, W. D. Taylor, C. J. Evans, O. H. Ramírez-Agudelo, H. Sana, L. A. Almeida, S. E. de Mink, P. L. Dufton, N. Langer

Abstract: We present results from the B-type Binaries Characterisation (BBC) programme, a multi-epoch spectroscopic study of 88 early B-type binary candidates in the 30 Doradus region of the Large Magellanic Cloud (LMC). From radial-velocity analysis of 29 observational epochs we confirm the binary status of 64 of our targets, comprising 50 SB1 and 14 SB2 B-type binaries. A further 20 systems (classified as… ▽ More We present results from the B-type Binaries Characterisation (BBC) programme, a multi-epoch spectroscopic study of 88 early B-type binary candidates in the 30 Doradus region of the Large Magellanic Cloud (LMC). From radial-velocity analysis of 29 observational epochs we confirm the binary status of 64 of our targets, comprising 50 SB1 and 14 SB2 B-type binaries. A further 20 systems (classified as SB1*) show clear signs of periodicity but with more tentative periods. Orbital solutions are presented for these 84 systems, providing the largest homogeneous sample to date of the binary properties of early B-type stars. Our derived orbital-period distribution is generally similar to those for samples of more massive (O-type) binaries in both the LMC and the Galaxy. This similarity with the properties of the more massive O-type binaries is important as early B-type stars are expected to account for the majority of core-collapse supernovae. Differences in the period distributions of the different samples start to increase above 4 d, and are also present between the earliest (B0-0.7) and later-type (B1-2.5) systems within the BBC sample, although further study is required to understand if this is an observational bias or a real physical effect. We have examined the semi-amplitude velocities and orbital periods of our sample to identify potential candidates that could hide compact companions. Comparing with probability distributions of finding black hole companions to OB-type stars from a recent theoretical study, we have found 16 binaries in the higher probability region that warrant further study. △ Less

Submitted 21 July, 2021; originally announced July 2021.

Comments: 42 pages (main body: 28 pages, supplementary material: 14 pages), 20 figures. Accepted for publication in MNRAS

Induced differential rotation and mixing in asynchronous binary stars

Authors: Gloria Koenigsberger, Edmundo Moreno, Norbert Langer

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

Submitted 19 July, 2021; originally announced July 2021.

Comments: 21 pages, 18 figures

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

Fallback supernova assembly of heavy binary neutron stars and light black hole-neutron star pairs and the common stellar ancestry of GW190425 and GW200115

Authors: Alejandro Vigna-Gómez, Sophie L. Schrøder, Enrico Ramirez-Ruiz, David R. Aguilera-Dena, Aldo Batta, Norbert Langer, Reinhold Willcox

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

Submitted 22 September, 2021; v1 submitted 23 June, 2021; originally announced June 2021.

Comments: 7 main pages, 3 main figures, plus appendices. Accepted for publication in ApJ Letters

A stringent upper limit on Be star fractions produced by binary interaction

Authors: B. Hastings, N. Langer, C. Wang, A. Schootemeijer, A. P. Milone

Abstract: Context. Binary evolution can result in fast-rotating stars, predicted to be observable as Be stars, through accretion of angular momentum during mass-transfer phases. Despite numerous observational evidence pointing to this possibly being the dominant Be formation channel, current models struggle to produce a satisfactory description of Be star populations. Aims. Given distinct uncertainties in… ▽ More Context. Binary evolution can result in fast-rotating stars, predicted to be observable as Be stars, through accretion of angular momentum during mass-transfer phases. Despite numerous observational evidence pointing to this possibly being the dominant Be formation channel, current models struggle to produce a satisfactory description of Be star populations. Aims. Given distinct uncertainties in detailed binary evolution calculations, we investigate a rigorous and model independent upper limit for the production of Be stars through binary interaction and aim to confront this limit with observations of Be stars in young star clusters. Methods. Using extreme assumptions, we calculate the number ratio of post-interaction to pre-interaction binary systems in a coeval population, which describes an upper limit to Be star formation through mass-transfer. A detailed comparison is made between our derived upper limit and relevant observations of Be stars, which allows us to probe several aspects of binary star physics. Results. We find that in coeval populations, binary interaction can at most account for one third of all main-sequence stars being Be stars. Near the cluster turn-off region, this limit appears to be realised in the clusters studied. Away from the turn-off, applying simple assumptions about which systems undergo unstable mass-transfer produces a good fit to the observed Be fraction as a function of mass. Conclusions. We find that assuming distinct physics, binary evolution alone can in principle match the high numbers of Be stars observed in open clusters. Whether the required binary physics is realised in nature remains to be investigated. △ Less

Submitted 23 June, 2021; originally announced June 2021.

Comments: Accepted for publication in A&A

A three-dimensional hydrodynamics simulation of oxygen-shell burning in the final evolution of a fast-rotating massive star

Authors: Takashi Yoshida, Tomoya Takiwaki, David R. Aguilera-Dena, Kei Kotake, Koh Takahashi, Ko Nakamura, Hideyuki Umeda, Norbert Langer

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

Submitted 30 June, 2021; v1 submitted 18 June, 2021; originally announced June 2021.

Comments: 6 pages, 6 figures, accepted for publication in MNRAS Letters

X-ray emission from BH+O star binaries expected to descend from the observed galactic WR+O binaries

Authors: K. Sen, X. -T. Xu, N. Langer, I. El Mellah, C. Schurmann, M. Quast

Abstract: In the Milky Way, $\sim$18 Wolf-Rayet+O (WR+O) binaries are known with estimates of their stellar and orbital parameters. Whereas black hole+O (BH+O) binaries are thought to evolve from the former, only one such system is known in the Milky Way. To resolve this disparity, it was suggested that upon core collapse, the WR stars receive large kicks such that most of the binaries are disrupted. We rea… ▽ More In the Milky Way, $\sim$18 Wolf-Rayet+O (WR+O) binaries are known with estimates of their stellar and orbital parameters. Whereas black hole+O (BH+O) binaries are thought to evolve from the former, only one such system is known in the Milky Way. To resolve this disparity, it was suggested that upon core collapse, the WR stars receive large kicks such that most of the binaries are disrupted. We reassess this issue, with emphasis on the uncertainty in the formation of an accretion disk around wind-accreting BHs in BH+O binaries, which is key to identifying such systems. We follow the methodology of previous work and apply an improved analytic criterion for the formation of an accretion disk around wind accreting BHs. We then use stellar models to predict the properties of the BH+O binaries which are expected to descend from the observed WR+O binaries, if the WR stars would form BHs without a natal kick. We find that disk formation depends sensitively on the O stars' wind velocity, the specific angular momentum carried by the wind, the efficiency of angular momentum accretion by the BH, and the spin of the BH. We show that the assumption of a low wind velocity may lead to predicting that most of the BH+O star binaries will have an extended X-ray bright period. However, this is not the case when typical wind velocities of O stars are considered. We find that a high spin of the BH can boost the duration of the X-ray active phase as well as the X-ray brightness during this phase, producing a strong bias for detecting high mass BH binaries in X-rays with high BH spin parameters. We conclude that large BH formation kicks are not required to understand the sparsity of X-ray bright BH+O stars in the Milky Way. Probing for a population of X-ray silent BH+O systems with alternative methods can inform us about BH kicks and the conditions for high energy emission from high mass BH binaries. (Abridged) △ Less

Submitted 2 June, 2021; originally announced June 2021.

Comments: 14 pages, 9 figures, accepted for publication in A&A

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

The young massive SMC cluster NGC 330 seen by MUSE II. Multiplicity properties of the massive-star population

Authors: J. Bodensteiner, H. Sana, C. Wang, N. Langer, L. Mahy, G. Banyard, A. de Koter, S. E. de Mink, C. J. Evans, Y. Götberg, L. R. Patrick, F. R. N. Schneider, F. Tramper

Abstract: Observations of massive stars in young open clusters (< ~8 Myr) have shown that a majority of them are in binary systems, most of which will interact during their life. Populations of massive stars older than ~20 Myr allow us to probe the outcome of such interactions after many systems have experienced mass and angular momentum transfer. Using multi-epoch integral-field spectroscopy, we investigat… ▽ More Observations of massive stars in young open clusters (< ~8 Myr) have shown that a majority of them are in binary systems, most of which will interact during their life. Populations of massive stars older than ~20 Myr allow us to probe the outcome of such interactions after many systems have experienced mass and angular momentum transfer. Using multi-epoch integral-field spectroscopy, we investigate the multiplicity properties of the massive-star population in NGC 330 (~40 Myr) in the Small Magellanic Cloud to search for imprints of stellar evolution on the multiplicity properties. From six epochs of VLT/MUSE observations supported by adaptive optics we extract spectra and measure radial velocities for stars brighter than F814W = 19. We identify single-lined spectroscopic binaries through significant RV variability as well as double-lined spectroscopic binaries, and quantify the observational biases for binary detection. The observed spectroscopic binary fraction is 13.2+/-2.0 %. Considering period and mass ratio ranges from log(P)=0.15-3.5, and q = 0.1-1.0, and a representative set of orbital parameter distributions, we find a bias-corrected close binary fraction of 34 +8 -7 %. This seems to decline for the fainter stars, which indicates either that the close binary fraction drops in the B-type domain, or that the period distribution becomes more heavily weighted towards longer orbital periods. Both fractions vary strongly in different regions of the color-magnitude diagram which probably reveals the imprint of the binary history of different groups of stars. We provide the first homogeneous RV study of a large sample of B-type stars at a low metallicity. The overall bias-corrected close binary fraction of B stars in NGC 330 is lower than the one reported for younger Galactic and LMC clusters. More data are needed to establish whether this result from an age or a metallicty effect. △ Less

Submitted 27 April, 2021; originally announced April 2021.

Comments: 19 page (incl. appendix), 10 figures, 4 tables, accepted for publication in A&A

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