Constraints On Short, Hard Gamma-Ray Burst Beaming Angles From Gravitational Wave Observations

The first detection of a binary neutron star merger, GW170817, and an associated short gamma-ray burst confirmed that neutron star mergers are responsible for at least some of these bursts. The prompt gamma ray emission from these events is thought to be highly relativistically beamed. We present a method for inferring limits on the extent of this beaming by comparing the number of short gamma-ray bursts observed electromagnetically to the number of neutron star binary mergers detected in gravitational waves. We demonstrate that an observing run comparable to the expected Advanced LIGO 2016--2017 run would be capable of placing limits on the beaming angle of approximately $\theta \in (2.88^\circ,14.15^\circ)$, given one binary neutron star detection. We anticipate that after a year of observations with Advanced LIGO at design sensitivity in 2020 these constraints would improve to $\theta \in (8.10^\circ,14.95^\circ)$.

The first observation of a NS-NS coalescence event, GW170817 Abbott et al. (2017a), its association with GRB 170817A Abbott et al. (2017b); Goldstein et al. (2017); Savchenko et al. (2017) and, later, multiwavelength electromagnetic emission, including a kilonova Abbott et al. (2017), confirmed that compact binary mergers are the engines of at least some SGRBs. The gravitational wave (GW) observation placed only weak constraints on the viewing angle due to a degeneracy between distance and inclination of the binary to the line of sight Abbott et al. (2017a). However, GRB 170817A was not typical of SGRBs, being around 10 4 times less energetic Goldstein et al. (2017). This, in addition to other aspects of the electromagnetic emission, has been widely interpreted as indicating that GRB 170817A was not viewed from within the cone of a canonical jet with a top-hat profile (see e.g. Fong et al. (2017); Kasliwal et al. (2017); Gottlieb et al. (2017); Haggard et al. (2017)).
A population of GW-SGRB observations could also allow us to measure the fraction of SGRBs associated with each progenitor type, and associated redshifts will enable a relatively systematics-free measurement of the Hubble parameter at low redshift, which would provide constraints on cosmological models Schutz (1986); Nissanke et al. (2010); Chen & Holz (2013); Abbott et al. (2017).
In this work we consider a population of binary merger sources, with and without SGRB counterparts, assuming that the vast majority of these counterparts would be viewed from within the cone of a standard jet. This is motivated by the fact that most mergers would be expected to occur at distances much greater than GW170817 and that weak, off-axis gamma ray emission would in all likelihood go undetected. With such a population we can constrain the average opening angle Chen & Holz (2013); Clark et al. (2015); Abbott et al. (2016a).
We investigate what statements can currently be made on the beaming angle itself using the bounds placed on the binary merger rate R from all-sky, all-time GW searches and explore the potential for direct inference of SGRB beaming angles in the advanced detector era. We first discuss the relationship between SGRBs and compact binary coalescences. In particular, we will focus on NS-NS inspirals as the progenitors of SGRBs. We then present our method for robustly inferring the jet opening angles using only GW observations. We demonstrate our method assuming the nominal number of GW signals observed from NS-NS inspirals expected for Advanced LIGO (aLIGO) and Advanced Virgo in planned observing scenarios, as defined in Abbott et al. (2013). Finally, we conclude with a discussion on the implications of our work as well as possible avenues for further extension of the work presented here.

SHORT GAMMA-RAY BURSTS AND COMPACT BINARY COALESCENCES
At their design sensitivities, the current generation of advanced GW detectors could observe NS-NS mergers out to distances of ∼400 Mpc at a rate of 0.1-200 yr −1 Abbott et al. (2013). It is worth noting that at galactic or near-galactic distances, soft gammaray repeater (SGR) hyperflares can resemble SGRBs. These SGR hyperflares are the likely explanations for GRB 070201 and GRB 051103, since compact binary coalescences at the distance of their probable host galaxies were excluded with greater than 90% confidence Abbott et al. (2008); Abadie et al. (2012a).
Given the link between SGRBs and compact binary coalescences, it is interesting to ask whether the SGRB beaming angle can be inferred from GW observations. As discussed in Clark et al. (2015), a comparison of the populations of observed SGRBs and NS-NS mergers may be the most promising avenue for this. Motivated by the study in Chen & Holz (2013), we note that if the SGRB population posseses a distribution of beaming angles then the observed rate of SGRBs is related to the rate of NS-NS coalescences R via, where angled brackets indicate the population mean and is the probability that a binary coalescence results in an observed SGRB. In this work, we assume an illustrative R grb = 10 Gpc −3 yr −1 Nakar (2007); Dietz (2011) and we shall refer to as the SGRB efficiency. The method we present, however, is amenable to using alternative values for R grb , or indeed to being extended to sampling values from a prior distribution on the SGRB rate. Generally, the efficiency with which NS-NS mergers produce SGRBs is unknown but will depend on a variety of progenitor physics. In particular, a significant fraction of NS-BH systems may be incapable of powering an SGRB Pannarale & Ohme (2014). Combining this knowledge with measurements of the binary parameters of of a population of GW-SGRB observations could be used to constrain . In this work, we will make no attempt to characterise and we simply aim to provide a framework which allows one to incorporate various levels of assumptions (or ignorance) regarding its value.
If the SGRB population has a distribution of beaming angles, as would seem likely from electromagnetic (EM) observations Fong et al. (2015), characterising the relative rates of SGRB and NS-NS coalescence will inform us as to the mean of that population, θ . To explore this point further we construct a simple Monte Carlo simulation to study the effect on the relative rates of SGRBs and NS-NS mergers. We arrange the following toy problem: 1. Set the number of 'observed' SGRBs to zero: N GRB = 0.
2. Draw N NS−NS values of orbital inclination ι from a distribution which is uniform in cos ι in the range [0, 1].
3. For each value of ι, draw a value for the beaming angle θ, from some distribution with finite width and limited to the range (0, 90] • .
4. If ι < θ then this combination of orbital inclination and beaming angle would result in an observable SGRB, so increment N GRB .
Such a simulation allows us to study the ratio of the number of observed SGRB to the total number of NS-NS mergers N GRB /N NS−NS . Since it is the comparison of the rates of these events that informs our inference on θ, studying the ratio N GRB /N NS−NS provides some intuition as to the effect and features of various θ distributions. Figure 1 plots this ratio as a function of various truncated normal distributions to demonstrate the effect of shifting the mean and scaling the width of the distribution. Points along the x-axis correspond to different choices of the distribution width σ θ , and the separate curves correspond to different choices of the distribution mean θ . Let us denote this truncated normal distribution N ( θ , σ θ ). We stress here that such θ distributions are not intended to represent the true distribution; they are merely intended to easily demonstrate the qualitative effects of different θ distributions on the ratio N GRB /N NS−NS . Figure 1 reveals that a population of SGRB beaming angles with a large mean but narrow width is, on the basis of rate measurements, indistinguishable from a population of SGRB beaming angles with a small mean and large width. For example, for the N (15, 9) and N (10, 13) beaming angle populations, the ratios of N GRB /N NS−NS are almost equal (∼ 4.8%). Thus, a sufficiently wide spread of SGRB beaming angles will yield relatively high rates for NS-NS and SGRBs that could lead to an overestimate of the mean beaming angle. The population-based constraints on θ must, therefore, be regarded as upper bounds on the mean of a distribution of beaming angles. Having said this, for a given mean value θ , the ratio is rather insensitive to the width.

FROM RATES TO BEAMING ANGLES
In this section, we discuss our approach to estimating the SGRB beaming angle based on the binary neutron star inspiral rate, estimated through a number of GW observations of NS-NS coalescence. We demonstrate the approach by considering plausible detection scenarios for aLIGO Abbott et al. (2013). Our ultimate goal is to develop a generic approach that folds in uncertainties in the NS-NS merger rate and our ignorance about the probability with which such mergers actually result in SGRBs. An overview of the general method is as follows: 1. Estimate the posterior probability distribution on the NS-NS merger rate in the local universe from a number of observed gravitational wave signals and our knowledge of the sensitivity of the detectors. We construct a joint posterior distribution on the NS-NS rate and the (unknown) probability that a given merger results in a SGRB.
2. Use equation 1, which relates the NS-NS merger and SGRB rates via the geometry of the beaming angle, to transform the rate posterior probability to a posterior probability on the mean SGRB beaming angle.
3. Marginalize over . We choose to consider a nuisance parameter because, to date, there is no accurate estimate of this parameter and it is not the main focus of our analysis.

Constructing The Rate Posterior
Our goal is to infer the posterior probability distribution for the mean SGRB beaming angle θ from GW constraints on the rate of NS-NS coalescence R. The core ingredient to the analysis is the posterior probability distribution on the coalescence rate p(R|D, I), where D represents some GW observation and I denotes other unenumerated prior information. We will first demonstrate how p(R|D, I) may be constructed for a few projected observing scenarios from Abbott et al. (2013). Later, in section 5, we will extend the analysis to place limits on θ based upon the lack of detection during O1. Previously, a comparison of rates was used to place a lower limit on the beaming angle in Abbott et al. (2016a).
To form the posterior on the coalescence rate, we begin by constructing the posterior on the signal rate. Note that these are not identical since only those NS-NS mergers which occur within a certain range yield a detectable signal. GW data analysis pipelines (e.g.  2017)) identify discrete 'candidate events' which are characterised by network signal-to-noise ratios (SNRs), ρ c , which, for the case of NS-NS searches, indicate the similarity between the detector data and a set of template NS-NS coalescence waveforms. The measured rate r of these events consists of two components: a population of true GW signals, s; and a background rate, b, due to noise fluctuations due to instrumental and environmental disturbances.
Typically for an all-sky, all-time analysis, like that described in Usman et al. (2016), the significance of a candidate event is empirically measured against 'background' data representative of the detector noise, which naturally varies from candidate to candidate. A detection requires this significance to be above some pre-determined threshold (e.g. 5σ for GW150914 and GW151226 Abbott et al. (2016b,c)). We follow the method in Abbott et al. (2013), which defines a detection as a candidate with ρ c ≥ 12, corresponding approximately to b = 10 −2 yr −1 . Since the background rate b is known, we are just left with the problem of inferring the signal rate s. Assuming a uniform prior on s and a Poisson process underlying the events, it may be shown (e.g., Gregory (2010)) that the posterior for the signal rate, given a known background rate b and n events observed over a time period T is, where, Finally, we can transform the posterior on the signal rate to the underlying coalescence rate via our knowledge of the sensitivity of the GW analysis. In particular, the signal detection rate is simply the product of the intrinsic coalescence rate R and the number of NS-NS mergers which would result in a GW signal with ρ c ≥ 12. Expressing the binary coalescence rate in terms of the number of mergers per Milky Way Equivalent Galaxy (MWEG), per year then we require the number of galaxies N G which may be probed by the GW analysis. At large distances, this is well approximated by LIGO Scientific Collaboration & Virgo Collaboration (2010), where D hor is the horizon distance (defined as the distance at which an optimally-oriented NS-NS merger yields ρ c ≥ 12), the factor of 2.26 results from averaging over sky-locations and orientations, and 1.16 × 10 −2 Mpc −3 is the extrapolated density of MWEG in space.
Finally, the posterior on the binary coalescence rate R is obtained from a trivial transformation of the posterior on the signal rate s, We see that in this approach, the rate posterior depends only on the number of signal detections n, the observation time T , the background rate b, and the horizon distance of the search D hor . It is precisely these quantities that comprise the detection scenarios outlined in Abbott et al. (2013). Before constructing expected rate posteriors, we outline the transformation from rate to beaming angle.

Constructing the beaming angle posterior
Inferences of the SGRB beaming angle are made from the posterior probability density on the beaming angle p(θ|D, I) where, as usual, D indicates some set of observations and I unenumerated prior knowledge. Our goal is to transform the measured posterior probability density on the rate R to a posterior on the beaming angle. First, note that we can express the joint distribution p(θ, |D, I) as a Jacobian transformation of the joint distribution p(R, |D, I): where we have dropped conditioning statements for notational convenience. The Jacobian determinant can be computed from equation 1. It is then straightforward to marginalize over to yield the posterior on θ itself: where we have assumed and R are logically independent such that, It is important to note that the entire procedure of deriving the jet angle posterior is completely independent of the approach used to derive the rate posterior. In the preceding section we adopted a straightforward Bayesian analysis of a Poisson rate which is amenable to a simple application of plausible future detection scenarios; there is no inherent requirement to use that method to derive the rate posterior. Given the posterior on the rate, p(R), the final ingredient in this approach is the specification of some prior distribution for . Given the lack of information on the value and distribution of , we choose three plausible priors and study their effects on our beaming angle inference. Our choice of priors are: Delta-function: p( ) = δ( = 0.5); the probability that NS-NS mergers yield SGRBs is known to be 50% exactly.
Uniform: p( ) = U (0, 1); the probability that NS-NS mergers yield SGRBs may lie anywhere ∈ (0, 1] with equal support in that range. Jeffreys: p( ) = β( 1 2 , 1 2 ); treating the outcome of a NS-NS merger as a Bernoulli trial in which a SGRB constitutes 'success' and is the probability of that success, the least informative prior, as derived from the square root of the determinant of the Fisher information for the Bernoulli distribution, is a β-distribution with shape parameters α = β = 1 2 .

PROSPECTS FOR BEAMING ANGLE CONSTRAINTS WITH ADVANCED LIGO
We now demonstrate the derivation of the rate posterior p(R) and the subsequent transformation to the beaming angle posterior p(θ). We consider four GW observation scenarios with aLIGO based on the work in Abbott et al. (2013). An observing scenario essentially consists of an epoch of aLIGO operation, which defines an expected search sensitivity (i.e., NS-NS horizon distance D hor ) and observation time T ; as well as an assumption on the rate of NS-NS coalescence in the local universe R. Each observing scenario ultimately results in an expectation for the number of observed GWs from NS-NS coalescences. For this study, we assume the 'realistic rate' for R as described in LIGO Scientific Collaboration & Virgo Collaboration (2010).
Our first goal is to establish the expected number of detections in each scenario. Given the observation time and horizon distance of the observation epoch we first compute the 4-volume accessible to the analysis, where the factor 2.26 arises from averaging over source sky location and orientation, T is the observation time and γ is the duty cycle for the science run. Following Abbott et al. (2013), we take γ = 0.5. For comparison, during the first observing run of aLIGO, the two interferometers observed in coincidence achieving γ coinc = 0.41.
Where there is a range in the horizon distances quoted in Abbott et al. (2013) to account for uncertainty in the sensitivity of the early configuration of the detectors, we use the arithmetic mean of the lower and upper bounds when computing the search volume.  Abbott et al. (2013). While the 2020+ and 2024+ scenarios appear identical in terms of the sensitivity of the detectors, the 2024+ scenario includes a third advanced LIGO detector in India. This expansion of the network is expected to lead to an increase in the network duty cycle, and an increase in the area of the sky which the network is sensitive to, resulting in a greater volume being searched per year. Figure 2 shows the NS-NS rate posteriors resulting from the observations in the scenarios in table 1 generated using the procedure described in section 3.1. Where a range of potential inspiral distances is given for a scenario we choose the median value, so for the 2015-2016 scenario we take D insp to be 60 Mpc, for example. Likewise we choose an illustrative value of n, the number of expected GW detections, from each range; these are listed in table 2.

Posterior Results
We now use these posteriors together with the prior distributions described in section 3.1 and the observed rate of SGRBs (as described in section 2, we use R grb = 10 Gpc −3 yr −1 Nakar (2007); Dietz (2011)) to derive the corresponding beaming angle posteriors.

Validation
Before we derive beaming angle posteriors corresponding to the aforementioned observing scenarios, it is useful to establish some form of validation for our procedure. This validation is performed by first selecting values of the beaming angle, the SGRB efficiency, and the rate of NS-NS coalescence. We choose θ = 10 • , and the 'realistic' NS-NS rate R = 10 −6 Mpc −3 yr −1 . We then compute the value of the SGRB rate that would correspond to these parameter choices. Finally, we simply   Table 2. Summary of the NS-NS rate posteriors for each of the observing scenarios which are considered in this work; these posteriors are plotted in figure 2. Here n is the number of GW events which were assumed to be observed in each scenario, chosen from the ranges in table 1.
use this artificial value for R grb in equation 10 when we compute the posterior on the beaming angle, with the understanding that the resulting posterior should yield an inference consistent with the 'true' value θ = 10 • . Figures 3 and 4 show the beaming angle posteriors which result from this analysis for the 2015-2016 and 2016-2017 scenarios respectively for each choice of prior distribution on the efficiency parameter. Unsurprisingly, the most accurate constraints arise when we already have the tightest possible constraints on the SGRB efficiency, . That is, the beaming angle posterior arising from the δ-function prior on is the narrowest, yielding the shortest possible credible interval. It is well worth remembering, however, that had we been incorrect regarding the value of when using the δ-function prior,  Figure 3. In order to validate the algorithm an artificial scenario was constructed with a known beaming angle by artificially setting a GRB rate of 36.7 yr −1 to induce a beaming angle of θ ≈ 10 • . The algorithm was then tested with the various priors used in the analysis, using the same horizon distance, observing time, and duty cycle as the 2015-2016 observing scenario. to ensure that the correct beaming angle was inferred. These posteriors are based on the simulated 2015-2016 observing scenario (see table 1).

Prior
Lower MAP Median Upper       the result would be significantly biased and our inference on the beaming angle would be incorrect. This highlights the necessity of building a suitable representation of our ignorance into the analysis. Finally, we note that the results from the uniform and Jeffreys distribution priors are broadly equivalent. Figures 5 and 6 show the beaming angle posteriors obtained for two of the detection scenarios. 12 Since it is a common assumption in related literature, we also now include a prior on the SGRB efficiency which dictates that all NS-NS produce a SGRB, p( |I) = δ( = 1), as well as our previous strong δ-function prior. For the 2016-2017 scenario where inferences are somewhat weak (i.e., broad posteriors) due to the sparsity of GW detections, the uncertainties are large enough that the results from each prior are broadly consistent. In the 2024+ scenario, where the posterior is more peaked, it is clear that the strong δ-function priors lead to inconsistent inferences on the SGRB beaming angle. The much weaker uniform and β distributions, by contrast, are again largely consistent with each other yielding more conservative and robust results, as well as being a more representative expression of our state of knowledge. The inferences drawn from each scenario and each prior are summarised in terms of the maximum a posteriori measurement and the 95% credible interval around the maximum in table 5.

BEAMING ANGLE CONSTRAINTS WITH NO GW DETECTIONS
While GW170817 provided a situation where GW signals from a NS-NS coalescence event were observed, our proposed approach is also valid in the regime where no GW signals from NS-NS coalescence have been observed, as was true during the first observing run of the advanced LIGO detectors when upper limits on binary merger rates were used to place lower limits on the beaming angle Abbott et al. (2016a).
In this scenario, our procedure is identical to before: construct the posterior probability density function on the NS-NS coalescence rate, transform to the joint posterior on the beaming angle and SGRB efficiency, , and marginalise over the nuisance parameter to yield the posterior on the beaming angle. Now, however, rather than quoting the maximum a posteriori estimate, together with some credible interval, we simply integrate the beaming angle posterior from θ = 0 until we reach that value which contains some desired confidence. Thus, we obtain an upper limit on the beaming angle, analogous to the rate upper limits set by past LIGO observations Abadie et al. (2012b). Figure 5 shows the four posteriors on the beaming angle, corresponding to the four priors on the SGRB efficiency, , using the observing 2015-2016 observing scenario from table 1, which corresponds closely to the conditions of the first science run of the advanced generation of ground based GW detectors. We define the upper limit on the beaming angle as the upper limit of the 95% credible interval where the credible interval is defined as the narrowest interval (θ ll , θ ul ) which satisfies the expression with p(θ|D, I) the posterior over which the interval is computed.
Similarly we define the lower limit as the lower limit (2.5 percentile) of the same credible interval. In this non-detection scenario, we choose to compute upper limit on the 95% credible interval on the beaming angle.
We see that here, where the rate posterior is rather uninformative, the results are dominated by the uncertainty in : there are substantive differences in the beaming angle upper limits yielded by the uniform (U (0, 1)) and β-distribution priors, while the δ-function priors yield dramatically different upper limits. Indeed, the most stringent (and mutually incompatible) upper limits are obtained using the strong δ-function priors. In fact, these beaming angle upper limits are also incompatible with the values of 3 • -8 • that have been inferred from observations of jet breaks in SGRB afterglows Fong et al. (2014);Panaitescu (2006);Nicuesa Guelbenzu et al. (2012). Recall, however, from the discussion in section 2 that we interpret the beaming angle inference from our rate measurements as the upper bound on the mean of a population of beaming angles. It would, therefore, seem premature to conclude that there is tension in these results; instead, we can only state that either the population of SGRBs have a distribution of beaming angles with some finite width or that the fraction of NS-NS mergers which yield a SGRB is smaller than 0.5.
It is also interesting to compare these upper limits on the beaming angle with those in Chen & Holz (2013), where the upper limit on the rate itself is used as a constraint (rather than transforming the posterior). This has the important implication that the constraint thus obtained is the smallest angle consistent with the rate: where R ul is the upper limit on the NS-NS rate. The same idea is used in Clark et al. (2015) to estimate beaming constraints in the advanced detector era. Thus, when comparing the constraints in e.g., Chen & Holz (2013) and the upper limits obtained from the transformed posterior (i.e., equation 10 and figure 7), one should remember that they are quite different quantities. There are two other noteworthy differences between Chen & Holz (2013) and this work: (i) the rate upper limit is computed based on the sensitivity of the initial LIGO-Virgo network (see e.g., Brady & Fairhurst (2008)), which gives R = 4.5×10 −4 Mpc −3 yr −1 (as compared with R = 1.3×10 −4 Mpc −3 yr −1 from the analysis in Abadie et al. (2012b)); and (ii) it is implicitly assumed that all NS-NS mergers yield an SGRB. That is, there is no factor or to account for the unknown fraction of mergers which successfully launch an SGRB jet. With these differences noted, the lower bound on the beaming angle is found to be θ ≥ 0.8 • , as compared with the  Figure 7. The upper-bound of the 95% credible interval on the beaming angle as a function of the rate of observed gravitational wave NS-NS events and the observed search 4-volume, taking a Jeffreys prior on the efficiency of GRB production from NS-NS events. The search volumes corresponding to observing scenarios are marked as vertical lines on the plot, with each line assuming that observations are carried out over the period of one year, achieving the search volume outlined in table 1.

BEAMING ANGLE CONSTRAINTS IN FUTURE SCENARIOS
With the advent of gravitational wave astronomy, and with the expectation of the detection of NS-NS gravitational wave signals during the lifetime of the advanced detectors it will become possible to place further constraints on the 95% credible interval of the SGRB beaming angle, as both the searched 4-volume of space increases, and the observed rate of gravitational wave NS-NS events is established.
In figure 7 we present the inferred upper-limit on the 95% credible interval for a range of search 4-volumes and gravitational wave event rates; overlayed on this plot are indications of the anticipated annual search volume for the advanced LIGO detectors in each of the observing scenarios detailed in table 1. These limits were determined by assuming a Jeffreys prior on the efficiency parameter of the model, and following the same procedure used to produce the posteriors in figures 5 and 6. In figure 8 we present a similar plot, showing the upper limits on the beaming angle under the stronger assumption that every NS-NS event also produces a GRB.
The lower limit (the 2.5% of the posterior) for the same range of scenarios is plotted in figure 9, with the same anticipated detector search volumes plotted, again assuming a Jeffreys prior on the efficiency, and in figure 10 we present those lower limits under the assumption that every NS-NS event produces a GRB.  Figure 9. The lower-bound of the 95% credible limit on the beaming angle as a function of the observed number of events and the observed search 4-volume, taking a Jeffreys prior on the efficiency of GRB production from NS-NS events. The search volumes corresponding to observing scenarios are marked as vertical lines on the plot.

CONCLUSION
We have presented a Bayesian analysis which demonstrates the ability of the current generation of advanced GW detectors to make observations that allow for the inference of SGRB jet beaming angles. In doing so we have made minimal assumptions about the processes which produce the jet, other than that NS-NS mergers are the progenitors and that, other than for rare nearby cases like GW170817, SGRBs are observed only by observers within the cone of the jet.
We demonstrate that with a year's worth of gravitational wave observations by the 2-detector aLIGO net-  Figure 10.
The lower-bound of the 95% credible limit on the beaming angle as a function of the observed number of events and the observed search 4-volume, assuming that every gravitational wave NS-NS event produces a GRB. The search volumes corresponding to observing scenarios are marked as vertical lines on the plot. work during its 2016-2017 observing run, and assuming a single NS-NS detection, that we can place a lower limit of 2.85 • , and an upper limit of 46.93 • on the jet beaming angle, given an uninformative prior on the efficiency at which NS-NS events produce observable SGRBs. Assuming that all NS-NS produce an observable SGRBs we can narrow these limits to between 2.88 • and 14.15 • .
When the advanced LIGO design sensitivity is acheived in 2020 the observation of 10 NS-NS events in gravitational waves is sufficient to place an upper limit of 56.99 • on the jet beaming angle, and can establishing the limit on the beaming angle to be between 7.82 • and 56.99 • , assuming an uninformative prior on the SGRBs production efficiency. These limits narrow between 8.10 • and 14.95 • if perfect efficiency is assumed.
The authors wish to thank Martin Hendry for many insightful and valuable discussions during the development of this technique. DW is supported by the Science and Technology Research Council (STFC) grant ST/N504075/1. JC acknowledges support from NSF awards PHYS-1505824 and PHYS-1505524SH. ISH is supported by STFC grant ST/L000946/1. This document has been assigned the control number LIGO-P1600102 by the LIGO document control centre.
The python code used to produce this analysis is available as both a Jupyter (iPython) notebook, and plain python scripts via Zonodo (doi: 10.5281/zenodo.1066019), along with the data used to produce figures 7 and 9.
Software: The analysis presented in this manuscript made extensive use of the numpy van der Walt et al.