No Arabic abstract
We present a detailed evaluation of the expected rate of joint gravitational-wave and short gamma-ray burst (GRB) observations over the coming years. We begin by evaluating the improvement in distance sensitivity of the gravitational wave search that arises from using the GRB observation to restrict the time and sky location of the source. We argue that this gives a 25% increase in sensitivity when compared to an all-sky, all-time search, corresponding to more than doubling the number of detectable gravitational wave signals associated with GRBs. Using this, we present the expected rate of joint observations with the advanced LIGO and Virgo instruments, taking into account the expected evolution of the gravitational wave detector network. We show that in the early advanced gravitational wave detector observing runs, from 2015-2017, there is only a small chance of a joint observation. However, as the detectors approach their design sensitivities, there is a good chance of joint observations provided wide field GRB satellites, such as Fermi and the Interplanetary Network, continue operation. The rate will also depend critically upon the nature of the progenitor, with neutron star--black hole systems observable to greater distances than double neutron star systems. The relative rate of binary mergers and GRBs will depend upon the jet opening angle of GRBs. Consequently, joint observations, as well as accurate measurement of both the GRB rate and binary merger rates, will allow for an improved estimation of the opening angle of GRBs.
Coalescing neutron star (NS)-black hole (BH) binaries are promising sources of gravitational-waves (GWs) to be detected within the next few years by current GW observatories. If the NS is tidally disrupted outside the BH innermost stable circular orbit, an accretion torus may form, and this could eventually power a short gamma-ray burst (SGRB). The observation of an SGRB in coincidence with gravitational radiation from an NS-BH coalescence would confirm the association between the two phenomena and also give us new insights on NS physics. We present here a new method to measure NS radii and thus constrain the NS equation of state using joint SGRB and GW observations of NS-BH mergers. We show that in the event of a joint detection with realistic GW signal-to-noise ratio (S/N) of 10, the NS radius can be constrained to $lesssim,$20% accuracy at 90% confidence.
The combined detection of a gravitational-wave signal, kilonova, and short gamma-ray burst (sGRB) from GW170817 marked a scientific breakthrough in the field of multi-messenger astronomy. But even before GW170817, there have been a number of sGRBs with possible associated kilonova detections. In this work, we re-examine these historical sGRB afterglows with a combination of state-of-the-art afterglow and kilonova models. This allows us to include optical/near-infrared synchrotron emission produced by the sGRB as well as ultraviolet/optical/near-infrared emission powered by the radioactive decay of $r$-process elements (i.e., the kilonova). Fitting the lightcurves, we derive the velocity and the mass distribution as well as the composition of the ejected material. The posteriors on kilonova parameters obtained from the fit were turned into distributions for the peak magnitude of the kilonova emission in different bands and the time at which this peak occurs. From the sGRB with an associated kilonova, we found that the peak magnitude in H bands falls in the range [-16.2, -13.1] ($95%$ of confidence) and occurs within $0.8-3.6,rm days$ after the sGRB prompt emission. In g band instead we obtain a peak magnitude in range [-16.8, -12.3] occurring within the first $18,rm hr$ after the sGRB prompt. From the luminosity distributions of GW170817/AT2017gfo, kilonova candidates GRB130603B, GRB050709 and GRB060614 (with the possible inclusion of GRB150101B) and the upper limits from all the other sGRBs not associated with any kilonova detection we obtain for the first time a kilonova luminosity function in different bands.
Using 72 Short Gamma Ray Bursts (SGRBs) with well determined spectral data observed by BATSE, we determine their redshift and the luminosity by applying $E_p$--$L_p$ correlation for SGRBs found by cite{tsutsui13}. For 53 SGRBs with the observed flux brighter than $4 times 10^{-6}~{rm erg~cm^{-2}s^{-1}}$, the cumulative redshift distribution up to $z=1$ agrees well with that of 22 {it Swift}~SGRBs. This suggests that the redshift determination by the $E_p$--$L_p$ correlation for SGRBs works well. The minimum event rate at $z=0$ is estimated as $rho_{SGRB}(0) = 6.3_{-3.9}^{+3.1} times 10^{-10}~{rm events~Mpc^{-3}yr^{-1}}$ so that the minimum beaming angle is $0.6^circ-7.8^circ$ assuming the merging rate of $10^{-7}-4times 10^{-6}~{rm events~Mpc^{-3}yr^{-1}}$ suggested from the binary pulsar data. Interestingly, this angle is consistent with that for SGRB130603B of $sim 4^circ-8^circ$citep{fong13b}. On the other hand, if we assume the beaming angle of $sim 6^circ$ suggested from four SGRBs with the observed value of beaming angle, the minimum event rate including off-axis SGRBs is estimated as $rho_{SGRB,all}^{min}(0)=1.15_{-0.71}^{+0.57}times 10^{-7}~{rm events~Mpc^{-3}yr^{-1}}$. If SGRBs are induced by coalescence of binary neutron stars (NSs) and/or black holes (BHs), this event rate leads to the minimum gravitational-wave detection rate of $rm 3.9_{-2.4}^{+1.9} (152_{-94}^{+75})~events~y^{-1}$ for NS-NS (NS-BH) binary, respectively, by a worldwide network with KAGRA, advanced-LIGO, advanced-Virgo, and GEO.
The equation of state plays a critical role in the physics of the merger of two neutron stars. Recent numerical simulations with microphysical equation of state suggest the outcome of such events depends on the mass of the neutron stars. For less massive systems, simulations favor the formation of a hypermassive, quasi-stable neutron star, whose oscillations produce a short, high frequency burst of gravitational radiation. Its dominant frequency content is tightly correlated with the radius of the neutron star, and its measurement can be used to constrain the supranuclear equation of state. In contrast, the merger of higher mass systems results in prompt gravitational collapse to a black hole. We have developed an algorithm which combines waveform reconstruction from a morphology-independent search for gravitational wave transients with Bayesian model selection, to discriminate between post-merger scenarios and accurately measure the dominant oscillation frequency. We demonstrate the efficacy of the method using a catalogue of simulated binary merger signals in data from LIGO and Virgo, and we discuss the prospects for this analysis in advanced ground-based gravitational wave detectors. From the waveforms considered in this work and assuming an optimally oriented source, we find that the post-merger neutron star signal may be detectable by this technique to $sim 10text{--}25$,Mpc. We also find that we successfully discriminate between the post-merger scenarios with $sim 95%$ accuracy and determine the dominant oscillation frequency of surviving post-merger neutron stars to within $sim 10$,Hz, averaged over all detected signals. This leads to an uncertainty in the estimated radius of a non-rotating 1.6,M$_{odot}$ reference neutron star of $sim 100,$m.
We use realistic Monte-Carlo simulations including both gravitational-wave and short gamma-ray burst selection effects to revisit the coincident rate of binary systems composed of two neutron stars or a neutron star and a black hole. We show that the fraction of GW triggers that can be observed in coincidence with sGRBs is proportional to the beaming factor at $z=0$, but increases with the distance, until it reaches 100 % at the GW detector horizon distance. When this is taken into account the rate is improved by a factor of $~3$ compared to the simple beaming factor correction. We provide an estimate of the performance future GRB detectors should achieve in order to fully exploit the potentiality of the planned third generation GW antenna Einstein Telescope, and we propose a simple method to constrain the beaming angle of sGRBs.