Aims. The mass discrepancy between the observed population of double neutron star binaries by radio pulsar observations and gravitational-wave observation requires an explanation. Methods. Binary population synthesis calculations are performed, and their results are compared with the radio and the gravitational-wave observations simultaneously. Results. Simulations of binary evolution are used to link different observations of double neutron star binaries with each other. The progenitor of GW190425 is investigated in more detail. A distribution of masses and merger times of the possible progenitors is presented. Conclusions. A mass discrepancy between the radio pulsars in the Milky Way with another neutron star companion and the inferred masses from gravitational-wave observations of those kind of merging systems is naturally found in binary evolution.

Tertiary tides (TTs), or the continuous tidal distortion of the tertiary in a hierarchical triple system, can extract energy from the inner binary, inducing within it a proclivity to merge. Despite previous work on the subject, which established that it is significant for certain close triple systems, it is still not a well-understood process. A portion of our ignorance in this regard stems from our inability to integrate a simulation of this phenomenon into conventional stellar evolution codes, since full calculations of these tidal interactions are computationally expensive on stellar evolution timescales. Thus, to attain a better understanding of how these TTs act on longer timescales, an empirical expression of its effects as a function of parameters of the triple system involved is required. In our work, we evaluate the rate at which TTs extract energy from the inner binary within a series of constructed hierarchical triple systems under varying parameters, and study the rate at which the inner binary orbital separation shrinks as a function of those parameters. We find that this rate varies little with the absolute values of the masses of the three component objects, but is very sensitive to the mass ratio of the inner binary $q$, the tertiary radius $R_{\rm 3}$, the inner binary orbital separation $a_{\rm 1}$, the outer orbital separation $a_{\rm 2}$, and the viscoelastic relaxation time of the tertiary $\tau$. More specifically, we find that the percentage by which $a_{\rm 1}$ shrinks per unit time can be reasonably approximated by (1/$a_{\rm 1}$)(d$a_{\rm 1}$/d$t$)=$\left(2.22{\times}10^{-8}{\rm yrs}^{-1}\right)4q\left(1+q\right)^{-2}(R_{\rm 3}/100{\rm R}_{\odot})^{5.2}(a_{\rm 1}/0.2{\rm AU})^{4.8}(a_{\rm 2}/2{\rm AU})^{-10.2}$ $({\tau}/0.534{\rm yrs})^{-1.0}$. We also provide tests of how precise this fitting function is.

The progenitors of many core-collapse supernovae (CCSNe) are expected to be in binary systems. By performing a series of three-dimensional hydrodynamical simulations, we investigate how CCSN explosions affect their binary companion. 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. Also, it is foud that the impact effects of CCSN ejecta on the structure of main-sequence (MS) companions, and thus their long term post-explosion evolution, is in general not be dramatic.

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 about 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).

Radio pulsars found in binary systems with short orbital periods are usually fast spinning as a consequence of recycling via mass transfer from their companion stars; this process is also thought to decrease the magnetic field of the neutron star being recycled. Here, we report on timing observations of the recently discovered binary PSR J1755$-$2550 and find that this pulsar is an exception: with a characteristic age of 2.1 Myr, it is relatively young; furthermore, with a spin period of 315 ms and a surface magnetic field strength at its poles of 0.88$\times$10$^{12}$ G the pulsar shows no sign of having been recycled. Based on its timing and orbital characteristics, the pulsar either has a massive white dwarf (WD) or a neutron star (NS) companion. To distinguish between these two cases, we searched radio observations for a potential recycled pulsar companion and analysed archival optical data for a potential WD companion. Neither work returned conclusive detections. We apply population synthesis modelling and find that both solutions are roughly equally probable. Our population synthesis also predicts a minimum mass of 0.90 M$_{\odot}$ for the companion star to PSR J1755$-$2550 and we simulate the systemic runaway velocities for the resulting WDNS systems which may merge and possibly produce Ca-rich supernovae. Whether PSR J1755$-$2550 hosts a WD or a NS companion star, it is certainly a member of a rare subpopulation of binary radio pulsars.

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 $\sim 0.02\;M_{\odot}$. 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 $\gamma$-ray burst.

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. 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 that process. We calculate envelope binding energies of donor stars with initial masses between 4 and 115 Msun for metallicities of Z=Zsun/2 and Z=Zsun/50, and derive minimum masses of in-spiralling objects needed to eject these envelopes. We find that CE evolution, besides from producing WD-WD and NS-NS binaries, may, in principle, also produce massive BH-BH systems with individual BH component masses up to ~50-60 Msun, in particular for donor stars evolved to giants. 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. Massive stars in a wide range of metallicities and with initial masses up to at least 115 Msun may possibly shed their envelopes and survive CE evolution, depending on their initial orbital parameters, similarly to the situation for intermediate mass and low-mass stars with degenerate cores. We conclude that based on stellar structure calculations, and in the view of the usual simple energy budget analysis, events like GW150914 and GW151226 could possibly be produced from the CE channel. Calculations of post-CE orbital separations, however, and thus the estimated LIGO detection rates, remain highly uncertain. [Abridged]

The progenitors of many CCSNe are expected to be in binary systems. After the SN explosion, the companion 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, the companion is ejected as a runaway and hypervelocity star. By performing a series of 3D hydrodynamical simulations of the collision of SN ejecta with the companion star, we investigate how CCSN explosions affect their companions. We use the BEC code to construct the detailed companion structure at the time of SN explosion. The impact of the SN blast wave on the companion is followed by means of 3D SPH simulations using the Stellar GADGET code. For main-sequence (MS) companions, we find that the amount of removed mass, impact velocity, and 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 SNe~Ia. 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 companions to 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. Our simulations show that the impact effects of Ib/c SN ejecta on the structure of MS companions, 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.