On January 14, 2019, the Major Atmospheric Gamma Imaging Cherenkov telescopes detected GRB 190114C above 0.2 TeV, recording the most energetic photons ever observed from a gamma-ray burst. We use this unique observation to probe an energy dependence of the speed of light in vacuo for photons as predicted by several quantum gravity models. Based on a set of assumptions on the possible intrinsic spectral and temporal evolution, we obtain competitive lower limits on the quadratic leading order of speed of light modification.

Gamma-ray bursts (GRBs) of the long-duration class are the most luminous sources of electromagnetic radiation known in the Universe. They are generated by outflows of plasma ejected at near the speed of light by newly formed neutron stars or black holes of stellar mass at cosmological distances. Prompt flashes of MeV gamma rays are followed by longer-lasting afterglow emission from radio waves to GeV gamma rays, due to synchrotron radiation by energetic electrons in accompanying shock waves. Although emission of gamma rays at even higher, TeV energies by other radiation mechanisms had been theoretically predicted, it had never been detected previously. Here we report the clear detection of GRB 190114C in the TeV band, achieved after many years of dedicated searches for TeV emission from GRBs. Gamma rays in the energy range 0.2--1 TeV are observed from about 1 minute after the burst (at more than 50 standard deviations in the first 20 minutes). This unambiguously reveals a new emission component in the afterglow of a GRB, whose power is comparable to that of the synchrotron component. The observed similarity in the radiated power and temporal behaviour of the TeV and X-ray bands points to processes such as inverse Compton radiation as the mechanism of the TeV emission, while processes such as synchrotron emission by ultrahigh-energy protons are disfavoured due to their low radiative efficiency.

Long-duration gamma-ray bursts (GRBs) originate from ultra-relativistic jets launched from the collapsing cores of dying massive stars. They are characterised by an initial phase of bright and highly variable radiation in the keV-MeV band that is likely produced within the jet and lasts from milliseconds to minutes, known as the prompt emission. Subsequently, the interaction of the jet with the external medium generates external shock waves, responsible for the afterglow emission, which lasts from days to months, and occurs over a broad energy range, from the radio to the GeV bands. The afterglow emission is generally well explained as synchrotron radiation by electrons accelerated at the external shock. Recently, an intense, long-lasting emission between 0.2 and 1 TeV was observed from the GRB 190114C. Here we present the results of our multi-frequency observational campaign of GRB~190114C, and study the evolution in time of the GRB emission across 17 orders of magnitude in energy, from $5\times10^{-6}$ up to $10^{12}$\,eV. We find that the broadband spectral energy distribution is double-peaked, with the TeV emission constituting a distinct spectral component that has power comparable to the synchrotron component. This component is associated with the afterglow, and is satisfactorily explained by inverse Compton upscattering of synchrotron photons by high-energy electrons. We find that the conditions required to account for the observed TeV component are not atypical, supporting the possibility that inverse Compton emission is commonly produced in GRBs.

We discuss the new surprising observational results that indicate quite convincingly that the prompt emission of Gamma-Ray Bursts (GRBs) is due to synchrotron radiation produced by a particle distribution that has a low energy cut-off. The evidence of this is provided by the low energy part of the spectrum of the prompt emission, that shows the characteristic F(nu) \propto nu^(1/3) shape followed by F(nu) \propto nu^(-1/2) up to the peak frequency. This implies that although the emitting particles are in fast cooling, they do not cool completely. This poses a severe challenge to the basic ideas about how and where the emission is produced, because the incomplete cooling requires a small value of the magnetic field, to limit synchrotron cooling, and a large emitting region, to limit the self-Compton cooling, even considering Klein-Nishina scattering effects. Some new and fundamental ingredient is required for understanding the GRBs prompt emission. We propose proton-synchrotron as a promising mechanism to solve the incomplete cooling puzzle.

The Gamma Ray Burst (GRB) 180720B is one of the brightest events detected by the Fermi satellite and the first GRB detected by the H.E.S.S. telescope above 100 GeV. We analyse the Fermi (GBM and LAT) and Swift (XRT and BAT) data and describe the evolution of the burst spectral energy distribution in the 0.5 keV - 10 GeV energy range over the first 500 seconds of emission. We reveal a smooth transition from the prompt phase, dominated by synchrotron emission in a moderately fast cooling regime, to the afterglow phase whose emission has been observed from the radio to the GeV energy range. The LAT (0.1 - 100 GeV) light curve initially rises ($F_{\rm LAT}\propto t^{2.4}$), peaks at $\sim$78 s, and falls steeply ($F_{\rm LAT}\propto t^{-2.2}$) afterwards. The peak, which we interpret as the onset of the fireball deceleration, allows us to estimate the bulk Lorentz factor $\Gamma_{0}\sim 150 \ (300)$ under the assumption of a wind-like (homogeneous) circum-burst medium density. We derive a flux upper limit in the LAT energy range at the time of H.E.S.S. detection, but this does not allow us to unveil the nature of the high energy component observed by H.E.S.S. We fit the prompt spectrum with a physical model of synchrotron emission from a non-thermal population of electrons. The 0 - 35 s spectrum after its $E F(E)$ peak (at 1 - 2 MeV) is a steep power law extending to hundreds of MeV. We derive a steep slope of the injected electron energy distribution $N(\gamma)\propto \gamma^{-5}$. Our fit parameters point towards a very low magnetic field ($B'\sim 1 $ G) in the emission region.

The first two LIGO and Virgo observation runs have been important milestones in the gravitational wave (GW) field, thanks to the detection of GW signals from ten binary black hole systems and a binary neutron star system. In order to fully characterize the emitting source, the remnant object and its environment, electromagnetic follow-up observations at different wavelengths are essential, as learned from the GW170817/GRB170817A case. Given the quite large localization uncertainties provided by interferometers, the main challenge faced by facilities with a narrow field of view (e.g. Imaging Atmospheric Cherenkov Telescopes, IACTs) is to setup a suitable follow-up strategy in order to observe sky regions with the highest probability to host the electromagnetic (EM) counterpart of the GW signal. As member of the EM follow-up community, the MAGIC collaboration joined this effort in 2014. As the third observation run (O3) is currently ongoing, where both LIGO and Virgo are expected to have much better sensitivities, MAGIC is refining its follow-up strategy to maximize the chances of observing the EM counterparts as soon as possible. In this contribution we will describe this strategy, focusing on the different observation cases, which mainly depends on the information available from both GW and EM partner facilities.

Gamma-Ray Bursts (GRBs) are one of the main targets for current and next generation Imaging Atmospheric Cherenkov Telescopes (IACTs). Given their transient behavior, especially in the case of their prompt emission phase, performing fast follow-up observations is challenging for IACTs, which have a narrow field of view and limited duty cycle. Despite this, MAGIC plays a major role in the search for Very High Energy (VHE, E>100 GeV) gamma-ray emission from GRBs: this is possible thanks to its fast repositioning speed, low energy threshold and high sensitivity at the lowest energies. In 2013 the MAGIC GRB automatic procedure was upgraded, increasing the number of GRBs followed in the prompt and early afterglow phases and decreasing dramatically hardware failures during fast repositioning. Currently, only GRB 190114C was firmly detected in the VHE band, while for other GRBs no significant detection was achieved. In such a case, upper limits (ULs) can give insight into the physics driving such eluding sources, especially on their emission mechanisms. In this contribution we report on the status of the GRB follow-up with MAGIC and focus on the ULs and results obtained from a sample of GRBs observed between 2013 and 2018. This GRB catalog is the result of the MAGIC well-designed and tested follow-up procedure, and it serves as a precursor of GRBs observation with the next generation IACT system, the Cherenkov Telescope Array (CTA).

The birth of gravitational-wave / electromagnetic astronomy was heralded by the joint observation of gravitational waves (GWs) from a binary neutron star (BNS) merger by Advanced LIGO and Advanced Virgo, GW170817, and of gamma-rays from the short gamma-ray burst GRB170817A by the Fermi Gamma-ray Burst Monitor (GBM) and INTEGRAL. This detection provided the first direct evidence that at least a fraction of BNSs are progenitors of short GRBs. GRBs are now also known to emit very-high-energy (VHE, > 100 GeV) photons as has been shown by recent independent detections of the GRBs 1901114C and 180720B by the ground-based gamma-ray detectors MAGIC and H.E.S.S. In the next years, the Cherenkov Telescope Array (CTA) will boost the searches for VHE counterparts thanks to its unprecedented sensitivity, rapid response and capability to monitor large sky areas via survey-mode operation. In this contribution, we present the CTA program of observations following the detection of GW events. We discuss various follow-up strategies and links to multi-wavelength and multi-messenger observations. Finally we outline the capabilities and prospects of detecting VHE emission from GW counterparts.

One of the central scientific goals of the next-generation Cherenkov Telescope Array (CTA) is the detection and characterization of gamma-ray bursts (GRBs). CTA will be sensitive to gamma rays with energies from about 20 GeV, up to a few hundred TeV. The energy range below 1 TeV is particularly important for GRBs. CTA will allow exploration of this regime with a ground-based gamma-ray facility with unprecedented sensitivity. As such, it will be able to probe radiation and particle acceleration mechanisms at work in GRBs. In this contribution, we describe POSyTIVE, the POpulation Synthesis Theory Integrated project for very high-energy emission. The purpose of the project is to make realistic predictions for the detection rates of GRBs with CTA, to enable studies of individual simulated GRBs, and to perform preparatory studies for time-resolved spectral analyses. The mock GRB population used by POSyTIVE is calibrated using the entire 40-year dataset of multi-wavelength GRB observations. As part of this project we explore theoretical models for prompt and afterglow emission of long and short GRBs, and predict the expected radiative output. Subsequent analyses are performed in order to simulate the observations with CTA, using the publicly available ctools and Gammapy frameworks. We present preliminary results of the design and implementation of this project.

Information on the spectral shape of prompt emission in gamma-ray bursts (GRB) is mostly available only at energies $\gtrsim10$ keV, where the main instruments for GRB detection are sensitive. The origin of this emission is still very uncertain because of the apparent inconsistency with synchrotron radiation, which is the most obvious candidate, and the resulting need for considering less straightforward scenarios. The inclusion of data down to soft X-rays ($\sim$ 0.5 keV), which are available only in a small fraction of GRBs, has firmly established the common presence of a spectral break in the low-energy part of prompt spectra, and the consistency of the overall spectral shape with synchrotron radiation in the moderately fast-cooling regime, the low-energy break being identified with the cooling frequency. In this work we further extend the range of investigation down to the optical band. In particular, we test the synchrotron interpretation by directly fitting a theoretically derived synchrotron spectrum and making use of optical to gamma-ray data. Secondly, we test an alternative model that considers the presence of a black-body component at $\sim$keV energies, in addition to a non-thermal component that is responsible for the emission at the spectral peak (100 keV-1 MeV). We find that synchrotron radiation provides a good description of the broadband data, while models composed of a thermal and a non-thermal component require the introduction of a low-energy break in the non-thermal component in order to be consistent with optical observations. Motivated by the good quality of the synchrotron fits, we explore the physical parameter space of the emitting region. In a basic prompt emission scenario we find quite contrived solutions for the magnetic field strength (5 G $<B^\prime<40$ G) and for the location of the region where the radiation is produced ($R_\gamma>10^{16}$ cm).

The long-lasting tension between the observed spectra of gamma ray bursts (GRBs) and the predicted synchrotron emission spectrum might be solved if electrons do not completely cool. Evidence for incomplete cooling was recently found in Swift GRBs with prompt observations down to 0.1 keV and in one bright Fermi burst, GRB 160625B. Here we systematically search for evidence of incomplete cooling in the spectra of the ten brightest short and long GRBs observed by Fermi. We find that in 8/10 long GRBs there is compelling evidence of a low energy break (below the peak energy) and good agreement with the photon indices of the synchrotron spectrum (respectively -2/3 and -3/2 below the break and between the break and the peak energy). Interestingly, none of the ten short GRBs analysed shows a break but the low energy spectral slope is consistent with -2/3. In a standard scenario, these results imply a very low magnetic field in the emission region (B' ~ 10 G in the comoving frame), at odd with expectations.

We update a flux-limited complete sample of Swift-based SGRBs (SBAT4, D'Avanzo et al. 2014), bringing it to 25 events and doubling its previous redshift range. We then evaluate the column densities of the events in the updated sample, in order to compare them with the NH distribution of LGRBs, using the sample BAT6ext (Arcodia et al. 2016). We rely on Monte Carlo simulations of the two populations and compare the computed NH distributions with a two sample Kolmogorov Smirnov (K-S) test. We then study how the K-S probability varies with respect to the redshift range we consider. We find that the K-S probability keeps decreasing as redshift increases until at z$\sim$1.8 the probability that short and long GRBs come from the same parent distribution drops below 1$\%$. This testifies for an observational difference among the two populations. This difference may be due to the presence of highly absorbed LGRBs above z$\sim$1.3, which have not been observed in the SGRB sample yet, although this may be due to our inability to detect them, or to the relatively small sample size.

We study the problem of the escape and transport of Cosmic-Rays (CR) from a source embedded in a fully ionised, hot phase of the interstellar medium (HIM). In particular, we model the CR escape and their propagation in the source vicinity taking into account excitation of Alfv\'enic turbulence by CR streaming and mechanisms damping the self-excited turbulence itself. Our estimates of escape radii and times result in large values (100 pc, $2\times10^5$ yr) for particle energies $\lesssim20$ GeV and smaller values for particles with increasing energies (35 pc and 14 kyr at 1 TeV). These escape times and radii, when used as initial conditions for the CR propagation outside the source, result in relevant suppression of the diffusion coefficient (by a factor 5-10) on time-scales comparable with their (energy dependent) escape time-scale. The damping mechanisms are fast enough that even on shorter time scales the Alfv\'enic turbulence is efficiently damped, and the ratio between random and ordered component of the magnetic field is $\delta B/B_0\ll 1$, justifying the use of quasi-linear theory. In spite of the suppressed diffusion coefficient, and then the increased residence time in the vicinity (<200 pc) of their source, the grammage accumulated by CRs after their escape is found to be negligible (at all energies) as compared to the one accumulated while diffusing in the whole Galaxy, due to the low density of the HIM.

The binary neutron star merger event GW170817 was detected through both electromagnetic radiation and gravitational waves. Its afterglow emission may have been produced by either a narrow relativistic jet or an isotropic outflow. High spatial resolution measurements of the source size and displacement can discriminate between these scenarios. We present Very Long Baseline Interferometry observations, performed 207.4 days after the merger, using a global network of 32 radio telescopes. The apparent source size is constrained to be smaller than 2.5 milliarcseconds at the 90% confidence level. This excludes the isotropic outflow scenario, which would have produced a larger apparent size, indicating that GW170817 produced a structured relativistic jet. Our rate calculations show that at least 10% of neutron star mergers produce such a jet.

e-ASTROGAM (enhanced ASTROGAM) is a breakthrough Observatory space mission, with a detector composed by a Silicon tracker, a calorimeter, and an anticoincidence system, dedicated to the study of the non-thermal Universe in the photon energy range from 0.3 MeV to 3 GeV - the lower energy limit can be pushed to energies as low as 150 keV for the tracker, and to 30 keV for calorimetric detection. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LIGO-Virgo-GEO600-KAGRA, SKA, ALMA, E-ELT, TMT, LSST, JWST, Athena, CTA, IceCube, KM3NeT, and LISA.

Detection of prompt emission by Swift-XRT provides a unique tool to study how the prompt spectrum of gamma-ray bursts (GRBs) extends down to the soft X-ray band. This energy band is particularly important for prompt emission studies, since it is towards low energies that the observed spectral shape is in disagreement with the synchrotron predictions. Unfortunately, the number of cases where XRT started observing the GRB location during the prompt phase is very limited. In this work we collect a sample of 34 GRBs and perform joint XRT+BAT spectral analysis of prompt radiation, extending a previous study focused on the 14 brightest cases. Fermi-GBM observations are included in the analysis when available (13 cases), allowing the characterization of prompt spectra from soft X-rays to MeV energies. In $62\%$ of the spectra, the XRT data reveal a hardening of the spectrum, well described by introducing an additional, low-energy power-law segment (with index $\alpha_1$) into the empirical fitting function. The break energy below which the spectrum hardens takes values between $\rm 3\,keV$ and $\rm 30\,keV$. A second power-law ($\alpha_2$) describes the spectrum between the break energy and the peak energy. The mean values of the photon indices are $\langle\alpha_1\rangle=-0.51$ ($\sigma=0.29$) and $\langle\alpha_2\rangle=-1.54$ ($\sigma=0.26$). These are consistent, within one $\sigma$, with the synchrotron values in fast cooling regime. As a test, if we exclude XRT data from the fits we find typical results: the spectrum below the peak energy is described by a power law with $\langle\alpha\rangle=-1.15$. This shows the relevance of soft X-ray data in revealing prompt emission spectra consistent with synchrotron spectra. Finally, we do not find any correlation between the presence of the X-ray break energy and the flux, fluence or duration of the prompt emission.

We report our observation of the short GRB 170817A, associated to the binary neutron star merger event GW 170817, perfomed in the X-ray band with XMM-Newton 135 d after the event (on the 29th December 2017). We find evidence for a flattening of the X-ray light curve with respect to the previously observed brightening. This is also supported by a nearly simultaneous optical Hubble Space Telescope and successive X-ray Chandra and low-frequency radio observations recently reported in the literature. Since the optical-to-X-ray spectral slope did not change with respect to previous observations, we exclude that the change in the temporal evolution of the light curve is due to the passage of the cooling frequency: its origin must be geometric or dynamical. We interpret all the existing afterglow data with two models: i) a structured jet and ii) a jet-less isotropic fireball with some stratification in its radial velocity structure. Both models fit the data and predict that the radio flux must decrease simultaneously with the optical and the X-ray one, making hard to distinguish between them at the present stage. Polarimetric measures and the rate of short GRB-GW association in future LIGO/Virgo runs will be key to disentangle these two geometrically different scenarios.

The number of Gamma-Ray Bursts (GRBs) detected at high energies ($\sim\,0.1-100$ GeV) has seen a rapid increase over the last decade, thanks to observations from the Fermi-Large Area Telescope. The improved statistics and quality of data resulted in a better characterisation of the high-energy emission properties and in stronger constraints on theoretical models. In spite of the many achievements and progresses, several observational properties still represent a challenge for theoretical models, revealing how our understanding is far from being complete. This paper reviews the main spectral and temporal properties of $\sim\,0.1-100$ GeV emission from GRBs and summarises the most promising theoretical models proposed to interpret the observations. Since a boost for the understanding of GeV radiation might come from observations at even higher energies, the present status and future prospects for observations at very-high energies (above $\sim$ 100 GeV) are also discussed. The improved sensitivity of upcoming facilities, coupled to theoretical predictions, supports the concrete possibility for future ground GRB detections in the high/very-high energy domain.

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.

Recent advances in fitting prompt emission spectra in gamma-ray bursts (GRBs) are boosting our understanding of the still elusive origin of this radiation. These progresses have been possible thanks to a more detailed analysis of the low-energy part ($<$\,100\,keV) of the prompt spectrum, where the spectral shape is sometimes found to deviate from a simple power-law shape. This deviation is well described by a spectral break or, alternatively by the addition of a thermal component. Spectral data extending down to less than 1\,keV are extremely relevant for these studies, but presently they are available only for a small subsample of {\it Swift} GRBs observed by XRT (the X-ray telescope, 0.3-10\,keV) during the prompt emission. The space mission \th\ will allow a systematic study of prompt spectra from 0.3\,keV to several MeV. We show that observations performed by \th\ will allow us to discriminate between different models presently considered for GRB prompt studies, solving the long-standing open issue about the nature of the prompt radiation, with relevant consequences on the location of the emitting region, magnetic field strength and presence of thermal components.

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

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.

The radiative process responsible for gamma-Ray Burst (GRB) prompt emission has not been identified yet. If dominated by fast-cooling synchrotron radiation, the part of the spectrum immediately below the $\nu F_\nu$ peak energy should display a power-law behavior with slope $\alpha_2=-3/2$, which breaks to a higher value $\alpha_1=-2/3$ (i.e. to a harder spectral shape) at lower energies. Prompt emission spectral data (usually available down to $\sim10-20\,$keV) are consistent with one single power-law behavior below the peak, with typical slope $\langle\alpha\rangle=-1$, higher than (and then inconsistent with) the expected value $\alpha_2=-3/2$. To better characterize the spectral shape at low energy, we analyzed 14 GRBs for which the Swift X-ray Telescope started observations during the prompt. When available, Fermi-GBM observations have been included in the analysis. For 67% of the spectra, models that usually give a satisfactory description of the prompt (e.g., the Band model) fail in reproducing the $0.5-1000\,$keV spectra: low-energy data outline the presence of a spectral break around a few keV.We then introduce an empirical fitting function that includes a low-energy power law $\alpha_1$, a break energy $E_{\rm break}$, a second power law $\alpha_2$, and a peak energy $E_{\rm peak}$. We find $\langle\alpha_1\rangle=-0.66$ ($ \rm \sigma=0.35$), $\langle \log (E_{\rm break}/\rm keV)\rangle=0.63$ ($ \rm \sigma=0.20$), $\langle\alpha_2\rangle=-1.46$ ($\rm \sigma=0.31$), and $\langle \log (E_{\rm peak}/\rm keV)\rangle=2.1$ ($ \rm \sigma=0.56$).The values $\langle\alpha_1\rangle$ and $\langle\alpha_2\rangle$ are very close to expectations from synchrotron radiation. In this context, $E_{\rm break}$ corresponds to the cooling break frequency.

Most of the diffuse Galactic GeV gamma-ray emission is produced via collisions of cosmic ray (CR) protons with ISM protons. As such the observed spectra of the gamma-rays and the CRs should be strongly linked. Recent observations of Fermi-LAT exhibit a hardening of the gamma-ray spectrum at around a hundred GeV, between the Sagittarius and Carina tangents, and a further hardening at a few degrees above and below the Galactic plane. However, standard CR propagation models that assume a time independent source distribution and a location independent diffusion cannot give rise to a spatially dependent CR (and hence gamma-ray) spectral slopes. Here we consider a dynamic spiral arm model in which the distribution of CR sources is concentrated in the (dynamic) spiral arms, and we study the effects of this model on the $\pi^0$-decay produced gamma-ray spectra. Within this model, near the Galactic arms the observed gamma-ray spectral slope is not trivially related to the CR injection spectrum and energy dependence of the diffusion coefficient. We find unique signatures that agree with the Fermi-LAT observations. This model also provides a physical explanation for the difference between the local CR spectral slope and the CR slope inferred from the average gamma-ray spectrum.

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 is largely recognized that Gamma-Ray Burst (GRB) jets involve ultra-relativistic motion. However, the value of the Lorentz factor Gamma_0 is still not clear and only lower limits are known for most bursts. We suggest here a new method to obtain upper limits on Gamma_0. The early high-energy synchrotron afterglow flux depends strongly on Gamma_0. Upper limits on GeV emission therefore provide uppers limit on Gamma_0. Applying this method to 190 Fermi GRBs that have not been detected by the Fermi-LAT we place upper limits on the high-energy afterglow flux, and in turn on Gamma_0. For bursts at a typical redshift z=2, we find values of the order of 200 (and above) for a homogeneous density medium, and in the range 100-400 for a wind-like medium. These upper limits are consistent with (and are very close to) lower limits and direct estimates inferred using other methods, suggesting that the typical Lorentz factors of GRB jets are of order a few hundred.

Supernova remnants are believed to be the main sources of galactic Cosmic Rays (CR). Within this framework, particles are accelerated at supernova remnant shocks and then released in the interstellar medium. The mechanism through which CRs are released and the way in which they propagate still remain open issues. The main difficulty is the high non-linearity of the problem: CRs themselves excite the magnetic turbulence that confines them close to their sources. We solve numerically the coupled differential equations describing the evolution in space and time of the escaping particles and of the waves generated through the CR streaming instability. The warm ionized and warm neutral phases of the interstellar medium are considered. These phases occupy the largest fraction of the disc volume, where most supernovae explode, and are characterised by the significant presence of neutral particles. The friction between those neutrals and ions results in a very effective wave damping mechanism. It is found that streaming instability affects the propagation of CRs even in the presence of ion-neutral friction. The diffusion coefficient can be suppressed by more than a factor of $\sim 2$ over a region of few tens of pc around the remnant. The suppression increases for smaller distances. The propagation of $\approx 10$ GeV particles is affected for several tens of kiloyears after escape, while $\approx 1$ TeV particles are affected for few kiloyears. This might have a great impact on the interpretation of gamma-ray observations of molecular clouds located in the vicinity of supernova remnants.

The prompt Gamma-Ray Bursts' (GRBs) efficiency is an important clue on the emission mechanism producing the $\gamma$-rays. Previous estimates of the kinetic energy of the blast waves, based on the X-ray afterglow luminosity $L_X$, suggested that this efficiency is large, with values above 90\% in some cases. This poses a problem to emission mechanisms and in particular to the internal shocks model. These estimates are based, however, on the assumption that the X-ray emitting electrons are fast cooling and that their Inverse Compton (IC) losses are negligible. The observed correlations between $L_X$ (and hence the blast wave energy) and $E_{\gamma\rm ,iso}$, the isotropic equivalent energy in the prompt emission, has been considered as observational evidence supporting this analysis. It is reasonable that the prompt gamma-ray energy and the blast wave kinetic energy are correlated and the observed correlation corroborates, therefore, the notion $L_X$ is indeed a valid proxy for the latter. Recent findings suggest that the magnetic field in the afterglow shocks is significantly weaker than was earlier thought and its equipartition fraction, $\epsilon_B$, could be as low as $10^{-4}$ or even lower. Motivated by these findings we reconsider the problem, taking now IC cooling into account. We find that the observed $L_X-E_{\gamma\rm ,iso}$ correlation is recovered also when IC losses are significant. For small $\epsilon_B$ values the blast wave must be more energetic and we find that the corresponding prompt efficiency is significantly smaller than previously thought. For example, for $\epsilon_B\sim10^{-4}$ we infer a typical prompt efficiency of $\sim15\%$.

Polarization measurements from relativistic outflows are a valuable tool to probe the geometry of the emission region and the microphysics of the particle distribution. Indeed, the polarization level depends on: (i) the local magnetic field orientation, (ii) the geometry of the emitting region with respect to the line of sight, and (iii) the electron pitch-angle distribution. Here we consider optically thin synchrotron emission and we extend the theory of circular polarization from a point source to an extended radially expanding relativistic jet. We present numerical estimates for both linear and circular polarization in such systems. We consider different configurations of the magnetic field, spherical and jetted outflows, isotropic and anisotropic pitch-angle distributions, and outline the difficulty in obtaining the reported high level of circular polarization observed in the afterglow of GRB 121024A. We conclude that the origin of the observed polarization cannot be intrinsic to an optically thin synchrotron process, even when the electron pitch-angle distribution is extremely anisotropic.

Recent advances on gamma-ray observations from SuperNova Remnants and Molecular Clouds offer the possibility to study in detail the properties of the propagation of escaping Cosmic Rays (CR). However, a complete theory for CR transport outside the acceleration site has not been developed yet. Two physical processes are thought to be relevant to regulate the transport: the growth of waves caused by streaming instability, and possible wave damping mechanisms that reduce the growth of the turbulence. Only a few attempts have been made so far to incorporate these mechanisms in the theory of CR diffusion. In this work we present recent advances in this subject. In particular, we show results obtained by solving the coupled equations for the diffusion of CRs and the evolution of Alfven waves. We discuss the importance of streaming instabilities and wave damping in different ISM phases.

We consider a sample of ten GRBs with long lasting ($\gtrsim10^2\rm\,sec$) emission detected by Fermi/LAT and for which X-ray data around $1\,$day are also available. We assume that both the X-rays and the GeV emission are produced by electrons accelerated at the external forward shock, and show that the X-ray and the GeV fluxes lead to very different estimates of the initial kinetic energy of the blast wave. The energy estimated from GeV is on average $\sim50$ times larger than the one estimated from X-rays. We model the data (accounting also for optical detections around $1\,$day, if available) to unveil the reason for this discrepancy and find that good modelling within the forward shock model is always possible and leads to two possibilities: either the X-ray emitting electrons (unlike the GeV emitting electrons) are in the slow cooling regime or ii) the X-ray synchrotron flux is strongly suppressed by Compton cooling, whereas, due to the Klein-Nishina suppression, this effect is much smaller at GeV energies. In both cases the X-ray flux is no longer a robust proxy for the blast wave kinetic energy. On average, both cases require weak magnetic fields ($10^{-6}\lesssim \epsilon_B \lesssim 10^{-3}$) and relatively large isotropic kinetic blast wave energies $10^{53}\rm\,erg<E_{0,kin}<10^{55}\rm\,erg$ corresponding to large lower limits on the collimated energies, in the range $10^{52}\rm\,erg<E_{\theta,kin}<5\times10^{52}\rm\,erg$ for an ISM environment with $n\sim 1\mbox{cm}^{-3}$ and $10^{52}\rm\,erg<E_{\theta,kin}<10^{53}\rm\,erg$ for a wind environment with $A_* \sim 1$. These energies are larger than those estimated from the X-ray flux alone, and imply smaller inferred values of the prompt efficiency mechanism, reducing the efficiency requirements on the still uncertain mechanism responsible for prompt emission.

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.

The gamma-ray observations of molecular clouds associated with supernova remnants are considered one of the most promising ways to search for a solution of the problem of cosmic ray origin. Here we briefly review the status of the field, with particular emphasis on the theoretical and phenomenological aspects of the problem.

A certain degree of linear polarization has been measured in several GRB afterglows. More surprisingly, circular polarization has been recently measured in GRB121024A. For synchrotron emission, the polarization level depends on: (i) the local magnetic field orientation (ii) the geometry of the emitting region with respect to the line of sight and (iii) the electron pitch-angle distribution. For this reason, polarization measurements are a valuable tool to probe afterglow micro-physics. We present numerical estimates of linear and circular polarization for different configurations (i.e., magnetic fields, geometries and pitch-angle distributions). For each different scenario, we study the conditions for reaching the maximum and minimum linear and circular polarization and provide their values. We discuss the implication of our results to the micro-physics of GRB afterglows in view of recent polarization measurements.

We examine GRBs with both Fermi-LAT and X-ray afterglow data. Assuming that the 100MeV (LAT) emission is radiation from cooled electrons accelerated by external shocks, we show that the kinetic energy of the blast wave estimated from the 100MeV flux is 50 times larger than the one estimated from the X-ray flux. This can be explained if either: i) electrons radiating at X-rays are significantly cooled by SSC (suppressing the synchrotron flux above the cooling frequency) or ii) if the X-ray emitting electrons, unlike those emitting at 100MeV energies, are in the slow cooling regime. In both cases the X-ray flux is no longer an immediate proxy of the blast wave kinetic energy. We model the LAT, X-ray and optical data and show that in general these possibilities are consistent with the data, and explain the apparent disagreement between X-ray and LAT observations. All possible solutions require weak magnetic fields: $10^{-6}< \epsilon_B < 10^{-3}$ (where $\epsilon_B$ is the fraction of shocked plasma energy in magnetic fields). Using the LAT emission as a proxy for the blast wave kinetic energy we find that the derived prompt efficiencies are of order 15%. This is considerably lower compared with previous estimates (87% and higher for the same bursts). This provides at least a partial solution to the "prompt high efficiency paradox".

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

The validity of the $E_{\rm p}$ - $E_{\rm iso}$ correlation in gamma-ray bursts and the possibility of explaining the prompt emission with internal shocks are highly debated questions. We study whether the $E_{\rm p}$ - $E_{\rm iso}$ correlation can be reproduced if internal shocks are indeed responsible for the prompt emission, or conversely, if the correlation can be used to constrain the internal shock scenario. We developed a toy model where internal shocks are limited to the collision of only two shells. Synthetic burst populations were constructed for various distributions of the model parameters, such as the injected power in the relativistic outflow, the average Lorentz factor, and its typical contrast between the shells. These parameters can be independent or linked by various relations. Synthetic $E_{\rm p}$ - $E_{\rm iso}$ diagrams are obtained in the different cases and compared with the observed correlation. The reference observed correlation is the one defined by the BAT6 sample, a sample of Swift bursts almost complete in redshift and affected by well-known and reproducible instrumental selection effects. The comparison is then performed with a subsample of synthetic bursts that satisfy the same selection criteria as were imposed on the BAT6 sample. A satisfactory agreement between model and data can often be achieved, but only if several strong constraints are satisfied on both the dynamics of the flow and the microphysics that governs the redistribution of the shock-dissipated energy.

We investigated the rest frame spectral lags of two complete samples of bright long (50) and short (6) gamma-ray bursts (GRB) detected by Swift. We analysed the Swift/BAT data through a discrete cross-correlation function (CCF) fitted with an asymmetric Gaussian function to estimate the lag and the associated uncertainty. We find that half of the long GRBs have a positive lag and half a lag consistent with zero. All short GRBs have lags consistent with zero. The distributions of the spectral lags for short and long GRBs have different average values. Limited by the small number of short GRBs, we cannot exclude at more than 2 sigma significance level that the two distributions of lags are drawn from the same parent population. If we consider the entire sample of long GRBs, we do not find evidence for a lag-luminosity correlation, rather the lag-luminosity plane appears filled on the left hand side, thus suggesting that the lag-luminosity correlation could be a boundary. Short GRBs are consistent with the long ones in the lag-luminosity plane.

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 physical origin of the >0.1 GeV emission detected from Gamma-Ray Bursts (GRBs) by the Fermi satellite has not yet been completely understood. In this work we consider the GeV light curves of ten GRBs with measured redshift detected by the Fermi-LAT. These light curves are characterised by a long-lived ($\gtrsim10^2$ seconds) emission, whose luminosity decays in time as a power-law. While the decay rate is similar for all GRBs (i.e. $L_{LAT}\propto t^{-1.2}$), the normalisation spans about two orders of magnitude in luminosity. However, after re-normalising the luminosities to the prompt energetics $E_{iso}$ the light curves overlap. We consider the scenario in which the temporally extended LAT emission is dominated by synchrotron radiation from electrons accelerated at the forward external shock. According to this model, at high-energies (i.e. above the typical synchrotron frequencies) a small dispersion of the $E_{iso}$-normalised light curves is expected. The fact that the LAT temporally extended emission follows this behaviour reinforces its interpretation in terms of afterglow radiation from external shocks. Assuming this scenario, we argue that the parameters $\epsilon_e$ and $\eta_\gamma$ (i.e., the fraction of shock-dissipated energy gained by the electrons, and the efficiency of the mechanism producing the prompt radiation, respectively) must be narrowly distributed.

We present a carefully selected sample of short gamma-ray bursts (SGRBs) observed by the Swift satellite up to June 2013. Inspired by the criteria we used to build a similar sample of bright long GRBs (the BAT6 sample), we selected SGRBs with favorable observing conditions for the redshift determination on ground, ending up with a sample of 36 events, almost half of which with a redshift measure. The redshift completeness increases up to about 70% (with an average redshift value of z = 0.85) by restricting to those events that are bright in the 15-150 keV Swift Burst Alert Telescope energy band. Such flux-limited sample minimizes any redshift-related selection effects, and can provide a robust base for the study of the energetics, redshift distribution and environment of the Swift bright population of SGRBs. For all the events of the sample we derived the prompt and afterglow emission in both the observer and (when possible) rest frame and tested the consistency with the correlations valid for long GRBs. The redshift and intrinsic X-ray absorbing column density distributions we obtain are consistent with the scenario of SGRBs originated by the coalescence of compact objects in primordial binaries, with a possible minor contribution (~10%-25%) of binaries formed by dynamical capture (or experiencing large natal kicks). This sample is expected to significantly increase with further years of Swift activity.

We present the rest-frame light curves in the optical and X-ray bands of an unbiased and complete sample of Swift long Gamma-Ray Bursts (GRBs), namely the BAT6 sample. The unbiased BAT6 sample (consisting of 58 events) has the highest level of completeness in redshift ($\sim$ 95%), allowing us to compute the rest-frame X-ray and optical light curves for 55 and 47 objects, respectively. We compute the X-ray and optical luminosities accounting for any possible source of absorption (Galactic and intrinsic) that could affect the observed fluxes in these two bands. We compare the behaviour observed in the X-ray and in the optical bands to assess the relative contribution of the emission during the prompt and afterglow phases. We unarguably demonstrate that the GRBs rest-frame optical luminosity distribution is not bimodal, being rather clustered around the mean value Log(L$_{R}$) = 29.9 $\pm$ 0.8 when estimated at a rest frame time of 12 hr. This is in contrast with what found in previous works and confirms that the GRB population has an intrinsic unimodal luminosity distribution. For more than 70% of the events the rest-frame light curves in the X-ray and optical bands have a different evolution, indicating distinct emitting regions and/or mechanisms. The X-ray light curves normalised to the GRB isotropic energy (E$_{\rm iso}$), provide evidence for X-ray emission still powered by the prompt emission until late times ($\sim$ hours after the burst event). On the other hand, the same test performed for the E$_{\rm iso}$-normalised optical light curves shows that the optical emission is a better proxy of the afterglow emission from early to late times.

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

Long-duration Gamma-Ray Bursts (GRBs) are an extremely rare outcome of the collapse of massive stars, and are typically found in the distant Universe. Because of its intrinsic luminosity ($L\sim 3 \times 10^{53}$ erg s$^{-1}$) and its relative proximity ($z=0.34$), GRB 130427A was a unique event that reached the highest fluence observed in the gamma-ray band. Here we present a comprehensive multiwavelength view of GRB 130427A with Swift, the 2-m Liverpool and Faulkes telescopes and by other ground-based facilities, highlighting the evolution of the burst emission from the prompt to the afterglow phase. The properties of GRB 130427A are similar to those of the most luminous, high-redshift GRBs, suggesting that a common central engine is responsible for producing GRBs in both the contemporary and the early Universe and over the full range of GRB isotropic energies.

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

Starting from the Swift sample we define a complete sub-sample of 58 bright long Gamma Ray Bursts (GRB), 55 of them (95%) with a redshift determination, in order to characterize their properties. Our sample (BAT6) allows us to study the properties of the long GRB population and their evolution with cosmic time. We focus in particular on the GRB luminosity function, on the spectral-energy correlations of their prompt emission, on the nature of dark bursts, on possible correlations between the prompt and the X-ray afterglow properties, and on the dust extinction.

We consider the 0.1-10 GeV rest frame light curves of 10 GRBs detected by LAT and with known redshift. In all cases the emission persists after the prompt has faded away. This extended emission decays in time as a power-law. The decay slope is similar among different bursts, while the normalization spans more than 2 decades. However, when the LAT luminosity is normalized to the 1 keV-10 MeV energetics of the prompt emission E_iso all light curves become consistent with having the same normalization, i.e. they cluster. At each given time the ratio between the LAT luminosity and the prompt energetics is narrowly distributed. We argue that this result is expected in the external shock scenario and it strengthens the interpretation of the GeV emission in terms of radiation from external shocks. In this context, we derive limits on the distribution of epsilon_e (the fraction of the shock energy that goes into electrons) and eta (the efficiency of the mechanism producing the prompt).

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

Forward shocks caused by the interaction between a relativistic blast wave and the circum-burst medium are thought to be responsible for the afterglow emission in Gamma-Ray Bursts (GRBs). We consider the hydrodynamics of a spherical relativistic blast wave expanding into the surrounding medium and we generalize the standard theory in order to account for several effects that are generally ignored. In particular, we consider the role of adiabatic and radiative losses on the hydrodynamical evolution of the shock, under the assumption that the cooling losses are fast. Our model can describe adiabatic, fully radiative and semi-radiative blast waves, and can describe the effects of a time-varying radiative efficiency. The equations we present are valid for arbitrary density profiles, and also for a circum-burst medium enriched with electron-positron pairs. The presence of pairs enhances the fraction of shock energy gained by the leptons, thus increasing the importance of radiative losses. Our model allows to study whether the high-energy (>0.1 GeV) emission in GRBs may originate from afterglow radiation. In particular, it is suitable to test whether the fast decay of the high-energy light curve observed in several Fermi LAT GRBs can be ascribed to an initial radiative phase, followed by the standard adiabatic evolution.

In this paper we compute rest-frame extinctions for the afterglows of a sample of gamma-ray bursts complete in redshift. The selection criteria of the sample are based on observational high-energy parameters of the prompt emission and therefore our sample should not be biased against dusty sight-lines. It is therefore expected that our inferences hold for the general population of gamma-ray bursts. Our main result is that the optical/near-infrared extinction of gamma-ray burst afterglows in our sample does not follow a single distribution. 87% of the events are absorbed by less than 2 mag, and 50% suffer from less than 0.3-0.4 mag extinction. The remaining 13% of the afterglows are highly absorbed. The true percentage of gamma-ray burst afterglows showing high absorption could be even higher since a fair fraction of the events without reliable redshift measurement are probably part of this class. These events may be due to highly dusty molecular clouds/star forming regions associated with the gamma-ray burst progenitor or along the afterglow line of sight, and/or to massive dusty host galaxies. No clear evolution in the dust extinction properties is evident within the redshift range of our sample, although the largest extinctions are at z~1.5-2, close to the expected peak of the star formation rate. Those events classified as dark are characterized, on average, by a higher extinction than typical events in the sample. A correlation between optical/near-infrared extinction and hydrogen-equivalent column density based on X-ray studies is shown although the observed NH appears to be well in excess compared to those observed in the Local Group. Dust extinction does not seem to correlate with GRB energetics or luminosity.

We recently found that Gamma Ray Burst energies and luminosities, in their comoving frame, are remarkably similar. This, coupled with the clustering of energetics once corrected for the collimation factor, suggests the possibility that all bursts, in their comoving frame, have the same peak energy E'peak (of the order of a few keV) and the same energetics of the prompt emission E'gamma (of the order of 2e48 erg). The large diversity of bursts energies is then due to the different bulk Lorentz factor Gamma and jet aperture angle theta_jet. We investigated, through a population synthesis code, what are the distributions of Gamma and theta_jet compatible with the observations. Both quantities must have preferred values, with log-normal best fitting distributions and <Gamma0> ~ 275 and <\theta_jet> ~ 8.7 degree. Moreover, the peak values of the Gamma and theta_jet distributions must be related - theta_jet^2.5 Gamma =const: the narrower the jet angle, the larger the bulk Lorentz factor. 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. Finally, we estimate that the local rate of GRBs is ~0.3% of all local SNIb/c and ~2.5% of local hypernovae, i.e. SNIb/c with broad absorption lines.

Cosmic rays that escape their acceleration site interact with the ambient medium and produce gamma rays as the result of inelastic proton-proton collisions. The detection of such diffuse emission may reveal the presence of an accelerator of cosmic rays, and also constrain the cosmic ray diffusion coefficient in its vicinity. Preliminary results in this direction have been obtained in the last years from studies of the gamma-ray emission from molecular clouds located in the vicinity of supernova remnants, which are the prime candidate for cosmic ray production. Hints have been found for a significant suppression of the diffusion coefficient with respect to the average one in the Galaxy. However, most of these studies rely on the assumption of isotropic diffusion, which may not be very well justified. Here, we extend this study to the case in which cosmic rays that escape an accelerator diffuse preferentially along the magnetic field lines. As a first approximation, we further assume that particles are strongly magnetized and that their transport perpendicular to the magnetic field is mainly due to the wandering of the field lines. The resulting spatial distribution of runaway cosmic rays around the accelerator is, in this case, strongly anisotropic. An application of the model to the case of the supernova remnant W28 demonstrates how the estimates of the diffusion coefficient from gamma-ray observations strongly depend on the assumptions made on the isotropy (or anisotropy) of diffusion. For higher levels of anisotropy of the diffusion, larger values of the diffusion coefficient are found to provide a good fit to data. Thus, detailed models for the propagation of cosmic rays are needed in order to interpret in a correct way the gamma-ray observations.

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 investigate wheter there is any correlation between the X-ray afterglow luminosity and the prompt emission properties of a carefully selected sub-sample of bright Swift long Gamma-Ray Bursts (GRBs) nearly complete in redshift (~90%). Being free of selection effects (except flux limit), this sample provides the possibility to compare the rest frame physical properties of GRB prompt and afterglow emission in an unbiased way. The afterglow X-ray luminosities are computed at four different rest frame times (5 min, 1 hr, 11 hr and 24 hr after trigger) and compared with the prompt emission isotropic energy E_iso, the isotropic peak luminosity L_iso and the rest frame peak energy E_peak. We find that the rest frame afterglow X-ray luminosity do correlate with these prompt emission quantities, but the significance of each correlation decreases over time. This result is in agreement with the idea that the GRB X-ray light curve can be described as the result of a combination of different components whose relative contribution and weight change with time, with the prompt and afterglow emission dominating at early and late time, respectively. In particular, we found evidence that the plateau and the shallow decay phase often observed in GRB X-ray light curves are powered by activity from the central engine. The existence of the L_X-E_iso correlation at late times (t_rf > 11 hr) suggests a similar radiative efficiency among different bursts with on average about 6% of the total kinetic energy powering the prompt emission.

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 Monte 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 present a carefully selected sub-sample of Swift Long Gamma-ray Bursts (GRBs), that is complete in redshift. The sample is constructed by considering only bursts with favorable observing conditions for ground-based follow-up searches, that are bright in the 15-150 keV Swift/BAT band, i.e. with 1-s peak photon fluxes in excess to 2.6 ph s^-1 cm^-2. The sample is composed by 58 bursts, 52 of them with redshift for a completeness level of 90%, while another two have a redshift constraint, reaching a completeness level of 95%. For only three bursts we have no constraint on the redshift. The high level of redshift completeness allows us for the first time to constrain the GRB luminosity function and its evolution with cosmic times in a unbiased way. We find that strong evolution in luminosity (d_l=2.3\pm 0.6) or in density (d_d=1.7\pm 0.5) is required in order to account for the observations. The derived redshift distribution in the two scenarios are consistent with each other, in spite of their different intrinsic redshift distribution. This calls for other indicators to distinguish among different evolution models. Complete samples are at the base of any population studies. In future works we will use this unique sample of Swift bright GRBs to study the properties of the population of long GRBs.

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 properties of the population of optically dark events present in a carefully selected complete sample of bright Swift long gamma-ray bursts. The high level of completeness in redshift of our sample (52 objects out of 58) allow us to establish the existence of a genuine dark population and we are able to estimate the maximum fraction of dark burst events (~30%) expected for the whole class of long gamma-ray burst. The redshift distribution of this population of dark bursts is similar to the one of the whole sample. Interestingly, the rest-frame X-ray luminosity (and the de-absorbed X-ray flux) of the sub-class of dark bursts is slightly higher than the average luminosity of the non-dark events. At the same time the prompt properties do not differ and the optical flux of dark events is at the lower tail of the optical flux distribution, corrected for Galactic absorption. All these properties suggest that dark bursts events generate in much denser environments with respect to normal bright events. We can therefore exclude the high-z and the low-density scenarios and conclude that the major cause of the origin of optically dark events is the dust extinction.

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

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 present the results of the spectral analysis of the public data of 438 Gamma Ray Bursts (GRBs) detected by the Fermi Gamma ray Burst Monitor (GBM) up to March 2010. For 432 bursts we could fit the time integrated spectrum. In 318 cases we can reliably constrain the peak energy Epeak of their \nu F_\nu spectrum by analyzing their time integrated spectrum between 8 keV and 35 MeV. 80% of these spectra are fitted by a power law with an exponential cutoff, and the remaining with the Band function. Among these 318 GRBs, 274 and 44 belong to the long and short GRB class, respectively. Long GRBs have a typical peak energy Epeak=160 keV and low energy spectral index alpha=-0.92. Short GRBs have harder peak energy (Epeak=490 keV) and harder low energy spectral index (alpha=-0.50) than long bursts. For each Fermi GRB we analyzed also the spectrum corresponding to the peak flux of the burst. On average, the peak spectrum has harder low energy spectral index but similar Epeak than the corresponding time-integrated spectrum for the same burst. The spectral parameters derived in our analysis of Fermi/GBM bursts are globally consistent with those reported in the GRB Cicular Network (GCN) archive after December 2008, while we found systematic differences, concerning the low energy power law index, for earlier bursts.

The Gamma-ray Burst Monitor (GBM) on board Fermi allows to study the spectra of Gamma Ray Bursts (GRBs) over an unprecedented wide energy range (8 keV - 35 MeV). We compare the spectral properties of short and long GRBs detected by the GBM (up to March 2010) with those of GRBs detected by the BATSE instrument on board the CGRO. GBM and BATSE long bursts have similar distributions of fluence (F), Epeak and peak flux (P) but GBM bursts have a slightly harder low-energy spectral index \alpha with respect to BATSE GRBs. GBM and BATSE short bursts have similar distributions of fluence, \alpha and peak flux, with GBM bursts having slightly larger Epeak. We discuss these properties in light of the found correlations between Epeak and the fluence and the peak flux. GBM bursts confirm that these correlations are not determined by instrumental selection effects. Indeed, GBM bursts extend the Epeak-Fluence and Epeak-P correlations both in fluence/peak flux and in peak energy. No GBM long burst with Epeak exceeding a few MeV is found, despite the possibility of detecting it. Similarly to what found with BATSE, there are 3% of GBM long bursts (and almost all short ones) that are outliers at more than 3\sigma of the Epeak-Eiso correlation. Instead there is no outlier of the Epeak-Lp,iso correlation, for both long and short GBM bursts.

We study the spectral evolution of 13 short duration Gamma Ray Bursts (GRBs) detected by the Gamma Burst Monitor (GBM) on board Fermi. We study spectra resolved in time at the level of 2-512 ms in the 8 keV-35 MeV energy range. We find a strong correlation between the observed peak energy Ep and the flux P within individual short GRBs. The slope of the Ep P^s correlation for individual bursts ranges between ~0.4 and ~1. There is no correlation between the low energy spectral index and the peak energy or the flux. Our results show that in our 13 short GRBs Ep evolves in time tracking the flux. This behavior is similar to what found in the population of long GRBs and it is in agreement with the evidence that long GRBs and (the still few) short GRBs with measured redshifts follow the same rest frame Ep-Liso correlation. Its origin is most likely to be found in the radiative mechanism that has to be the same in both classes of GRBs.

We compare the spectral properties of 227 Gamma Ray Bursts (GRBs) detected by the Fermi Gamma Ray Burst Monitor (GBM) up to February 2010 with those of bursts detected by the CGRO/BATSE instrument. Out of 227 Fermi GRBs, 166 have a measured peak energy E_peak_obs of their \nuF(\nu) spectrum: of these 146 and 20 belong the long and short class, respectively. Fermi long bursts follow the correlations defined by BATSE bursts between their E_peak_obs vs fluence and peak flux: as already shown for the latter ones, these correlations and their slopes do not originate from instrumental selection effects. Fermi/GBM bursts extend such correlations toward lower fluence/peak energy values with respect to BATSE ones whereas no GBM long burst with E_peak_obs exceeding a few MeV is found, despite the possibility of detecting them. Again as for BATSE, $\sim$ 5% of long and almost all short GRBs detected by Fermi/GBM are outliers of the E_peak-isotropic equivalent energy ("Amati") correlation while no outlier (neither long nor short) of the E_peak-isotropic equivalent luminosity ("Yonetoku") correlation is found. Fermi long bursts have similar typical values of E_peak_obs but a harder low energy spectral index with respect to all BATSE events, exacerbating the inconsistency with the limiting slopes of the simplest synchrotron emission models. Although the short GRBs detected by Fermi are still only a few, we confirm that their E_peak_obs is greater and the low energy spectrum is harder than those of long ones. We discuss the robustness of these results with respect to observational biases induced by the differences between the GBM and BATSE instruments.

We study the emission observed at energies greater than 100 MeV of 11 Gamma Ray Bursts (GRBs) detected by the Fermi/Large Area Telescope (LAT) until October 2009. The GeV emission has three main properties: (i) its duration is often longer than the duration of the softer emission detected by the Gamma Burst Monitor (GBM) onboard Fermi [this confirms earlier results from the Energetic Gamma-Ray Experiment Telescope (EGRET)]; (ii) its spectrum is consistent with F(v) propto v^(-1) and does not show strong spectral evolution; (iii) for the brightest bursts, the flux detected by the LAT decays as a power law with a typical slope: t^(-1.5). We argue that the observed >0.1 GeV flux can be interpreted as afterglow emission shortly following the start of the prompt phase emission as seen at smaller frequencies. The decay slope is what expected if the fireball emission is produced in the radiative regime, i.e. all dissipated energy is radiated away. We also argue that the detectability in the GeV energy range depends on the bulk Lorentz factor Gamma of the bursts, being strongly favoured in the case of large Gamma. This implies that the fraction of bursts detected at high energies corresponds to the fraction of bursts having the largest Gamma. The radiative interpretation can help to explain why the observed X-ray and optical afterglow energetics are much smaller than the energetics emitted during the prompt phase, despite the fact that the collision with the external medium should be more efficient than internal shocks in producing the radiation we see.

We study the emission of the short/hard GRB 090510 at energies > 0.1 GeV as observed by the Large Area Telescope (LAT) onboard the Fermi satellite. The GeV flux rises in time as t^2 and decays as t^-1.5 up to 200 s. The peak of the high energy flux is delayed by 0.2 s with respect to the main ~MeV pulse detected by the Fermi Gamma Burst Monitor (GBM). Its energy spectrum is consistent with F(E)=E^-1. The time behavior and the spectrum of the high energy LAT flux are strong evidences of an afterglow origin. We then interpret it as synchrotron radiation produced by the forward shock of a fireball having a bulk Lorentz factor Gamma ~ 2000. The afterglow peak time is independent of energy in the 0.1-30 GeV range and coincides with the arrival time of the highest energy photon (~ 30 GeV). Since the flux detected by the GBM and the LAT have different origins, the delay between these two components is not entirely due to possible violation of the Lorentz invariance. It is the LAT component by itself that allows to set a stringent lower limit on the quantum-gravity mass of 4.7 times the Planck mass.

We study the spectra of all long Gamma Ray Bursts (GRBs) of known redshift detected by the Fermi satellite. Their fluxes and fluences are large enough to allow a time dependent study of their spectral characteristics in the 8 keV-1 MeV energy range. We find that the peak energy Ep of their EL(E) spectrum correlates with the luminosity in a remarkable tight way within individual bursts. This time resolved Ep-Liso correlation is very similar for all the considered bursts, and has a slope and normalisation similar to the analogous Ep-Liso correlation defined by the time integrated spectra of different bursts detected by several different satellites. For a few of the considered GRBs, we could also study the behaviour of the Ep-Liso correlation during the rising and decaying phases of individual pulses within each burst, finding no differences. Our results indicate the presence of a similar physical mechanism, operating for the duration of different GRBs, linking tightly the burst luminosity with the peak energy of the spectrum emitted at different times. Such a physical mechanism is the same during the rise and decay phase of individual pulses composing a GRB. These results, while calling for a robust physical interpretation, strongly indicate that the Ep-Liso spectral energy correlation found considering the time integrated spectra of different bursts is real, and not the result of instrumental selection effects.

Instrumental selection effects can act upon the estimates of the peak energy Ep, the fluence F and the peak flux P of GRBs. If this were the case, then the correlations involving the corresponding rest frame quantities (i.e. Ep, Eiso and the peak luminosity Liso) would be questioned. We estimated, as a function of Ep, the minimum peak flux necessary to trigger a GRB and the minimum fluence a burst must have to determine the value of Ep by considering different instruments (BATSE, Swift, BeppoSAX). We find that the latter dominates over the former. We then study the Ep-fluence (and flux) correlation in the observer plane. GRBs with redshift show well defined Ep-F and Ep-P correlations: in this planes the selection effects are present, but do not determine the found orrelations. This is not true for Swift GRBs with redshift, for which the spectral analysis threshold does affect their distribution in the observer planes. Extending the sample to GRBs without z, we still find a significant Ep-F correlation, although with a larger scatter than that defined by GRBs with redshift. We find that 6% are outliers of the Amati correlation. The Ep-P correlation of GRBs with or without redshift is the same and no outlier is found among bursts without redshift.

We compare the spectral properties of 79 short and 79 long Gamma-Ray Bursts (GRBs) detected by BATSE and selected with the same limiting peak flux. Short GRBs have a low-energy spectral component harder and a peak energy slightly higher than long GRBs, but no difference is found when comparing short GRB spectra with those of the first 1-2 sec emission of long GRBs. These results confirm earlier findings for brighter GRBs. The bolometric peak flux of short GRBs correlates with their peak energy in a similar way to long bursts. Short and long GRBs populate different regions of the bolometric fluence-peak energy plane, short bursts being less energetic by a factor similar to the ratio of their durations. If short and long GRBs had similar redshift distributions, they would have similar luminosities yet different energies, which correlate with the peak energy E_peak for the population of long GRBs. We also test whether short GRBs are consistent with the E_peak-E_iso and E_peak-L_iso correlations for the available sample of short (6 events) and long (92 events) GRBs with measured redshifts and E_peak,obs: while short GRBs are inconsistent with the E_peak-E_iso correlation of long GRBs, they could follow the E_peak-L_iso correlation of long bursts. All the above indications point to short GRBs being similar to the first phases of long bursts. This suggests that a similar central engine (except for its duration) operates in GRBs of different durations.

For a sample of long GRBs with known redshift, we study the distribution of the evolutionary tracks on the rest-frame luminosity-peak energy Liso-Ep' diagram. We are interested in exploring the extension of the `Yonetoku' correlation to any phase of the prompt light curve, and in verifying how the high-signal prompt duration time, Tf, in the rest frame correlates with the residuals of such correlation (Firmani et al. 2006). For our purpose, we analyse separately two samples of time-resolved spectra corresponding to 32 GRBs with peak fluxes >1.8 phot cm^-2 s^-1 from the Swift-BAT detector, and 7 bright GRBs from the CGRO-BATSE detector previously processed by Kaneko et al. (2006). After constructing the Liso-Ep' diagram, we discuss the relevance of selection effects, finding that they could affect significantly the correlation. However, we find that these effects are much less significant in the Liso x Tf-Ep' diagram, where the intrinsic scatter reduces significantly. We apply further corrections for reducing the intrinsic scatter even more. For the sub-samples of GRBs (7 from Swift and 5 from CGRO) with measured jet break time, we analyse the effects of correcting Liso by jet collimation. We find that (i) the scatter around the correlation is reduced, and (ii) this scatter is dominated by the internal scatter of the individual evolutionary tracks. These results suggest that the time, integrated `Amati' and `Ghirlanda' correlations are consequences of the time resolved features, not of selection effects, and therefore call for a physical origin. We finally remark the relevance of looking inside the nature of the evolutionary tracks.

We study a sample of Gamma-Ray Bursts detected by the Swift satellite with known redshift which show a precursor in the Swift-BAT light curve. We analyze the spectra of the precursors and compare them with the time integrated spectra of the prompt emission. We find neither a correlation between the two slopes nor a tendency for the precursors spectra to be systematically harder or softer than the prompt ones. The energetics of the precursors are large: on average, they are just a factor of a few less energetic (in the source rest frame energy range 15-150 keV) than the entire bursts. These properties do not depend upon the quiescent time between the end of the precursor and the start of the main event. These results suggest that what has been called a "precursor" is not a phenomenon distinct from the main event, but is tightly connected with it, even if, in some case, the quiescent time intervals can be longer than 100 seconds.

The spectral-energy and (luminosity) correlations in long GRBs are being hotly debated to establish, first of all, their reality against possible selection effects. These are best studied in the observer planes, namely the peak energy E_peak_obs vs the fluence F or the peak flux P. In a recent paper we started to attack this problem considering all GRBs with known z and spectral properties. Here we consider instead all bursts with known E_peak_obs, irrespective of z, adding to those a sample of 100 faint BATSE bursts representative of a larger population. This allows us to construct a complete, fluence limited, sample, to study the selection/instrumental effects. We found that fainter bursts have smaller E_peak_obs than those of bright events. As a consequence, the E_peak_obs of these bursts is correlated with the fluence, though with a slope flatter than that defined by bursts with z. Selection effects, which are present, are shown not to be responsible for the existence of such a correlation. About 6% of these bursts are surely outliers of the E_peak-E_iso correlation (updated to include 83 bursts), since they are inconsistent with it for any z. E_peak_obs correlates also with P, with a slope similar to the E_peak-L_iso correlation.In this case there is only one sure outlier.The scatter of the E_peak_obs-P correlation defined by the BATSE bursts of our sample is smaller than the E_peak_obs-F correlation of the same bursts, while for the bursts with known z the E_peak-E_iso correlation is tighter than the E_peak-L_iso one. Once a very large number of bursts with E_peak_obs and z will be available, we thus expect that the E_peak-L_iso correlation will be similar to that currently found, whereas it is likely that the E_peak-E_iso correlation will become flatter and with a larger scatter.

We study the distribution of long Gamma Ray Bursts in the Ep-Eiso and in the Ep,obs-Fluence planes through an updated sample of 76 bursts, with measured redshift and spectral parameters, detected up to September 2007. We confirm the existence of a strong rest frame correlation Ep ~ Eiso^0.54+-0.01. Contrary to previous studies, no sign of evolution with redshift of the Ep-Eiso correlation (either its slope and normalisation) is found. The 76 bursts define a strong Ep,obs-Fluence correlation in the observer frame (Ep,obs ~ F^0.32+-0.05) with redshifts evenly distributed along this correlation. We study possible instrumental selection effects in the observer frame Ep,obs-Fluence plane. In particular, we concentrate on the minimum peak flux necessary to trigger a given GRB detector (trigger threshold) and the minimum fluence a burst must have to determine the value of Ep,obs (spectral analysis threshold). We find that the latter dominates in the Ep,obs-Fluence plane over the former. Our analysis shows, however, that these instrumental selection effects do not dominate for bursts detected before the launch of the Swift satellite, while the spectral analysis threshold is the dominant truncation effect of the Swift GRB sample (27 out of 76 events). This suggests that the Ep,obs-Fluence correlation defined by the pre--Swift sample could be affected by other, still not understood, selection effects. Besides we caution about the conclusions on the existence of the Ep,obs-Fluence correlation based on our Swift sample alone.

We consider the interaction between a relativistic fireball and material assumed to be still located just outside the progenitor star. Only a small fraction of the expected mass is sufficient to efficiently decelerate the fireball, leading to dissipation of most of its kinetic energy. Since the scattering optical depths are still large at distances comparable to the progenitor radius, the dissipated energy is trapped in the system, accelerating it to relativistic velocities. The process resembles the birth of another fireball at radii R~1e11 cm, not far from the transparency radius, and with a starting bulk Lorentz factors Gamma_c~10. As seen in the observer frame, this "re--generated" fireball appears collimated within an angle theta_j=1/Gamma_c. If the central engine works intermittently, the funnel can, at least partially, refill and the process can repeat itself. We discuss how this idea can help solving some open issues of the more conventional internal shock scenario for interpreting the Gamma-Ray Burst properties.

We study the spectral and energetics properties of 47 long-duration gamma-ray bursts (GRBs) with known redshift, all of them detected by the Swift satellite. Due to the narrow energy range (15-150 keV) of the Swift-BAT detector, the spectral fitting is reliable only for fitting models with 2 or 3 parameters. As high uncertainty and correlation among the errors is expected, a careful analysis of the errors is necessary. We fit both the power law (PL, 2 parameters) and cut--off power law (CPL, 3 parameters) models to the time-integrated spectra of the 47 bursts, and present the corresponding parameters, their uncertainties, and the correlations among the uncertainties. The CPL model is reliable only for 29 bursts for which we estimate the nuf_nu peak energy Epk. For these GRBs, we calculate the energy fluence and the rest- frame isotropic-equivalent radiated energy, Eiso, as well as the propagated uncertainties and correlations among them. We explore the distribution of our homogeneous sample of GRBs on the rest-frame diagram E'pk vs Eiso. We confirm a significant correlation between these two quantities (the "Amati" relation) and we verify that, within the uncertainty limits, no outliers are present. We also fit the spectra to a Band model with the high energy power law index frozen to -2.3, obtaining a rather good agreement with the "Amati" relation of non-Swift GRBs.

In their recent paper, Campana et al. (2007) found that 5 bursts, among those detected by Swift, are outliers with respect to the E_peak-E_gamma ("Ghirlanda") correlation. We instead argue that they are not.

The early X-ray light curve of Gamma Ray Bursts (GRBs) is complex, and shows a typical steep-flat-steep behaviour. The time T_a at which the flat (plateau) part ends may bear some important physical information, especially if it plays the same role of the so called jet break time t_jet. To this aim, stimulated by the recent analysis of Willingale et al., we have assembled a sample of GRBs of known redshifts, spectral parameters of the prompt emission, and T_a. By using T_a as a jet angle indicator, and then estimating the collimation corrected prompt energetics, we find a correlation between the latter quantity and the peak energy of the prompt emission. However, this correlation has a large dispersion, similar to the dispersion of the Amati correlation and it is not parallel to the Ghirlanda correlation. Furthermore, we show that the correlation itself results mainly from the dependence of the jet opening angle on the isotropic prompt energy, with the time T_a playing no role, contrary to what we find for the jet break time t_jet. We also find that for the bursts in our sample T_a weakly correlates with E_iso of the prompt emission, but that this correlation disappears when considering all bursts of known redshift and T_a. There is no correlation between T_a and the isotropic energy of the plateau phase.

The flat decay phase in the first 1e2-1e4 seconds of the X-ray light curve of Gamma Ray Bursts (GRBs) has not yet found a convincing explanation. The fact that the optical and X-ray lightcurves are often different, with breaks at different times, makes contrived any explanation based on the same origin for both the X-ray and optical fluxes. We here assume that the central engine can be active for a long time, producing shells of decreasing bulk Lorentz factors Gamma. We also assume that the internal dissipation of these late shells produces a continuous and smooth emission (power-law in time), usually dominant in X-rays and sometimes in the optical. When Gamma of the late shells is larger than 1/theta_j, where theta_j is the jet opening angle, we see only a portion of the emitting surface. Eventually, Gamma becomes smaller than 1/theta_j, and the entire emitting surface is visible. Thus there is a break in the light curve when Gamma=1/theta_j, which we associate to the time at which the plateau ends. After the steeply decaying phase which follows the early prompt, we see the sum of two emission components: the "late-prompt" emission (due to late internal dissipation), and the "real afterglow" emission (due to external shocks). A variety of different optical and X-ray light curves are then possible, explaining why the X-ray and the optical light curves often do not track each other (but sometimes do), and often they do not have simultaneous breaks.

We test the spectral-energy correlation including the new bursts detected (mostly) by Swift with firm measurements of their redshifts and peak energy. The problem of identifying the jet breaks is discussed in the complex and multibreak/flaring X-ray light curves observed by Swift. We use the optical data as the most reliable source for the identification of the jet break, since the X-ray flux may be produced by a mechanism different from the external shocks between the fireball and the circumburst medium, which are responsible for the optical afterglow. We show that the presence of an underlying SN event in XRF 050416A requires a break to occur in the afterglow optical light curve at around the expected jet break time. The possible presence of a jet break in the optical light curve of GRB 050401 is also discussed. We point out that, for measuring the jet break, it is mandatory that the optical light curve extends after the epoch where the jet break is expected. The interpretation of the early optical breaks in GRB 050922C and GRB 060206 as jet breaks is controversial because they might instead correspond to the flat-to-steep decay transition common in the early X-ray light curves. All the 16 bursts coming from Swift are consistent with the E_p-E_gamma and E_p-E_iso-t_jet correlation. No outlier is found to date. Moreover, the small dispersion of these correlations, confirmed also by the Swift bursts, strengthens the case of using GRBs as standard candles.

Ghirlanda et al. (2004) derived the collimation-corrected energy E_gamma for a sample of 15 bursts under the assumption of a homogeneous circumburst density profile. They found a correlation (the so-called Ghirlanda correlation) between E_gamma and the rest frame peak energy of the nuF_nu prompt spectrum (E_p). Nava et al. (2006) showed that, assuming a circumburst density distribuited with a r^-2 wind profile, the Ghirlanda correlation remains tight and becomes linear. This implies that: i) it remains linear also in the comoving frame, no matter the distribution of bulk Lorentz factors, ii) it entails that different bursts have the same number of relevant photons. We have updated these findings including recently detected bursts (21 in total), stressing the two important implications.

As soon as it was realized that long GRBs lie at cosmological distances, attempts have been made to use them as cosmological probes. Besides their use as lighthouses, a task that presents mainly the technological challenge of a rapid deep high resolution follow-up, researchers attempted to find the Holy Grail: a way to create a standard candle from GRB observables. We discuss here the attempts and the discovery of the Ghirlanda correlation, to date the best method to standardize the GRB candle. Together with discussing the promises of this method, we will underline the open issues, the required calibrations and how to understand them and keep them under control. Even though GRB cosmology is a field in its infancy, ongoing work and studies will clarify soon if and how GRBs will be able to keep up to the promises.

The use of Gamma Ray Bursts as ``standard candles'' has been made possible by the recent discovery of a very tight correlation between their rest frame intrinsic properties. This correlation relates the GRB prompt emission peak spectral energy E_peak to the energy E_gamma corrected for the collimation angle theta_jet of these sources. The possibility to use GRBs to constrain the cosmological parameters and to study the nature of Dark Energy are very promising.

We present the constraints on the cosmological parameters obtained with the $E_{\rm peak}$--$E_{\gamma}$ correlation found with the most recent sample of 19 GRBs with spectroscopically measured redshift and well determined prompt emission spectral and afterglow parameters. We compare our results obtained in the two possible uniform jet scenarios, i.e. assuming a homogeneous density profile (HM) or a wind density profile (WM) for the circumburst medium. Better constraints on $\Omega_{M}$ and $\Omega_{\Lambda}$ are obtained with the (tighter) $E_{\rm peak}$--$E_{\gamma}$ correlation derived in the wind density scenario. We explore the improvements to the constraints of the cosmological parameters that could be reached with a large sample, $\sim$ 150 GRBs, in the future. We study the possibility to calibrate the slope of these correlations. Our optimization analysis suggests that $\sim 12$ GRBs with redshift $z\in(0.9,1.1)$ can be used to calibrate the $E_{\rm peak}$--$E_{\gamma}$ with a precision better than 1%. The same precision is expected for the same number of bursts with $z\in(0.45,0.75)$. This result suggests that we do not necessarily need a large sample of low z GRBs for calibrating the slope of these correlations.

Recently, Liang & Zhang (2005) found a tight correlation involving only observable quantities, namely the isotropic emitted energy $E_{\gamma,iso}$, the energy of the peak of the prompt spectrum $E^\prime_{p}$, and the jet break time $t^\prime_{j}$. This phenomenological correlation can have a first explanation in the framework of jetted fireballs, whose semiaperture angle $\theta_{j}$ is measured by the jet break time $t^\prime_{j}$. By correcting $E_{\gamma, iso}$ for the angle $\theta_{j}$ one obtains the so called Ghirlanda correlation linking the collimation corrected energy $E_\gamma$ and $E^\prime_{p}$. There are two ways to derive $\theta_{j}$ from $t^\prime_{j}$ in the standard scenario, corresponding to an homogeneous or to a wind-like circumburst medium. We show that the Ghirlanda correlation with a wind-like medium is as tight as (if not tighter) than the Ghirlanda correlation found in the case of an homogeneous medium. There are hence two Ghirlanda correlations, both entirely consistent with the phenomenological Liang & Zhang relation. We consider the difference between the observed correlations and the ones one would see in the comoving frame (i.e. moving with the same bulk Lorentz factor of the fireball). Since both $E_{p}$ and $E_\gamma$ transform in the same way, the wind-like Ghirlanda relation, which is linear, remains linear also in the comoving frame, no matter the distribution of bulk Lorentz factors. Instead, in the homogeneous density case, one is forced to assume the existence of a strict relation between the bulk Lorentz factor and the total energy, which in turn put constraints on the radiation mechanisms of the prompt emission. The wind-like Ghirlanda correlation, being linear, corresponds to different bursts having the same number of photons.