We present time resolved spectral analysis of prompt emission from GRB 160625B, one of the brightest bursts ever detected by Fermi in its nine years of operations. Standard empirical functions fail to provide an acceptable fit to the GBM spectral data, which instead require the addition of a low-energy break to the fitting function. We introduce a new fitting function, called 2SBPL, consisting of three smoothly connected power laws. Fitting this model to the data, the goodness of the fits significantly improves and the spectral parameters are well constrained. We also test a spectral model that combines non-thermal and thermal (black body) components, but find that the 2SBPL model is systematically favoured. The spectral evolution shows that the spectral break is located around $E_{\rm break}\sim$ 100 keV, while the usual $\nu F_{\nu}$ peak energy feature $E_{\rm peak}$ evolves in the 0.5-6 MeV energy range. The slopes below and above $E_{\rm break}$ are consistent with the values -0.67 and -1.5, respectively, expected from synchrotron emission produced by a relativistic electron population with a low energy cut-off. If $E_{\rm break}$ is interpreted as the synchrotron cooling frequency, the implied magnetic field in the emitting region is $\sim$ 10 Gauss, i.e. orders of magnitudes smaller than the value expected for a dissipation region located at $\sim 10^{13-14}$ cm from the central engine. The low ratio between $E_{\rm peak}$ and $E_{\rm break}$ implies that the radiative cooling is incomplete, contrary to what is expected in strongly magnetized and compact emitting regions.

The prompt emission of Gamma Ray Bursts extends from the early pulses observed in gamma-rays (>15 keV) to very late flares of X-ray photons (0.3-10 keV). The duration of prompt gamma-ray pulses is rather constant while the width of X-ray flares correlates with their peak time suggesting a possible different origin. However, pulses and flares have similar spectral properties. Considering internal and external shock scenarios, we derive how the energy and duration of pulses scale with their time of occurrence and we compare with observations. The absence of an observed correlation between prompt emission pulse duration and its time of occurrence favours an "internal" origin and confirms the earlier results of Ramirez-Ruiz & Fenimore. We show that also the energetic and temporal properties of X-ray flares are consistent with being produced by internal shocks between slow fireballs with a small contrast between their bulk Lorentz factors. These results relax the requirement of a long lasting central engine to explain the latest X-ray flares.

The merger of two neutron stars is predicted to give rise to three major detectable phenomena: a short burst of gamma-rays, a gravitational wave signal, and a transient optical/near-infrared source powered by the synthesis of large amounts of very heavy elements via rapid neutron capture (the r-process). Such transients, named "macronovae" or "kilonovae", are believed to be centres of production of rare elements such as gold and platinum. The most compelling evidence so far for a kilonova was a very faint near-infrared rebrightening in the afterglow of a short gamma-ray burst at z = 0.356, although findings indicating bluer events have been reported. Here we report the spectral identification and describe the physical properties of a bright kilonova associated with the gravitational wave source GW 170817 and gamma-ray burst GRB 170817A associated with a galaxy at a distance of 40 Mpc from Earth. Using a series of spectra from ground-based observatories covering the wavelength range from the ultraviolet to the near-infrared, we find that the kilonova is characterized by rapidly expanding ejecta with spectral features similar to those predicted by current models. The ejecta is optically thick early on, with a velocity of about 0.2 times light speed, and reaches a radius of about 50 astronomical units in only 1.5 days. As the ejecta expands, broad absorption-like lines appear on the spectral continuum indicating atomic species produced by nucleosynthesis that occurs in the post-merger fast-moving dynamical ejecta and in two slower (0.05 times light speed) wind regions. Comparison with spectral models suggests that the merger ejected 0.03-0.05 solar masses of material, including high-opacity lanthanides.

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.

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 present the results of a spectroscopic search for narrow emission and absorption features in the X-ray spectra of long gamma-ray burst (GRB) afterglows. Using XMM-Newton data, both EPIC and RGS spectra, of six bright (fluence >10^{-7} erg cm^{-2}) and relatively nearby (z=0.54-1.41) GRBs, we performed a blind search for emission or absorption lines that could be related to a high cloud density or metal-rich gas in the environ close to the GRBs. We detected five emission features in four of the six GRBs with an overall statistical significance, assessed through Monte Carlo simulations, of <3.0 sigma. Most of the lines are detected around the observed energy of the oxygen edge at ~0.5 keV, suggesting that they are not related to the GRB environment but are most likely of Galactic origin. No significant absorption features were detected. A spectral fitting with a free Galactic column density (N_H) testing different models for the Galactic absorption confirms this origin because we found an indication of an excess of Galactic N_H in these four GRBs with respect to the tabulated values.

Long gamma-ray bursts (LGRBs) are associated with the deaths of massive stars and could thus be a potentially powerful tool to trace cosmic star formation. However, especially at low redshifts (z < 1.5) LGRBs seem to prefer particular types of environment. Our aim is to study the host galaxies of a complete sample of bright LGRBs to investigate the impact of the environment on GRB formation. We study host galaxy spectra of the Swift/BAT6 complete sample of 14 z < 1 bright LGRBs. We use the detected nebular emission lines to measure the dust extinction, star formation rate (SFR) and nebular metallicity (Z) of the hosts and supplement the data set with previously measured stellar masses M$_{\star}$. The distributions of the obtained properties and their interrelations (e.g. mass-metallicity and SFR-M$_{\star}$ relations) are compared to samples of field star-forming galaxies.We find that LGRB hosts at z < 1 have on average lower SFRs than if they were direct star-formation tracers. By directly comparing metallicity distributions of LGRB hosts and star-forming galaxies, we find a good match between the two populations up to Z $\sim 8.4-8.5$, after which the paucity of metal-rich LGRB hosts becomes apparent. The LGRB host galaxies of our complete sample are not inconsistent with the mass-metallicity relation at similar mean redshift and stellar masses. The cutoff against high metallicities (and high masses) can explain the low SFR values of LGRB hosts. We find a hint of increased incidence of starburst galaxies in the Swift/BAT6 z < 1 sample with respect to that of a field star-forming population. Given that the SFRs are low on average, the latter is ascribed to low stellar masses. Nevertheless the limits on the completeness and metallicity availability of current surveys, coupled with the limited number of LGRB host galaxies, prevent us from reaching more quantitative conclusions.

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.

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.

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.

Gamma ray bursts with blackbody spectra are only a few and in most cases this spectral component is accompanied by a dominating non-thermal one. Only four bursts detected by Batse have a pure blackbody spectrum throughout their duration. We present the new case of GRB 100507 detected by the Gamma Burst Monitor on board the Fermi satellite. GRB 100507 has a blackbody spectrum for the entire duration (~30 s) of the prompt emission. The blackbody temperature varies between 25 and 40 keV. The flux varies between 1e-7 and 4e-7 erg/cm2 s. There is no clear evidence of a correlation between the temperature and the blackbody flux. If the thermal emission in GRB 100507 is due to the fireballs becoming transparent, we can estimate the radius R_T and bulk Lorentz factor Gamma_T corresponding to this transition and the radius R_0 where the fireballs are created. We compare these parameters with those derived for the other four bursts with a pure blackbody spectrum. In all but one burst, for fiducial assumptions on the radiative efficiency and distance of the sources, R_0~1e9-1e10 cm, i.e. much larger than the gravitational radius of a few solar mass black hole. Possible solutions of this apparent inconsistency are tentatively discussed considering the dependence of R_0 on the unknown parameters. Alternatively, such a large R_0 could be where the fireball, still opaque, converts most of its kinetic energy into internal energy (due to the impact with some material left over by the progenitor star) and starts to re-accelerate.