We estimate the rate and the luminosity function of short (hard) Gamma-Ray Bursts (sGRBs) that are non-Collapsars, using the peak fluxes and redshifts of BATSE, Swift and Fermi GRBs. Following Bromberg2013 we select a sub-sample of Swift bursts which are most likely non-Collapsars. We find that these sGRBs are delayed relative to the global star formation rate (SFR) with a typical delay time of a 3-4 Gyr (depending on the SFR model). However, if two or three sGRB at high redshifts have been missed because of selection effects, a distribution of delay times of ~1/t would be also compatible. The current event rate of these non-Collapsar sGRBs with L_iso > 5*10^49 erg/s is 4.1(-1.9,+2.3)Gpc^-3 yr^-1. The rate was significantly larger around z ~ 1 and it declines since that time. The luminosity function we find is a broken power law with a break at 2.0(-0.4,+1.4) * 10^52~erg/s and power-law indices 0.95(-0.12,+0.12) and 2.0(-0.8,+1.0). When considering the whole Swift sGRB sample we find that it is composed of two populations: One group (~ 60%-80% of Swift sGRBs) with the above rate and time delay and a second group (~ 20%-40% of Swift sGRBs) of potential "impostors" that follow the SFR with no delay. These two populations are in very good agreement with the division of sGRBs to non-Collapsars and Collapsars suggested recently by Bromberg2013. If non-Collapsar sGRBs arise from neutron star merger this rate suggest a detection rate of 3-100 yr^-1 by a future gravitational wave detectors (e.g. Advanced Ligo/Virgo with detection horizon on 300 Mpc), and a co-detection with Fermi (Swift) rate of 0.1-1 yr^-1 (0.02-0.14 yr^-1). We estimate that about 4 * 10^5 (f_b^-1 / 30) mergers took place in the Milky Way. If $0.025 m_\odot$ were ejected in each event this would have been sufficient to produce all the heavy r-process material in the Galaxy.

The new and extreme population of GRBs detected by Fermi-LAT shows several new features in high energy gamma-rays that are providing interesting and unexpected clues into GRB prompt and afterglow emission mechanisms. Over the last 6 years, it has been Swift that has provided the robust dataset of UV/optical and X-ray afterglow observations that opened many windows into components of GRB emission structure. The relationship between the LAT detected GRBs and the well studied, fainter, less energetic GRBs detected by Swift-BAT is only beginning to be explored by multi-wavelength studies. We explore the large sample of GRBs detected by BAT only, BAT and Fermi-GBM, and GBM and LAT, focusing on these samples separately in order to search for statistically significant differences between the populations, using only those GRBs with measured redshifts in order to physically characterize these objects. We disentangle which differences are instrumental selection effects versus intrinsic properties, in order to better understand the nature of the special characteristics of the LAT bursts.

We invert directly the redshift - luminosity distribution of observed long Swift GRBs to obtain their rate and luminosity function. Our best fit rate is described by a broken power law that rises like (1+z)^2.1{+0.5-0.6} for 0<z<3 and decrease like (1+z)^-1.4{+2.4-1.0} for z>3. The local rate is 1.3^{+0.6-0.7} [Gpc^-3 yr^-1]. The luminosity function is well described by a broken power law with a break at L* = 10^52.5{+-0.2}[erg/sec] and with indices alpha = 0.2^{+0.2-0.1} and beta = 1.4^{+0.3-0.6}. The recently detected GRB 090423, with redshift ~8, fits nicely into the model's prediction, verifying that we are allowed to extend our results to high redshifts. While there is a possible agreement with the star formation rate (SFR) for z<3, the high redshift slope is shallower than the steep decline in the SFR for 4<z. However we cannot rule out a GRB rate that follows one of the recent SFR models.