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.

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.

We compute the local rate of events similar to GRB 170817A, which has been recently found to be associated with a kilonova (KN) outburst. Our analysis finds an observed rate of such events of R$_{KN}\sim 352^{+810}_{-281}$ Gpc$^{-3}$yr$^{-1}$. After comparing at their face values this density of sGRB outbursts with the much higher density of Binary Neutron Star (BNS) mergers of 1540$^{+3200}_{-1220}$ Gpc$^{-3}$yr$^{-1}$, estimated by LIGO-Virgo collaboration, one can conclude, admittedly with large uncertainty that either only a minor fraction of BNS mergers produces sGRB/KN events or the sGRBs associated with BNS mergers are beamed and observable under viewing angles as large as $\theta$ $\leq$ $40^{\circ}$. Finally we provide preliminary estimates of the number of sGRB/KN events detected by future surveys carried out with present/future ground-based/space facilities, such as LSST, VST, ZTF, SKA and THESEUS.

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.

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.

We derive the contribution to the extragalactic gamma-ray background (EGB) from AGN winds and star-forming galaxies by including a physical model for the gamma-ray emission produced by relativistic protons accelerated by AGN-driven and supernova-driven shocks into a state-of-the-art semi-analytic model of galaxy formation. This is based on galaxy interactions as triggers of AGN accretion and starburst activity and on expanding blast wave as the mechanism to communicate outwards the energy injected into the interstellar medium by the active nucleus. We compare the model predictions with the latest measurement of the EGB spectrum performed by the Fermi-LAT in the range between 100 MeV and 820 GeV. We find that AGN winds can provide ~35$\pm$15% of the observed EGB in the energy interval E_{\gamma}=0.1-1 GeV, for ~73$\pm$15% at E_{\gamma}=1-10 GeV, and for ~60$\pm$20% at E_{\gamma}>10 GeV. The AGN wind contribution to the EGB is predicted to be larger by a factor of 3-5 than that provided by star-forming galaxies (quiescent plus starburst) in the hierarchical clustering scenario. The cumulative gamma-ray emission from AGN winds and blazars can account for the amplitude and spectral shape of the EGB, assuming the standard acceleration theory, and AGN wind parameters that agree with observations. We also compare the model prediction for the cumulative neutrino background from AGN winds with the most recent IceCube data. We find that for AGN winds with accelerated proton spectral index p=2.2-2.3, and taking into account internal absorption of gamma-rays, the Fermi-LAT and IceCube data could be reproduced simultaneously.

We report on a search for ultra-high-energy (UHE) neutrinos from gamma-ray bursts (GRBs) in the data set collected by the Testbed station of the Askaryan Radio Array (ARA) in 2011 and 2012. From 57 selected GRBs, we observed no events that survive our cuts, which is consistent with 0.12 expected background events. Using NeuCosmA as a numerical GRB reference emission model, we estimate upper limits on the prompt UHE GRB neutrino fluence and quasi-diffuse flux from $10^{7}$ to $10^{10}$ GeV. This is the first limit on the prompt UHE GRB neutrino quasi-diffuse flux above $10^{7}$ GeV.

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.

We compute the non-thermal emissions produced by relativistic particles accelerated by the AGN-driven shocks in NGC 1068, and we compare the model predictions with the observed gamma-ray and radio spectra . The former is contributed by pion decay, inverse Compton scattering, and bremsstrahlung, while the latter is produced by synchrotron radiation. We derive the gamma-ray and radio emissions by assuming the standard acceleration theory, and we discuss how our results compare with those corresponding to other commonly assumed sources of gamma-ray and radio emissions, like Supernova remnants (SNR) or AGN jets. We find that the AGN-driven shocks observed in the circumnuclear molecular disk of such a galaxy provide a contribution to the gamma-ray emission comparable to that provided by the starburst activity when standard particle acceleration efficiencies are assumed, while they can yield the whole gamma-ray emission only when the parameters describing the acceleration efficiency and the proton coupling with the molecular gas are tuned to values larger than those assumed in standard, SNR-driven shocks. We discuss the range of acceleration efficiencies (for protons and electrons) and of proton calorimetric fractions required to account for the observed gamma-ray emission in the AGN outflow model. We further compare the neutrino flux expected in our model with constraints from current experiments, and we provide predictions for the detections by the upcoming KM3NeT neutrino telescope. This analysis strongly motivates observations of NGC 1068 at >TeV energies with current and future Cherenkov telescopes in order to gain insight into the nature of the gamma-rays source.

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.

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.

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]

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.

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.

This is a White Paper in support of the mission concept of the Large Observatory for X-ray Timing (LOFT), proposed as a medium-sized ESA mission. We discuss the potential of LOFT for the study of gamma-ray bursts. For a summary, we refer to the paper.

The origin of the high-energy neutrinos (0.1$-$1 PeV range) detected by IceCube remains a mystery. In this work we explore the possibility that efficient neutrino production can occur in structured jets of BL Lac objects, characterized by a fast inner spine surrounded by a slower layer. This scenario has been widely discussed in the framework of the high-energy emission models for BL Lacs and radiogalaxies. One of the relevant consequences of a velocity structure is the enhancement of the inverse Compton emission caused by the radiative coupling of the two zones. We show that a similar boosting could occur for the neutrino output of the spine through the photo-meson reaction of high-energy protons scattering off the amplified soft target photon field of the layer. Assuming the local density and the cosmological evolution of $\gamma$-ray BL Lac derived from {\it Fermi}-LAT data, we calculate the expected diffuse neutrino intensity, that can match the IceCube data for a reasonable choice of the parameters.

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.

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.

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.

Many of the astrophysical sources and violent phenomena observed in our Universe are potential emitters of gravitational waves (GW) and high-energy neutrinos (HEN). Both GWs and HENs may escape very dense media and travel unaffected over cosmological distances, carrying information from the innermost regions of the astrophysical engines. Such messengers could also reveal new, hidden sources that have not been observed by conventional photon-based astronomy. Coincident observation of GWs and HENs may thus play a critical role in multimessenger astronomy. This is particularly true at the present time owing to the advent of a new generation of dedicated detectors: IceCube, ANTARES, VIRGO and LIGO. Given the complexity of the instruments, a successful joint analysis of this data set will be possible only if the expertise and knowledge of the data is shared between the two communities. This review aims at providing an overview of both theoretical and experimental state-of-the-art and perspectives for such a GW+HEN multimessenger astronomy.

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.

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.

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$.

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.

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.

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.

The Large Area Telescope (LAT) on board of the Fermi satellite detected emission above 30 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 - 10 MeV. Those bursts that were detected by the LAT were among the brightest GBM bursts. We examine a sample of the most luminous GBM bursts with no LAT detection and obtain upper limits on their high energy fluence. We find an average upper limit of LAT/GBM fluence ratio of 0.13 for GeV fluence during $T_{90}$ and an average upper limit ratio of 0.45 for GeV fluence during the first 600 seconds after the trigger. These ratios strongly constrain various emission models and in particular rule out SSC models for the prompt emission. In about a third of both LAT detected and LAT non-detected bursts, we find that the extrapolation of the MeV range Band spectrum to the GeV range is larger than the observed GeV fluence (or its upper limit). While this excess is not highly significant for any specific burst, the overall excess in a large fraction of the bursts suggests a decline in the high energy spectral slope in at least some of these bursts. Possibly an evidence for the long sought after pair creation limit.

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.

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.

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.

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.

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).

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.

Gamma-ray bursts (GRBSs) are produced by rare types of massive stellar explosions. Their rapidly fading afterglows are often bright enough at optical wavelengths, that they are detectable up to cosmological distances. Hirtheto, the highest known redshift for a GRB was z=6.7, for GRB 080913, and for a galaxy was z=6.96. Here we report observations of GRB 090423 and the near-infrared spectroscopic measurement of its redshift z=8.1^{+0.1}_{-0.3}. This burst happened when the Universe was only ~4% of its current age. Its properties are similar to those of GRBs observed at low/intermediate redshifts, suggesting that the mechanisms and progenitors that gave rise to this burst about 600 million years after the Big Bang are not markedly different from those producing GRBs ~10 billion years later.

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.

One of the most intriguing features discovered by Swift is a plateau phase in the X-ray flux decay of about 70% of the afterglows of gamma-ray bursts (GRBs). The physical origin of this feature is still being debated. We constrain the proposed interpretations, based on the intrinsic temporal properties of the plateau phase. We selected and analyzed all the Swift/XRT GRB afterglows at known redshift observed between March 2005 and June 2008 featuring a shallow decay phase in their X-ray lightcurves. For our sample of 21 GRBs we find an anticorrelation of the logarithm of the duration of the shallow phase with re dshift, with a Spearman rank-order correlation coefficient of r=-0.4 and a null hypothesis probability of 5%. When we correct the durations for cosmological dilation, the anticorrelation strenghtens, with r=-0.6 and a null hypothesis probability of 0.4%. Considering only those GRBs in our sample that have a well-measured burst peak energy (8 out of 21), we find an anticorrelation between the energy of the burst and the shallow phase duration, with r=-0.80 and a null hypothesis probability of 1.8%. If the burst energy anticorrelation with the shallow phase duration is real, then the dependence of the shallow phase on redshift could be the result of a selection effect, since on average high-redshift bursts with lower energies and longer plateaus would be missed. A burst energy anticorrelation with the shallow phase duration would be expected if the end of the plateau arises from a collimated outflow. Alternative scenarios are briefly discussed involving a possible cosmological evolution of the mechanism responsible for the X-ray shallow decay.

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.

GRB080319B reached 5th optical magnitude during the burst prompt emission. Thanks to the VLT/UVES rapid response mode, we observed its afterglow just 8m:30s after the GRB onset when the magnitude was R ~ 12. This allowed us to obtain the best signal-to-noise, high resolution spectrum of a GRB afterglow ever (S/N per resolution element ~ 50). The spectrum is rich of absorption features belonging to the main system at z=0.937, divided in at least six components spanning a total velocity range of 100 km/s. The VLT/UVES observations caught the absorbing gas in a highly excited state, producing the strongest Fe II fine structure lines ever observed in a GRB. A few hours later the optical depth of these lines was reduced by a factor of 4-20, and the optical/UV flux by a factor of ~ 60. This proves that the excitation of the observed fine structure lines is due to "pumping" by the GRB UV photons. A comparison of the observed ratio between the number of photons absorbed by the excited state and those in the Fe II ground state suggests that the six absorbers are ~ 2-6 kpc from the GRB site, with component I ~ 3 times closer to the GRB site than components III to VI. Component I is characterized also by the lack of Mg I absorption, unlike all other components. This may be due both to a closer distance and a lower density, suggesting a structured ISM in this galaxy complex.

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.

We present the results from an ESO/VLT campaign aimed at studying the afterglow properties of the short/hard gamma ray burst GRB 070707. Observations were carried out at ten different epochs from ~0.5 to ~80 days after the event. The optical flux decayed steeply with a power-law decay index greater than 3, later levelling off at R~27.3 mag; this is likely the emission level of the host galaxy, the faintest yet detected for a short GRB. Spectroscopic observations did not reveal any line features/edges that could unambiguously pinpoint the GRB redshift, but set a limit z < 3.6. In the range of allowed redshifts, the host has a low luminosity, comparable to that of long-duration GRBs. The existence of such faint host galaxies suggests caution when associating short GRBs with bright, offset galaxies, where the true host might just be too dim for detection. The steepness of the decay of the optical afterglow of GRB 070707 challenges external shock models for the optical afterglow of short/hard GRBs. We argue that this behaviour might results from prolonged activity of the central engine or require alternative scenarios.

We present and discuss the results of an extensive observational campaign devoted to GRB071010A, a long-duration gamma-ray burst detected by the Swift satellite. This event was followed for almost a month in the optical/near-infrared (NIR) with various telescopes starting from about 2min after the high-energy event. Swift-XRT observations started only later at about 0.4d. The light-curve evolution allows us to single out an initial rising phase with a maximum at about 7min, possibly the afterglow onset in the context of the standard fireball model, which is then followed by a smooth decay interrupted by a sharp rebrightening at about 0.6d. The rebrightening was visible in both the optical/NIR and X-rays and can be interpreted as an episode of discrete energy injection, although various alternatives are possible. A steepening of the afterglow light curve is recorded at about 1d. The entire evolution of the optical/NIR afterglow is consistent with being achromatic. This could be one of the few identified GRB afterglows with an achromatic break in the X-ray through the optical/NIR bands. Polarimetry was also obtained at about 1d, just after the rebrightening and almost coincident with the steepening. This provided a fairly tight upper limit of 0.9% for the polarized-flux fraction.

In the redshift range z = 0-1, the gamma ray burst (GRB) redshift distribution should increase rapidly because of increasing differential volume sizes and strong evolution in the star formation rate. This feature is not observed in the Swift redshift distribution and to account for this discrepancy, a dominant bias, independent of the Swift sensitivity, is required. Furthermore, despite rapid localization, about 40-50% of Swift and pre-Swift GRBs do not have a measured redshift. We employ a heuristic technique to extract this redshift bias using 66 GRBs localized by Swift with redshifts determined from absorption or emission spectroscopy. For the Swift and HETE+BeppoSAX redshift distributions, the best model fit to the bias in z < 1 implies that if GRB rate evolution follows the SFR, the bias cancels this rate increase. We find that the same bias is affecting both Swift and HETE+BeppoSAX measurements similarly in z < 1. Using a bias model constrained at a 98% KS probability, we find that 72% of GRBs in z < 2 will not have measurable redshifts and about 55% in z > 2. To achieve this high KS probability requires increasing the GRB rate density in small z compared to the high-z rate. This provides further evidence for a low-luminosity population of GRBs that are observed in only a small volume because of their faintness.

In this paper we study synchrotron and synchrotron self Compton (SSC) emission from internal shocks (IS) during the prompt and X-ray flare phases of Gamma-Ray Bursts (GRBs). The aim is to test the IS model for the flare emission and for whether GRBs can be GeV sources. We determine the parameters for which the IS model can account for the observed prompt and X-ray flares emission, and study the detectability of the high energy SSC emission by the AGILE and GLAST satellites. We find that the detectability of the SSC emission during the prompt phase of GRBs improves for higher values of the fireball Lorentz factor and of the temporal variability. If IS is the mechanism responsible of the flare emission, and the Lorentz factor of the shells producing the flare is of the order of 100, the flare light curves are expected to present some substructures with temporal variability of 10-100 ms which are much smaller than the average duration of flares, and similar to those observed during the prompt phase of GRBs. If one assumes lower Lorentz factors, such as 10-25, then a larger temporal variaibility of 40 s can also account for the observed flare properties. However in this case we predict that X-ray flares do not have a counterpart at very high energies (MeV-GeV). An investigation on the substructures of the X-ray flares light curves, and simultaneous X-ray and high energy observations, will allow us to corroborate the hypothesis that late IS are responsible of the X-ray flares.

Gamma-ray burst are thought to be produced by highly relativistic outflows. Although upper and lower limits for the outflow initial Lorentz factor $\Gamma_0$ are available, observational efforts to derive a direct determination of $\Gamma_0$ have so far failed or provided ambiguous results. As a matter of fact, the shape of the early-time afterglow light curve is strongly sensitive on $\Gamma_0$ which determines the time of the afterglow peak, i.e. when the outflow and the shocked circumburst material share a comparable amount of energy. We now comment early-time observations of the near-infrared afterglows of GRB 060418 and GRB 060607A performed by the REM robotic telescope. For both events, the afterglow peak was singled out and allowed us to determine the initial fireball Lorentz, $\Gamma_0\sim 400$.

Our knowledge of the intrinsic properties of short duration Gamma-Ray Bursts has relied, so far, only upon a few cases for which the estimate of the distance and an extended, multiwavelength monitoring of the afterglow have been obtained. We carried out multiwavelength observations of the short GRB 061201 aimed at estimating its distance and studying its properties. We performed a spectral and timing analysis of the prompt and afterglow emission and discuss the results in the context of the standard fireball model. A clear temporal break was observed in the X-ray light curve about 40 minutes after the burst trigger. We find that the spectral and timing behaviour of the X-ray afterglow is consistent with a jet origin of the observed break, although the optical data can not definitively confirm this and other scenarios are possible. No underlying host galaxy down to R~26 mag was found after fading of the optical afterglow. Thus, no secure redshift could be measured for this burst. The nearest galaxy is at z=0.111 and shows evidence of star formation activity. We discuss the association of GRB 061201 with this galaxy and with the ACO S 995 galaxy cluster, from which the source is at an angular distance of 17'' and 8.5', respectively. We also test the association with a possible undetected, positionally consistent galaxy at z~1. In all these cases, in the jet interpretation, we find a jet opening angle of 1-2 degrees.

We present observations of GRB 060124, the first event for which both the prompt and the afterglow emission could be observed simultaneously and in their entirety by the three Swift instruments and by Konus-Wind. Thanks to these exceptional circumstances, the temporal and spectral properties of the prompt emission could be studied in the optical, X-ray and gamma-ray ranges (up to 2 MeV). While the X-ray emission (0.2-10 keV) clearly tracks the gamma-ray burst, the optical component follows a different pattern, likely indicating a different origin, possibly the onset of external shocks. The prompt GRB spectrum shows significant spectral evolution, with both the peak energy and the spectral index varying. As observed in several long GRBs, significant lags are measured between the hard- and low-energy components, showing that this behaviour extends over 3 decades in energy. The GRB peaks are also much broader at soft energies. This is related to the temporal evolution of the spectrum, and can be accounted for by the softening of the electron spectral index with time. The burst energy (E_iso~5x10^{53} erg at z=2.297) and average peak energy (E_p~300 keV) make GRB 060124 consistent with the Amati relation. The X-ray afterglow is characterized by a decay which presents a break at t_b~10^5s.

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.

Gamma-ray burst (GRB) emission is believed to originate in highly relativistic fireballs. Currently, only lower limits were securely set to the initial fireball Lorentz factor Gamma_0. We aim to provide a direct measure of Gamma_0. The early-time afterglow light curve carries information about Gamma_0, which determines the time of the afterglow peak. We have obtained early observations of the near-infrared afterglows of GRB 060418 and GRB 060607A with the REM robotic telescope. For both events, the afterglow peak could be clearly singled out, allowing a firm determination of the fireball Lorentz of Gamma_0 ~ 400, fully confirming the highly relativistic nature of GRB fireballs. The deceleration radius was inferred to be R_dec ~ 10^17 cm. This is much larger than the internal shocks radius (believed to power the prompt emission), thus providing further evidence for a different origin of the prompt and afterglow stages of the GRB.

Long Gamma Ray Bursts hold the promise of probing star-formation and metal enrichment up to very high redshifts. The present GRB samples with redshift determination are largely incomplete and therefore a careful analysis of selection effects is mandatory before any conclusion can be drawn from the observed GRB redshift distribution. We study and compare three well defined samples of long GRBs detected by Swift, HETE2 and BeppoSAX. We find that Swift GRBs are slighly fainter and harder than BeppoSAX and HETE2 GRBs, as expected due to the higher energy range in which Swift GRBs are detected and localized, compared to BeppoSAX and HETE2. Gas and dust obscuration plays a role in shaping the GRB samples and the present samples of GRBs with redshift. We argue that the majority of the bright Swift GRBs without redshift might actually be z<~2 events therefore the present Swift GRB sample with redshift is biased against low-z GRBs. On the other hand, the detection of bright UV rest-frame afterglows from high-z GRBs, and even from those with large X-ray obscuration, implies a dust amount lower than in nearby GRBs,and/or a different dust composition. If this is the case, the Swift sample of GRBs with redshifts is probably a fair sample of the real high-z GRB population. The absence of high-z GRBs in the BeppoSAX and HETE2 samples of GRBs with redshifts is probably due to the fact at the time of BeppoSAX and HETE2 follow-up faint afterglows of high redshift GRBs will have weaken below the spectroscopic capabilities of even 10m class telescopes. The redshift distribution of a subsample of Swift GRBs with distributions of peak-fluxes, X-ray obscuration and optical magnitude at a fixed observing time similar to those of the BeppoSAX and HETE2 samples, is roughly consistent with BeppoSAX+HETE2 redshift distribution.

We present the results of a detailed spectral and temporal analysis of Swift and XMM-Newton observations of the high redshift (z=3.969) GRB 050730. The X-ray afterglow of GRB 050730 was found to decline with time with superimposed intense flaring activity that extended over more than two orders of magnitude in time. Seven distinct re-brightening events starting from 236 s up to 41.2 ks after the burst were observed. The underlying decay of the afterglow was well described by a double broken power-law model with breaks at t_1= 237 +/- 20 s and t_2 = 10.1 (-2.2) (+4.6) ks. The temporal decay slopes before, between and after these breaks were alpha_1 = 2.1 +/- 0.3, alpha_2 = 0.44 (-0.08) (+0.14) and alpha_3 = 2.40 (+0.07) (-0.09), respectively. The spectrum of the X-ray afterglow was well described by a photoelectrically absorbed power-law with an absorbing column density N_H=(1.28 +/- 0.26) 10^22 cm^-2 in the host galaxy. Strong X-ray spectral evolution during the flaring activity was present. In the majority of the flares (6/7) the ratio Delta_t/t_p between the duration of the event and the time when the flare peaks was nearly constant and about 0.6-0.7. We showed that the observed spectral and temporal properties of the first three flares are consistent with being due both to high-latitude emission, as expected if the flares were produced by late internal shocks, or to refreshed shocks, i.e. late time energy injections into the main afterglow shock by slow moving shells ejected from the central engine during the prompt phase. The event fully satisfies the E_p-E_iso Amati relation while is not consistent with the E_p-E_jet Ghirlanda relation.

We analyze high resolution spectroscopic observations of the optical afterglow of GRB050730, obtained with UVES@VLT about hours after the GRB trigger. The spectrum shows that the ISM of the GRB host galaxy at z = 3.967 is complex, with at least five components contributing to the main absorption system. We detect strong CII*, SiII*, OI* and FeII* fine structure absorption lines associated to the second and third component. For the first three components we derive information on the relative distance from the site of the GRB explosion. Component 1, which has the highest redshift, does not present any fine structure nor low ionization lines; it only shows very high ionization features, such as CIV and OVI, suggesting that this component is very close to the GRB site. From the analysis of low and high ionization lines and fine structure lines, we find evidences that the distance of component 2 from the site of the GRB explosion is 10-100 times smaller than that of component 3. We evaluated the mean metallicity of the z=3.967 system obtaining values about 0.01 of the solar metallicity or less. However, this should not be taken as representative of the circumburst medium, since the main contribution to the hydrogen column density comes from the outer regions of the galaxy while that of the other elements presumably comes from the ISM closer to the GRB site. Furthermore, difficulties in evaluating dust depletion correction can modify significantly these values. The mean [C/Fe] ratio agrees well with that expected by single star-formation event models. Interestingly the [C/Fe] of component 2 is smaller than that of component 3, in agreement with GRB dust destruction scenarios, if component 2 is closer than component 3 to the GRB site.

We present a detailed analysis of the selection effects that plague GRB observations. We find that these effects may partially explain the different redshift distributions between BeppoSax/HETE2 and Swift bursts. It is mandatory to consider these effects to determine the redshift evolution of GRBs.

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.

We present a detailed study of the spectral and temporal properties of the X-ray flares emission of several GRBs. We select a sample of GRBs which X-ray light curve exhibits large amplitude variations with several rebrightenings superposed on the underlying three-segment broken powerlaw that is often seen in Swift GRBs. We try to understand the origin of these fluctuations giving some diagnostic in order to discriminate between refreshed shocks and late internal shocks. For some bursts our time-resolved spectral analysis supports the interpretation of a long-lived central engine, with rebrightenings consistent with energy injection in refreshed shocks as slower shells generated in the central engine prompt phase catch up with the afterglow shock at later times.

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.

We present a detailed study of the spectral and temporal properties of the X-ray and optical emission of GRB050713a up to 0.5 day after the main GRB event. The X-ray light curve exhibits large amplitude variations with several rebrightenings superposed on the underlying three-segment broken powerlaw that is often seen in Swift GRBs. Our time-resolved spectral analysis supports the interpretation of a long-lived central engine, with rebrightenings consistent with energy injection in refreshed shocks as slower shells generated in the central engine prompt phase catch up with the afterglow shock at later times. Our sparsely-sampled light curve of the optical afterglow can be fitted with a single power law without large flares. The optical decay index appears flatter than the X-ray one, especially at later times.

Context. GRB050721 was detected by Swift and promptly followed-up, in the X-ray by Swift itself and, in the optical band, by the VLT operated, for the first time, in rapid response mode. A multiwavelength monitoring campaign was performed in order to study its afterglow behavior. Aims. We present the analysis of the early and late afterglow emission in both the X-ray and optical bands, as observed by Swift, a robotic telescope, and the VLT. We compare early observations with late afterglow observations obtained with Swift and the VLT in different bands in order to constrain the density of the medium in which the fireball is expanding. Methods. We have analyzed both the X-ray and the optical light curves and compared the spectral energy distribution of the afterglow at two different epochs. Results. We observed an intense rebrightening in the optical band at about one day after the burst which was not seen in the X-ray band. This is the first observation of a GRB afterglow in which a rebrightening is observed in the optical but not in the X-ray band. The lack of detection in X-ray of such a strong rebrightening at lower energies can be described with a variable external density profile. In such a scenario, the combined X-ray and optical observations allow us to derive the matter density at 10^15 cm from the burst. This is about a factor of 10 higher than in the inner region.

We present observations of GRB 060124, the first event for which both the prompt and the afterglow emission could be observed simultaneously and in their entirety by the three Swift instruments. Indeed, Swift-BAT triggered on a precursor ~570s before the main burst peak, and this allowed Swift to repoint the narrow field instruments to the burst position ~350s before the main burst occurred. GRB 060124 also triggered Konus-Wind, which observed the prompt emission in a harder gamma-ray band (up to 2MeV). Thanks to these exceptional circumstances, the temporal and spectral properties of the prompt emission can be studied in the optical, X-ray and gamma-ray ranges. While the X-ray emission (0.2-10keV) clearly tracks the gamma-ray burst, the optical component follows a different pattern, likely indicating a different origin, possibly the onset of external shocks. The prompt GRB spectrum shows significant spectral evolution, with both the peak energy and the spectral index varying. As observed in several long GRBs, significant lags are measured between the hard- and low-energy components, showing that this behaviour extends over 3 decades in energy. The GRB peaks are also much broader at soft energies. This is related to the temporal evolution of the spectrum, and can be accounted for by assuming that the electron spectral index softened with time. The burst energy (E_iso~5x10^{53} erg) and average peak energy (E_p~300keV) make GRB 060124 consistent with the Amati relation. The X-ray afterglow is characterized by a decay which presents a break at t_b~10^5s.

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.

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.

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.

We present optical observations of the short/hard gamma-ray burst GRB 050709, the first such event with an identified optical counterpart. The object is coincident with a weak X-ray source and is located inside a galaxy at redshift z = 0.1606 +- 0.0002. Multiband photometry allowed us to study the broad-band spectral energy distribution. Late-time monitoring places strong limits on any supernova simultaneous with the GRB. The host galaxy is not of early type. Spectra show that the dominant stellar population is relatively young (~ 1 Gyr), and that ongoing star formation is present at a level of 2-3 (L/L*) Msun/yr. This is at least 2 orders of magnitude larger than that observed in the elliptical hosts of the short GRB 050509B and GRB 050724. This shows that at least some short GRBs originate in a young population. Short/hard GRB models based on the merger of a binary degenerate system are compatible with the host galaxy characteristics, although there is still the possibility of a connection between young stars and at least a fraction of such events.

Using the BATSE peak flux distribution we rederive the short GRBs luminosity function and compare it with the observed redshift distribution of long bursts. We show that both distributions are compatible with the assumption that short as well as long bursts follow the star formation rate. In this case the difference between the two observed distributions can be interpreted as arising mostly from differences in the detector's sensitivity to long and short bursts, while the local rate of short bursts is 0.1 h_70^3 Gpc^-3 yr^-1 . We also consider the possibility that short GRBs may be associated with binary neutron star mergers and estimate the effect of the merging time delay on the luminosity function and redshift distribution. We find that in this case the local rate of short GRBs is ~0.8 h_70^3 Gpc^-3 yr^-1. Assuming that all binary merging systems lead to short GRBs, we find a typical jet opening angle of 1.6^o

A general formalism is developed for calculating the luminosity function and the expected number $N$ of observed GRBs above a peak photon flux $S$ for any GRB jet structure. This new formalism directly provides the true GRB rate without the need for a `correction factor'. We apply it to the uniform jet (UJ) and universal structured jet (USJ) models for the structure of GRB jets and perform fits to the observed log(N)-log(S) distribution from the GUSBAD catalog which contains 2204 BATSE bursts. A core angle $\theta_c$ and an outer edge at $\theta_{max}$ are introduced for the structured jet, and a finite range of half-opening angles $\theta_{min}\leq\theta_j\leq\theta_{max}$ is assumed for the uniform jets. The efficiency $\epsilon_\gamma$ for producing gamma-rays, and the energy per solid angle $\epsilon$ in the jet are allowed to vary with $\theta_j$ (the viewing angle $\theta_{obs}$) in the UJ (USJ) model, $\epsilon_\gamma\propto\theta^{-b}$ and $\epsilon\propto\theta^{-a}$. We find that a single power-law luminosity function provides a good fit to the data. Such a luminosity function arises naturally in the USJ model, while in the UJ model it implies a power-law probability distribution for $\theta_j$, $P(\theta_j)\propto\theta_j^{-q}$. The value of $q$ cannot be directly determined from the fit to the observed log(N)-log(S) distribution, and an additional assumption on the value of $a$ or $b$ is required. Alternatively, an independent estimate of the true GRB rate would enable one to determine $a$, $b$ and $q$. The implied values of $\theta_c$ (or $\theta_{min}$) and $\theta_{max}$ are close to the current observational limits. The true GRB rate for the USJ model is found to be $R_{GRB}(z=0)=0.86^{+0.14}_{-0.05} Gpc^{-3} yr^{-1}$.

We study the probability that three GRBs (980425, 030329, 031203) are found within z=0.17, given the luminosity functions consistent with the log N-log S relationship for classical cosmological bursts (i.e., those observed by BATSE). We show that, in order for the probability of these three low-z events to be non-negligible (thus making it more likely that they belong to the same class of the classical cosmological bursts), the bursts' luminosity function must be a broken powerlaw. By reasoning in analogy with beamed AGNs, we show that observations are consistent with the expectations if GRB 980425 and GRB 031203 are indeed normal bursts seen sideways. Within this model, no bright burst within z=0.17 should be observed by a HETE--like instrument within the next \sim 20 yr.

The realization that the total energy of GRBs is correlated with their jet break angles motivates the search for a similar relation between the peak luminosity, $L$, and the jet break angles, $L\propto\theta^{-2}$. Such a relation implies that the GRB luminosity function determines the angular distribution. We re-derive the GRB luminosity function using the BATSE peak flux distribution and compare the predicted distribution with the observed redshift distribution. The luminosity function can be approximated by a broken power law with a break peak luminosity of $4.4 \times 10^{51}$ erg/sec, a typical jet angle of 0.12 rad and a local GRB rate of $0.44 h_{65}^{3}$ Gpc$^{-3}$yr$^{-1}$. The angular distribution implied by $L\propto\theta^{-2}$ agrees well with the observed one, and implies a correction factor to the local rate due to beaming of $75\pm25$ (instead of 500 as commonly used). The inferred overall local GRB rate is $33\pm11 h_{65}^{3}$ Gpc$^{-3}$yr$^{-1}$. The luminosity function and angle distribution obtained within the universal structured jet model, where the angular distribution is essentially $\propto \theta$ and hence the luminosity function must be $\propto L^{-2}$, deviate from the observations at low peak fluxes and, correspondingly, at large angles. The corresponding correction factor for the universal structure jet is $\sim 20 \pm 10$.

The two leading models for the structure of GRB jets are the uniform jet model and the universal structured jet (USJ) model. In the latter, all GRB jets are intrinsically identical and the energy per solid angle drops as the inverse square of the angle from the jet axis. The simplicity of the USJ model gives it a strong predictive power, including a specific prediction for the observed GRB distribution as a function of both the redshift $z$ and the viewing angle $\theta$. We show that the current sample of GRBs with known $z$ and estimated $\theta$ does not agree with the predictions of the USJ model. This can be best seen for a relatively narrow range in $z$, in which the USJ model predicts that most GRBs should be near the upper end of the observed range in $\theta$, while in the observed sample most GRBs are near the lower end of that range. Since the current sample is very inhomogeneous (i.e. involves many different detectors), it should be taken with care and cannot be used to rule out the USJ model. Nevertheless, this sample strongly disfavors the USJ model. Comparing the prediction for the observed GRB distribution both in $\theta$ and in $z$, with a larger and more homogeneous GRB sample, like the one expected from Swift, would either clearly rule out the USJ model, or alternatively, provide a strong support for it. The test presented here is general, and can be used to test any model that predicts both a luminosity function and a luminosity-angle relation.

We consider internal shocks as the main dissipation mechanism responsible for the emission in blazars and show that it can satisfactorily account for the properties of all blazars. In particular, we extend previous work (Spada et al. 2001) on powerful objects, to intermediate (BL Lac) and low power sources (Mkn 421), in order to reproduce the whole of the blazar sequence. The model self-consistently treats the dynamics, spectral emission and its variability. The key parameters driving the phenomenological sequence are the jet power and the properties of the broad line region, which regulate the cooling efficiency of the emitting particles and thus the shape of the spectral energy distribution. By assuming that the remaining parameters are similar for all objects it has been possible to reproduce the full range of the observed spectral "states". A more detailed comparison of the variability properties shows (for Mkn 421) a good agreement in the X-ray band, while in the optical the simulated flux appears to be too variable. For BL Lac lags (~10 days) are predicted between the gamma--rays and the infrared emission

The gamma-ray burst (GRB) of October 17, 1994 (941017), showed a distinct high energy spectral component extending from $\lesssim$ a few MeV to $\gtrsim$ 200 MeV, in addition to the typical GRB emission which peaked at $\lesssim$ a few hundred keV. The high energy component carried at least $\sim$ 3 times more energy than the lower energy component. It displayed an almost constant flux with a rather hard spectrum ($F_\nu \propto \nu^{-\alpha}$ with $\alpha \sim 0$) from $\lesssim$ 20 s into the burst up to $\sim$ 200 s, while the duration of the GRB, where 90% of the energy in the lower energy component was emitted, was only 77 s. Such a high energy component was seen in only one out of $\sim$ 30 GRBs in which a similar component could have been detected, and thus appears to be quite rare. We examine possible explanations for this high energy spectral component and find that most models fail. The only emission region that provides the right temporal behavior is the reverse shock that goes into the GRB ejecta as it is decelerated by the ambient medium, or possibly the very early forward shock while the reverse shock is still going on. The best candidate for the emission mechanism is synchrotron self-Compton emission from the reverse shock. Even in this model the most natural spectral slope is only marginally consistent with the observed value, and some degree of fine tuning is required in order to improve the agreement. This might suggest that an additional or alternative emission mechanism is at work here. A prediction of this interpretation is that such a high energy component should be accompanied by a bright optical transient, similar to the one observed in GRB 990123.

We calculate the neutrino emission from individual gamma-ray bursts observed by the BATSE detector on the Compton Gamma-Ray Observatory. Neutrinos are produced by photoproduction of pions when protons interact with photons in the region where the kinetic energy of the relativistic fireball is dissipated allowing the acceleration of electrons and protons. We also consider models where neutrinos are predominantly produced on the radiation surrounding the newly formed black hole. From the observed redshift and photon flux of each individual burst, we compute the neutrino flux in a variety of models based on the assumption that equal kinetic energy is dissipated into electrons and protons. Where not measured, the redshift is estimated by other methods. Unlike previous calculations of the universal diffuse neutrino flux produced by all gamma-ray bursts, the individual fluxes (compiled at http://www.arcetri.astro.it/~dafne/grb/) can be directly compared with coincident observations by the AMANDA telescope at the South Pole. Because of its large statistics, our predictions are likely to be representative for future observations with larger neutrino telescopes.

The supranova model for Gamma-Ray Bursts (GRBs) is becoming increasingly more popular. In this scenario the GRB occurs weeks to years after a supernova explosion, and is located inside a pulsar wind bubble (PWB). Protons accelerated in the internal shocks that emit the GRB may interact with the external PWB photons producing pions which decay into \sim 10^{16} eV neutrinos. A km^2 neutrino detector would observe several events per year correlated with the GRBs.

The supranova model for gamma-ray bursts (GRBs) has recently gained popularity. In this scenario the GRB occurs weeks to years after a supernova explosion, and is located inside a pulsar wind bubble (PWB). High energy protons from the PWB can interact with photons from the rich radiation field inside the PWB or collide with cold protons from the supernova remnant, producing pions which decay into \sim 10-10^3 TeV neutrinos. The predicted neutrino flux from the PWBs that host the GRBs should be easily detectable by planned 1 km^2 detectors.

We re-examine critically the arguments raised against the theory that Ultra High Energy Cosmic Rays observed at Earth are produced in Gamma Ray Bursts. These include the limitations to the highest energy attainable by protons around the bursts' shocks, the spectral slope at the highest energies, the total energy released in non--thermal particles, the occurrence of doublets and triplets in the data reported by AGASA. We show that, to within the uncertainties in our current knowledge of GRBs, none of these objections is really fatal to the scenario. In particular, we show that the total energy budget of GRBs easily accounts for the energy injection rate necessary to account for UHECRs as observed at Earth. We also compute the expected particle spectrum at Earth, showing that it fits the HiRes and AGASA data to within statistical uncertainties. We consider the existence of multiplets in AGASA' data. To this end, we present a Langevin--like treatment for the motion of a charged particle in the IGM magnetic field, which allows us to estimate both the average and the rms timedelay for particles of given energy; we discuss when particles of identical energies reach the Earth in bunches, or spread over the rms timedelay, showing that multiplets pose no problem for an explosive model for the sources of UHECRs. We compare our model with a scenario where the particles are accelerated at internal shocks, underlining differences and advantages of particle acceleration at external shocks.

We present the main observational features expected for Gamma-Ray Bursts (GRBs) that occur inside pulsar wind bubbles (PWBs). This is the most natural outcome of the supranova model, where initially a supernova (SN) explosion takes place, leaving behind a supra-massive neutron star, which loses its rotational energy over a time t_sd and collapses to a black hole, triggering a GRB. We find that the time delay, t_sd, between the SN and GRB events is the most important parameter that determines the behavior of the system. We consider the afterglow, prompt GRB and PWB emission. Constraints on the model are derived for a spherical PWB, from current afterglow observations and the lack of direct detection of the PWB emission. We find that a simple spherical model cannot account for the X-ray features detected in some afterglows, together with the typical afterglow emission that was observed for the same GRBs. The discrepancies with observations may be reconciled by resorting to a non-spherical geometry, where the PWB is elongated along the polar axis. Finally, we predict that the inverse Compton upscattering of PWB photons by the relativistic electrons of the afterglow (external Compton) should lead to high energy emission during the early afterglow that may explain the GeV photons detected by EGRET for a few GRBs, and should be detectable by future missions such as GLAST.

The resonant cyclotron absorption model is very successful in describing particle acceleration in plerions. A sensible prediction of this model is the presence of a substantial amount of relativistic protons in pulsar winds. Although difficult to detect, these protons may show up through their interactions either with the photons in the plerion environment or with the thermal gas in the supernova ejecta. Inelastic proton-proton (p-p) collisions are expected to be very effective in young objects, resulting in a copious production of neutral and charged pions. Charged pions produced during the first few hundred years after the supernova explosion may have time to decay into muons, whose subsequent decay may provide an additional source of electrons and positrons in these nebulae, that sums up to the pulsar input. These secondary leptons evolve just as the pairs in the pulsar wind, and signatures of their presence could be found, in principle, even in the synchrotron spectrum of older objects. p-p collisions may remain fairly efficient even in moderately old objects resulting in the production of TeV gamma-rays and neutrinos. We apply our calculations to the case of the Crab Nebula, the best studied plerion insofar, and find that existing data already allow to infer interesting constraints on the physical properties of the Crab pulsar wind.

We consider GRBs that occur inside pulsar wind bubbles (PWBs), in the context of the supranova model, where initially a supernova explosion takes place, leaving behind a supra-massive neutron star (SMNS), which eventually loses its rotational energy and collapses to a black hole, triggering a GRB explosion. The most natural mechanism by which the SMNS can lose its rotational energy is through a strong pulsar type wind, between the supernova and the GRB events, which is expected to create a PWB. We analyze in some detail the observational implications of such a plerionic environment on the afterglow (AG) and prompt GRB emissions, and the prospect for direct detection of the plerion emission. We find that for a simple spherical model, GRBs with iron lines detected in the X-ray AG should not have a detectable radio afterglow, and should have small jet break times and non-relativistic transition times, in disagreement with observations for some of the GRBs with X-ray lines. These discrepancies with the observations may be reconciled by resorting to a non-spherical geometry, where the PWB is elongated along the polar axis. We find that lack of detection of the plerion emission provides interesting constraints on the model parameters. Finally, we predict that the inverse Compton upscattering of the PWB photons by the relativistic electrons of the AG (external Compton, EC) should lead to high energy emission during the early AG that might explain the GeV photons detected by EGRET for a few GRBs, and should be detectable by future missions such as GLAST.

The supranova model of gamma-ray bursts (GRBs), in which the GRB event is preceded by a supernova (SN) explosion by a few months to years, has recently gained support from Fe line detections in X-ray afterglows. A crucial ingredient of this model yet to be studied is the fast-rotating pulsar that should be active during the time interval between the SN and the GRB, driving a powerful wind and a luminous plerionic nebula. We discuss some observational consequences of this precursor plerion, which should provide important tests for the supranova model: 1) the fragmentation of the outlying SN ejecta material by the plerion and its implications for Fe line emission; and 2) the effect of inverse Compton cooling and emission in the GRB external shock due to the plerion radiation field. The plerion-induced inverse Compton emission can dominate in the GeV-TeV energy range during the afterglow, being detectable by GLAST from redshifts $z \lesssim 1.5$ and distinguishable from self-Compton emission by its spectrum and light curve. The prospects for direct detection and identification of the precursor plerion emission are also briefly considered.

We present the main observational features expected for GRBs that occur inside pulsar wind bubbles (PWBs). This is the most natural outcome of supranova model where initially a supernova explosion takes place, leaving a supra-massive NS, which loses its rotational energy over a time t_sd and collapses to a BH, triggering a GRB explosion. We find that the time delay t_sd between the supernova and GRB events is the most important parameter that determines the behavior of the system. We consider the afterglow, prompt GRB and direct emission from the PWB. The observational signatures for different t_sd are described and joined together into one coherent framework. Constraints on the model are derived for a spherical PWB, from the lack of detection of the PWB emission and afterglow observations. For very low values of t_sd<1 hr the supranova model reduces to the collapsar model; 0.4<t_sd/1yr<1 is required to produce the iron lines seen in some X-ray afterglows. However, we find that for a spherical PWB, this implies no detectable radio afterglow, a small jet break time and non-relativistic transition time, in disagreement with observations for some of these GRBs. These discrepancies with the observations may be reconciled by resorting to a non-spherical geometry. We find that light element lines, recently detected in a few X-ray afterglows, are expected to dominate over iron lines for small t_sd, while for large t_sd the situation is reversed. Finally, we predict that external Compton scattering of the PWB photons by the afterglow relativistic electrons, should lead to high energy emission during the early afterglow that may explain the GeV photons detected by EGRET in some GRBs, and should be detectable by future missions (GLAST).

We consider the internal shock model as the dissipation mechanism responsible for the emission in blazars. It has been shown that this model is successful in reproducing the observed spectral energy distribution and the variability properties of a powerful blazar like 3C 279. However, the blazar family covers a wide range of spectral characteristics which appear to be correlated and the whole class can be seen as a sequence: the frequency and the intensity of the low energy versus high energy peak intensity increase with decreasing luminosity. We show that the internal shock model can satisfactorily account also for the properties of the low power blazars like BL Lac and Mkn 421 and it is successful in reproducing the blazar sequence.

We study the synchrotron and synchrotron self-Compton (SSC) emission from internal shocks that are responsible for the prompt gamma-ray emission in GRBs, and consider the relation between these two components, taking into account the high energy (HE) cutoff due to pair production and Thomson scattering. We find that in order for the peak energy of the synchrotron to be E_p\sim 300 keV with a variability time t_v>1 ms, a Lorentz factor \Gamma<350 is needed, implying no HE emission above \sim 30 MeV and the synchrotron component would dominate at all energies. If we want both E_p\sim 300 keV and prompt HE emission up to 2 GeV, as detected by EGRET for GRB 940217, we need \Gamma\sim 600 and t_v\sim 0.1 ms, which might be resolved by super AGILE. If such prompt HE emission is common in GRBs, as may be tested by GLAST, then for t_v\gtrsim 1 ms we need \Gamma\gtrsim 350, which implies E_p\lesssim 100 keV. Therefore if X-ray flashes are GRBs with high values of t_v and \Gamma, they should produce \gtrsim 1 GeV emission. For an electron power law index p>2, the SSC component dominates the emission above 100 MeV. Future observations by GLAST may help determine the value of p and whether the HE emission is consistent with a single power law (one component--the synchrotron, dominates) or has a break where the \nuF_\nu slope turns from negative to positive, implying that the SSC component becomes dominant above \sim 100 MeV. The HE emission is expected to show similar variability and time structure to that of the soft gamma-ray emission. Finally, we find that in order to see delayed HE emission from the prompt GRB due to pair production with the cosmic IR background, extremely small intergalactic magnetic fields (\lessim 10^{-22} G) are required.

We consider a class of plerions that have been detected in the TeV range, and investigate the possibility that the emission is due to \pi^0 decay. From the TeV flux we derive what is the expected \nu flux at Earth and find that this is likely to exceed the detection threshold of the upcoming km^2 neutrino detectors.

We discuss the possibility that microquasar jets may be powerful emitters of TeV neutrinos. We estimate the neutrino fluxes produced by photopion production in the jets of a sample of identified microquasars and microquasar candidates, for which available data enables rough determination of the jet parameters. We demonstrate that in several of the sources considered, the neutrino flux at Earth, produced in events similar to those observed, can exceed the detection threshold of a km^2 neutrino detector. Sources with bulk Lorentz factors larger than those characteristic of the sample considered here, directed along our line of sight may be very difficult to resolve at radio wavelengths and hence may be difficult to identify as microqusar candidates. However these sources can be identified through their neutrino and gamma-ray emission.

It has been proposed recently that Galactic microquasars may be prodigious emitters of TeV neutrinos that can be detected by upcoming km^2 neutrino telescopes. In this paper we consider a sample of identified microquasars and microquasar candiates, for which available data enables rough determination of the jet parameters. By employing the parameters inferred from radio observations of various jet ejection events, we determine the neutrino fluxes that should have been produced during these events by photopion production in the jet. Despite the large uncertainties in our analysis, we demonstrate that in several of the sources considered, the neutrino flux at Earth, produced in events similar to those observed, would exceed the detection threshold of a km^2 neutrino detector. The class of microquasars may contain also sources with bulk Lorentz factors larger than those characteristic of the sample considered here, directed along our line of sight. Such sources, which may be very difficult to resolve at radio wavelengths and hence may be difficult to identify as microqusar candidates, may emit neutrinos with fluxes significantly larger than typically obtained in the present analysis. These sources may eventually be identified through their neutrino and gamma-ray emission.

If the first (PopIII) stars were very massive, their final fate is to collapse into very massive black holes. Once a proto-black hole has formed into the stellar core, accretion continues through a disk. It is widely accepted, although not confirmed, that magnetic fields drive an energetic jet which produces a burst of TeV neutrinos by photon-meson interaction, and eventually breaks out of the stellar envelope appearing as a Gamma Ray Burst (GRB). Based on recent numerical simulations and neutrino emission models, we predict the expected neutrino diffuse flux from these PopIII GRBs and compare it with the capabilities of present and planned detectors as AMANDA and IceCube. If beamed into 1% of the sky, we find that the rate of PopIII GRBs is $\le 4 \times 10^6$ yr$^{-1}$. High energy neutrinos from PopIII GRBs could dominate the overall flux in two energy bands [$10^4 - 10^5$] GeV and [$10^5 - 10^6$] GeV of neutrino telescopes. The enhanced sensitivities of forthcoming detectors in the high-energy band (AMANDA-II, IceCube) will provide a fundamental insight on the characteristic explosion energies of PopIII GRBs and will constitute a unique probe of the the Initial Mass Function (IMF) of the first stars and of the redshift $z_f$ marking the metallicity-driven transition from a top-heavy to a normal IMF. The current upper limit set by AMANDA-B10 implies that such transition must have occurred not later than $z_f =9.2$ for the most plausible jet energies. Based on such results, we speculate that PopIII GRBs, if not chocked, could be associated with a new class of events detected by BeppoSax, the Fast X-ray Transient (FXTs), which are bright X-ray sources, with peak energies in the 2-10 keV band and durations between 10-200 s.

Observations imply that gamma-ray bursts (GRBs) are produced by the dissipation of the kinetic energy of a highly relativistic fireball. Photo-meson interactions of protons with gamma-rays within the fireball dissipation region are expected to convert a significant fraction of the fireball energy into high energy neutrinos. In this talk we summarize the results of an analysis of the internal shock model of GRBs, where production of synchrotron photons and photo-meson neutrinos with energies about 10^14 eV are self-consistently calculated. These neutrinos will be coincident with the GRB. We show that the fraction of fireball energy converted into high energy neutrinos is not sensitive to uncertainties in the fireball model parameters, such as the expansion Lorentz factor and characteristic variability time. Other processes of neutrinos emission from GRBs are also reviewed. Photomeson interactions within the plasma region shocked by the reverse shock, may produce a burst of 10^18 eV neutrinos following the GRB on the time scale of 10 s. Inelastic p-n nuclear collisions result in the production of a burst of around 10 GeV neutrinos in coincidence with the GRB. Planned 1 km^3 neutrino telescopes are expected to detect ten 100 TeV neutrino events and several 10^18 eV events, correlated with GRBs, per year. A suitably densely spaced detector may allow the detection of several 10 GeV events per year.

We present an analysis of the Internal Shock Model of GRBs, where gamma-rays are produced by internal shocks within a relativistic wind. We show that observed GRB characteristics impose stringent constraints on wind and source parameters. We find that a significant fraction, of order 20 %, of the wind kinetic energy can be converted to radiation, provided the distribution of Lorentz factors within the wind has a large variance and provided the minimum Lorentz factor is higher than 10^(2.5)L_(52)^(2/9), where L=10^(52)L_(52)erg/s is the wind luminosity. For a high, >10 %, efficiency wind, spectral energy breaks in the 0.1 to 1 MeV range are obtained for sources with dynamical time R/c < 1 ms, suggesting a possible explanation for the observed clustering of spectral break energies in this range. The lower limit to wind Lorenz factor and the upper limit, around (R/10^7 cm)^(-5/6) MeV to observed break energies are set by Thomson optical depth due to electron positron pairs produced by synchrotron photons. Natural consequences of the model are absence of bursts with peak emission energy significantly exceeding 1 MeV, and existence of low luminosity bursts with low, 1 keV to 10 keV, break energies.

Observations imply that gamma-ray bursts (GRBs) are produced by the dissipation of the kinetic energy of a highly relativistic fireball. Photo-meson interactions of protons with gamma-rays within the fireball dissipation region are expected to convert a significant fraction of fireball energy to >10^14 eV neutrinos. We present an analysis of the internal shock model of GRBs, where production of synchrotron photons and photo-meson neutrinos are self-consistently calculated, and show that the fraction of fireball energy converted to high energy neutrinos is not sensitive to uncertainties in fireball model parameters, such as the expansion Lorentz factor and characteristic variability time. This is due in part to the constraints imposed on fireball parameters by observed GRB characteristics, and in part to the fact that for parameter values for which the photo-meson optical depth is high (implying high proton energy loss to pion production) neutrino production is suppressed by pion and muon synchrotron losses. The neutrino flux is therefore expected to be correlated mainly with the observed gamma-ray flux. The time averaged neutrino intensity predicted by the model, ~10^-8.5 GeV/cm^2 s sr, is consistent with the flux predicted by the assumption that GRBs are the sources of >10^19 eV cosmic-rays.