No Arabic abstract
We summarize our results on the Galactic merger rate of double neutron stars (DNS) in view of the recent discovery of PSR J0737-3039. We also present previously unpublished results for the global probability distribution of merger rate values that incorporate the presently known systematics from the radio pulsar luminosity function. The most likely value obtained from the global distribution is only ~15 per Myr, but a re-analysis of the current pulsar sample and radio luminosities is needed for a reliable assessment of the best fitting distribution. Finally, we use our theoretical understanding of DNS formation to calculate a possible upper limit on the DNS merger rate from current Type Ib/c supernova rate estimates.
We revisit the merger rate for Galactic double neutron star (DNS) systems in light of recent observational insight into the longitudinal and latitudinal beam shape of the relativistic DNS PSR J1906$+$0746. Due to its young age and its relativistic orbit, the pulsar contributes significantly to the estimate of the joint Galactic merger rate. We follow previous analyses by modelling the underlying pulsar population of nine merging DNS systems and study the impact and resulting uncertainties when replacing simplifying assumptions made in the past with actual knowledge of the beam shape, its extent and the viewing geometry. We find that the individual contribution of PSR J1906$+$0746 increases to $R = 6^{+28}_{-5}$ Myr$^{-1}$ although the values is still consistent with previous estimates given the uncertainties. We also compute contributions to the merger rates from the other DNS systems by applying a generic beam shape derived from that of PSR J1906+0746, evaluating the impact of previous assumptions. We derive a joint Galactic DNS merger rate of $R^{rm{gen}}_{rm{MW}} = 32^{+19}_{-9}$Myr$^{-1}$, leading to a LIGO detection rate of ${R}^{rm{gen}}_{rm{LIGO}} = 3.5^{+2.1}_{-1.0}$Myr$^{-1}$ (90% conf. limit), considering the upcoming O3 sensitivity of LIGO. As these values are in good agreement with previous estimates, we conclude that the method of estimating the DNS merger and LIGO detection rates via the study of the radio pulsar DNS population is less prone to systematic uncertainties than previously thought.
We estimate the coalescence rate of close binaries with two neutron stars (NS) and discuss the prospects for the detection of NS-NS inspiral events by ground-based gravitational-wave observatories, such as LIGO. We derive the Galactic coalescence rate using the observed sample of close NS-NS binaries (PSR B1913+16 and PSR B1534+12) and examine in detail each of the sources of uncertainty associated with the estimate. Specifically, we investigate (i) the dynamical evolution of NS-NS binaries in the Galactic potential and the vertical scale height of the population, (ii) the pulsar lifetimes, (iii) the effects of the faint end of the radio pulsar luminosity function and their dependence on the small number of observed objects, (iv) the beaming fraction, and (v) the extrapolation of the Galactic rate to extragalactic distances expected to be reachable by LIGO. We find that the dominant source of uncertainty is the correction factor (up to about 200) for faint (undetectable) pulsars. All other sources are much less important, each with uncertainty factors smaller than 2. Despite the relatively large uncertainty, the derived coalescence rate is approximately consistent with previously derived upper limits, and is more accurate than rates obtained from population studies. We obtain a most conservative lower limit for the LIGO II detection rate of 2 events per year. Our upper limit on the detection rate lies between 300 to more than 1000 events per year.
LIGO/Virgo Collaboration reported the detection of the most massive black hole - black hole (BH-BH) merger up to date with component masses of 85 Msun and 66 Msun (GW190521). Motivated by recent observations of massive stars in the 30 Doradus cluster in the Large Magellanic Cloud (>200 Msun; e.g. R136a) and employing newly estimated uncertainties on pulsational pair-instability mass-loss (that allow for possibility of forming BHs with mass up to 90Msun) we show that it is trivial to form such massive BH-BH mergers through the classical isolated binary evolution (with no assistance from either dynamical interactions or exotica). A binary consisting of two massive (180+150 Msun) Population II stars (Z=0.0001) evolves through a stable Roche lobe overflow and common envelope episode. Both exposed stellar cores undergo direct core-collapse and form massive BHs while avoiding pair-instability pulsation mass-loss or total disruption. LIGO/Virgo observations show that the merger rate density of light BH-BH mergers (both components: <50 Msun) is of the order of 10-100 Gpc^-3 yr^-1, while GW190521 indicates that the rate of heavier mergers is 0.02-0.43 Gpc^-3 yr^-1. Our model (with standard assumptions about input physics) but extended to include 200 Msun stars and allowing for the possibility of stellar cores collapsing to 90 Msun BHs produces the following rates: 63 Gpc^-3 yr^-1 for light BH-BH mergers and 0.04 Gpc^-3 yr^-1 for heavy BH-BH mergers. We do not claim that GW190521 was formed by an isolated binary, but it appears that such a possibility can not be excluded.
The astrophysical r-process site where about half of the elements heavier than iron are produced has been a puzzle for several decades. Here we discuss the role of neutron star mergers (NSMs) in the light of the first direct detection of such an event in both gravitational (GW) and electromagnetic (EM) waves. We analyse bolometric and NIR lightcurves of the first detected double neutron star merger and compare them to nuclear reaction network-based macronova models. The slope of the bolometric lightcurve is consistent with the radioactive decay of neutron star ejecta with $Y_e lesssim 0.3$ (but not larger), which provides strong evidence for an r-process origin of the electromagnetic emission. This rules out in particular nickel winds as major source of the emission. We find that the NIR lightcurves can be well fitted either with or without lanthanide-rich ejecta. Our limits on the ejecta mass together with estimated rates directly confirm earlier purely theoretical or indirect observational conclusions that double neutron star mergers are indeed a major site of cosmic nucleosynthesis. If the ejecta mass was {em typical}, NSMs can easily produce {em all} of the estimated Galactic r-process matter, and --depending on the real rate-- potentially even more. This could be a hint that the event ejected a particularly large amount of mass, maybe due to a substantial difference between the component masses. This would be compatible with the mass limits obtained from the GW-observation. The recent observations suggests that NSMs are responsible for a broad range of r-process nuclei and that they are at least a major, but likely the dominant r-process site in the Universe.
The first neutron star-neutron star (NS-NS) merger was discovered on August 17, 2017 through gravitational waves (GW170817) and followed with electromagnetic observations. This merger was detected in an old elliptical galaxy with no recent star formation. We perform a suite of numerical calculations to understand the formation mechanism of this merger. We probe three leading formation mechanisms of double compact objects: classical isolated binary star evolution, dynamical evolution in globular clusters and nuclear cluster formation to test whether they are likely to produce NS-NS mergers in old host galaxies. Our simulations with optimistic assumptions show current NS-NS merger rates at the level of 10^-2 yr^-1 from binary stars, 5 x 10^-5 yr^-1 from globular clusters and 10^-5 yr^-1 from nuclear clusters for all local elliptical galaxies (within 100 Mpc^3). These models are thus in tension with the detection of GW170817 with an observed rate 1.5 yr^-1 (per 100 Mpc^3; LIGO/Virgo estimate). Our results imply that either (i) the detection of GW170817 by LIGO/Virgo at their current sensitivity in an elliptical galaxy is a statistical coincidence; or that (ii) physics in at least one of our three models is incomplete in the context of the evolution of stars that can form NS-NS mergers; or that (iii) another very efficient (unknown) formation channel with a long delay time between star formation and merger is at play.