skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Host galaxies and electromagnetic counterparts to binary neutron star mergers across the cosmic time: detectability of GW170817-like events
ABSTRACT The association of GRB170817A with a binary neutron star (BNS) merger has revealed that BNSs produce at least a fraction of short gamma-ray bursts (SGRBs). As gravitational wave (GW) detectors push their horizons, it is important to assess coupled electromagnetic (EM)/GW probabilities and maximize observational prospects. Here, we perform BNS population synthesis calculations with the code mobse, seeding the binaries in galaxies at three representative redshifts, $$z$$ = 0.01, 0.1, and 1 of the Illustris TNG50 simulation. The binaries are evolved and their locations numerically tracked in the host galactic potentials until merger. Adopting the microphysics parameters of GRB170817A, we numerically compute the broad-band light curves of jets from BNS mergers, with the afterglow brightness dependent on the local medium density at the merger site. We perform Monte Carlo simulations of the resulting EM population assuming either a random viewing angle with respect to the jet, or a jet aligned with the orbital angular momentum of the binary, which biases the viewing angle probability for GW-triggered events. We find a gamma-ray detection probability of $$\sim\!2{{\rm per\ cent}},10{{\rm per\ cent}},\mathrm{and}\ 40{{\rm per\ cent}}$$ for BNSs at $$z$$ = 1, 0.1, and 0.01, respectively, for the random case, rising to $$\sim\!75{{\rm per\ cent}}$$ for the $$z$$ = 0.01, GW-triggered aligned case. Afterglow detection probabilities of GW-triggered BNS mergers vary in the range of $$\sim \! 0.3 \!-\! 0.5{{\rm per\ cent}}$$, with higher values for aligned jets, and are comparable across the high- and low-energy bands, unlike gamma-ray-triggered events (cosmological SGRBs) which are significantly brighter at higher energies. We further quantify observational biases with respect to host galaxy masses.  more » « less
Award ID(s):
2006839 1907955
PAR ID:
10364458
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
512
Issue:
2
ISSN:
0035-8711
Page Range / eLocation ID:
p. 2654-2668
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT We examine massive black hole (MBH) mergers and their associated gravitational wave signals from the large-volume cosmological simulation Astrid . Astrid includes galaxy formation and black hole models recently updated with an MBH seed population between 3 × 104h−1M⊙ and 3 × 105h−1M⊙ and a sub-grid dynamical friction (DF) model to follow the MBH dynamics down to 1.5 ckpc h−1. We calculate the initial eccentricities of MBH orbits directly from the simulation at kpc-scales, and find orbital eccentricities above 0.7 for most MBH pairs before the numerical merger. After approximating unresolved evolution on scales below $${\sim 200\, \text{pc}}$$, we find that the in-simulation DF on large scales accounts for more than half of the total orbital decay time ($$\sim 500\, \text{Myr}$$) due to DF. The binary hardening time is an order of magnitude longer than the DF time, especially for the seed-mass binaries (MBH < 2Mseed). As a result, only $$\lesssim 20{{\rm per \,cent}}$$ of seed MBH pairs merge at z > 3 after considering both unresolved DF evolution and binary hardening. These z > 3 seed-mass mergers are hosted in a biased population of galaxies with the highest stellar masses of $$\gt 10^9\, {\rm M}_\odot$$. With the higher initial eccentricity prediction from Astrid , we estimate an expected merger rate of 0.3−0.7 per year from the z > 3 MBH population. This is a factor of ∼7 higher than the prediction using the circular orbit assumption. The Laser Interferometer Space Antenna events are expected at a similar rate, and comprise $$\gtrsim 60\,{\rm{per\,cent}}$$ seed-seed mergers, $$\sim 30\,{\rm{per\,cent}}$$ involving only one seed-mass MBH, and $$\sim 10\,{\rm{per\,cent}}$$ mergers of non-seed MBHs. 
    more » « less
  2. One of the most surprising gamma-ray burst (GRB) features discovered with theSwiftX-ray telescope (XRT) is a plateau phase in the early X-ray afterglow light curves. These plateaus are observed in the majority of long GRBs, while their incidence in short GRBs (SGRBs) is still uncertain due to their fainter X-ray afterglow luminosity with respect to long GRBs. An accurate estimate of the fraction of SGRBs with plateaus is of utmost relevance given the implications that the plateau may have for our understanding of the jet structure and possibly of the nature of the binary neutron star (BNS) merger remnant. This work presents the results of an extensive data analysis of the largest and most up-to-date sample of SGRBs observed with the XRT, and for which the redshift has been measured. We find a plateau incidence of 18–37% in SGRBs, which is a significantly lower fraction than that measured in long GRBs (> 50%). Although still debated, the plateau phase could be explained as energy injection from the spin-down power of a newly born magnetized neutron star (NS; magnetar). We show that this scenario can nicely reproduce the observed short GRB (SGRBs) plateaus, while at the same time providing a natural explanation for the different plateau fractions between short and long GRBs. In particular, our findings may imply that only a minority of BNS mergers generating SGRBs leave behind a sufficiently stable or long-lived NS to form a plateau. From the probability distribution of the BNS remnant mass, a fraction 18–37% of short GRB plateaus implies a maximum NS mass in the range ∼2.3 − 2.35 M
    more » « less
  3. ABSTRACT Some short gamma-ray bursts (SGRBs) show a longer lasting emission phase, called extended emission (EE) lasting $${\sim}10^{2\!-\!3}\, \rm s$$, as well as a plateau emission (PE) lasting $${\sim}10^{4\!-\!5}\, \rm s$$. Although a long-lasting activity of the central engines is a promising explanation for powering both emissions, their physical origin and their emission mechanisms are still uncertain. In this work, we study the properties of the EEs and their connection with the PEs. First, we constrain the minimal Lorentz factor Γ of the outflows powering EEs, using compactness arguments and find that the outflows should be relativistic, Γ ≳ 10. We propose a consistent scenario for the PEs, where the outflow eventually catches up with the jet responsible for the prompt emission, injecting energy into the forward shock formed by the prior jet, which naturally results in a PE. We also derive the radiation efficiency of EEs and the Lorentz factor of the outflow within our scenario for 10 well-observed SGRBs accompanied by both EE and PE. The efficiency has an average value of $${\sim}3\, {{\ \rm per\ cent}}$$ but shows a broad distribution ranging from ∼0.01 to $${\sim}100{{\ \rm per\ cent}}$$. The Lorentz factor is ∼20–30, consistent with the compactness arguments. These results suggest that EEs are produced by a slower outflow via more inefficient emission than the faster outflow that causes the prompt emission with a high radiation efficiency. 
    more » « less
  4. Abstract The association of GRB170817A with GW170817 has confirmed the long-standing hypothesis that binary neutron star (BNS) mergers are the progenitors of at least some short gamma-ray bursts (SGRBs). This connection has ushered in an era in which broadband observations of SGRBs, together with measurements of the time delay between the gravitational waves and the electromagnetic radiation, allow for probing the properties of the emitting outflow and its engine to an unprecedented detail. Because the structure of the radiating outflow is molded by the interaction of a relativistic jet with the binary ejecta, it is of paramount importance to study the system in a realistic setting. Here we present a three-dimensional hydrodynamic simulation of a relativistic jet propagating in the ejecta of a BNS merger, which were computed with a general relativistic magnetohydrodynamic simulation. We find that the jet’s centroid oscillates around the axis of the system, due to inhomogeneities encountered in the propagation. These oscillations allow the jet to find the path of least resistance and travel faster than an identical jet in smooth ejecta. In our setup the breakout time is ∼0.6 s, which is comparable to the expected central engine duration in SGRBs and possibly a non-negligible fraction of the total delay between the gravitational and gamma-ray signals. Our simulation also shows that energy is carried in roughly equal amounts by the jet and by the cocoon, and that about 20% of the injected energy is transferred to the ejecta via mechanical work. 
    more » « less
  5. ABSTRACT We perform a comparative analysis of nucleosynthesis yields from binary neutron star (BNS) mergers, black hole-neutron star (BHNS) mergers, and core-collapse supernovae (CCSNe) with the goal of determining which are the most dominant sources of r-process enrichment observed in stars. We find that BNS and BHNS binaries may eject similar mass distributions of robust r-process nuclei post-merger (up to third peak and actinides, A ∼ 200−240), after accounting for the volumetric event rates. Magnetorotational (MR) CCSNe likely undergo a weak r-process (up to A ∼ 140) and contribute to the production of light element primary process (LEPP) nuclei, whereas typical thermal, neutrino-driven CCSNe only synthesize up to first r-process peak nuclei (A ∼ 80−90). We also find that the upper limit to the rate of MR CCSNe is $$\lesssim 1~{{\ \rm per\ cent}}$$ the rate of typical thermal CCSNe; if the rate was higher, then weak r-process nuclei would be overproduced. Although the largest uncertainty is from the volumetric event rate, the prospects are encouraging for confirming these rates in the next few years with upcoming surveys. Using a simple model to estimate the resulting kilonova light curve from mergers and our set of fiducial merger parameters, we predict that ∼7 BNS and ∼2 BHNS events will be detectable per year by the Vera C. Rubin Observatory (LSST), with prior gravitational wave (GW) triggers. 
    more » « less