Search

Some previous studies based on IceCube neutrinos had found intriguing preliminary evidence that some of them might be GRB neutrinos with travel times affected by quantum properties of spacetime delaying them proportionally to their energy, an effect often labeled as "quantum-spacetime-induced in-vacuo dispersion". Those previous studies looked for candidate GRB neutrinos in a fixed (neutrino-energy-independent) time window after the GRB onset and relied rather crucially on crude estimates of the redshift of GRBs whose redshift has not been measured. We here introduce a complementary approach to the search of quantum-spacetime-affected GRB neutrinos which restricts the analysis to GRBs of sharply known redshift, and, in a way that we argue is synergistic with having sharp information on redshift, adopts a neutrino-energy-dependent time window. We find that knowing the redshift of the GRBs strengthens the analysis enough to compensate for the fact that of course the restriction to GRBs of known redshift reduces the number of candidate GRB neutrinos. And rather remarkably our estimate of the magnitude of the in-vacuo-dispersion effects is fully consistent with what had been found using the previous approach. Our findings are still inconclusive, since their significance is quantified by a $p$-value of little less than $0.01$, but provide motivation for monitoring the accrual of neutrino observations by IceCube and KM3NeT as well as for further refinements of the strategy of analysis here proposed.

We estimate the neutrino flux from different kinds of galactic sources and compare it with the recently diffuse neutrino flux detected by IceCube. We find that the flux from these sources may contribute to ~ 20% of the IceCube neutrino flux. Most of the sources selected in this work populate the southern hemisphere, therefore a detector like KM3NeT could help in resolving the sources out of the observed diffused galactic neutrino flux., Comment: 14 pages, 4 figures, accepted for publication in ApJ

Supernova (SN) 2023ixf in M101 is the closest SN explosion observed in the last decade. Therefore it is a suitable test bed to study the role of jets in powering the SN ejecta. With this aim, we explored the idea that high-energy neutrinos could be produced during the interaction between the jets and the intense radiation field produced in the SN explosion and eventually be observed by the IceCube neutrino telescope. The lack of detection of such neutrinos has significantly constrained both the fraction of stellar collapses that produce jets and/or the theoretical models for neutrino production. Finally, we investigated the possibility of detecting low-energy neutrinos from SN 2023ixf with the Super- and Hyper-Kamiokande experiments, obtaining in both cases sub-threshold estimates., Comment: 9 pages, 2 figures. Accepted for publication in Astrophysical Journal Letters

The emission from shock breakouts (SBOs) represents the earliest electromagnetic (EM) signal emitted by cataclysmic events involving the formation or the merger of neutron stars (NSs). As such, SBOs carry unique information on the structure of their progenitors and on the explosion energy. The characteristic~SBO emission is expected in the UV range, and its detection is one of the key targets of~the ULTRASAT satellite. Among SBO sources, we focus on a specific class involving the formation of fast spinning magnetars in the core-collapse (CC) of massive stars. Fast spinning magnetars are expected to produce a specific signature in the early UV supernova light curve, powered by the extra spin energy quickly released by the NS. Moreover, they are considered as optimal candidates for the emission of long-transient gravitational wave (GW) signals, the detection of which requires early EM triggers to boost the sensitivity of dedicated GW search pipelines. We calculate early supernova UV light curves in the presence of a magnetar central engine, as a function of the explosion energy, ejecta mass and magnetar parameters. We then estimate the ULTRASAT detection horizon (z < 0.15) as a function of the same physical parameters, and the overall expected detection rate finding that magnetar-powered SBOs may represent up to 1/5 of the total events detected by ULTRASAT. Moreover, at the expected sensitivity of the LIGO/Virgo/Kagra O5 science run, one such event occurring within 5 Mpc will providean ideal trigger for a GW long transient search. Future GW detectors like the Einstein Telescope will push the horizon for joint EM-GW detections to 35-40 Mpc., Comment: 7 pages, 3 figures. Accepted for publication in ApJ

The recently discovered high energy emission from the recurrent nova RS Ophiuchi by Fermi-LAT ($>$ 100 MeV), H.E.S.S. and MAGIC ($>$ 100 GeV), hints towards a possible hadronic origin of this radiation component. From the observed high energy photon flux we derive the expected number of neutrino events that could be detected by present and future neutrino telescopes in the different energy ranges. We find that both hadronic and leptonic processes remain valid interpretations for this $ \ gamma $ emission. Preliminary estimates indicate that with the "next-generation" instrument IceCube-Gen2, the expected number of neutrino detections from Galactic novae is of the order of $\sim $ once per decade. Given the current uncertainties in the frequency of the occurrence of nova outbursts, the detection rate may possibly increase to up to once every $\sim $ three years., Comment: 14 pages, 2 figures, accepted for publication in JCAP

There exists an inevitable scatter in intrinsic luminosity of Gamma Ray Bursts(GRBs). If there is relativistic beaming in the source, viewing angle variation necessarily introduces variation in the intrinsic luminosity function (ILF). Scatter in the ILF can cause a selection bias where distant sources that are detected have a larger median luminosity than those detected close by. Median luminosity divides any given population into equal halves. When the functional form of a distribution is unknown, it can be a more robust diagnostic than any that use trial functional forms. In this work we employ a statistical test based on median luminosity and apply it to test a class of models for GRBs. We assume that the GRB jet has a finite opening angle and that the orientation of the GRB jet is random relative to the observer.We calculate $L_{median}$ as a function of redshift by simulating GRBs empirically, theoretically and use the luminosity vs redshift {\it Swift} data in order to compare the theoretical results with the observed ones. The method accounts for the fact that at some redshifts there may be some GRBs that go undetected. We find that $L_{median}$ is extremely insensitive to the on-axis (i.e. maximal) luminosity of the jet., Comment: 9 pages, 9 figures, 3 tables, accepted for publication in A&A. arXiv admin note: substantial text overlap with arXiv:2010.04810

The joint detection of gravitational waves and the gamma-ray counterpart of a binary neutron star merger event, GW170817, unambiguously validates the connection between short gamma-ray bursts and compact binary object (CBO) mergers. We focus on a special scenario where short gamma-ray bursts produced by CBO mergers are embedded in disks of active galactic nuclei (AGN), and we investigate the $\gamma$-ray emission produced in the internal dissipation region via synchrotron, synchrotron self-Compton and external inverse-Compton (EIC) processes. In this scenario, isotropic thermal photons from the AGN disks contribute to the EIC component. We show that a low-density cavity can be formed in the migration traps, leading to the embedded mergers producing successful GRB jets. We find that the EIC component would dominate the GeV emission for typical CBO mergers with an isotropic-equivalent luminosity of $L_{j,\rm iso}=10^{48.5}~\rm erg~s^{-1}$ which are located close to the central supermassive black hole. Considering a long-lasting jet of duration $T_{\rm dur}\sim10^2-10^3$ s, we find that the future CTA will be able to detect its $25-100$ GeV emission out to a redshift $z=1.0$. {In the optimistic case, it is possible to detect the on-axis extended emission simultaneously with GWs within one decade using MAGIC, H.E.S.S., VERITAS, CTA, and LHAASO-WCDA. Early diagnosis of prompt emissions with Fermi-GBM and HAWC can provide valuable directional information for the follow-up observations., Comment: 11 pages, 5 figures. Accepted for publication in The Astrophysical Journal

Multi-messenger astrophysics is becoming a major avenue to explore the Universe, with the potential to span a vast range of redshifts. The growing synergies between different probes is opening new frontiers, which promise profound insights into several aspects of fundamental physics and cosmology. In this context, THESEUS will play a central role during the 2030s in detecting and localizing the electromagnetic counterparts of gravitational wave and neutrino sources that the unprecedented sensitivity of next generation detectors will discover at much higher rates than the present. Here, we review the most important target signals from multi-messenger sources that THESEUS will be able to detect and characterize, discussing detection rate expectations and scientific impact., Comment: Submitted to Experimental Astronomy

The strong constraints from the Fermi-LAT data on the isotropic gamma-ray background suggest that the neutrinos observed by IceCube might possibly come from sources that are hidden to gamma-ray observations. A possibility recently discussed in the literature is that neutrinos may come from jets of collapsing massive stars which fail to break out of the stellar envelope, and for this reason they are known as choked jets, or choked Gamma-Ray Bursts (GRBs). In this paper, we estimate the neutrino flux and spectrum expected from these sources, focusing on Type II SNe. We perform detailed calculations of pg interactions, accounting for all the neutrino production channels and scattering angles. We provide predictions of expected event rates for operating neutrino telescopes, such as ANTARES and IceCube, as well as for the future generation telescope KM3NeT. We find that for GRB energies channeled into protons spanning between 10^51 - 10^53 erg, choked GRBs may substantially contribute to the observed astrophysical neutrino flux, if their local rate is 80 - 1 Gpc^-3 yr^-1 respectively., Comment: 24 pages, 5 figures, published on JCAP

Different forms of long gamma-ray bursts (GRBs) Luminosity Functions are considered on the basis of an explicit physical model. The inferred flux distributions are compared with the observed ones from two samples of GRBs, Swift and Fermi GBM. The best fit parameters of the Luminosity functions are found and the physical interpretations are discussed. The results are consistent with the observation of a comparable number of flat phase afterglows and monotonic decreasing ones., Comment: 12 pages, 4 figures, 2 tables, current version accepted by ApJ

There exists an inevitable scatter in intrinsic luminosity of Gamma Ray Bursts(GRBs). If there is relativistic beaming in the source, viewing angle variation necessarily introduces variation in the intrinsic luminosity function(ILF). Scatter in the ILF can cause a selection bias where distant sources that are detected have a larger median luminosity than those detected close by. Median luminosity, as we know, divides any given population into equal halves. When the functional form of a distribution is unknown, it can be a more robust diagnostic than any that use trial functional forms. In this work we employ a statistical test based on median luminosity and apply it to test a class of models for GRBs. We assume that the GRB jet has a finite opening angle and that the orientation of the GRB jet is random relative to the observer. We parameterize the jet with constant Lorentz factor $\Gamma$ and opening angle $\theta_0$. We calculate $L_{median}$ as a function of redshift with an average of 17 grbs in each redshift bin($dz=0.01$) empirically, theoretically and use Fermi GBM data, noting that SWIFT data is problematic as it is biased, specially at high redshifts. We find that $L_{median}$ is close to $L_{max}$ for sufficiently extended GRB jet and does not fit the data. We find an acceptable fit with the data when $\Gamma$ is between $100$ and $200$, $\theta_0\leq 0.1$, provided that the jet material along the line of sight to the on axis observer is optically thick, such that the shielded maximum luminosity is well below the bare $L_{max}$. If we associate an on-axis observer with a classically projected monotonically decreasing afterglow, we find that their ILF is similar to those of off-jet observer which we associate with flat phase afterglows., Comment: Withdrawn.We have completely changed the paper-calculations, figures and our analysis. The authors have also changed as one of the authors, Prof. David EIchler, recently passed away

The joint observation of core-collapse supernovae with gamma-ray bursts shows that jets can be launched in the aftermath of stellar core collapse, likely by a newly formed black hole that accretes matter from the star. Such gamma-ray bursts have only been observed accompanying Type Ibc supernovae, indicating a stellar progenitor that lost its Hydrogen envelope before collapse. According to recent hypothesis it is possible that jets are launched in core-collapse events even when the progenitors still retain their Hydrogen envelopes, however, such jets are not able to burrow through the star and will be stalled into the interior of the progenitor star before escaping. These jets are called chocked jets. High-energy neutrinos produced by such choked jets could escape the stellar envelope and could be observed. Here we examine how multi-messenger searches for high-energy neutrinos and core-collapse supernovae can detect or limit the fraction of stellar collapses that produce jets. We find that a high fraction of jet production is already limited by previous observational campaigns. We explore possibilities with future observations using LSST, IceCube and Km3NET., Comment: 5 pages, accepted for publication in MNRAS

Multimessenger observations may hold the key to learn about the most energetic sources in the universe. The recent construction of large scale observatories opened new possibilities in testing non thermal cosmic processes with alternative probes, such as high energy neutrinos and gravitational waves. We propose to combine information from gravitational wave detections, neutrino observations and electromagnetic signals to obtain a comprehensive picture of some of the most extreme cosmic processes. Gravitational waves are indicative of source dynamics, such as the formation, evolution and interaction of compact objects. These compact objects can play an important role in astrophysical particle acceleration, and are interesting candidates for neutrino and in general high-energy astroparticle studies. In particular we will concentrate on the most promising gravitational wave emitter sources: compact stellar remnants. The merger of binary black holes, binary neutron stars or black hole-neutron star binaries are abundant gravitational wave sources and will likely make up the majority of detections. However, stellar core collapse with rapidly rotating core may also be significant gravitational wave emitter, while slower rotating cores may be detectable only at closer distances. The joint detection of gravitational waves and neutrinos from these sources will probe the physics of the sources and will be a smoking gun of the presence of hadrons in these objects which is still an open question. Conversely, the non-detection of neutrinos or gravitational waves from these sources will be fundamental to constrain the hadronic content., Comment: Published in EPJ Web Conf. Volume 209, 01036, pg.9, 2019; RICAP18, 7th Roma International Conference on Astroparticle Physics

The first, long awaited, detection of a gravitational wave (GW) signal from the merger of a binary neutron-star (NS-NS) system was finally achieved (GW$\,$170817), and was also accompanied by an electromagnetic counterpart -- the short-duration GRB 170817A. It occurred in the nearby ($D\approx40\;$Mpc) elliptical galaxy NGC$\,$4993, and showed optical, IR and UV emission from half a day up to weeks after the event, as well as late time X-ray (at $\geq 8.9\;$days) and radio (at $\geq 16.4\;$days) emission. There was a delay of $\Delta t \approx 1.74\;$s between the GW merger chirp signal and the prompt-GRB emission onset, and an upper limit of $\theta_{\rm obs}<28^\circ$ was set on the viewing angle w.r.t the jet's symmetry axis from the GW signal. In this letter we examine some of the implications of these groundbreaking observations. The delay $\Delta t$ sets an upper limit on the prompt-GRB emission radius, $R_\gamma\lesssim 2c\Delta t/(\theta_{\rm obs}-\theta_0)^2$, for a jet with sharp edges at an angle $\theta_0<\theta_{\rm obs}$. GRB 170817A's relatively low isotropic equivalent $\gamma$-ray energy-output may suggest a viewing angle slightly outside the jet's sharp edge, $\theta_{\rm{}obs}-\theta_0\sim(0.05-0.1)(\Gamma/100)^{-1}$, but its peak $\nu{}F_\nu$ photon energy and afterglow emission suggest instead that the jet does not have sharp edges and the prompt emission was dominated by less energetic material along our line of sight, at $\theta_{\rm{}obs}\gtrsim 2\theta_0$. Finally, we consider the type of remnant that is produced by the NS-NS merger and find that a relatively long-lived ($>2\;$s) massive NS is strongly disfavored, while a hyper-massive NS of lifetime $\sim1\;$s appears to be somewhat favored over the direct formation of a black hole., Comment: 5 pages, 3 figures, minor changes; accepted to ApJL

The short-duration ($\lesssim2\;$s) GRB 170817A in the nearby ($D=40\;$Mpc) elliptical galaxy NGC 4993 is the first electromagnetic counterpart of the first gravitational wave (GW) detection of a binary neutron-star (NS-NS) merger. It was followed by optical, IR, and UV emission from half a day up to weeks after the event, as well as late time X-ray and radio emission. The early UV, optical, and IR emission showed a quasi-thermal spectrum suggestive of radioactive-decay powered kilonova-like emission. Comparison to kilonova models favors the formation of a short-lived ($\sim1\;$s) hypermassive NS, which is also supported by the $\Delta t\approx1.74\;$s delay between the GW chirp signal and the prompt GRB onset. However, the late onset of the X-ray (8.9$\;$days) and radio (16.4$\;$days) emission, together with the low isotropic equivalent $\gamma$-ray energy output ($E_{\rm\gamma,iso}\approx5\times10^{46}\;$erg), strongly suggest emission from a narrow relativistic jet viewed off-axis. Here we set up a general framework for off-axis GRB jet afterglow emission, comparing analytic and numerical approaches, and showing their general predictions for short-hard GRBs that accompany binary NS mergers. The prompt GRB emission suggests a viewing angle well outside the jet's core, and we compare the afterglow lightcurves expected in such a case to the X-ray to radio emission from GRB 170817A. We fit an afterglow off-axis jet model to the X-ray and radio data and find that the observations are explained by a viewing angle $\theta_{\rm obs}\approx16^\circ-26^\circ$, GRB jet energy $E\sim10^{48.5}-10^{49.5}~{\rm erg}$, and external density $n\sim10^{-5}-10^{-1}~{\rm cm}^{-3}$ for a $\xi_e\sim 0.1$ non-thermal electron acceleration efficiency., Comment: 14 pages, 10 figures; Minor changes made to text; Fig. 10 now shows results from least-squares fit instead of MCMC simulation; updated references; Submitted to MNRAS

The origin of the high-energy flux of neutrinos detected by IceCube is still unknown. Recent works report the evidence for a possible positional correlation between the reconstructed neutrino arrival directions and the positions in the sky of low power, high-energy emitting BL Lac objects (HBL). Assuming that $\gamma$-ray emitting HBL form the bulk of the sources of high-energy neutrinos above 100 TeV, we intend to calculate the number of events expected to be detected for each source by IceCube and KM3NeT. Based on a simple theoretically-motivated framework inspired by the structured jet scenario for these sources, we postulate a direct proportionality between high-energy $\gamma$-ray and neutrino fluxes. We calculate the expected neutrino event rate for the HBL sources of the Second Fermi-LAT Catalog of High-Energy Sources (2FHL) for IceCube and the presently under construction KM3NeT using declination-dependent and exposure-weighted effective areas. We provide a list of 2FHL HBL with the calculated number of events. For IceCube, the derived count rate for several sources is relatively high, of the order of $\lesssim$1 yr$^{-1}$, consistently with the recent findings of a possible positional correlation. For KM3NeT the calculated rates are higher, with several sources with expected rate exceeding 1 yr$^{-1}$. This, coupled with the improved angular resolution, implies that the HBL origin can be effectively tested with few years of observation of KM3NeT (and IceCube Gen2, for which similar performances are foreseen) through the direct association of neutrinos and single HBL. Our results show that if -- as hinted by recent works -- HBL represent a possible population of high-energy neutrino emitters, several single sources should be identified in few years of exposure of KM3NeT, highlighting the importance of the improved angular resolution anticipated for KM3NeT and IceCube Gen2., Comment: 8 pages, 3 figures, accepted for publication in A&A

Several Pulsar Wind Nebulae (PWNe) have been detected in the TeV band in the last decade.The TeV emission is typically interpreted in a purely leptonic scenario, but this usually requires that the magnetic field in the Nebula be much lower than the equipartition value and the assumption of an enhanced target radiation at IR frequencies. In this work we consider the possibility that, in addition to the relativistic electrons, also relativistic hadrons are present in these nebulae. Assuming that part of the emitted TeV photons are of hadronic origin, we compute the associated flux of $\sim 1-100$ TeV neutrinos. We use the IceCube non detection to put constraints on the fraction of TeV photons that might be contributed by hadrons and estimate the number of neutrino events that can be expected from these sources in ANTARES and in KM3Net., Comment: Submitted to ApJ, 12 pages, 3 figures

We consider a sample of 107 Gamma Ray Bursts (GRBs) for which early UV emission was measured by Swift, and extrapolate the photon intensity to lower energies. Protons accelerated in the GRB jet may interact with such photons to produce charged pions and subsequently ultra high energy neutrinos $\varepsilon_\nu\geq 10^{16}$ eV. We use simple energy conversion efficiency arguments to predict the maximal neutrino flux expected from each GRB. We estimate the neutrino detection rate at large area radio based neutrino detectors and conclude that the early afterglow neutrino emission is too weak to be detected even by next generation neutrino observatories.

The IceCube telescope has detected diffuse neutrino emission, 20 events of which were reported to be above 60~TeV. In this paper, we fit the diffuse neutrino spectrum using Poisson statistics, which are the most appropriate for the low counts per energy bin. We extend the fitted energy range and exploit the fact that no neutrinos were detected above 2~PeV, despite the high detector sensitivity around the Glashow resonance at 6.3\,PeV and beyond. A best-fit power-law slope of $\alpha=2.9\pm 0.3$ is found with no evidence for a high-energy cutoff. This slope is steeper than $\alpha=2.3\pm 0.3$ found by the IceCube team using a different fitting method. Such a steep spectrum facilitates the identification of high energy ($\gg$ PeV) neutrinos, if detected, to be due to the GZK effect of cosmic-ray protons interacting with the Extragalactic Background Light. We use the ratio of EeV to PeV neutrinos in GZK models to show that the currently detected PeV neutrinos could not be due to the GZK effect, because this would imply many more higher-energy neutrinos that should have been detected, but were not. The non-detection of GZK neutrinos by IceCube despite more than essentially 1200 observing days, has already ruled out (at 95\% confidence) models that predict rates of $\sim1$ neutrino/yr or more. We use this non-detection to quantify the confidence at which GZK models are ruled out, and compute the additional IceCube and (in the future) ARA observing time that would rule them out with 95\% confidence if no detection is made., Comment: 6 pages, 4 figures. Submitted to Astrophysical Journal

The latest IceCube data suggest that the all-flavor cosmic neutrino flux may be as large as 10^-7 GeV/cm^2/s/sr around 30 TeV. We show that, if sources of the TeV-PeV neutrinos are transparent to gamma rays with respect to two-photon annihilation, strong tensions with the isotropic diffuse gamma-ray background measured by Fermi are unavoidable, independently of the production mechanism. We further show that, if the IceCube neutrinos have a photohadronic (pgamma) origin, the sources are expected to be opaque to 1-100 GeV gamma rays. With these general multimessenger arguments, we find that the latest data suggest a population of cosmic-ray accelerators hidden in GeV-TeV gamma rays as a neutrino origin. Searches for x-ray and MeV gamma-ray counterparts are encouraged, and TeV-PeV neutrinos themselves will serve as special probes of dense source environments., Comment: 6+3 pages, 2+1 figures, version to appear in PRL, with expanded Supplementary Material

We compare the isotropic equivalent 15-2000 keV gamma-ray energy, E_gamma, emitted by a sample of 91 swift Gamma-Ray Bursts (GRBs) with known redshifts, with the isotropic equivalent fireball energy, E_fb, as estimated within the fireball model framework from X-ray afterglow observations of these bursts. The uncertainty in E_gamma, which spans the range of ~10^51 erg to ~10^53.5 erg, is approximately 25% on average, due mainly to the extrapolation from the BAT detector band to the 15-2000 keV band. The uncertainty in E_fb is approximately a factor of 2, due mainly to the X-ray measurements' scatter. We find E_gamma and E_fb to be tightly correlated. The average(std) of {\eta}^11hr_gamma is approximately log_10(E_gamma/(3{\epsilon} _eE^11hr_fb)) are -0.34(0.60), and the upper limit on the intrinsic spread of {\eta}_gamma is approximately 0.5 ({\epsilon}_e is the fraction of shocked plasma energy carried by electrons and E^x hr_fb is inferred from the X-ray flux at x hours). We also find that E_fb inferred from X-ray observations at 3 and 11 hours are similar, with an average(std) of log_10(E^3hr_fb/E^11hr_fb) of 0.04(0.28). The small variance of {\eta}_gamma implies that burst-to-burst variations in {\epsilon}_e and in the efficiency of fireball energy conversion to gamma-rays are small, and suggests that both are of order unity. The small variance of {\eta}_gamma and the similarity of E^3hr_fb and E^11hr_fb further imply that {\epsilon}_e does not vary significantly with shock Lorentz factor, and that for most bursts the modification of fireball energy during the afterglow phase, by processes such as radiative losses or extended duration energy injection, are not significant. Finally, our results imply that if fireballs are indeed jets, then the jet opening angle satisfies {\theta}>0.1 for most cases. [abridged], Comment: An expanded discussion of the absence of energy injection and 'flat X-ray phases' in the majority of bursts

In this review I discuss the ultra-high energy neutrinos (UHEN) originated from Cosmic-Rays propogation (GZK neutrinos) and from Gamma Ray Bursts (GRBs), and discuss their detectability in kilometers scale detectors like ARA and IceCube. While GZK neutrinos are expected from cosmic ray interactions on the CMB, the GRB neutrinos depend on the physics inside the sources. GRBs are predicted to emit UHEN in the prompt and in the later 'after-glow' phase. I discuss the constraints on the hadronic component of GRBs derived from the search of four years of IceCube data for a prompt neutrino fux from gamma-ray bursts (GRBs) and more in general I present the results of the search for high-energy neutrinos interacting within the IceCube detector between 2010 and 2013., Comment: Accepted for publication in Journal of High Energy Astrophysics. arXiv admin note: text overlap with arXiv:1407.0155

The IceCube high-energy neutrino telescope has been collecting data since 2006. Conversely, hundreds of Gamma Ray Bursts (GRBs) have been detected by the GBM on board Fermi, since its launch in 2008. So far no neutrino event has been associated with a GRB, despite many models predicting the generation of high energy neutrinos through GRB photon interaction with PeV protons in the GRB jet. We use the non-detection of neutrinos to constrain the hadronic content of GRB jets independent of jet model parameters. Assuming a generic particle spectrum of $E^{-\alpha}$ with $\alpha = 2$, we find that the ratio of the energy carried by pions to that in electrons has to be small $f_\pi / f_e \lesssim 0.24$ at 95\% confidence level. A distribution of spectral slopes can lower $f_\pi / f_e$ by orders of magnitude. Another limit, independent of neutrinos, is obtained if one ascribes the measured Fermi/LAT GeV gamma-ray emission to pair-photon cascades of high-energy photons resulting from (the same photon-hadronic interactions and subsequent) neutral pion decays. Based on the generally observed MeV to GeV GRB fluence ratio of $\approx 10$, we show that $f_\pi / f_e \lesssim 0.3$. In some bursts, where this ratio is as low as unity, $f_\pi / f_e \lesssim 0.03$. These findings add to the mounting doubts regarding the presence of PeV protons in GRB jets., Comment: 6 pages, 5 figures. Submitted to Astrophysical Journal

Most Gamma-Ray Bursts (GRBs) detected by the Fermi Gamma-ray Space Telescope exhibit a delay of up to about 10 seconds between the trigger time of the hard X-ray signal as measured by the Fermi GBM and the onset of the MeV-GeV counterpart detected by the LAT. This delay may hint at important physics, whether it is due to the intrinsic variability of the inner engine or it is related to quantum dispersion effects in the velocity of light propagation from the sources to the observer. It is critical to have a proper assessment of how these time delays affect the overall properties of the light curves. We cross-correlated the 5 brightest GRBs of the 1st Fermi LAT Catalog by means of the continuous correlation function (CCF) and of the Discrete Correlation Function (DCF). A maximum in the DCF suggests the presence of a time lag between the curves, whose value and uncertainty are estimated through a Gaussian fitting of the DCF profile and light curve simulation via a Monte Carlo approach. The cross-correlation of the observed LAT and GBM light curves yields time lags that are mostly similar to those reported in the literature, but they are formally consistent with zero. The cross-correlation of the simulated light curves yields smaller errors on the time lags and more than one time lag for GRBs 090902B and 090926A; for all 5 GRBs, the time lags are significantly different from zero and consistent with those reported in the literature, when only the secondary maxima are considered for those two GRBs. The DCF method evidences the presence of time lags between the LAT and GBM light curves and underlines their complexity. While this suggests that the delays should be ascribed to intrinsic physical mechanisms, more sensitivity and larger statistics are needed to assess whether time lags are universally present in the early GRB emission and which dynamical time scales they trace., Comment: 9 pages, 3 figures, accepted for publication in Astronomy & Astrophysics

We process the Fermi LAT data on GRB130427A using the Fermi Science Tools, and we summarize some of the key facts that render this observation truly remarkable, especially concerning the quality of information on high-energy emission by GRBs. We then exploit this richness for a search of spectral lags, of the type that has been recently of interest for its relevance in quantum-spacetime research. We do find some evidence of systematic soft spectral lags: when confining the analysis to photons of energies greater than 5 GeV there is an early hard development of minibursts within this long burst. The effect turns out to be well characterized quantitatively by a linear dependence, within such a miniburst, of the detection time on energy. With the guidance of our findings for GRB130427A we can then recognize that some support for these features is noticeable also in earlier Fermi-LAT GRBs, particularly for the presence of hard minibursts whose onset is marked by the highest-energy photon observed for the GRB. A comparison of these features for GRBs at different redshifts provides some encouragement for a redshift dependence of the effects of the type expected for a quantum-spacetime interpretation, but other aspects of the analysis appear to invite the interpretation as intrinsic properties of GRBs.

One of the primary goals of neutrino telescopes, such as IceCube, is the discovery of neutrinos emitted by gamma-ray bursts (GRBs). Another source of interest in the results obtained by these telescopes is their possible use for tests of the applicability of Einstein's Special Relativity to neutrinos, particularly with respect to modifications that lead to Lorentz invariance violation that have been conjectured by some models of quantum space-time. We examine here the fascinating scenario in which these two aspects of neutrino-telescope physics require a combined analysis. We discuss how neutrinos that one would not associate to a GRB, when assuming a classical spacetime picture, may well be GRB neutrinos if the possibility that Lorentz invariance is broken at very high energies is taken into account. As an illustrative example we examine three IceCube high energy neutrinos that arrived hours before GRBs (but from the same direction) and we find that the available, IceCube data, while inconclusive, is compatible with a scenario in which one or two of these neutrinos were GRB neutrinos and their earlier arrival reflects Lorentz invariance violation. We outline how future analyses of neutrino data should be done in order to systematically test this possibility., Comment: 6 pages, 2 figures

The Large Area Telescope (LAT) on-board the Fermi satellite detected emission above 20 MeV only in a small fraction of the long gamma-ray bursts (GRBs) detected by the Fermi Gamma-ray Burst Monitor (GBM) at 8 keV-40 MeV. Those bursts that were detected by the LAT were among the brightest GBM bursts. The emission detected by LAT seerms to be delayed respect to the one detected by the GBM. In this review I will show the main implications of these LAT observations on the GRB models and discuss the importance of a synergy between Fermi and other telescopes., Comment: Invited Review, 2012 Fermi Symposium proceedings - eConf C121028

In this work we give an estimate of the neutrino flux that can be expected from the microquasar Cyg X-3. We calculate the muon neutrino flux expected here on Earth as well as the corresponding number of neutrino events in the IceCube telescope based on the so-called hypersoft X-ray state of Cyg~X-3. If the average emission from Cyg~X-3 over a period of 5 yr were as high as during the used X-ray state, a total of 0.8 events should be observed by the full IceCube telescope. We also show that this conclusion holds by a factor of a few when we consider the other measured X-ray states. Using the correlation of AGILE data on the flaring episodes in 2009 June and July to the hypersoft X-ray state we calculate that the upper limits on the neutrino flux given by IceCube are starting to constrain the hadronic models, which have been introduced to interpret the high-energy emission detected by AGILE., Comment: 12 pages, 2 figures; Updated to match version published in ApJ, contains color versions of figures

We employ realistic constraints on astrophysical and instrumental selection effects to model the Gamma-Ray Burst (GRB) redshift distribution using {\it Swift} triggered redshift samples acquired from optical afterglows (OA) and the TOUGH survey. Models for the Malmquist bias, redshift desert, and the fraction of afterglows missing because of host galaxy dust extinction, are used to show how the "true" GRB redshift distribution is distorted to its presently observed biased distribution. We also investigate another selection effect arising from a correlation between $E_{{\rm iso}}$ and $L_{{\rm opt}}$. The analysis, which accounts for the missing fraction of redshifts in the two data subsets, shows that a combination of selection effects (both instrumental and astrophysical) can describe the observed GRB redshift distribution. Furthermore, the observed distribution is compatible with a GRB rate evolution that tracks the global SFR, although the rate at high-$z$ cannot be constrained with confidence. Taking selection effects into account, it is not necessary to invoke high-energy GRB luminosity evolution with redshift to explain the observed GRB rate at high-$z$., Comment: Version 2. Includes new data, figures and refined analysis

It is clear that optical selection effects have distorted the "true" GRB redshift distribution to its presently observed biased distribution. We constrain a statistically optimal model that implies GRB host galaxy dust extinction could account for up to 40% of missing optical afterglows and redshifts in $z = 0-3$, but the bias is negligible at very high-$z$. The limiting sensitivity of the telescopes, and the time to acquire spectroscopic/photometric redshifts, are significant sources of bias for the very high-$z$ sample. We caution on constraining star formation rate and luminosity evolution using the GRB redshift distribution without accounting for these selection effects., Comment: 4 pages, 1 figure, Poster: to appear in: Proceedings of Science, Gamma-Ray Bursts 2012 Conference - GRB2012, May 07-11, 2012, Munich Germany

Presently only 30% of short gamma ray bursts (SGRBs) have accurate redshifts, and this sample is highly biased by the limited sensitivity of {\it Swift} to detect SGRBs. We account for the dominant biases to calculate a realistic SGRB rate density out to $z = 0.5$ using the {\it Swift} sample of peak fluxes, redshifts, and those SGRBs with a beaming angle constraint from X-ray/optical observations. Assuming a significant fraction of binary neutron star mergers produce SGRBs, we calculate lower and upper detection rate limits of (1-180) per Yr by an advanced LIGO and Virgo coincidence search. Our detection rate is compatible with extrapolations using Galactic pulsar observations and population synthesis., Comment: 6 pages, 1 figure, to appear in: Proceedings of Science, Gamma-Ray Bursts 2012 Conference - GRB2012, May 07-11, 2012, Munich Germany

Short gamma-ray bursts (SGRBs) observed by {\it Swift} are potentially revealing the first insight into cataclysmic compact object mergers. To ultimately acquire a fundamental understanding of these events requires pan-spectral observations and knowledge of their spatial distribution to differentiate between proposed progenitor populations. Up to April 2012 there are only some 30% of SGRBs with reasonably firm redshifts, and this sample is highly biased by the limited sensitivity of {\it Swift} to detect SGRBs. We account for the dominant biases to calculate a realistic SGRB rate density out to $z\approx0.5$ using the {\it Swift} sample of peak fluxes, redshifts, and those SGRBs with a beaming angle constraint from X-ray/optical observations. We find an SGRB lower rate density of $8^{+5}_{-3}$ $\mathrm{Gpc}^{-3}\mathrm{yr}^{-1}$ (assuming isotropic emission), and a beaming corrected upper limit of $1100^{+700}_{-470} $ $\mathrm{Gpc}^{-3}\mathrm{yr}^{-1}$. Assuming a significant fraction of binary neutron star mergers produce SGRBs, we calculate lower and upper detection rate limits of $(1-180)$ yr$^{-1}$ by an aLIGO and Virgo coincidence search. Our detection rate is similar to the lower and realistic rates inferred from extrapolations using Galactic pulsar observations and population synthesis., Comment: version 2: 7 pages, 1 figure, updated table & references, submitted to MNRAS

The statistical distribution of energies among particles responsible for long Gamma Ray Burst (GRB) emission is analyzed in light of recent results of the Fermi Observatory. The allsky flux, $F_{\gamma}$, recorded by the Gamma Ray Burst Monitor (GBM) is shown, despite its larger energy range, to be not significantly larger than that reported by the Burst and Transient Explorer (BATSE), suggesting a relatively small flux in the 3 - 30 MeV energy range. The present-day energy input rate in $\gamma$-rays recorded by the GBM from long GRB is found, assuming star-formation rates in the literature, to be $\dot W(0)=0.5 F_{\gamma} H/c = 5 \times 10^{42}\ \rm{erg/Mpc^3 yr}$. The Large Area Telescope (LAT) fluence, when observed, is about 5-10\% per decade of the total, in good agreement with the predictions of saturated, non-linear shock acceleration. The high-energy component of long GRBs, as measured by Fermi, is found to contain only $\sim 10^{-2.5}$ of the energy needed to produce ultrahigh-energy cosmic rays (UHECR) above 4 Eev, assuming the latter to be extragalactic, when various numerical factors are carefully included, if the cosmic ray source spectrum has a spectral index of -2. The observed $\gamma$-ray fraction of the required UHECR energy is even smaller if the source spectrum is softer than $E^{-2}$. The AMANDA II limits rule out such a GRB origin for UHECR if much more than $10^{-2}$ of the cosmic ray energy goes into neutrinos that are within, and simultaneous with, the $\gamma$-ray beam. It is suggested that "orphan" neutrinos out of the $\gamma$-ray beam might be identifiable via orphan afterglow { or other wide angle signatures of GRB in lieu of coincidence with prompt $\gamma$-rays}, and it is recommended that feasible single neutrino trigger criteria be established to search for such coincidences., Comment: to appear in The Astrophysical Journal

We investigate a model for the shallow decay phases of Gamma-ray Burst (GRB) afterglows discovered by Swift/XRT in the first hours following a GRB event. In the context of the fireball scenario, we consider the possibility that long-lived energy injection from a millisecond spinning, ultramagnetic neutron star (magnetar) powers afterglow emission during this phase. We consider the energy evolution in a relativistic shock subject to both radiative losses and energy injection from a spinning down magnetar in spherical symmetry. We model the energy injection term through magnetic dipole losses and discuss an approximate treatment for the dynamical evolution of the blastwave. We obtain an analytic solution for the energy evolution in the shock and associated lightcurves. To fully illustrate the potential of our solution we calculate lightcurves for a few selected X-ray afterglows observed by Swift and fit them using our theoretical lightcurves. Our solution naturally describes in a single picture the properties of the shallow decay phase and the transition to the so-called normal decay phase. In particular, we obtain remarkably good fits to X-ray afterglows for plausible parameters of the magnetar. Even though approximate, our treatment provides a step forward with respect to previously adopted approximations and provides additional support to the idea that a millisecond spinning (1-3 ms), ultramagnetic (B$\sim 10^{14}-10^{15}$ G) neutron star loosing spin energy through magnetic dipole radiation can explain the luminosity, durations and shapes of X-ray GRB afterglows., Comment: 7 pages, 2 figures, submitted to Astronomy & Astrophysics - referee's comments included

We investigate the constraints imposed on the luminosity function (LF) of long duration Gamma Ray Bursts (LGRBs) by the flux distribution of bursts detected by the GBM at ~1 MeV, and the implications of the non detection of the vast majority, ~95%, of the LGRBs at higher energy, ~1 GeV, by the LAT detector. We find a LF that is consistent with those determined by BATSE and Swift. The non detections by LAT set upper limits on the ratio R of the prompt fluence at ~1 GeV to that at ~1 MeV. The upper limits are more stringent for brighter bursts, with R<{0.1,0.3,1} for {5,30,60}% of the bursts. This implies that for most bursts the prompt ~1 GeV emission may be comparable to the ~1 MeV emission, but can not dominate it. The value of R is not universal, with a spread of (at least) an order of magnitude around R~10^(-1). For several bright bursts with reliable determination of the photon spectral index at ~1 MeV, the LAT non detection implies an upper limit to the ~100 MeV flux which is <0.1 of the flux obtained by extrapolating the ~1 MeV flux to high energy. For the widely accepted models, in which the ~1 MeV power-law photon spectrum reflects the power-law energy distribution of fast cooling electrons, this suggests that either the electron energy distribution does not follow a power-law over a wide energy range, or that the high energy photons are absorbed. Requiring an order unity pair production optical depth at ~100 MeV sets an upper limit for the Lorentz factor, Gamma<=10^(2.5)., Comment: 12 pages, 6 figures. Submitted to A&A

Gravitational waves (GW) can be emitted from coalescing neutron star (NS) and black hole-neutron star (BH-NS) binaries, which are thought to be the sources of short hard gamma ray bursts (SHBs). The gamma ray fireballs seem to be beamed into a small solid angle and therefore only a fraction of detectable GW events is expected to be observationally coincident with SHBs. Similarly ultrahigh energy (UHE) neutrino signals associated with gamma ray bursts (GRBs) could fail to be corroborated by prompt gamma-ray emission if the latter is beamed in a narrower cone than the neutrinos. Alternative ways to corroborate non-electromagnetic signals from coalescing neutron stars are therefore all the more desirable. It is noted here that the extended X-ray tails (XRT) of SHBs are similar to X-ray flashes (XRFs), and that both can be attributed to an off-axis line of sight and thus span a larger solid angle than the hard emission. It is proposed that a higher fraction of detectable GW events may be coincident with XRF/XRT than with hard gamma-rays, thereby enhancing the possibility to detect it as a GW or neutrino source. Scattered gamma-rays, which may subtend a much larger solid angle that the primary gamma ray jet, are also candidates for corroborating non-electromagnetic signals., Comment: 13 pages, accepted for publication in Astrophysical Journal Letters

We investigate the origin of the prompt and delayed emission observed in the short GRB 090510. We use the broad-band data to test whether the most popular theoretical models for gamma-ray burst emission can accommodate the observations for this burst. We first attempt to explain the soft-to-hard spectral evolution associated with the delayed onset of a GeV tail with the hypothesis that the prompt burst and the high energy tail both originate from a single process, namely synchrotron emission from internal shocks. Considerations on the compactness of the source imply that the high-energy tail should be produced in a late-emitted shell, characterized by a Lorentz factor greater than the one generating the prompt burst. However, in this hypothesis, the predicted evolution of the synchrotron peak frequency does not agree with the observed soft-to-hard evolution. Given the difficulties of a single-mechanism hypothesis, we test two alternative double-component scenarios. In the first, the prompt burst is explained as synchrotron radiation from internal shocks, and the high energy emission (up to about 1 s following the trigger) as internal shock synchrotron-self-Compton. In the second scenario, in view of its long duration (\sim 100 s), the high energy tail is decoupled from the prompt burst and has an external shock origin. In this case, we show that a reasonable choice of parameters does indeed exist to accommodate the optical-to-GeV data, provided the Lorentz factor of the shocked shell is sufficiently high. Finally, we attempt to explain the chromatic break observed around \sim 1e3 s with a structured jet model. We find that this might be a viable explanation, and that it lowers the high value of the burst energy derived assuming isotropy, \sim 1e53 erg, below \sim 1e49 erg, more compatible with the energetics from a binary merger progenitor., Comment: 19 pages, 3 color figures, accepted for publication in the Astrophysical Journal. New version addresses Referee's comments

We compare the luminosity function and rate inferred from the GBM long bursts peak flux distribution with those inferred from the Swift and BATSE peak flux distribution. We find that the GBM, BATSE and the Swift peak fluxes can be fitted by the same luminosity function implying the consistency of these three samples. Using the trigger algorithm of the LAT instrument we derive important information on the flux at 100 MeV compared to lower energy detected by the GBM. We find that the simple extension of the synchrotron emission to high energy cannot justify the low rate of GRBs detected by LAT and for several GRBs detected by the GBM, the flux at >100 MeV should be suppressed. Two bursts, GRB090217 and GRB 090202B, detected by LAT have very soft spectra in the GBM and therefore their high energy emission cannot be due to an extension of the synchrotron., Comment: 12 pages, 7 figures. Proceeding for the workshop "7th Agile Meeting & The Bright Gamma-Ray Sky" Frascati, Italy, 29 September- 1 October 2009

We investigate whether the high energy tail detected by the Fermi/LAT for the short GRB 081024B can be caused by synchrotron and self-Compton emission in the context of either the internal or external shock models. For the internal shock scenario, we explore the possibility of generating the high energy photons directly by means of the synchrotron process, or inverse Compton emission in which target photons are synchrotron photons produced in internal shocks taking place in a lately emitted shell (delayed internal shocks). In the external shock scenario, we test whether the high energy tail can be an extension of the afterglow synchrotron emission, or alternatively the inverse Compton component associated with the afterglow synchrotron photons. For the internal shock scenario, we conclude that only an inverse Compton component from delayed internal shocks can explain the high energy tail that extends to the GeV range. In the external shock scenario, we show that the high energy tail may be interpreted as synchrotron afterglow emission, if the slow cooling phase starts as early as a few seconds after the trigger. On the other hand, the observed high energy tail is consistent with an inverse Compton component of the afterglow in the fast cooling regime., Comment: 8 pages, 3 figures; new version addresses Referee's comments; accepted for publication in Astronomy and Astrophysics

Merging binary systems consisting of two collapsed objects are among the most promising sources of high frequency gravitational wave, GW, signals for ground based interferometers. Double neutron star or black hole/neutron star mergers are also believed to give rise to short hard bursts, SHBs, a subclass of gamma ray bursts. SHBs might thus provide a powerful way to infer the merger rate of two-collapsed object binaries. Under the hypothesis that most SHBs originate from double neutron star or black hole/neutron star mergers, we outline here a method to estimate the incidence of merging events from dynamically formed binaries in globular clusters and infer the corresponding GW event rate that can be detected with Advanced LIGO/Virgo. In particular a sizeable fraction of detectable GW events is expected to be coincident with SHBs. The beaming and redshift distribution of SHBs are reassessed and their luminosity function constrained by using the results from recent SHBs observations. We confirm that a substantial fraction of SHBs goes off at low redshifts, where the merging of systems formed in globular clusters through dynamical interactions is expected., Comment: 6 pages, 1 figure. Accepted for publication in Astronomy and Astrophysics. Some quotations and references changed

It is noted that X-ray tails (XRTs) of short, hard $\gamma$-ray bursts (SHBs) are similar to X-ray flashes (XRFs). We suggest a universal central engine hypothesis, as a way of accounting for this curiosity, in which SHBs differ from long $\gamma$-ray bursts (GRBs) in prompt emission because of the differences in the host star and attendant differences in the environment they present to the compact central engine (as opposed to differences in the central engine itself). Observational constraints and implications are discussed, especially for confirming putative detections of gravitational waves from merging compact objects., Comment: 3 figures, one with two panels

We compare the luminosity function and rate inferred from the BATSE long bursts peak flux distribution with those inferred from the Swift peak flux distribution. We find that both the BATSE and the Swift peak fluxes can be fitted by the same luminosity function and the two samples are compatible with a population that follows the star formation rate. The estimated local long GRB rate (without beaming corrections) varies by a factor of five from 0.05 Gpc^(-3)yr^(-1) for a rate function that has a large fraction of high redshift bursts to 0.27 Gpc^(-3)yr^(-1) for a rate function that has many local ones. We then turn to compare the BeppoSax/HETE2 and the Swift observed redshift distributions and compare them with the predictions of the luminosity function found. We find that the discrepancy between the BeppoSax/HETE2 and Swift observed redshift distributions is only partially explained by the different thresholds of the detectors and it may indicate strong selection effects. After trying different forms of the star formation rate (SFR) we find that the observed Swift redshift distribution, with more observed high redshift bursts than expected, is inconsistent with a GRB rate that simply follows current models for the SFR. We show that this can be explained by GRB evolution beyond the SFR (more high redshift bursts). Alternatively this can also arise if the luminosity function evolves and earlier bursts were more luminous or if strong selection effects affect the redshift determination., Comment: 15 pages, 8 figures, accepted for publication in JCAP

We measure the local rates of ``low-luminosity'' (LL-GRBs, i.e. L<10^{48--49}erg/sec) and ``high-luminosity'' Gamma-ray Bursts (HL-GRBs). The values are in the range n0=100--1800 Gpc^(-3)yr^(-1) and n0=100--550 Gpc^(-3) yr^(-1), respectively, and the ratios to SNe-Ibc \sim 1%-9% and 0.4% -3%. These data may suggest the existence of two physically distinct classes of GRBs in which LL-GRBs are (intrinsically) more frequent events than HL-GRBs. However, with the present data we cannot exclude the possibility of a single population of GRBs which give rise to both an isotropic low-luminous emission (LL-GRBs: detectable only in nearby GRBs) and to a highly collimated high-luminous emission (HL-GRBs: detectable preferentially at high-z). We compute also the rate of SNe-Ibc characterized by broad-lined spectra (Hypernovae) and found it to be about 1.5\times 10^(-4)HNe yr^(-1) 10^(10) L_(B\odot) (i.e less than 10% of SNe-Ibc occurring in Spirals). This result implies that the ratio HL-GRBs/HNe is smaller than 1, possibly in the range 0.04--0.3. We have used the ratio between Hypernovae and LL-GRBs to constrain their beaming factor to f_b^(-1)\sim 10 or less., Comment: 13 pages, 2 figures, ApJ Letters in press

We compare the luminosity function and rate inferred from the BATSE short hard bursts (SHBs) peak flux distribution with the redshift and luminosity distributions of SHBs observed by Swift/HETE II. The Swift/HETE II SHB sample is incompatible with SHB population that follows the star formation rate. However, it is compatible with a distribution of delay times after the SFR. This would be the case if SHBs are associated with the mergers of double neutron star (DNS) systems. DNS may be ``primordial'' or can form dynamically by binary exchange interaction in globular clusters during core collapse. The implied SHB rates that we find range from \sim 8 to \sim 30h_(70)^3 Gpc^(-3)yr^(-1). This rate is a much higher than what was previously estimated and, when beaming is taken into account, it is comparable to the rate of neutron star mergers estimated from statistics of binary pulsars. If GRBs are produced in mergers the implied rate practically guarantees detection by LIGO II and possibly even by LIGO I., Comment: 6 pages, 2figures. Submitted to the proceedings of the Swift-Venice 2006 meeting to be published by "Il Nuovo Cimento"

Short-hard gamma-ray bursts (SHBs) may arise from gravitational wave (GW) driven mergers of double neutron star (DNS) systems. DNSs may be "primordial" or can form dynamically by binary exchange interactions in globular clusters during core-collapse. For primordial binaries, the time delay between formation and merger is expected to be short, tau~0.1 Gyr, implying that the redshift distribution of merger events should follow that of star-formation. We point out here that for dynamically formed DNSs, the time delay between star-formation and merger is dominated by the cluster core-collapse time, rather than by the GW inspiral time, yielding delays comparable to the Hubble time. We derive the redshift distribution of merger events of dynamically formed DNSs, and find it to differ significantly from that typically expected for primordial binaries. The observed redshift distribution of SHBs favors dynamical formation, although a primordial origin cannot be ruled out due to possible detection biases. Future red-shift observations of SHBs may allow to determine whether they are dominated by primordial or dynamically formed DNSs., Comment: Accepted to ApJL; Minor revisions

We compare the luminosity function and rate inferred from the BATSE short hard bursts (SHBs) peak flux distribution with the redshift and luminosity distributions of SHBs observed by Swift/HETE II. The Swift/HETE II SHB sample is incompatible with SHB population that follows the star formation rate. However, it is compatible with a distribution of delay times after the SFR. This would be the case if SHBs are associated with binary neutron star mergers. The implied SHB rates that we find range from \sim 8 to \sim 30h_{70}^3Gpc^{-3}yr^{-1}. This rate is a much higher than what was previously estimated and, when beaming is taken into account, it is comparable to the rate of neutron star mergers estimated from statistics of binary pulsars. If GRBs are produced in mergers the implied rate practically guarantees detection by LIGO II and possibly even by LIGO I, if we are lucky. Our analysis, which is based on observed short hard burst is limited to bursts with luminosities above 10^{49}erg/sec. Weaker bursts may exist but if so they are hardly detected by BATSE or Swift and hence their rate is very weakly constrained by current observations. Thus the rate of mergers that lead to a detection of a gravitational radiation signal might be even higher., Comment: Proceedings of the conference on "Gamma Ray Bursts in the Swift Era", November 29,-December 2, Washington, DC

We compare the luminosity function and rate inferred from the BATSE short hard bursts (SHBs) peak flux distribution with the redshift and luminosity distributions of SHBs observed by Swift/HETE II. While the Swift/HETE II SHB sample is incompatible with SHB population that follows the star formation rate, it is compatible with a SHB rate that reflect a distribution of delay times after the SFR. This would be the case if SHBs are associated with binary neutron star mergers. The available data allows, however, different interpretations. For example, a population whose rate is independent of the redshift fits the data very well. The implied SHB rates that we find range from $\sim 8$ to $\sim 30h_{70}^3$Gpc$^{-3}$yr$^{-1}$. This is a much higher rate than what was previously estimated. A detailed (2 dimensional) look at the best fit models shows, however, some discrepancy between the four Swift/HETE II SHBs and the models based on BATSE SHBs. This could be a statistical fluke. It could also arise from wrong estimates of the triggering criteria or from selection effects. If real it may indicate the existence of two SHB populations with different luminosity functions and redshift distributions., Comment: Revised version includes additional bursts and revised redshift of older bursts