Double neutron star (DNS) systems represent extreme physical objects and the endpoint of an exotic journey of stellar evolution and binary interactions. Large numbers of DNS systems and their mergers are anticipated to be discovered using the Square-Kilometre-Array searching for radio pulsars and high-frequency gravitational wave detectors (LIGO/VIRGO), respectively. Here we discuss all key properties of DNS systems, as well as selection effects, and combine the latest observational data with new theoretical progress on various physical processes with the aim of advancing our knowledge on their formation. We examine key interactions of their progenitor systems and evaluate their accretion history during the high-mass X-ray binary stage, the common envelope phase and the subsequent Case BB mass transfer, and argue that the first-formed NSs have accreted at most ∼ 0.02 M ⊙ . We investigate DNS masses, spins and velocities, and in particular correlations between spin period, orbital period and eccentricity. Numerous Monte Carlo simulations of the second supernova (SN) events are performed to extrapolate pre-SN stellar properties and probe the explosions. All known close-orbit DNS systems are consistent with ultra-stripped exploding stars. Although their resulting NS kicks are often small, we demonstrate a large spread in kick magnitudes which may, in general, depend on the past interaction history of the exploding star and thus correlate with the NS mass. We analyze and discuss NS kick directions based on our SN simulations. Finally, we discuss the terminal evolution of close-orbit DNS systems until they merge and possibly produce a short γ-ray burst.

The first gravitational wave detections of mergers between black holes and neutron stars represent a remarkable new regime of high-energy transient astrophysics. The signals observed with LIGO-Virgo detectors come from mergers of extreme physical objects which are the end products of stellar evolution in close binary systems. To better understand their origin and merger rates, we have performed binary population syntheses at different metallicities using the new grid-based binary population synthesis code ComBinE. Starting from newborn pairs of stars, we follow their evolution including mass loss, mass transfer and accretion, common envelopes and supernova explosions. We apply the binding energies of common envelopes based on dense grids of detailed stellar structure models, make use of improved investigations of the subsequent Case BB Roche-lobe overflow and scale supernova kicks according to the stripping of the exploding stars. We demonstrate that all the double black hole mergers, GW150914, LVT151012, GW151226, GW170104, GW170608 and GW170814, as well as the double neutron star merger GW170817, are accounted for in our models in the appropriate metallicity regime. Our binary interaction parameters are calibrated to match the accurately determined properties of Galactic double neutron star systems, and we discuss their masses and types of supernova origin. Using our default values for the input physics parameters, we find a double neutron star merger rate of 3.0 Myr −1 for Milky-Way equivalent galaxies. Our upper limit to the merger-rate density of double neutron stars is R 400 yr −1 Gpc −3 in the local Universe (z = 0).

Context. The recently detected gravitational wave signals (GW150914 and GW151226) of the merger event of a pair of relatively massive stellar-mass black holes (BHs) calls for an investigation of the formation of such progenitor systems in general. Aims. We analyse the common-envelope (CE) stage of the traditional formation channel in binaries where the first-formed compact object undergoes an in-spiral inside the envelope of its evolved companion star and ejects the envelope in this process. Methods. We calculated envelope binding energies of donor stars with initial masses between 4 and 115 M ⊙ for metallicities of Z = Z Milky Way ≃ Z ⊙ /2 and Z = Z ⊙ /50, and derived minimum masses of in-spiralling objects needed to eject these envelopes.Results. In addition to producing double white dwarf and double neutron star binaries, CE evolution may also produce massive BH-BH systems with individual BH component masses of up to ∼ 50 − 60 M ⊙ , in particular for donor stars evolved to giants beyond the Hertzsprung gap. However, the physics of envelope ejection of massive stars remains uncertain. We discuss the applicability of the energy-budget formalism, the location of the bifurcation point, the recombination energy, and the accretion energy during in-spiral as possible energy sources, and also comment on the effect of inflated helium cores. Conclusions. Massive stars in a wide range of metallicities and with initial masses of up to at least 115 M ⊙ may shed their envelopes and survive CE evolution, depending on their initial orbital parameters, similarly to the situation for intermediate-and low-mass stars with degenerate cores. In addition to being dependent on stellar radius, the envelope binding energies and λ-values also depend on the applied convective core-overshooting parameter, whereas these structure parameters are basically independent of metallicity for stars with initial masses below 60 M ⊙ . Metal-rich stars > ∼ 60 M ⊙ become luminous blue variables and do not evolve to reach the red giant stage. We conclude that based on stellar structure calculations, and in the view of the usual simple energy budget analysis, events like GW150914 and GW151226 might be produced by the CE channel. Calculations of post-CE orbital separations, however, and thus the estimated LIGO detection rates, remain highly uncertain.

Context. The progenitors of many core-collapse supernovae (CCSNe) are expected to be in binary systems. After the SN explosion in a binary, the companion star may suffer from mass stripping and be shock heated as a result of the impact of the SN ejecta. If the binary system is disrupted by the SN explosion, the companion star is ejected as a runaway star, and in some cases as a hypervelocity star. Aims. By performing a series of three-dimensional (3D) hydrodynamical simulations of the collision of SN ejecta with the companion star, we investigate how CCSN explosions affect their binary companion. Methods. We use the BEC stellar evolution code to construct the detailed companion structure at the moment of SN explosion. The impact of the SN blast wave on the companion star is followed by means of 3D smoothed particle hydrodynamics (SPH) simulations using the Stellar GADGET code. Results. For main-sequence (MS) companion stars, we find that the amount of removed stellar mass, the resulting impact velocity, and the chemical contamination of the companion that results from the impact of the SN ejecta strongly increases with decreasing binary separation and increasing explosion energy. Their relationship can be approximately fitted by power laws, which is consistent with the results obtained from impact simulations of Type Ia SNe. However, we find that the impact velocity is sensitive to the momentum profile of the outer SN ejecta and, in fact, may decrease with increasing ejecta mass, depending on the modeling of the ejecta. Because most companion stars to Type Ib/c CCSNe are in their MS phase at the moment of the explosion, combined with the strongly decaying impact effects with increasing binary separation, we argue that the majority of these SNe lead to inefficient mass stripping and shock heating of the companion star following the impact of the ejecta. Conclusions. Our simulations show that the impact effects of Type Ib/c SN ejecta on the structure of MS companion stars, and thus their long-term post-explosion evolution, is in general not dramatic. We find that at most 10% of their mass is lost and their resulting impact velocities are less than 100 km s −1 .

We present a catalog containing 839 candidate post–common envelope systems. Common envelope evolution is very important in stellar astrophysics, particularly in the context of very compact and short-period binaries, including cataclysmic variables, as progenitors of, e.g., supernovae Type Ia or mergers of black holes and/or neutron stars. At the same time, it is a barely understood process in binary evolution. Due to limitations, since partially remedied, on direct simulation, early investigations were mainly focused on providing analytic prescriptions of the outcome of common envelope evolution. In recent years, detailed hydrodynamical calculations have produced deeper insight into the previously elusive process of envelope ejection. However, a direct link between the observations and theory of this relatively short-lived phase in binary evolution has not been forthcoming. Therefore, the main insight to be gained from observations has to be derived from the current state of systems likely to have gone through a common envelope. Here we present an extensive catalog of such observations as found in the literature. The aim of this paper is to provide a reliable set of data, obtained from observations, to be used in the theoretical modeling of common envelope evolution. In this catalog, the former common envelope donor star is commonly observed as a white dwarf or hot subdwarf star. This catalog includes period and mass estimates wherever obtainable. Some binaries are borderline cases to allow an investigation of the transition between a common envelope formation and other mass-transfer processes.