Knowledge of the bulk Lorentz factor $\Gamma_{0}$ of GRBs allows us to compute their comoving frame properties shedding light on their physics. Upon collisions with the circumburst matter, the fireball of a GRB starts to decelerate, producing a peak or a break (depending on the circumburst density profile) in the light curve of the afterglow. Considering all bursts with known redshift and with an early coverage of their emission, we find 67 GRBs with a peak in their optical or GeV light curves at a time $t_{\rm p}$. For another 106 GRBs we set an upper limit $t_{\rm p}^{\rm UL}$. We show that $t_{\rm p}$ is due to the dynamics of the fireball deceleration and not to the passage of a characteristic frequency of the synchrotron spectrum across the optical band. Considering the $t_{\rm p}$ of 66 long GRBs and the 85 most constraining upper limits, using censored data analysis methods, we reconstruct the most likely distribution of $t_{\rm p}$. All $t_{\rm p}$ are larger than the time $t_{\rm p,g}$ when the prompt emission peaks, and are much larger than the time $t_{\rm ph}$ when the fireball becomes transparent. The reconstructed distribution of $\Gamma_0$ has median value $\sim$300 (150) for a uniform (wind) circumburst density profile. In the comoving frame, long GRBs have typical isotropic energy, luminosity, and peak energy $\langle E_{\rm iso}\rangle=3(8)\times 10^{50}$ erg, $\langle L_{\rm iso}\rangle=3(15) \times 10^{47}$ erg s$^{-1}$ , and $\langle E_{\rm peak}\rangle =1(2)$ keV in the homogeneous (wind) case. We confirm that the significant correlations between $\Gamma$ and the rest frame isotropic energy ($E_{\rm iso}$), luminosity ($L_{\rm iso}$) and peak energy ($E_{\rm peak}$) are not due to selection effects. Assuming a typical opening angle of 5 degrees, we derive the distribution of the jet baryon loading which is centered around a few $10^{-6} {\rm M_{\odot}}$.

We present a comprehensive temporal and spectral analysis of the long Swift GRB 120327A afterglow data to investigate the possible causes of the observed early time colour variations. We collected data from various instruments/telescopes in different bands (X-rays, ultra- violet, optical and near-infrared) and determined the shapes of the afterglow early-time light curves. We studied the overall temporal behaviour and the spectral energy distributions from early to late times. The ultra-violet, optical, and near-infrared light curves can be modelled with a single power-law component between 200 and 2e4 s after the burst event. The X-ray light curve shows a canonical steep-shallow-steep behaviour, typical of long gamma-ray bursts. At early times a colour variation is observed in the ultra-violet/optical bands, while at very late times a hint of a re-brightening is visible. The observed early time colour change can be explained as a variation in the intrinsic optical spectral index, rather than an evolution of the optical extinction.

We present a multiwavelength study of GRB 151027A. This is the 999th GRB detected by the Swift satellite and it has a densely sampled emission in the X-ray and optical band and has been observed and detected in the radio up to 140 days after the prompt. The multiwavelength light curve from 500 s to 140 days can be modelled through a standard forward shock afterglow but requires an additional component to reproduce the early X-ray and optical emission. We present TNG and LBT optical observations performed 19.6, 33.9 and 92.3 days after the trigger which show a bump with respect to a standard afterglow flux decay and are possibly interpreted as due to the underlying SN and host galaxy (of 0.4 uJy in the R band). Radio observations, performed with SRT, Medicina, EVN and VLBA between day 4 and 140, suggest that the burst exploded in an environment characterised by a density profile scaling with the distance from the source (wind profile). A remarkable feature of the prompt emission is the presence of a bright flare 100 s after the trigger, lasting 70 seconds in the soft X-ray band, which was simultaneously detected from the optical band up to the MeV energy range. By combining Swift-BAT/XRT and Fermi-GBM data, the broadband (0.3-1000 keV) time resolved spectral analysis of the flare reveals the coexistence of a non-thermal (power law) and thermal blackbody components. The BB component contributes up to 35% of the luminosity in the 0.3-1000 keV band. The gamma-ray emission observed in Swift-BAT and Fermi-GBM anticipates and lasts less than the soft X-ray emission as observed by Swift-XRT, arguing against a Comptonization origin. The BB component could either be produced by an outflow becoming transparent or by the collision of a fast shell with a slow, heavy and optically thick fireball ejected during the quiescent time interval between the initial and later flares of the burst.

It has been proposed that blazar jets are structured, with a fast spine surrounded by a slower sheath or layer. This structured jet model explains some properties of their emission and morphology. Because of their relative motion, the radiation produced by one component is seen amplified by the other, thus enhancing the inverse Compton emission of both. Radiation is emitted anisotropically in the comoving frames, and causes the emitting plasma to recoil. As seen in the observer frame, this corresponds to a deceleration of the fastest component (the spine) and an acceleration of the slower one (the layer). While the deceleration of the spine has already been investigated, here we study for the first time the acceleration of the sheath and find self-consistent velocity profile solutions for both the spine and the sheath while accounting for radiative cooling. We find that the sheath can be accelerated to the velocities required by the observations if its leptons remain energetic in the acceleration region, assumed to be of the order of 100 Schwarzschild radii, demanding continuous injection of energetic particles in that region.

We derive the luminosity function and redshift distribution of short Gamma Ray Bursts (SGRBs) using (i) all the available observer-frame constraints (i.e. peak flux, fluence, peak energy and duration distributions) of the large population of Fermi SGRBs and (ii) the rest-frame properties of a complete sample of Swift SGRBs. We show that a steep $\phi(L)\propto L^{-a}$ with a>2.0 is excluded if the full set of constraints is considered. We implement a Monte Carlo Markov Chain method to derive the $\phi(L)$ and $\psi(z)$ functions assuming intrinsic Ep-Liso and Ep-Eiso correlations or independent distributions of intrinsic peak energy, luminosity and duration. To make our results independent from assumptions on the progenitor (NS-NS binary mergers or other channels) and from uncertainties on the star formation history, we assume a parametric form for the redshift distribution of SGRBs. We find that a relatively flat luminosity function with slope ~0.5 below a characteristic break luminosity ~3$\times10^{52}$ erg/s and a redshift distribution of SGRBs peaking at z~1.5-2 satisfy all our constraints. These results hold also if no Ep-Liso and Ep-Eiso correlations are assumed. We estimate that, within ~200 Mpc (i.e. the design aLIGO range for the detection of GW produced by NS-NS merger events), 0.007-0.03 SGRBs yr$^{-1}$ should be detectable as gamma-ray events. Assuming current estimates of NS-NS merger rates and that all NS-NS mergers lead to a SGRB event, we derive a conservative estimate of the average opening angle of SGRBs: $\theta_{jet}$~3-6 deg. Our luminosity function implies an average luminosity L~1.5$\times 10^{52}$ erg/s, nearly two orders of magnitude higher than previous findings, which greatly enhances the chance of observing SGRB "orphan" afterglows. Efforts should go in the direction of finding and identifying such orphan afterglows as counterparts of GW events.

If gamma-ray burst prompt emission originates at a typical radius, and if material producing the emission moves at relativistic speed, then the variability of the resulting light curve depends on the viewing angle. This is due to the fact that the pulse evolution time scale is Doppler contracted, while the pulse separation is not. For off-axis viewing angles $\theta_{\rm view} \gtrsim \theta_{\rm jet} + \Gamma^{-1}$, the pulse broadening significantly smears out the light curve variability. This is largely independent of geometry and emission processes. To explore a specific case, we set up a simple model of a single pulse under the assumption that the pulse rise and decay are dominated by the shell curvature effect. We show that such a pulse observed off-axis is (i) broader, (ii) softer and (iii) displays a different hardness-intensity correlation with respect to the same pulse seen on-axis. For each of these effects, we provide an intuitive physical explanation. We then show how a synthetic light curve made by a superposition of pulses changes with increasing viewing angle. We find that a highly variable light curve, (as seen on-axis) becomes smooth and apparently single-pulsed (when seen off-axis) because of pulse overlap. To test the relevance of this fact, we estimate the fraction of off-axis gamma-ray bursts detectable by \textit{Swift} as a function of redshift, finding that a sizable fraction (between 10\% and 80\%) of nearby ($z<0.1$) bursts are observed with $\theta_{\rm view} \gtrsim \theta_{\rm jet} + \Gamma^{-1}$. Based on these results, we argue that low luminosity gamma-ray bursts are consistent with being ordinary bursts seen off-axis.

We derive, adopting a direct method, the luminosity function and the formation rate of long Gamma Ray Bursts through a complete, flux-limited, sample of Swift bursts which has a high level of completeness in redshift z (~82%). We parametrise the redshift evolution of the GRB luminosity as L = L_0(1+ z)^k and we derive k = 2.5, consistently with recent estimates. The de-evolved luminosity function of GRBs can be represented by a broken power law with slopes a = -1.32 +- 0.21 and b = -1.84 +- 0.24 below and above, respectively, a characteristic break luminosity L_0,b = 10^51.45+-0.15 erg/s. Under the hypothesis of luminosity evolution we find that the GRB formation rate increases with redshift up to z~2, where it peaks, and then decreases in agreement with the shape of the cosmic star formation rate. We test the direct method through numerical simulations and we show that if it is applied to incomplete (both in redshift and/or flux) GRB samples it can misleadingly result in an excess of the GRB formation rate at low redshifts.

High-redshift gamma-ray bursts have several advantages for the study of the distant universe, providing unique information about the structure and properties of the galaxies in which they exploded. Spectroscopic identification with large ground-based telescopes has improved our knowledge of the class of such distant events. We present the multi-wavelength analysis of the high-$z$ Swift gamma-ray burst GRB140515A ($z = 6.327$). The best estimate of the neutral hydrogen fraction of the intergalactic medium (IGM) towards the burst is $x_{HI} \leq 0.002$. The spectral absorption lines detected for this event are the weakest lines ever observed in gamma-ray burst afterglows, suggesting that GRB140515A exploded in a very low density environment. Its circum-burst medium is characterised by an average extinction (A$_{\rm V} \sim 0.1$) that seems to be typical of $z \ge 6$ events. The observed multi-band light curves are explained either with a very flat injected spectrum ($p = 1.7$) or with a multi-component emission ($p = 2.1$). In the second case a long-lasting central engine activity is needed in order to explain the late time X-ray emission. The possible origin of GRB140515A from a Pop III (or from a Pop II stars with local environment enriched by Pop III) massive star is unlikely.

With this paper we introduce the concept of apparent structure of a GRB jet, as opposed to its intrinsic structure. The latter is customarily defined specifying the functions epsilon(theta) (the energy emitted per jet unit solid angle) and Gamma(theta) (the Lorentz factor of the emitting material); the apparent structure is instead defined by us as the isotropic equivalent energy E_iso(theta_v) as a function of the viewing angle theta_v. We show how to predict the apparent structure of a jet given its intrinsic structure. We find that a Gaussian intrinsic structure yields a power law apparent structure: this opens a new viewpoint on the Gaussian (which can be understood as a proxy for a realistic narrow, well collimated jet structure) as a possible candidate for a quasi-universal GRB jet structure. We show that such a model (a) is consistent with recent constraints on the observed luminosity function of GRBs; (b) implies fewer orphan afterglows with respect to the standard uniform model; (c) can break out the progenitor star (in the collapsar scenario) without wasting an unreasonable amount of energy; (d) is compatible with the explanation of the Amati correlation as a viewing angle effect; (e) can be very standard in energy content, and still yield a very wide range of observed isotropic equivalent energies.

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.

The structure of Gamma-Ray Burst (GRB) jets impacts on their prompt and afterglow emission properties. Insights into the still unknown structure of GRBs can be achieved by studying how different structures impact on the luminosity function (LF): i) we show that low ($10^{46} < L_{\rm iso} < 10^{48}$ erg/s) and high (i.e. with $L_{\rm iso} > 10^{50}$ erg/s) luminosity GRBs can be described by a unique LF; ii) we find that a uniform jet (seen on- and off-axis) as well as a very steep structured jet (i.e. $\epsilon(\theta) \propto \theta^{-s}$ with $s > 4$) can reproduce the current LF data; iii) taking into account the emission from the whole jet (i.e. including contributions from mildly relativistic, off-axis jet elements) we find that $E_{\rm iso}(\theta_{\rm v})$ (we dub this quantity "apparent structure") can be very different from the intrinsic structure $\epsilon(\theta)$: in particular, a jet with a Gaussian intrinsic structure has an apparent structure which is more similar to a power law. This opens a new viewpoint on the quasi-universal structured jet hypothesis.

About 15% of Gamma Ray Bursts have precursors, i.e. emission episodes preceding the main event, whose spectral and temporal properties are similar to the main emission. We propose that precursors have their own fireball, producing afterglow emission due to the dissipation of the kinetic energy via external shock. In the time lapse between the precursor and the main event, we assume that the central engine is not completely turned off, but it continues to eject relativistic material at a smaller rate, whose emission is below the background level. The precursor fireball generates a first afterglow by the interaction with the external circumburst medium. Matter injected by the central engine during the "quasi-quiescent" phase replenishes the external medium with material in relativistic motion. The fireball corresponding to the main prompt emission episode crashes with this moving material, producing a second afterglow, and finally catches up and merges with the first precursor fireball. We apply this new model to GRB 091024, an event with a precursor in the prompt light curve and two well defined bumps in the optical afterglow, obtaining an excellent agreement with the existing data.

The structure of Gamma Ray Burst (GRB) jets impacts on their prompt and afterglow emission properties. The jet of GRBs could be uniform, with constant energy per unit solid angle within the jet aperture, or it could instead be structured, namely with energy and velocity that depend on the angular distance from the axis of the jet. We try to get some insight about the still unknown structure of GRBs by studying their luminosity function. We show that low (1e46-1e48 erg/s) and high (i.e. with L > 1e50 erg/s) luminosity GRBs can be described by a unique luminosity function, which is also consistent with current lower limits in the intermediate luminosity range (1e48-1e50} erg/s). We derive analytical expressions for the luminosity function of GRBs in uniform and structured jet models and compare them with the data. Uniform jets can reproduce the entire luminosity function with reasonable values of the free parameters. A structured jet can also fit adequately the current data, provided that the energy within the jet is relatively strongly structured, i.e. E propto theta^{-k} with k > 4. The classical E propto theta^{-2} structured jet model is excluded by the current data.