It has to be cool: Supergiant progenitors of binary black hole mergers from common-envelope evolution

Common-envelope (CE) evolution in massive binary systems is thought to be one of the most promising channels for the formation of compact binary mergers. In the case of merging binary black holes (BBHs), the essential CE phase takes place at a stage when the first BH is already formed and the companion star expands as a supergiant. We aim to decipher the kinds of BH binaries with supergiant companions that could potentially evolve through and survive a CE phase. To this end, we compute envelope binding energies from detailed massive stellar models at different evolutionary stages and metallicities. We make multiple physically extreme choices of assumptions that favor easier CE ejection as well as account for recent advancements in mass-transfer stability criteria. We find that even with the most optimistic assumptions, a successful CE ejection in BH binaries is only possible if the donor is a massive convective-envelope giant, namely a red supergiant (RSG). The same is true for neutron-star binaries with massive companions. In other words, pre-CE progenitors of BBH mergers are BH binaries with RSG companions. We find that because of its influence on the radial expansion of massive giants, metallicity has an indirect but a very strong effect on the chemical profile, density structure, and the binding energies of RSG envelopes. Our results suggest that merger rates from population-synthesis models could be severely overestimated, especially at low metallicity. Additionally, the lack of observed RSGs with luminosities above log(L/L⊙) ≈ 5.6 − 5.8, corresponding to stars with M ≳ 40 M⊙, puts into question the viability of the CE channel for the formation of the most massive BBH mergers. Either such RSGs elude detection due to very short lifetimes, or they do not exist and the CE channel can only produce BBH systems with total mass ≲50 M⊙. Finally, we discuss an alternative CE scenario in which a partial envelope ejection is followed by a phase of possibly long and stable mass transfer.

The Tarantula Massive Binary Monitoring

Context. Accurate stellar parameters of individual objects in binary systems are essential to constrain the effects of binarity on stellar evolution. These parameters serve as a prerequisite to probing existing and future theoretical evolutionary models. Aims. We aim to derive the atmospheric parameters of the 31 double-lined spectroscopic binaries in the Tarantula Massive Binary Monitoring sample. This sample, composed of detached, semi-detached and contact systems with at least one of the components classified as an O-type star, is an excellent test-bed to study how binarity can impact our knowledge of the evolution of massive stars. Methods. In the present paper, 32 epochs of FLAMES/GIRAFFE spectra are analysed by using spectral disentangling to construct the individual spectra of 62 components. We then apply the CMFGEN atmosphere code to determine their stellar parameters and their helium, carbon, and nitrogen surface abundances. Results. Among the 31 systems that we study in the present paper, we identify between 48 and 77% of them as detached, likely pre-interacting systems, 16% as semi-detached systems, and between 5 and 35% as systems in or close to contact phase. Based on the properties of their components, we show that the effects of tides on chemical mixing are limited. Components on longer-period orbits show higher nitrogen enrichment at their surface than those on shorter-period orbits, in contrast to expectations of rotational or tidal mixing, implying that other mechanisms play a role in this process. For semi-detached systems, components that fill their Roche lobe are mass donors. They exhibit higher nitrogen content at their surface and rotate more slowly than their companions. By accreting new material, their companions spin faster and are likely rejuvenated. Their locations in the N − v sin i diagram tend to show that binary products are good candidates to populate the two groups of stars (slowly rotating, nitrogen-enriched objects and rapidly rotating non-enriched objects) that cannot be reproduced through single-star population synthesis. Finally, we find no peculiar surface abundances for the components in (over-)contact systems, as has been suggested by evolutionary models for tidal mixing. Conclusions. This sample, consisting of 31 massive binary systems, is the largest sample of binaries composed of at least one O-type star to be studied in such a homogeneous way by applying spectral disentangling and atmosphere modelling. The study of these objects gives us strong observational constraints to test theoretical binary evolutionary tracks.

The evolution of massive binary systems

The evolution of massive stars in close binary systems is significantly different from single star evolution due to a series of interactions between the two stellar components. Such massive close binary systems are linked to various astrophysical phenomena, for example Wolf-Rayet stars, supernova type Ib and Ic, X-ray binaries and gamma-ray bursts. Also, the emission of gravitational waves, recently observed by the LIGO-Virgo detectors, is associated with mergers in binary systems containing compact objects, relics of massive stars - black holes and neutron stars. Evolutionary calculations of massive close binary systems were performed by various authors, but many aspects are not yet fully understood. In this paper, the main concepts of massive close binary evolution are reviewed, together with the most important parameters that can influence the final outcome of the binary system evolution, such as rotation, magnetic fields, stellar wind mass loss and mass accretion efficiency during interactions. An extensive literature overview of massive close binary models in the light of exciting observations connected with those systems is presented.

The radial distribution of supernovae compared to star formation tracers

ABSTRACT Given the limited availability of direct evidence (pre-explosion observations) for supernova (SN) progenitors, the location of supernovae (SNe) within their host galaxies can be used to set limits on one of their most fundamental characteristics, their initial progenitor mass. We present our constraints on SN progenitors derived by comparing the radial distributions of 80 SNe in the Survey for Ionised Neutral Gas Galaxies (SINGG) and Survey of Ultraviolet emission in Neutral Gas Galaxies (SUNGG) to the R-band, H α, and UV light distributions of the 55 host galaxies. The strong correlation of Type Ia SNe with R-band light is consistent with models containing only low-mass progenitors, reflecting earlier findings. When we limit the analysis of Type II SNe to apertures containing 90 per cent of the total flux, the radial distribution of these SNe best traces far-ultraviolet emission, consistent with recent direct detections indicating Type II SNe have moderately massive red supergiant progenitors. Stripped-envelope (SE) SNe have the strongest correlation with H α fluxes, indicative of very massive progenitors (M* ≳ 20 M⊙). This result contradicts a small, but growing, number of direct detections of SE SN progenitors, indicating they are moderately massive binary systems. Our result is consistent, however, with a recent population analysis, suggesting binary SE SN progenitor masses are regularly underestimated. SE SNe are centralized with respect to Type II SNe and there are no SE SNe recorded beyond half the maximum disc radius in the optical and one third the disc radius in the ultraviolet. The absence of SE SNe beyond these distances is consistent with reduced massive star formation efficiencies in the outskirts of the host galaxies.

Massive runaway and walkaway stars

We perform an extensive numerical study of the evolution of massive binary systems to predict the peculiar velocities that stars obtain when their companion collapses and disrupts the system. Our aim is to (i) identify which predictions are robust against model uncertainties and assess their implications, (ii) investigate which physical processes leave a clear imprint and may therefore be constrained observationally, and (iii) provide a suite of publicly available model predictions to allow for the use of kinematic constraints from the Gaia mission. We find that 22+26−8% of all massive binary systems merge prior to the first core-collapse in the system. Of the remainder, 86+11−9% become unbound because of the core-collapse. Remarkably, this rarely produces runaway stars (observationally defined as stars with velocities above 30 km s−1). These are outnumbered by more than an order of magnitude by slower unbound companions, or “walkaway stars”. This is a robust outcome of our simulations and is due to the reversal of the mass ratio prior to the explosion and widening of the orbit, as we show analytically and numerically. For stars more massive than 15 M⊙, we estimate that 10+5−8% are walkaways and only 0.5+1.0−0.4% are runaways, nearly all of which have accreted mass from their companion. Our findings are consistent with earlier studies; however, the low runaway fraction we find is in tension with observed fractions of about 10%. Thus, astrometric data on presently single massive stars can potentially constrain the physics of massive binary evolution. Finally, we show that the high end of the mass distributions of runaway stars is very sensitive to the assumed black hole natal kicks, and we propose this as a potentially stringent test for the explosion mechanism. We also discuss companions remaining bound that can evolve into X-ray and gravitational wave sources.

Hints about the multiplicity of WR 133 based on multiepoch radio observations

Several tens of massive binary systems display indirect, or even strong evidence for non-thermal radio emission, hence their particle accelerator status. These objects are referred to as particle-accelerating colliding-wind binaries (PACWBs). WR 133 is one of the shortest period Wolf-Rayet + O systems in this category, and is therefore critical to characterize the boundaries of the parameter space adequate for particle acceleration in massive binaries. Our methodology consists in analyzing JVLA observations of WR 133 at different epochs to search for compelling evidence for a phase-locked variation attributable to synchrotron emission produced in the colliding-wind region. New data obtained during two orbits reveal a steady and thermal emission spectrum, in apparent contradiction with the previous detection of non-thermal emission. The thermal nature of the radio spectrum along the 112.4-d orbit is supported by the strong free–free absorption by the dense stellar winds, and shows that the simple binary scenario cannot explain the non-thermal emission reported previously. Alternatively, a triple system scenario with a wide, outer orbit would fit with the observational facts reported previously and in this paper, albeit no hint for the existence of a third component exists to date. The epoch-dependent nature of the identification of synchrotron radio emission in WR 133 emphasizes the issue of observational biases in the identification of PACWBs, that undoubtedly affect the present census of PACWB among colliding-wind binaries.

Resolving the MYSO binaries PDS 27 and PDS 37 with VLTI/PIONIER

Context. Binarity and multiplicity appear to be a common outcome in star formation. In particular, the binary fraction of massive (OB-type) stars can be very high. In many cases, the further stellar evolution of these stars is affected by binary interactions at some stage during their lifetime. The origin of this high binarity and the binary parameters are poorly understood because observational constraints are scarce, which is predominantly due to a dearth of known young massive binary systems. Aims. We aim to identify and describe massive young binary systems in order to fill in the gaps of our knowledge of primordial binarity of massive stars, which is crucial for our understanding of massive star formation. Methods. We observed the two massive young stellar objects (MYSOs) PDS 27 and PDS 37 at the highest spatial resolution provided by VLTI/PIONIER in the H-band (1.3 mas). We applied geometrical models to fit the observed squared visibilities and closure phases. In addition, we performed a radial velocity analysis using published VLT/FORS2 spectropolarimetric and VLT/X-shooter spectroscopic observations. Results. Our findings suggest binary companions for both objects at 12 mas (30 au) for PDS 27 and at 22–28 mas (42–54 au) for PDS 37. This means that they are among the closest MYSO binaries resolved to date. Conclusions. Our data spatially resolve PDS 27 and PDS 37 for the first time, revealing two of the closest and most massive (>8 M⊙) YSO binary candidates to date. PDS 27 and PDS 37 are rare but great laboratories to quantitatively inform and test the theories on formation of such systems.

Non-synchronous rotations in massive binary systems

Context. Most massive stars are in binary or multiple systems. Several massive stars have been detected as double-lined spectroscopic binaries and among these, the OWN Survey has detected a non-negligible number whose components show very different spectral line broadening (i.e., projected rotational velocities). This fact raises a discussion about the contributing processes, such as angular-momentum transfer and tidal forces. Aims. We seek to constrain the physical and evolutionary status of one of such systems, the O+O binary HD 93343. Methods. We analyzed a series of high-resolution multiepoch optical spectra to determine the orbital parameters, projected rotational velocities, and evolutionary status of the system. Results. HD 93343 is a binary system comprised of two O7.5 Vz stars that each have minimum masses of approximately 22 M⊙ in a wide and eccentric orbit (e = 0.398±0.004; P = 50.432±0.001 d). Both stars have very similar stellar parameters, and hence ages. As expected from the qualitative appearance of the combined spectrum of the system, however, these stars have very different projected rotational velocities (~65 and ~325 km s−1, respectively). Conclusions. The orbits and stellar parameters obtained for both components seem to indicate that their youth and relative separation is enough to discard the effects of mass transfer and tidal friction. Thus, non-synchronization should be intrinsic to their formation.

Properties of six short-period massive binaries: A study of the effects of binarity on surface chemical abundances

Context.A significant percentage of massive stars are found in multiple systems. The effect of binarity on stellar evolution is poorly constrained. In particular, the role of tides and mass transfer on surface chemical abundances is not constrained observationally.Aims.The aim of this study is to investigate the effect of binarity on the stellar properties and surface abundances of massive binaries.Methods.We performed a spectroscopic analysis of six Galactic massive binaries. We obtained the spectra of individual components via a spectral disentangling method and subsequently analyzed these spectra by means of atmosphere models. The stellar parameters and CNO surface abundances were determined.Results.Most of these six systems are comprised of main-sequence stars. Three systems are detached, two are in contact, and no information is available for the sixth system. For 11 out of the 12 stars studied, the surface abundances are only mildly affected by stellar evolution and mixing. The surface abundances are not different from those of single stars within the uncertainties. The secondary of XZ Cep is strongly chemically enriched. Considering previous determinations of surface abundances in massive binary systems suggests that the effect of tides on chemical mixing is limited, whereas the mass transfer and removal of outer layers of the mass donor leads to the appearance of chemically processed material at the surface, although this is not systematic. The evolutionary masses of the components of our six systems are on average 16.5% higher than the dynamical masses. Some systems seem to have reached synchronization, while others may still be in a transitory phase.