Gamma Ray Bursts are detectable in the gamma-ray band if their jets are oriented towards the observer. However, for each GRB with a typical theta_jet, there should be ~2/theta_jet^2 bursts whose emission cone is oriented elsewhere in space. These off
-axis bursts can be eventually detected when, due to the deceleration of their relativistic jets, the beaming angle becomes comparable to the viewing angle. Orphan Afterglows (OA) should outnumber the current population of bursts detected in the gamma-ray band even if they have not been conclusively observed so far at any frequency. We compute the expected flux of the population of orphan afterglows in the mm, optical and X-ray bands through a population synthesis code of GRBs and the standard afterglow emission model. We estimate the detection rate of OA by on-going and forthcoming surveys. The average duration of OA as transients above a given limiting flux is derived and described with analytical expressions: in general OA should appear as daily transients in optical surveys and as monthly/yearly transients in the mm/radio band. We find that ~ 2 OA yr^-1 could already be detected by Gaia and up to 20 OA yr^-1 could be observed by the ZTF survey. A larger number of 50 OA yr^-1 should be detected by LSST in the optical band. For the X-ray band, ~ 26 OA yr^-1 could be detected by the eROSITA. For the large population of OA detectable by LSST, the X-ray and optical follow up of the light curve (for the brightest cases) and/or the extensive follow up of their emission in the mm and radio band could be the key to disentangle their GRB nature from other extragalactic transients of comparable flux density.
Gamma Ray Bursts (GRBs) are a powerful probe of the high redshift Universe. We present a tool to estimate the detection rate of high-z GRBs by a generic detector with defined energy band and sensitivity. We base this on a population model that reprod
uces the observed properties of GRBs detected by Swift, Fermi and CGRO in the hard X-ray and gamma-ray bands. We provide the expected cumulative distributions of the flux and fluence of simulated GRBs in different energy bands. We show that scintillator detectors, operating at relatively high energies (e.g. tens of keV to the MeV), can detect only the most luminous GRBs at high redshifts due to the link between the peak spectral energy and the luminosity (Ep-Liso) of GRBs. We show that the best strategy for catching the largest number of high-z bursts is to go softer (e.g. in the soft X-ray band) but with a very high sensitivity. For instance, an imaging soft X-ray detector operating in the 0.2-5 keV energy band reaching a sensitivity, corresponding to a fluence of ~10^-8 erg cm^-2, is expected to detect ~40 GRBs yr^-1 sr^-1 at z>5 (~3 GRBs yr^-1 sr^-1 at z>10). Once high-z GRBs are detected the principal issue is to secure their redshift. To this aim we estimate their NIR afterglow flux at relatively early times and evaluate the effectiveness of following them up and construct usable samples of events with any forthcoming GRB mission dedicated to explore the high-z Universe.
Orphan Afterglows (OA) are slow transients produced by Gamma Ray Bursts seen off-axis that become visible on timescales of days/years at optical/NIR and radio frequencies, when the prompt emission at high energies (X and gamma rays) has already cease
d. Given the typically estimated jet opening angle of GRBs theta_jet ~ 3 deg, for each burst pointing to the Earth there should be a factor ~ 700 more GRBs pointing in other directions. Despite this, no secure OAs have been detected so far. Through a population synthesis code we study the emission properties of the population of OA at radio frequencies. OAs reach their emission peak on year-timescales and they last for a comparable amount of time. The typical peak fluxes (which depend on the observing frequency) are of few micro Jy in the radio band with only a few OA reaching the mJy level. These values are consistent with the upper limits on the radio flux of SN Ib/c observed at late times. We find that the OA radio number count distribution has a typical slope -1.7 at high fluxes and a flatter (-0.4) slope at low fluxes with a break at a frequency-dependent flux. Our predictions of the OA rates are consistent with the (upper) limits of recent radio surveys and archive searches for radio transients. Future radio surveys like VAST/ASKAP at 1.4 GHz should detect ~ 3x10^-3 OA deg^-2 yr-1, MeerKAT and EVLA at 8.4 GHz should see ~ 3x10^-1 OA deg-2 yr-1. The SKA, reaching the micro Jy flux limit, could see up to ~ 0.2-1.5 OA deg^-2 yr^-1. These rates also depend on the duration of the OA above a certain flux limit and we discuss this effect with respect to the survey cadence.
We estimate the initial bulk Lorentz factors Gamma_0 for GRBs that show the onset of the afterglow in their optical light curves. We find that Gamma_0 is strongly correlated with both the isotropic equivalent luminosity L_iso and energy E_iso and, wi
th a larger scatter, also with the rest frame peak energy E_peak. These new correlations allow us to interpret the spectral energy correlations E_peak-L_iso (-E_iso) as a sequence of Gamma_0 factors. By accounting for the beaming effects, we find that the comoving frame properties of GRBs result clustered around typical values (e.g. L_iso~5x10^48 erg/s). Moreover, it is theoretically predicted that there should be a link between the jet dynamics (Gamma_0) and its geometry (theta_jet). Through a population synthesis code we reconstruct the Gamma_0 and theta_jet distributions and search for a possible link between them. We find that Gamma_0 and theta_jet in GRBs should have log-normal distributions and they should be anti correlated (i.e. theta_jet^2*Gamma_0=const).
Radio observations of Gamma Ray Bursts afterglows are fundamental in providing insights into their physics and environment, and in constraining the true energetics of these sources. Nonetheless, radio observations of GRB afterglows are presently spar
se in the time/frequency domain. Starting from a complete sample of 58 bright Swift long bursts (BAT6), we constructed a homogeneous sub-sample of 38 radio detections/upper limits which preserves all the properties of the parent sample. One half of the bursts have detections between 1 and 5 days after the explosion with typical fluxes F>100 muJy at 8.4 GHz. Through a Population SYnthesis Code coupled with the standard afterglow Hydrodynamical Emission model (PSYCHE) we reproduce the radio flux distribution of the radio sub-sample. Based on these results we study the detectability in the time/frequency domain of the entire long GRB population by present and future radio facilities. We find that the GRBs that typically trigger Swift can be detected at 8.4 GHz by JVLA within few days with modest exposures even at high redshifts. The final SKA can potentially observe the whole GRB population provided that there will be a dedicated GRB gamma-ray detector more sensitive than Swift. For a sizable fraction (50%) of these bursts, SKA will allow us to perform radio-calorimetry, after the trans-relativistic transition (occurring ~100 d), providing an estimate of the true (collimation corrected) energetics of GRBs.
The jet opening angle theta_jet and the bulk Lorentz factor Gamma_0 are crucial parameters for the computation of the energetics of Gamma Ray Bursts (GRBs). From the ~30 GRBs with measured theta_jet or Gamma_0 it is known that: (i) the real energetic
E_gamma, obtained by correcting the isotropic equivalent energy E_iso for the collimation factor ~theta_jet^2, is clustered around 10^50-10^51 erg and it is correlated with the peak energy E_p of the prompt emission and (ii) the comoving frame E_p and E_gamma are clustered around typical values. Current estimates of Gamma_0 and theta_jet are based on incomplete data samples and their observed distributions could be subject to biases. Through a population synthesis code we investigate whether different assumed intrinsic distributions of Gamma_0 and theta_jet can reproduce a set of observational constraints. Assuming that all bursts have the same E_p and E_gamma in the comoving frame, we find that Gamma_0 and theta_jet cannot be distributed as single power-laws. The best agreement between our simulation and the available data is obtained assuming (a) log-normal distributions for theta_jet and Gamma_0 and (b) an intrinsic relation between the peak values of their distributions, i.e theta_jet^2.5*Gamma_0=const. On average, larger values of Gamma_0 (i.e. the faster bursts) correspond to smaller values of theta_jet (i.e. the narrower). We predict that ~6% of the bursts that point to us should not show any jet break in their afterglow light curve since they have sin(theta_jet)<1/Gamma_0. Finally, we estimate that the local rate of GRBs is ~0.3% of all local SNIb/c and ~4.3% of local hypernovae, i.e. SNIb/c with broad-lines.
We study the possible effects of selection biases on the Ep-Liso correlation caused by the unavoidable presence of flux-limits in the existing samples of Gamma Ray Bursts (GRBs). We consider a well defined complete sample of Swift GRBs and perform Mo
nte Carlo simulations of the GRB population under different assumptions for their luminosity functions. If we assume that there is no correlation between the peak energy Ep and the isotropic luminosity Liso, we are unable to reproduce it as due to the flux limit threshold of the Swift complete sample. We can reject the null hypothesis that there is no intrinsic correlation between Ep and Liso at more than 2.7 sigma level of confidence. This result is robust against the assumptions of our simulations and it is confirmed if we consider, instead of Swift, the trigger threshold of the Batse instrument. Therefore, there must be a physical relation between these two quantities. Our simulations seem to exclude, at a lower confidence level of 1.6 sigma, the possibility that the observed Ep-Liso correlation among different bursts is caused by a boundary, i.e. such that for any given Ep, we see only the largest Liso, which has a flux above the threshold of the current instruments.
We use a nearly complete sample of Gamma Ray Bursts (GRBs) detected by the Swift satellite to study the correlations between the spectral peak energy Ep of the prompt emission, the isotropic energetics Eiso and the isotropic luminosity Liso. This GRB
sample is characterized by a high level of completeness in redshift (90%). This allows us to probe in an unbiased way the issue related to the physical origin of these correlations against selection effects. We find that one burst, GRB 061021, is an outlier to the Ep-Eiso correlation. Despite this case, we find strong Ep-Eiso and Ep-Liso correlations for the bursts of the complete sample. Their slopes, normalisations and dispersions are consistent with those found with the whole sample of bursts with measured redshift and Ep. This means that the biases present in the total sample commonly used to study these correlations do not affect their properties. Finally, we also find no evolution with redshift of the Ep-Eiso and Ep-Liso correlations.
We study the spectral evolution on second and sub--second timescales in 11 long and 12 short Gamma Ray Bursts (GRBs) with peak flux >8.5e-6 erg/cm2 s (8 keV-35 MeV) detected by the Fermi satellite. The peak flux correlates with the time-averaged peak
energy in both classes of bursts. The peak energy evolution, as a function of time, tracks the evolution of the flux on short timescales in both short and long GRBs. We do not find evidence of an hard-to-soft spectral evolution. While short GRBs have observed peak energies larger than few MeV during most of their evolution, long GRBs can start with a softer peak energy (of few hundreds keV) and become as hard as short ones (i.e. with Ep,obs larger than few MeV) at the peak of their light curve. Six GRBs in our sample have a measured redshift. In these few cases we find that their correlations between the rest frame Ep and the luminosity Liso are less scattered than their correlations in the observer frame between the peak energy Ep,obs and the flux P. We find that the rest frame Ep of long bursts can be as high or even larger than that of short GRBs and that short and long GRBs follow the same Ep-Liso correlation, despite the fact that they likely have different progenitors.
We estimate the bulk Lorentz factor Gamma_0 of 31 GRBs using the measured peak time of their afterglow light curves. We consider two possible scenarios for the estimate of Gamma_0: the case of a homogeneous circumburst medium or a wind density profil
e. The values of Gamma_0 are broadly distributed between few tens and several hundreds with average values ~138 and ~66 for the homogeneous and wind density profile, respectively. We find that the isotropic energy and luminosity correlate in a similar way with Gamma_0, i.e. Eiso Gamma_0^2 and Liso Gamma_0^2, while the peak energy Epeak Gamma_0. These correlations are less scattered in the wind density profile than in the homogeneous case. We then study the energetics, luminosities and spectral properties of our bursts in their comoving frame. The distribution of Liso is very narrow with a dispersion of less than a decade in the wind case, clustering around Liso=5x10^48 erg/s. Peak photon energies cluster around Epeak=6 keV. The newly found correlations involving Gamma_0 offer a general interpretation scheme for the spectral-energy correlations of GRBs. The Epeak-Eiso and Epeak-Liso correlations are due to the different Gamma_0 factors and the collimation-corrected correlation, Epeak-Egamma (obtained by correcting the isotropic quantities for the jet opening angle theta_j), can be explained if theta_j^2*Gamma_0=constant. Assuming the Epeak-Egamma correlation as valid, we find a typical value of theta_j*Gamma_0 ~ 6-20, in agreement with the predictions of magnetically accelerated jet models.
