Double neutron stars (DNSs) have been observed as Galactic radio pulsars, and the recent discovery of gravitational waves from the DNS merger GW170817 adds to the known DNS population. We perform rapid population synthesis of massive binary stars and discuss model predictions, including DNS formation rates, mass distributions, and delay time distributions. We vary assumptions and parameters of physical processes such as mass transfer stability criteria, supernova natal kick distributions, remnant mass prescriptions and common-envelope energetics. We compute the likelihood of observing the orbital period-eccentricity distribution of the Galactic DNS population under each of our population synthesis models, allowing us to quantitatively compare the models. We find that mass transfer from a stripped post-helium-burning secondary (case BB) onto a neutron star is most likely dynamically stable. We also find that a natal kick distribution composed of both low (Maxwellian σ = 30 km s −1 ) and high (σ = 265 km s −1 ) components is preferred over a single high-kick component. We conclude that the observed DNS mass distribution can place strong constraints on model assumptions.

Detections of gravitational waves are now starting to probe the mass distribution of stellar mass black holes (BHs). Robust predictions from stellar models are needed to interpret these. Theory predicts the existence of a gap in the BH mass distribution because of pair-instability supernovae. The maximum BH mass below the gap is the result of pulsational mass loss. We evolve massive helium stars through their late hydrodynamical phases of evolution using the open-source MESA stellar evolution code. We find that the location of the lower edge of the mass gap at 45 M  is remarkably robust against variations in the metallicity (≈3 M ), the treatment of internal mixing (≈1 M ), and stellar wind mass loss (≈4 M ), making it the most robust predictor for the final stages of the evolution of massive stars. The reason is that the onset of the instability is dictated by the near-final core mass, which in turn sets the resulting BH mass. However, varying the a g C , O 12 16 () reaction rate within its 1σ uncertainties shifts the location of the gap between 40 M  and 56 M . We provide updated analytic fits for population synthesis simulations. Our results imply that the detection of merging BHs can provide constraints on nuclear astrophysics. Furthermore, the robustness against metallicity suggests that there is a universal maximum for the location of the lower edge of the gap, which is insensitive to the formation environment and redshift for first-generation BHs. This is promising for the possibility to use the location of the gap as a "standard siren" across the universe.

Pair-instability and pulsational pair-instability supernovae (PPISN) have not been unambiguously observed so far. They are, however, promising candidates for the progenitors of the heaviest binary black hole (BBH) mergers detected. If these BBHs are the product of binary evolution, then PPISNe could occur in very close binaries. Motivated by this, we discuss the implications of a PPISN happening with a close binary companion, and what impact these explosions have on the formation of merging BBHs through binary evolution. For this, we have computed a set of models of metal-poor (Z /10) helium stars using the MESA software instrument. For PPISN progenitors with pre-explosion masses > 50M we find that, after a pulse, heat deposited throughout the layers of the star that remain bound cause it to expand to more than 100R for periods of 10 2 − 10 4 yrs depending on the mass of the progenitor. This results in long-lived phases of Roche-lobe overflow or even common-envelope events if there is a close binary companion, leading to additional electromagnetic transients associated to PPISN eruptions. If we ignore the effect of these interactions, we find that mass loss from PPISNe reduces the final black hole spin by ∼ 30%, induces eccentricities that can be detected by the LISA observatory, and can produce a double-peaked distribution of measured chirp masses in BBH mergers observed by ground-based detectors.

We present the SuperNova Explosion Code (SNEC), an open-source Lagrangian code for the hydrodynamics and equilibrium-diffusion radiation transport in the expanding envelopes of supernovae. Given a model of a progenitor star, an explosion energy, and an amount and distribution of radioactive nickel, SNEC generates the bolometric light curve, as well as the light curves in different broad bands assuming black body emission. As a first application of SNEC, we consider the explosions of a grid of 15 M (at zero-age main sequence) stars whose hydrogen envelopes are stripped to different extents and at different points in their evolution. The resulting light curves exhibit plateaus with durations of ∼20 − 100 days if 1.5 − 2 M of hydrogen-rich material is left and no plateau if less hydrogen-rich material is left. If these shorter plateau lengths are not seen for Type IIP supernovae in nature, it suggests that, at least for zero-age main sequence masses 20 M , hydrogen mass loss occurs as an all or nothing process. This perhaps points to the important role binary interactions play in generating the observed mass-stripped supernovae (i.e., Type Ib/c events). These light curves are also unlike what is typically seen for Type IIL supernovae, arguing that simply varying the amount of mass loss cannot explain these events. The most stripped models begin to show double-peaked light curves similar to what is often seen for Type IIb supernovae, confirming previous work that these supernovae can come from progenitors that have a small amount of hydrogen and a radius of ∼ 500 R .

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.

Stars stripped of their hydrogen-rich envelope through interaction with a binary companion are generally not considered when accounting for ionizing radiation from stellar populations, despite the expectation that stripped stars emit hard ionizing radiation, form frequently and live 10 − 100 times longer than single massive stars. We compute the first grid of evolutionary and spectral models specially made for stars stripped in binaries for a range of progenitor masses (2-20 M ) and metallicities ranging from solar to values representative for pop II stars. For stripped stars with masses in the range 0.3-7 M , we find consistently high effective temperatures (20 000-100 000 K, increasing with mass), small radii (0.2-1 R ) and high bolometric luminosities, comparable to that of their progenitor before stripping. The spectra show a continuous sequence that naturally bridge subdwarf-type stars at the low mass end and Wolf-Rayet like spectra at the high mass end. For intermediate masses we find hybrid spectral classes showing a mixture of absorption and emission lines. These appear for stars with mass loss rates of 10 −8 − 10 −6 M yr −1 , which have semi-transparent atmospheres. At low metallicity, substantial hydrogen-rich layers are left at the surface and we predict spectra that resemble O-type stars instead. We obtain spectra undistinguishable from subdwarfs for stripped stars with masses up to 1.7 M , which questions whether the widely adopted canonical value of 0.47 M is uniformly valid. Only a handful of stripped stars of intermediate mass have currently been identified observationally. Increasing this sample will provide necessary tests for the physics of interaction, internal mixing and stellar winds. We use our model spectra to investigate the feasibility to detect stripped stars next to an optically bright companion and recommend systematic searches for their UV excess and possible emission lines, most notably HeII λ4686 in the optical and HeII λ1640 in the UV. Our models are publicly available for further investigations or inclusion in spectral synthesis simulations.

Gravitational-wave detections are starting to allow us to probe the physical processes in the evolution of very massive stars through the imprints they leave on their final remnants. Stellar evolution theory predicts the existence of a gap in the black hole mass distribution at high mass due to the effects of pair instability. Previously, we showed that the location of the gap is robust against model uncertainties, but it does depend sensitively on the uncertain a g C , O 12 16 300. Degeneracies with other model uncertainties need to be investigated further, but probing nuclear stellar astrophysics poses a promising science case for the future gravitational-wave detectors.

Most massive stars, the progenitors of core-collapse supernovae, are in close binary systems and may interact with their companion through mass transfer or merging. We undertake a population synthesis study to compute the delay-time distribution of core-collapse supernovae, that is, the supernova rate versus time following a starburst, taking into account binary interactions. We test the systematic robustness of our results by running various simulations to account for the uncertainties in our standard assumptions. We find that a significant fraction, 15 +9 −8 %, of core-collapse supernovae are 'late', that is, they occur 50-200 Myrs after birth, when all massive single stars have already exploded. These late events originate predominantly from binary systems with at least one, or, in most cases, with both stars initially being of intermediate mass (4 − 8M ). The main evolutionary channels that contribute often involve either the merging of the initially more massive primary star with its companion or the engulfment of the remaining core of the primary by the expanding secondary that has accreted mass at an earlier evolutionary stage. Also, the total number of core-collapse supernovae increases by 14 +15 −14 % because of binarity for the same initial stellar mass. The high rate implies that we should have already observed such late core-collapse supernovae, but have not recognized them as such. We argue that φ Persei is a likely progenitor and that eccentric neutron star -white dwarf systems are likely descendants. Late events can help explain the discrepancy in the delay-time distributions derived from supernova remnants in the Magellanic Clouds and extragalactic type Ia events, lowering the contribution of prompt Ia events. We discuss ways to test these predictions and speculate on the implications for supernova feedback in simulations of galaxy evolution.