On 2017 August 17 a binary neutron star coalescence candidate (later designated GW170817) with merger time 12:41:04 UTC was observed through gravitational waves by the Advanced LIGO and Advanced Virgo detectors. The Fermi Gamma-ray Burst Monitor independently detected a gamma-ray burst (GRB 170817A) with a time delay of ∼ 1.7 s with respect to the merger time. From the gravitational-wave signal, the source was initially localized to a sky region of 31 deg2 at a luminosity distance of 40 − 8 + 8 Mpc and with component masses consistent with neutron stars. The component masses were later measured to be in the range 0.86 to 2.26 M ⊙ . An extensive observing campaign was launched across the electromagnetic spectrum leading to the discovery of a bright optical transient (SSS17a, now with the IAU identification of AT 2017gfo) in NGC 4993 (at ∼ 40 Mpc ) less than 11 hours after the merger by the One-Meter, Two Hemisphere (1M2H) team using the 1 m Swope Telescope. The optical transient was independently detected by multiple teams within an hour. Subsequent observations targeted the object and its environment. Early ultraviolet observations revealed a blue transient that faded within 48 hours. Optical and infrared observations showed a redward evolution over ∼10 days. Following early non-detections, X-ray and radio emission were discovered at the transient’s position ∼ 9 and ∼ 16 days, respectively, after the merger. Both the X-ray and radio emission likely arise from a physical process that is distinct from the one that generates the UV/optical/near-infrared emission. No ultra-high-energy gamma-rays and no neutrino candidates consistent with the source were found in follow-up searches. These observations support the hypothesis that GW170817 was produced by the merger of two neutron stars in NGC 4993 followed by a short gamma-ray burst (GRB 170817A) and a kilonova/macronova powered by the radioactive decay of r-process nuclei synthesized in the ejecta.

We have carried out a detailed binary population synthesis (BPS) study of the formation of subdwarf B (sdB) stars and related objects (sdO, sdOB stars) using the latest version of the BPS code developed by Han and co-workers. We systematically investigate the importance of the five main evolutionary channels in which the sdB stars form after one or two common-envelope (CE) phases, one or two phases of stable Roche lobe overflow (RLOF) or as the result of the merger of two helium white dwarfs (WDs). Our best BPS model can satisfactorily explain the main observational characteristics of sdB stars, in particular their distributions in the orbital period-minimum companion mass (log P-M comp ) diagram and in the effective temperaturesurface gravity (T eff -log g) diagram, their distributions of orbital period, log (gθ 4 ) (θ = 5040 K /T eff ) and mass function, their binary fraction and the fraction of sdB binaries with WD companions, their birth rates and their space density. We obtain a Galactic formation rate for sdB stars of 0.014-0.063 yr −1 with a best estimate of ∼0.05 yr −1 and a total number in the Galaxy of 2.4-9.5 × 10 6 with a best estimate of ∼6 × 10 6 ; half of these may be missing in observational surveys owing to selection effects. The intrinsic binary fraction is 76-89 per cent, although the observed frequency may be substantially lower owing to the selection effects. The first CE ejection channel, the first stable RLOF channel and the merger channel are intrinsically the most important channels, although observational selection effects tend to increase the relative importance of the second CE ejection and merger channels. We also predict a distribution of masses for sdB stars that is wider than is commonly assumed and that some sdB stars have companions of spectral type as early as B. The percentage of A-type stars with sdB companions can in principle be used to constrain some of the important parameters in the binary evolution model. We conclude that (i) the first RLOF phase needs to be more stable than is commonly assumed, either because the critical mass ratio q crit for dynamical mass transfer is higher or because of tidally enhanced stellar wind mass loss; (ii) mass transfer in the first stable RLOF phase is non-conservative, and the mass lost from the system takes away a specific angular momentum similar to that of the system; and (iii) common-envelope ejection is very efficient.

No abstract

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.