We present multi-wavelength observations of SN 2014C during the first 500 days. These observations represent the first solid detection of a young extragalactic stripped-envelope SN out to high-energy X-rays ∼40 keV. SN 2014C shows ordinary explosion parameters (E k ∼1.8×10 51 erg and M ej ∼1.7 M e ). However, over an ∼1 year timescale, SN 2014C evolved from an ordinary hydrogen-poor supernova into a strongly interacting, hydrogen-rich supernova, violating the traditional classification scheme of type-I versus type-II SNe. Signatures of the SN shock interaction with a dense medium are observed across the spectrum, from radio to hard X-rays, and revealed the presence of a massive shell of ∼1 M e of hydrogen-rich material at ∼6×10 16 cm. The shell was ejected by the progenitor star in the decades to centuries before collapse. This result challenges current theories of massive star evolution, as it requires a physical mechanism responsible for the ejection of the deepest hydrogen layer of H-poor SN progenitors synchronized with the onset of stellar collapse. Theoretical investigations point at binary interactions and/or instabilities during the last nuclear burning stages as potential triggers of the highly time-dependent mass loss. We constrain these scenarios utilizing the sample of 183 SNe Ib/c with public radio observations. Our analysis identifies SN 2014C-like signatures in ∼10% of SNe. This fraction is reasonably consistent with the expectation from the theory of recent envelope ejection due to binary evolution if the ejected material can survive in the close environment for 10 3 -10 4 years. Alternatively, nuclear burning instabilities extending to core C-burning might play a critical role.

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.

Context. After years of scientific progress, the origin of stellar binary black holes is still a great mystery. Several formation channels for merging black holes have been proposed in the literature. As more merger detections are expected with future gravitational-wave observations, population synthesis studies can help to distinguish between them. Aims. We study the formation of coalescing binary black holes via the evolution of isolated field binaries that go through the common envelope phase in order to obtain the combined distributions of observables such as black-hole spins, masses and cosmological redshifts of mergers. Methods. To achieve this aim, we used a hybrid technique that combines the parametric binary population synthesis code COMPAS with detailed binary evolution simulations performed with the MESA code. We then convolved our binary evolution calculations with the redshift- and metallicity-dependent star-formation rate and the selection effects of gravitational-wave detectors to obtain predictions of observable properties. Results. By assuming efficient angular momentum transport, we are able to present a model that is capable of simultaneously predicting the following three main gravitational-wave observables: the effective inspiral spin parameter χeff, the chirp mass Mchirp and the cosmological redshift of merger zmerger. We find an excellent agreement between our model and the ten events from the first two advanced detector observing runs. We make predictions for the third observing run O3 and for Advanced LIGO design sensitivity. We expect approximately 80% of events with χeff <  0.1, while the remaining 20% of events with χeff ≥ 0.1 are split into ∼10% with Mchirp <  15 M⊙ and ∼10% with Mchirp ≥ 15 M⊙. Moreover, we find that Mchirp and χeff distributions are very weakly dependent on the detector sensitivity. Conclusions. The favorable comparison of the existing LIGO/Virgo observations with our model predictions gives support to the idea that the majority, if not all of the observed mergers, originate from the evolution of isolated binaries. The first-born black hole has negligible spin because it lost its envelope after it expanded to become a giant star, while the spin of the second-born black hole is determined by the tidal spin up of its naked helium star progenitor by the first-born black hole companion after the binary finished the common-envelope phase.