The Fermi-LAT collaboration presented the second gamma-ray burst (GRB) catalog covering its first 10 years of operations. A significant fraction of afterglow-phase light curves in this catalog cannot be explained by the closure relations of the standard synchrotron forward-shock model, suggesting that there could be an important contribution from another process. In view of the above, we derive the synchrotron self-Compton (SSC) light curves from the reverse shock in the thick- and thin-shell regime for a uniform-density medium. We show that this emission could explain the GeV flares exhibited in some LAT light curves. Additionally, we demonstrate that the passage of the forward shock synchrotron cooling break through the LAT band from jets expanding in a uniform-density environment may be responsible for the late time ($\approx10^2$ s) steepening of LAT GRB afterglow light curves. As a particular case, we model the LAT light curve of GRB 160509A that exhibited a GeV flare together with a break in the long-lasting emission, and also two very high energy photons with energies of 51.9 and 41.5 GeV observed 76.5 and 242 s after the onset of the burst, respectively. Constraining the microphysical parameters and the circumburst density from the afterglow observations, we show that the GeV flare is consistent with a SSC reverse-shock model, the break in the long-lasting emission with the passage of the synchrotron cooling break through the Fermi-LAT band and the very energetic photons with SSC emission from the forward shock when the outflow carries a significant magnetic field ($R_{\rm B} \simeq 30$) and it decelerates in a uniform-density medium with a very low density ($n=4.554^{+1.128}_{-1.121}\times 10^{-4}\,{\rm cm^{-3}}$).

We present the afterglow light curves produced by the deceleration of the non-relativistic ejecta mass in a stratified circumstellar medium with a density profile $n(r)\propto r^{-k}$ with $k=0$, $1$, $1.5$, $2$ and $2.5$. Once the ejecta mass is launched with equivalent kinetic energy parametrized by $E(>\beta)\propto \beta^{-\alpha}$ (where beta is the ejecta velocity) and propagates into the surrounding circumstellar medium, it first moves with constant velocity (the free-coasting phase), and later it decelerates (the Sedov-Taylor expansion). We present the predicted synchrotron and synchrotron-self Compton light curves during the free-coasting phase, and the subsequent Sedov-Taylor expansion. In particular cases, we show the corresponding light curves generated by the deceleration of several ejecta masses with different velocities launched during the coalescence of binary compact objects and the core-collapse of dying massive stars which will contribute at distinct timescales, frequencies, and intensities. We find that before the radiation escapes from the kilonova (KN) and produces a peak in the light curve, the KN afterglow ejecta is very weak and the corresponding ejecta velocity moves as $\propto t^{-\frac{2}{\alpha+1}}$. Finally, using the multi-wavelength observations and upper limits collected by a large campaign of orbiting satellites and ground telescopes, we constrain the parameter space of both the KN afterglow in GW170817 and the possibly generated KN afterglow in S190814bv. Further observations on timescales of years post-merger are needed to derive tighter constraints.

Gamma-ray bursts (GRBs) are spectacularly energetic events, with the potential to inform on the early universe and its evolution, once their redshifts are known. Unfortunately, determining redshifts is a painstaking procedure requiring detailed follow-up multi-wavelength observations often involving various astronomical facilities, which have to be rapidly pointed at these serendipitous events. Here we use Machine Learning algorithms to infer redshifts from a collection of observed temporal and spectral features of GRBs. We obtained a very high correlation coefficient ($0.96$) between the inferred and the observed redshifts, and a small dispersion (with a mean square error of $0.003$) in the test set. The addition of plateau afterglow parameters improves the predictions by $61.4\%$ compared to previous results. The GRB luminosity function and cumulative density rate evolutions, obtained from predicted and observed redshift are in excellent agreement indicating that GRBs are effective distance indicators and a reliable step for the cosmic distance ladder.

We extend our study of energetic radio loud and quiet gamma-ray bursts (GRBs), suggesting these GRBs potentially come from two separate progenitor systems. We expand the sample from our previous paper (Lloyd-Ronning & Fryer, 2017) and find our results are strengthened - radio quiet GRBs have significantly shorter intrinsic prompt duration, and are also less energetic on average. However, the tenuous correlation between isotropic energy and intrinsic duration in the radio dark sample remains tenuous and is slightly weakened by adding more bursts. Interestingly, we find an anti-correlation between the intrinsic duration and redshift in the radio bright sample but not the radio dark sample, further supporting that these two samples may come from separate progenitors. We also find that very high energy (0.1 - 100 GeV) extended emission is only present in the radio loud sample. There is no significant difference between the presence of X-ray/optical plateaus or the average jet opening angles between the two samples. We explore the interpretation of these results in the context of different progenitor models. The data are consistent with the radio loud GRBs coming from a Helium-merger system and the radio quiet GRBs coming from a collapsar system, but may also reflect other dichotomies in the inner engine such as a neutron star versus black hole core.

The X-ray afterglow plateau emission observed in many Gamma-ray Bursts (GRBs) has been interpreted as either being fueled by fallback onto a newly formed black hole, or by the spin-down luminosity of an ultra-magnetized millisecond neutron star. If the latter model is assumed, GRB X-ray afterglow light curves can be analytically reproduced. We fit a sample of GRB X-ray plateaus, interestingly yielding a distribution in the magnetic field versus spin period (B-P) diagram consistent with $B\propto P^{7/6}$. This is expected from the well-established physics of the spin-up line minimum period for Galactic millisecond pulsars. The normalisation of the relation we obtain perfectly matches spin-up line predictions for the expected masses ($\sim 1 M_{\odot}$) and radii ($\sim 10 {\rm ~km}$) of newly born magnetars, and mass accretion rates consistent with GRB expectations of $10^{-4} M_{\odot}/{\rm s} <\dot{M}< 10^{-1} M_{\odot}/{\rm s}$. Short GRBs with extended emission (SEE) appear towards the high period end of the distribution, while the long GRBs (LGRBs) towards the short period end. This result is consistent with spin-up limit expectations where the total accreted mass determines the position of the neutron star in the B-P diagram. The P-B distribution for LGRBs and SEE are statistically different, further supporting the idea that the fundamental plane relation \citep{dainotti16c,Dainotti2017} is a powerful discriminant among those populations. Our conclusions are robust against suppositions regarding the GRB collimation angle and magnetar breaking index, which shifts the resulting magnetar properties parallel to the spin-up line, and strongly support a magnetar origin for GRBs presenting X-ray plateaus.

The recent discovery of the electromagnetic counterpart of the gravitational wave source GW170817, has demonstrated the huge informative power of multi-messenger observations. During the next decade the nascent field of multi-messenger astronomy will mature significantly. Around 2030, third generation gravitational wave detectors will be roughly ten times more sensitive than the current ones. At the same time, neutrino detectors currently upgrading to multi km^3 telescopes, will include a 10 km^3 facility in the Southern hemisphere that is expected to be operational around 2030. In this review, we describe the most promising high frequency gravitational wave and neutrino sources that will be detected in the next two decades. In this context, we show the important role of the Transient High Energy Sky and Early Universe Surveyor (THESEUS), a mission concept proposed to ESA by a large international collaboration in response to the call for the Cosmic Vision Programme M5 missions. THESEUS aims at providing a substantial advancement in early Universe science as well as playing a fundamental role in multi-messenger and time-domain astrophysics, operating in strong synergy with future gravitational wave and neutrino detectors as well as major ground- and space-based telescopes. This review is an extension of the THESEUS white paper (Amati et al. 2017), also in light of the discovery of GW170817/GRB170817A that was announced on October 16th, 2017.

THESEUS is a space mission concept aimed at exploiting Gamma-Ray Bursts for investigating the early Universe and at providing a substantial advancement of multi-messenger and time-domain astrophysics. These goals will be achieved through a unique combination of instruments allowing GRB and X-ray transient detection over a broad field of view (more than 1sr) with 0.5-1 arcmin localization, an energy band extending from several MeV down to 0.3 keV and high sensitivity to transient sources in the soft X-ray domain, as well as on-board prompt (few minutes) follow-up with a 0.7 m class IR telescope with both imaging and spectroscopic capabilities. THESEUS will be perfectly suited for addressing the main open issues in cosmology such as, e.g., star formation rate and metallicity evolution of the inter-stellar and intra-galactic medium up to redshift $\sim$10, signatures of Pop III stars, sources and physics of re-ionization, and the faint end of the galaxy luminosity function. In addition, it will provide unprecedented capability to monitor the X-ray variable sky, thus detecting, localizing, and identifying the electromagnetic counterparts to sources of gravitational radiation, which may be routinely detected in the late '20s / early '30s by next generation facilities like aLIGO/ aVirgo, eLISA, KAGRA, and Einstein Telescope. THESEUS will also provide powerful synergies with the next generation of multi-wavelength observatories (e.g., LSST, ELT, SKA, CTA, ATHENA).

The mechanism responsible for the prompt emission of gamma-ray bursts (GRBs) is still a debated issue. The prompt phase-related GRB correlations can allow to discriminate among the most plausible theoretical models explaining this emission. We present an overview of the observational two-parameter correlations, their physical interpretations, their use as redshift estimators and possibly as cosmological tools. The nowadays challenge is to make GRBs, the farthest stellar-scaled objects observed (up to redshift $z=9.4$), standard candles through well established and robust correlations. However, GRBs spanning several orders of magnitude in their energetics are far from being standard candles. We describe the advances in the prompt correlation research in the past decades, with particular focus paid to the discoveries in the last 20 years.

Long gamma-ray bursts (GRBs) with a plateau phase in their X-ray afterglows obeys a three-dimensional (3D) relation (Dainotti et al. 2016), between the rest-frame time at the end of the plateau, Ta, its corresponding X-ray luminosity, La, and the peak luminosity in the prompt emission, Lpeak, an extension of the two-dimensional Dainotti relation. This 3D relation identifies a GRB fundamental plane whose existence we here confirm. We extend the original analysis with X-ray data untill July 2016 gathering 183 Swift GRBs with afterglow plateaus and known redshifts. We compare several GRB categories, such as shorts with extended emission, SEE, X-ray Flashes, GRBs associated with SNe, a sample of only long-duration GRBs (132), selected from the total sample by excluding GRBs of the previous categories, and the gold sample, composed by GRBs with light curves with good data coverage and relatively flat plateaus. The relation planes for each of these categories are not statistically different from the gold fundamental plane, with the exception of the SEE, which are hence identified as a physically distinct class of objects. The gold fundamental plane has an intrinsic scatter smaller than any plane derived from the other sample categories. Thus, the distance of any particular GRB category from this plane becomes a key parameter. Additionally, we tested this 3D relations by using GRBs observed at high energy, namely the peak luminosity values derived by the Fermi-Gamma Ray Burst Monitor (GBM). The 3D relation is also confirmed for GRBs observed by the GBM, thus showing its independence from the energy range. Furthermore, we computed the several category planes with T*a as a dependent parameter obtaining for each category smaller intrinsic scatters (reaching a reduction of $24\%$ for the long GRBs). The fundamental plane is independent from several prompt and afterglow parameters.

The mechanism responsible for the afterglow emission of Gamma Ray Bursts (GRBs) and its connection to the prompt $\gamma$-ray emission is still a debated issue. Relations between intrinsic properties of the prompt or afterglow emission can help to discriminate between plausible theoretical models of GRB production. Here we present an overview of the afterglow and prompt-afterglow two parameter relations, their physical interpretations, their use as redshift estimators and as possible cosmological tools. A similar task has already been correctly achieved for Supernovae (SNe) Ia by using the peak magnitude-stretch relation, known in the literature as the Phillips relation \citep{phillips93}. The challenge today is to make GRBs, which are amongst the farthest objects ever observed, standardizable candles as the SNe Ia through well established and robust relations. Thus, the study of relations amongst the observable and physical properties of GRBs is highly relevant together with selection biases in their physical quantities. Therefore, we describe the state of the art of the existing GRB relations, their possible and debated interpretations in view of the current theoretical models and how relations are corrected for selection biases. We conclude that only after an appropriate evaluation and correction for selection effects can GRB relations be used to discriminate among the theoretical models responsible for the prompt and afterglow emission and to estimate cosmological parameters.

The prompt emission mechanism of gamma-ray bursts (GRB) even after several decades remains a mystery. However, it is believed that correlations between observable GRB properties given their huge luminosity/radiated energy and redshift distribution extending up to at least $z \approx 9$, are promising possible cosmological tools. They also may help to discriminate among the most plausible theoretical models. Nowadays, the objective is to make GRBs standard candles, similar to supernovae (SNe) Ia, through well-established and robust correlations. However, differently from SNe Ia, GRBs span over several order of magnitude in their energetics and hence cannot yet be considered standard candles. Additionally, being observed at very large distances their physical properties are affected by selection biases, the so called Malmquist bias or Eddington effect. We describe the state of the art on how GRB prompt correlations are corrected for these selection biases in order to employ them as redshift estimators and cosmological tools. We stress that only after an appropriate evaluation and correction for these effects, the GRB correlations can be used to discriminate among the theoretical models of prompt emission, to estimate the cosmological parameters and to serve as distance indicators via redshift estimation.

The analysis of 176 gamma ray burst (GRB) afterglow plateaus observed by Swift from GRBs with known redshifts revealed that the subsample of long GRBs associated with supernovae (LONG-SNe) - 19 events - presents a very high correlation coefficient between the luminosity at the end of the plateau phase La and the end time of the plateau T*a, hereafter Dainotti relation. Furthermore, these SNe Ib/c associated with GRBs also obey the peak-magnitude stretch relation, similar to that used to standardize the SNe Ia. We here investigate a category of GRBs with plateau and associated with SNe to compare the Dainotti correlation for this sample with the correlation for long GRBs for which no associated SN has been observed (hereafter LONG-NO-SNe, 128 GRBs) and to check whether there is a difference among these sub-samples. We first adopted a non-parametric statistical method to take redshift evolution into account and to check if and how this effect may steepen the slope for the LONG-NO-SNe sample. Therefore, removing selection bias is the first step before any comparison among samples observed at different redshifts could be properly performed. Then, we applied the T-student test to evaluate a statistical difference between the slopes of the two samples. We demonstrate that there is no evolution for the slope of the LONG-NO-SNe sample and no evolution is expected for the LONG-SNe sample at small redshifts. The difference between the slope of the LONG-NO-SNe and the slope of LONG-SNe with firm spectral detection of SN components, is statistically significant. This possibly suggests that, unlike LONG-NO-SNe, LONG-SNe with firm spectroscopic features of the associated SNe might not require a standard energy reservoir in the plateau phase. Therefore, this analysis may open new perspectives in future theoretical investigations of the GRBs with plateau emission and that are associated with SNe.

A class of long gamma-ray bursts (GRBs) with a plateau phase in their X-ray afterglows obeys a three-dimensional (3D) relation (Dainotti et al. 2016), between the rest-frame time at the end of the plateau, Ta, its corresponding X-ray luminosity, La, and the peak luminosity in the prompt emission, Lpeak. We extended the original analysis with X-ray data from July 2014 to July 2016 achieving a total sample of 183 Swift GRBs with afterglow plateaus and known redshifts. We added the most recent GRBs to the previous 'gold sample' (now including 45 GRBs) and obtained a relation plane with intrinsic scatter compatible within one sigma with the previous result. We compared several GRB categories, such as short with extended emission, X-ray Flashes, GRBs associated with SNe, long-duration GRBs, and the gold sample, composed only by GRBs with light curves with good data coverage and relatively flat plateaus and evaluated their relation planes. We found that they are not statistically different from the fundamental plane derived from the gold sample and that the fundamental plane still has the smallest scatter. We compared the jet opening angles tabulated in literature with the angles derived using the Eiso-Egamma relation of the method in Pescalli et al. (2015) and calculated the relation plane for a sample of long GRBs accounting for the different jet opening angles. We observed that this correction does not significantly reduce the scatter. In an extended analysis, we found that the fundamental plane is independent from several prompt and afterglow parameters.

In this work we study the distribution of temporal power-law decay indices, $\alpha$, in the Gamma Ray Burst (GRB) afterglow phase, fitted for $176$ GRBs (139 long GRBs, 12 short GRBs {\it with extended emission} and 25 X-Ray Flashes (XRFs)) with known redshifts. These indices are compared with the temporal decay index, $\alpha_W$, derived with the light curve fitting using the \cite{willingale07} model. This model fitting yields similar distributions of $\alpha_W$ to the fitted $\alpha$, but for individual bursts a difference can be significant. Analysis of ($\alpha$, $L_a$) distribution, where $L_a$ is the characteristic luminosity at the end of the plateau, reveals only a weak correlation of these quantities. However, we discovered a significant regular trend when studying GRB $\alpha$ values along the \cite{dainotti2008} correlation between $L_a$ and the end time of the plateau emission in the rest frame, $T_a^*$, hereafter LT correlation. We note a systematic variation of the $\alpha$ parameter distribution with luminosity for any selected $T_a^*$. We analyze this systematics with respect to the fitted LT correlation line, expecting that the presented trend may allow to constrain the GRB physical models. We also attempted to use the derived correlation of $\alpha(T_a)$ versus $L_a(T_a)$ to diminish the luminosity scatter related to the variations of $\alpha$ along the LT distribution, a step forward in the effort of standardizing GRBs. A proposed toy model accounting for this systematics applied to the analyzed GRB distribution results in a slight increase of the LT correlation coefficient.

A class of long Gamma-Ray Bursts (GRBs) presenting light curves with an extended plateau phase in their X-ray afterglows obeys a correlation between the rest frame end time of the plateau, $T_a$, and its corresponding X-ray luminosity, $L_{a}$, Dainotti et al. (2008). In this work we perform an analysis of a total sample of 176 {\it Swift} GRBs with known redshifts, exhibiting afterglow plateaus. By adding a third parameter, that is the peak luminosity in the prompt emission, $L_{peak}$, we discover the existence of a new three parameter correlation. The scatter of data about this plane becomes smaller when a class-specific GRB sample is defined. This sample of 122 GRBs is selected from the total sample by excluding GRBs with associated Supernovae (SNe), X-ray flashes and short GRBs with extended emission. {\bf With this sample the three parameter correlation identifies a GRB `fundamental plane'}. Moreover, we further limit our analysis to GRBs with lightcurves having good data coverage and almost flat plateaus, 40 GRBs forming our `gold sample'. The intrinsic scatter, $\sigma_{int}=0.27 \pm 0.04$, for the three-parameter correlation for this last subclass is more than twice smaller than the value for the $L_{a}-T_a$ one, making this the tightest three parameter correlation involving the afterglow plateau phase. Finally, we also show that a slightly less tight correlation is present between $L_{peak}$ and a proxy for the total energy emitted during the plateau phase, $L_a T_a$, confirming the existence of an energy scaling between prompt and afterglow phases.

A large fraction of Gamma Ray Bursts (GRBs) displays an X-ray plateau phase within <10^{5} s from the prompt emission, proposed to be powered by the spin-down energy of a rapidly spinning newly born magnetar. In this work we use the properties of the Galactic neutron star population to constrain the GRB-magnetar scenario. We re-analyze the X-ray plateaus of all Swift GRBs with known redshift, between January 2005 and August 2014. From the derived initial magnetic field distribution for the possible magnetars left behind by the GRBs, we study the evolution and properties of a simulated GRB-magnetar population using numerical simulations of magnetic field evolution, coupled with Monte Carlo simulations of Pulsar Population Synthesis in our Galaxy. We find that if the GRB X-ray plateaus are powered by the rotational energy of a newly formed magnetar, the current observational properties of the Galactic magnetar population are not compatible with being formed within the GRB scenario (regardless of the GRB type or rate at z=0). Direct consequences would be that we should allow the existence of magnetars and "super-magnetars" having different progenitors, and that Type Ib/c SNe related to Long GRBs form systematically neutron stars with higher initial magnetic fields. We put an upper limit of <16 "super-magnetars" formed by a GRB in our Galaxy in the past Myr (at 99% c.l.). This limit is somewhat smaller than what roughly expected from Long GRB rates, although the very large uncertainties do not allow us to draw strong conclusion in this respect.

Matter mixing is one important topic in the study of core-collapse supernova (CCSN) explosions. In this paper, we perform two-dimensional hydrodynamic simulations to reproduce the high velocity $^{56}$Ni clumps observed in SN 1987A. This is the first time that large density perturbation is proposed in the CCSN progenitor to generate Rayleigh-Taylor (RT) instability and make the effective matter mixing. In the case of a spherical explosion, RT instability is efficient at both C+O/He and He/H interfaces of the SN progenitor. Radial coherent structures shown in perturbation patterns are important for obtaining high velocity $^{56}$Ni clumps. We can also obtain matter mixing features and high velocity $^{56}$Ni clumps in some cases of aspherical explosion. We find that one of the most favorable models in our work has a combination of bipolar and equatorially asymmetric explosions in which at least 25\% of density perturbation is introduced at different composition interfaces of the CCSN progenitor. These simulation results are comparable to the observational findings of SN 1987A.

We present an analysis of 123 Gamma-ray bursts (GRBs) with known redshifts possessing an afterglow plateau phase. We reveal that $L_a-T^{*}_a$ correlation between the X-ray luminosity $L_a$ at the end of the plateau phase and the plateau duration, $T^*_a$, in the GRB rest frame has a power law slope different, within more than 2 $\sigma$, from the slope of the prompt $L_{f}-T^{*}_{f}$ correlation between the isotropic pulse peak luminosity, $L_{f}$, and the pulse duration, $T^{*}_{f}$, from the time since the GRB ejection. Analogously, we show differences between the prompt and plateau phases in the energy-duration distributions with the afterglow emitted energy being on average $10\%$ of the prompt emission. Moreover, the distribution of prompt pulse versus afterglow spectral indexes do not show any correlation. In the further analysis we demonstrate that the $L_{peak}-L_a$ distribution, where $L_{peak}$ is the peak luminosity from the start of the burst, is characterized with a considerably higher Spearman correlation coefficient, $\rho=0.79$, than the one involving the averaged prompt luminosity, $L_{prompt}-L_a$, for the same GRB sample, yielding $\rho=0.60$. Since some of this correlation could result from the redshift dependences of the luminosities, namely from their cosmological evolution we use the Efron-Petrosian method to reveal the intrinsic nature of this correlation. We find that a substantial part of the correlation is intrinsic. We apply a partial correlation coefficient to the new de-evolved luminosities showing that the intrinsic correlation exists.

In this work we present a procedure to infer the mass of progenitors and remnants of Gamma Ray Bursts (GRB), starting from the observed energy $E_{iso}^{GRB}$ emitted isotropically and considering the associated emission of Gravitational Waves (GW) $ E_{iso}^{GW}$ in the different phases. We assume that the GW energy of the progenitor $E_{PROG}^{GW}$ is emitted partially during a star collapse, and the residual energy is related to the GW energy emitted by the remnant. We take a sample of $237$ Long GRB, and use an hybrid Montecarlo procedure to explore, for each of them, a region of possible solutions of $ E_{iso}^{GW}$ as a function of the masses, radii, oblateness, rotation frequencies of progenitor and remnant and the fraction of energy $k$ emitted as GW by the GRB. We discriminate between a Neutron Star (NS) or Black Hole (BH) for the remnant and obtain interesting values for the GW emitted by the remnant NS or BH, for the conversion factor $k$ of and for the masses and radii of GRB progenitor stars. We also obtain remnant populations with mean masses, mean GW frequencies and GRB frequency of GW emission in agreement with the most accepted models.

Gamma Ray Bursts (GRBs) visible up to very high redshift have become attractive targets as potential new distance indicators. It is still not clear whether the relations proposed so far originate from an unknown GRB physics or result from selection effects. We investigate this issue in the case of the $L_X-T^*_a$ correlation (hereafter LT) between the X-ray luminosity $L_X (T_a)$ at the end of the plateau phase, $T_a$, and the rest frame time $T^{*}_a$. We devise a general method to build mock data sets starting from a GRB world model and taking into account selection effects on both time and luminosity. This method shows how not knowing the efficiency function could influence the evaluation of the intrinsic slope of any correlation and the GRB density rate. We investigate biases (small offsets in slope or normalization) that would occur in the LT relation as a result of truncations, possibly present in the intrinsic distributions of $L_X$ and $T^*_a$. We compare these results with the ones in Dainotti et al. (2013) showing that in both cases the intrinsic slope of the LT correlation is $\approx -1.0$. This method is general, therefore relevant to investigate if any other GRB correlation is generated by the biases themselves. Moreover, because the farthest GRBs and star-forming galaxies probe the reionization epoch, we evaluate the redshift-dependent ratio $\Psi(z)=(1+z)^{\alpha}$ of the GRB rate to star formation rate (SFR). We found a modest evolution $-0.2\leq \alpha \leq 0.5$ consistent with Swift GRB afterglow plateau in the redshift range $0.99<z<9.4$.

An intrinsic correlation has been identified between the luminosity and duration of plateaus in the X-ray afterglows of Gamma-Ray Bursts (GRBs; Dainotti et al. 2008), suggesting a central engine origin. The magnetar central engine model predicts an observable plateau phase, with plateau durations and luminosities being determined by the magnetic fields and spin periods of the newly formed magnetar. This paper analytically shows that the magnetar central engine model can explain, within the 1$\sigma$ uncertainties, the correlation between plateau luminosity and duration. The observed scatter in the correlation most likely originates in the spread of initial spin periods of the newly formed magnetar and provides an estimate of the maximum spin period of ~35 ms (assuming a constant mass, efficiency and beaming across the GRB sample). Additionally, by combining the observed data and simulations, we show that the magnetar emission is most likely narrowly beamed and has $\lesssim$20% efficiency in conversion of rotational energy from the magnetar into the observed plateau luminosity. The beaming angles and efficiencies obtained by this method are fully consistent with both predicted and observed values. We find that Short GRBs and Short GRBs with Extended Emission lie on the same correlation but are statistically inconsistent with being drawn from the same distribution as Long GRBs, this is consistent with them having a wider beaming angle than Long GRBs.

We explore the spectral and polarization properties of photospheric emissions from stratified jets in which multiple components, separated by a sharp velocity shear regions, are distributed in lateral direction. Propagation of thermal photons injected at high optical depth region are calculated until they escape from the photosphere. It is found that presence of the lateral structure within the jet leads to non-thermal feature of the spectra and significant polarization signal in the resulting emission. The deviation from thermal spectra as well as the polarization degree tends to be enhanced as the velocity gradient in the shear region increases. In particular, we show that emissions from multi-component jet can reproduce the typical observed spectra of gamma-ray bursts (GRBs) irrespective to the position of the observer when a velocity shear region is closely spaced in various lateral ($\theta$) positions. The degree of polarization associated in the emission is significant (> few%) at wide range of observer angles and can be higher than 30%.

We study the dark energy equation of state as a function of redshift in a non-parametric way, without imposing any {\it a priori} $w(z)$ (ratio of pressure over energy density) functional form. As a check of the method, we test our scheme through the use of synthetic data sets produced from different input cosmological models which have the same relative errors and redshift distribution as the real data. Using the luminosity-time $L_{X}-T_{a}$ correlation for GRB X-ray afterglows (the Dainotti et al. correlation), we are able to utilize GRB sample from the {\it Swift} satellite as probes of the expansion history of the Universe out to $z \approx 10$. Within the assumption of a flat FLRW universe and combining SNeIa data with BAO constraints, the resulting maximum likelihood solutions are close to a constant $w=-1$. If one imposes the restriction of a constant $w$, we obtain $w=-0.99 \pm 0.06$ (consistent with a cosmological constant) with the present day Hubble constant as $H_{0}=70.0 \pm 0.6$ ${\rm km} \, {\rm s}^{-1} {\rm Mpc}^{-1}$ and density parameter as $\Omega_{\Lambda 0}=0.723 \pm 0.025$, while non-parametric $w(z)$ solutions give us a probability map which is centred at $H_{0}=70.04 \pm 1$ ${\rm km} \, {\rm s}^{-1} {\rm Mpc}^{-1}$ and $\Omega_{\Lambda 0}=0.724 \pm 0.03$. Our chosen GRB data sample with full correlation matrix allows us to estimate the amount, as well as quality (errors) of data, needed to constrain $w(z)$ in the redshift range extending an order of magnitude in beyond the farthest SNeIa measured.

Gamma -ray bursts (GRBs) observed up to redshifts $z>9.4$ can be used as possible probes to test cosmological models. Here we show how changes of the slope of the {\it luminosity $L^*_X$ -break time $T^*_a$} correlation in GRB afterglows, hereafter the LT correlation, affect the determination of the cosmological parameters. With a simulated data set of 101 GRBs with a central value of the correlation slope that differs on the intrinsic one by a $5\sigma$ factor, we find an overstimated value of the matter density parameter, $\Omega_M$, compared to the value obtained with SNe Ia, while the Hubble constant, $H_0$, best fit value is still compatible in 1$\sigma$ compared to other probes. We show that this compatibility of $H_0$ is due to the large intrinsic scatter associated with the simulated sample. Instead, if we consider a subsample of high luminous GRBs ($HighL$), we find that both the evaluation of $H_0$ and $\Omega_M$ are not more compatible in 1$\sigma$ and $\Omega_M$ is underestimated by the $13\%$. However, the $HighL$ sample choice reduces dramatically the intrinsic scatter of the correlation, thus possibly identifying this sample as the standard canonical `GRBs' confirming previous results presented in Dainotti et al. (2010,2011). Here, we consider the LT correlation as an example, but this reasoning can be extended also for all other GRB correlations. In literature so far GRB correlations are not corrected for redshift evolution and selection biases, therefore we are not aware of their intrinsic slopes and consequently how far the use of the observed correlations can influence the derived `best' cosmological settings. Therefore, we conclude that any approach that involves cosmology should take into consideration only intrinsic correlations not the observed ones.

Gamma-ray bursts (GRBs), which have been observed up to redshifts $z \approx 9.5$ can be good probes of the early universe and have the potential of testing cosmological models. The analysis by Dainotti of GRB Swift afterglow lightcurves with known redshifts and definite X-ray plateau shows an anti-correlation between the \underline{rest frame} time when the plateau ends (the plateau end time) and the calculated luminosity at that time (or approximately an anti-correlation between plateau duration and luminosity). We present here an update of this correlation with a larger data sample of 101 GRBs with good lightcurves. Since some of this correlation could result from the redshift dependences of these intrinsic parameters, namely their cosmological evolution we use the Efron-Petrosian method to estimate the luminosity and time evolution and to correct for this effects to determine the intrinsic nature of this correlation.

Gamma-ray bursts (GRBs), which have been observed up to redshifts z approx 9.5 can be good probes of the early universe and have the potential of testing cosmological models. The analysis by Dainotti of GRB Swift afterglow lightcurves with known redshifts and definite X-ray plateau shows an anti-correlation between the rest frame time when the plateau ends (the plateau end time) and the calculated luminosity at that time (or approximately an anti-correlation between plateau duration and luminosity). We present here an update of this correlation with a larger data sample of 101 GRBs with good lightcurves. Since some of this correlation could result from the redshift dependences of these intrinsic parameters, namely their cosmological evolution we use the Efron-Petrosian method to reveal the intrinsic nature of this correlation. We find that a substantial part of the correlation is intrinsic and describe how we recover it and how this can be used to constrain physical models of the plateau emission, whose origin is still unknown. The present result could help clarifing the debated issue about the nature of the plateau emission.

Relativistic jets propagating through an ambient medium must produce some observational effects along their side boundaries because of interactions across the large velocity gradient. One possible effect of such an interaction would be a sheared magnetic field structure at the jet boundaries, leading to a characteristic radio polarization pattern. As proposed by Ostrowski, another effect can come from the generation of a high energy cosmic ray component at the boundary, producing dynamic effects on the medium surrounding the jet and forming a cocoon dominated by cosmic rays with a decreased thermal gas emissivity. We selected this process for our first attempt to look for the effects of this type of interaction. We analyzed the Chandra X-ray data for the radio galaxy M87 in order to verify if the expected regions of diminished emissivity may be present near the spectacular X-ray jet in this source. The detailed analysis of the data, merged from 42 separate observations, shows signatures of lower emissivity surrounding the jet. In particular we detect an intensity dip along the part of the jet, which would be approximately 150 pc x 2 kpc in size, if situated along the jet which is inclined toward us. Due to a highly non-uniform X-ray background in the central region we are not able to claim the discovery of a cosmic ray cocoon around the M87 jet: we only have demonstrated that the data show morphological structures which could be accounted for if a cosmic ray cocoon exists.

Thanks to their enormous energy release, Gamma Rays Bursts (GRBs) have recently attracted a lot of interest to probe the Hubble diagram (HD) deep into the matter dominated era and hence complement Type Ia Supernovae (SNeIa). We consider here three different calibration methods based on the use of a fiducial LCDM model, on cosmographic parameters and on the local regression on SNeIa to calibrate the scaling relations proposed as an equivalent to the Phillips law to standardize GRBs finding any significant dependence. We then investigate the evolution of these parameters with the redshift to obtain any statistical improvement. Under this assumption, we then consider possible systematics effects on the HDs introduced by the calibration method, the averaging procedure and the homogeneity of the sample arguing against any significant bias.

To reveal and understand astrophysical processes responsible for the Gamma Ray Burst (GRB) phenomenon, it is crucial to discover and understand relations between their observational properties. The presented study is performed in the GRB rest frames and it uses a sample of 62 long GRBs from our sample of 77 Swift GRBs with known redshifts. Following the earlier analysis of the afterglow {\it characteristic luminosity $L^*_a$ -- break time $T^*_a$} correlation for a sample of long GRBs \citep{Dainotti2010} we extend it to correlations between the afterglow and the prompt emission GRB physical parameters. We reveal a tight physical scaling between the mentioned afterglow luminosity $ L^*_a$ and the prompt emission {\it mean} luminosity $<L^*_p>_{45} \equiv E_{iso}/T^*_{45}$. The distribution, with the Spearman correlation coefficient reaching 0.95 for the data subsample with most regular light curves, can be fitted with approximately $L^*_a \propto {<L^*_p>_{45}}^{0.7}$. We also analyzed correlations of $L^*_a$ with several other prompt emission parameters, including the isotropic energy $E_{iso}$, the peak energy in the $\nu F_{\nu}$ spectrum, $E_{peak}$, and the variability parameter, $V$, defined by \cite{N000}. As a result, we reveal significant correlations also between these quantities, with an exception of the variability parameter. The main result of the present study is the discovery that the highest correlated GRB subsample in the \citet{Dainotti2010} afterglow analysis, for the GRBs with canonical X\,-\,ray light curves, leads also to the highest {\it prompt-afterglow} correlations and such events can be considered to form a sample of standard GRBs for astrophysics and cosmology.

Gamma Ray Bursts (GRBs) are the most energetic sources in the universe and among the farthest known astrophysical sources. These features make them appealing candidates as standard candles for cosmological applications so that studying the physical mechanisms for the origin of the emission and correlations among their observable properties is an interesting task. We consider here the luminosity L*X - break time Ta* (hereafter LT) correlation and investigate whether there are systematics induced by selection effects or redshift dependent calibra- tion. We perform this analysis both for the full sample of 77 GRBs with known redshift and for the subsample of GRBs having canonical X-ray light curves, hereafter called U0095 sample. We do not find any systematic bias thus con- firming the existence of physical GRB subclasses revealed by tight correlations of their afterglow properties. Furthermore, we study the possibility of applying the LT correlation as a redshift estimator both for the full distribution and for the canonical lightcurves. The large uncertainties and the non negligible intrin- sic scatter make the results not so encouraging, but there are nevertheless some hints motivating a further analysis with an increased U0095 sample.

Gamma Ray Bursts (GRB) observed up to redshifts $z>8$ are fascinating objects to study due to their still unexplained relativistic outburst mechanisms and a possible use to test cosmological models. Our analysis of 77 GRB afterglows with known redshifts revealed a physical subsample of long GRBs with canonical {\it plateau breaking to power-law} light curves with a significant {\it luminosity $L^*_X$ - break time $T^*_a$} correlation in the GRB rest frame. This subsample forms approximately the {\it upper envelope} of the studied distribution. We have also found a similar relation for a small sample of GRB afterglows that belong to the intermediate class (IC) between the short and the long ones. It proves that within the full sample of afterglows there exist physical subclasses revealed here by tight correlations of their afterglow properties. The afterglows with regular (`canonical') light curves obey not only a mentioned tight physical scaling, but -- for a given $T^*_a$ -- the more regular progenitor explosions lead to preferentially brighter afterglows.

We present the Hubble diagram (HD) of 66 Gamma Ray Bursts (GRBs) derived using only data from their X - ray afterglow lightcurve. To this end, we use the recently updated L_X - T_a correlation between the break time T_a and the X - ray luminosity L_X measured at T_a calibrated from a sample of Swift GRBs with lightcurves well fitted by the Willingale et al. (2007) model. We then investigate the use of this HD to constrain cosmological parameters when used alone or in combination with other data showing that the use of GRBs leads to constraints in agreement with previous results in literature. We finally argue that a larger sample of high luminosity GRBs can provide a valuable information in the search for the correct cosmological model.

(Shortened) [...] After recalling the basic features of the "fireshell model", we emphasize the following novel results: 1) the interpretation of the X-ray flares in GRB afterglows as due to the interaction of the optically thin fireshell with isolated clouds in the CircumBurst Medium (CBM); 2) an interpretation as "fake - disguised" short GRBs of the GRBs belonging to the class identified by Norris & Bonnell [...] consistent with an origin from the final coalescence of a binary system in the halo of their host galaxies with particularly low CBM density [...]; 3) the first attempt to study a genuine short GRB with the analysis of GRB 050509B, that reveals indeed still an open question; 4) the interpretation of the GRB-SN association in the case of GRB 060218 via the "induced gravitational collapse" process; 5) a first attempt to understand the nature of the "Amati relation", a phenomenological correlation between the isotropic-equivalent radiated energy of the prompt emission E_{iso} with the cosmological rest-frame \nu F_{\nu} spectrum peak energy E_{p,i}. In addition, recent progress on the thermalization of the electron-positron plasma close to their formation phase, as well as the structure of the electrodynamics of Kerr-Newman Black Holes are presented. An outlook for possible explanation of high-energy phenomena in GRBs to be expected from the AGILE and the Fermi satellites are discussed. As an example of high energy process, the work by Enrico Fermi dealing with ultrarelativistic collisions is examined. It is clear that all the GRB physics points to the existence of overcritical electrodynamical fields. In this sense we present some progresses on a unified approach to heavy nuclei and neutron stars cores, which leads to the existence of overcritical fields under the neutron star crust.

Gamma ray bursts (GRBs) have recently attracted much attention as a possible way to extend the Hubble diagram to very high redshift. To this aim, the luminosity (or isotropic emitted energy) of a GRB at redshift z must be evaluated from a correlation with a distance independent quantity so that one can then solve for the luminosity distance D_L(z) and hence the distance modulus mu(z). Averaging over five different two parameters correlations and using a fiducial cosmological model to calibrate them, Schaefer (2007) has compiled a sample of 69 GRBs with measured mu(z) which has since then been widely used to constrain cosmological parameters. We update here that sample by many aspects. First, we add a recently found correlation for the X - ray afterglow and use a Bayesian inspired fitting method to calibrate the different GRBs correlations known insofar assuming a fiducial LCDM model in agreement with the recent WMAP5 data. Averaging over six correlations, we end with a new GRBs Hubble diagram comprising 83 objects. We also extensively explore the impact of varying the fiducial cosmological model considering how the estimated mu(z) change as a function of the $(\Omega_M, w_0, w_a)$ parameters of the Chevallier - Polarski - Linder phenomenological dark energy equation of state. In order to avoid the need of assuming an {\it a priori} cosmological model, we present a new calibration procedure based on a model independent local regression estimate of mu(z) using the Union SNeIa sample to calibrate the GRBs correlations. This finally gives us a GRBs Hubble diagram made out of 69 GRBs whose estimated distance modulus mu(z) is almost independent on the underlying cosmological model.

(Shortened) CONTEXT: [...] GRB060614 is the first nearby long duration GRB clearly not associated to a bright Ib/c supernova. Moreover, its duration (T_{90} ~ 100s) makes it hardly classifiable as a short GRB. It presents strong similarities with GRB970228, the prototype of the new class of "fake" short GRBs that appear to originate from the coalescence of binary neutron stars or white dwarfs spiraled out into the galactic halo. AIMS: Within the "canonical" GRB scenario based on the "fireshell" model, we test if GRB060614 can be a "fake" or "disguised" short GRB. [...] METHODS: We fit GRB060614 light curves in Swift's BAT (15-150keV) and XRT (0.2-10keV) energy bands. Within the fireshell model, light curves are formed by two well defined and different components: the Proper-GRB (P-GRB), emitted at the fireshell transparency, and the extended afterglow, due to the interaction between the leftover accelerated baryonic and leptonic shell and the CBM. RESULTS: We determine the two free parameters describing the GRB source within the fireshell model. [...] A small average CBM density [...] is inferred, typical of galactic halos. The first spikelike emission is identified with the P-GRB and the following prolonged emission with the extended afterglow peak.[...] CONCLUSIONS: The anomalous GRB060614 finds a natural interpretation within our canonical GRB scenario: it is a "disguised" short GRB. [...] This result points to an old binary system, likely formed by a white dwarf and a neutron star, as the progenitor of GRB060614 and well justify the absence of an associated SN Ib/c. Particularly important for further studies of the final merging process are the temporal structures in the P-GRB down to 0.1s.

The fireshell model for Gamma-Ray Bursts (GRBs) naturally leads to a canonical GRB composed of a proper-GRB (P-GRB) and an afterglow. P-GRBs, introduced by us in 2001, are sometimes considered "precursors" of the main GRB event in the current literature. We show in this paper how the fireshell model leads to the understanding of the structure of GRBs, with precise estimates of the time sequence and intensities of the P-GRB and the of the afterglow. It leads as well to a natural classification of the canonical GRBs which overcomes the traditional one in short and long GRBs.

The Gamma-Ray Bursts (GRBs) offer the unprecedented opportunity to observe for the first time the blackholic energy extracted by the vacuum polarization during the process of gravitational collapse to a black hole leading to the formation of an electron-positron plasma. The uniqueness of the Kerr-Newman black hole implies that very different processes originating from the gravitational collapse a) of a single star in a binary system induced by the companion, or b) of two neutron stars, or c) of a neutron star and a white dwarf, do lead to the same structure for the observed GRB. The recent progress of the numerical integration of the relativistic Boltzmann equations with collision integrals including 2-body and 3-body interactions between the particles offer a powerful conceptual tool in order to differentiate the traditional "fireball" picture, an expanding hot cavity considered by Cavallo and Rees, as opposed to the "fireshell" model, composed of an internally cold shell of relativistically expanding electron-positron-baryon plasma. The analysis of the fireshell naturally leads to a canonical GRB composed of a proper-GRB and an extended afterglow. By recalling the three interpretational paradigms for GRBs we show how the fireshell model leads to an understanding of the GRB structure and to an alternative classification of short and long GRBs.

Gamma ray bursts (GRBs) have recently attracted much attention as a possible way to extend the Hubble diagram to very high redshift. However, the large scatter in their intrinsic properties prevents directly using them as distance indicator so that the hunt is open for a relation involving an observable property to standardize GRBs in the same way as the Phillips law makes it possible to use Type Ia Supernovae (SNeIa) as standardizable candles. We use here the data on the X - ray decay curve and spectral index of a sample of GRBs observed with the Swift satellite. These data are used as input to a Bayesian statistical analysis looking for a correlation between the X - ray luminosity L_X(T_a) and the time constant T_a of the afterglow curve. We find a linear relation between \log{[L_X(T_a)]} and \log{[T_a/(1+z)]} with an intrinsic scatter sigma_{int} = 0.33 comparable to previously reported relations. Remarkably, both the slope and the intrinsic scatter are almost independent on the matter density Omega_M and the constant equation of state w of the dark energy component thus suggesting that the circularity problem is alleviated for the $L_X - T_a$ relation.

(Shortened) CONTEXT: [...] AIMS: Motivated by the relation proposed by Amati and collaborators, we look within the ``fireshell'' model for a relation between the peak energy E_p of the \nu F_\nu total time-integrated spectrum of the afterglow and the total energy of the afterglow E_{aft}, which in our model encompasses and extends the prompt emission. METODS: [...] Within the fireshell model [...] We can then build two sets of ``gedanken'' GRBs varying the total energy of the electron-positron plasma E^{e^\pm}_{tot} and keeping the same baryon loading B of GRB050315. The first set assumes for the effective CBM density the one obtained in the fit of GRB050315. The second set assumes instead a constant CBM density equal to the average value of the GRB050315 prompt phase. RESULTS: For the first set of ``gedanken'' GRBs we find a relation E_p\propto (E_{aft})^a, with a = 0.45 \pm 0.01, whose slope strictly agrees with the Amati one. Such a relation, in the limit B \to 10^{-2}, coincides with the Amati one. Instead, in the second set of ``gedanken'' GRBs no correlation is found. CONCLUSIONS: Our analysis excludes the Proper-GRB (P-GRB) from the prompt emission, extends all the way to the latest afterglow phases and is independent on the assumed cosmological model, since all ``gedanken'' GRBs are at the same redshift. The Amati relation, on the other hand, includes also the P-GRB, focuses on the prompt emission only, and is therefore influenced by the instrumental threshold which fixes the end of the prompt emission, and depends on the assumed cosmology. This may well explain the intrinsic scatter observed in the Amati relation.

We show the preliminary results of the application of our "fireshell" model to GRB060124. This source is very peculiar because it is the first event for which both the prompt and the afterglow emission were observed simultaneously by the three Swift instruments: BAT (15-350 keV), XRT (0.2-10 keV) and UVOT (170-650 nm), due to the presence of a precursor ~ 570 s before the main burst. We analyze GRB060124 within our "canonical" GRB scenario, identifying the precursor with the P-GRB and the prompt emission with the afterglow peak emission. In this way we reproduce correctly the energetics of both these two components. We reproduce also the observed time delay between the precursor (P-GRB) and the main burst. The effect of such a time delay in our model will be discussed.

Within the "fireshell" model for the Gamma-Ray Bursts (GRBs) we define a "canonical GRB" light curve with two sharply different components: the Proper-GRB (P-GRB), emitted when the optically thick fireshell of electron-positron plasma originating the phenomenon reaches transparency, and the afterglow, emitted due to the collision between the remaining optically thin fireshell and the CircumBurst Medium (CBM). We outline our "canonical GRB" scenario, with a special emphasis on the discrimination between "genuine" and "fake" short GRBs.

Within the "fireshell" model we define a "canonical GRB" light curve with two sharply different components: the Proper-GRB (P-GRB), emitted when the optically thick fireshell of electron-positron plasma originating the phenomenon reaches transparency, and the afterglow, emitted due to the collision between the remaining optically thin fireshell and the CircumBurst Medium (CBM). We here present the consequences of such a scenario on the theoretical interpretation of the nature of "long" and "short" GRBs.

(Shortened) We show by example how the uncoding of Gamma-Ray Bursts (GRBs) offers unprecedented possibilities to foster new knowledge in fundamental physics and in astrophysics. After recalling some of the classic work on vacuum polarization in uniform electric fields by Klein, Sauter, Heisenberg, Euler and Schwinger, we summarize some of the efforts to observe these effects in heavy ions and high energy ion collisions. We then turn to the theory of vacuum polarization around a Kerr-Newman black hole, leading to the extraction of the blackholic energy, to the concept of dyadosphere and dyadotorus, and to the creation of an electron-positron-photon plasma. We then present a new theoretical approach encompassing the physics of neutron stars and heavy nuclei. It is shown that configurations of nuclear matter in bulk with global charge neutrality can exist on macroscopic scales and with electric fields close to the critical value near their surfaces. These configurations may represent an initial condition for the process of gravitational collapse, leading to the creation of an electron-positron-photon plasma: the basic self-accelerating system explaining both the energetics and the high energy Lorentz factor observed in GRBs. We then turn to recall the two basic interpretational paradigms of our GRB model. [...] We then turn to the special role of the baryon loading in discriminating between "genuine" short and long or "fake" short GRBs [...] We finally turn to the GRB-Supernova Time Sequence (GSTS) paradigm: the concept of induced gravitational collapse. [...] We then present some general conclusions.

In the "fireshell" model we define a "canonical GRB" light curve with two sharply different components: the Proper-GRB (P-GRB), emitted when the optically thick fireshell of electron-positron plasma originating the phenomenon reaches transparency, and the afterglow, emitted due to the collision between the remaining optically thin fireshell and the CircumBurst Medium (CBM). We outline our "canonical GRB" scenario, originating from the gravitational collapse to a black hole, with a special emphasis on the discrimination between "genuine" and "fake" short GRBs.

(shortened) We study the Gamma-Ray Burst (GRB) 060218: a particularly close source at z=0.033 with an extremely long duration, namely T_{90} ~ 2000 s, related to SN 2006aj. [...] I present the fitting time consuming procedure. In order to show its sensitivity I also present two examples of fits with the same value of B and different value of E_{e^\pm}^{tot}. We fit the X- and \gamma-ray observations by Swift of GRB 060218 in the 0.1-150 keV energy band during the entire time of observations from 0 all the way to 10^6 s within a unified theoretical model. The free parameters of our theory are only three, namely the total energy E_{e\pm}^{tot} of the e^\pm plasma, its baryon loading B \equiv M_Bc^2/E_{e\pm}^{tot}, as well as the CircumBurst Medium (CBM) distribution. We justify the extremely long duration of this GRB by a total energy E_{e\pm}^{tot} = 2.32\times 10^{50} erg, a very high value of the baryon loading B=1.0\times 10^{-2} and the effective CircumBurst Medium (CBM) density which shows a radial dependence n_{cbm} \propto r^{-\alpha} with 1.0 \leq \alpha \leq 1.7 and monotonically decreases from 1 to $10^{-6}$ particles/cm$^3$. We recall that this value of the $B$ parameter is the highest among the sources we have analyzed and it is very close to its absolute upper limit expected. [...] We also think that the smallest possible black hole, formed by the gravitational collapse of a neutron star in a binary system, is consistent with the especially low energetics of the class of GRBs associated with SNe Ib/c.

On the basis of the recent understanding of GRB050315 and GRB060218, we return to GRB970228, the first Gamma-Ray Burst (GRB) with detected afterglow. We proposed it as the prototype for a new class of GRBs with "an occasional softer extended emission lasting tenths of seconds after an initial spikelike emission". Detailed theoretical computation of the GRB970228 light curves in selected energy bands for the prompt emission are presented and compared with observational BeppoSAX data. From our analysis we conclude that GRB970228 and likely the ones of the above mentioned new class of GRBs are "canonical GRBs" have only one peculiarity: they exploded in a galactic environment, possibly the halo, with a very low value of CBM density. Here we investigate how GRB970228 unveils another peculiarity of this class of GRBs: they do not fulfill the "Amati relation". We provide a theoretical explanation within the fireshell model for the apparent absence of such correlation for the GRBs belonging to this new class.

(Shortened) The Swift and HETE-2 discovery of an afterglow associated possibly with short GRBs opened the new problematic of their nature and classification. This has been further enhanced by the GRB060614 observation and by a re-analysis of the BATSE catalog leading to the identification of a new GRB class with "an occasional softer extended emission lasting tenths of seconds after an initial spikelike emission". We plan: a) to fit this new class of "hybrid" sources within our "canonical GRB" scenario, where all GRBs are generated by a "common engine" (i.e. the gravitational collapse to a black hole); b) to propose GRB970228 as the prototype of the such a class. We analyze BeppoSAX data on GRB970228 in the 40-700 keV and 2-26 keV energy bands within the "fireshell" model. We find that GRB970228 is a "canonical GRB", like e.g. GRB050315, with the main peculiarity of a particularly low CircumBurst Medium (CBM) average density n_{cbm}~10^{-3} #/cm^3. We also simulate the light curve corresponding to a rescaled CBM density profile with n_{cbm}=1 #/cm^3. From such a comparison it follows that the total time-integrated luminosity is a faithful indicator of the GRB nature, contrary to the peak luminosity which is merely a function of the CBM density. We call attention on discriminating the short GRBs between the "genuine" and the "fake" ones. The "genuine" ones are intrinsically short, with baryon loading B \la 10^{-5}, as stated in our original classification. The "fake" ones, characterized by an initial spikelike emission followed by an extended emission lasting tenths of seconds, have a baryon loading 10^{-4} \la B \leq 10^{-2}. They are observed as such only due to an underdense CBM consistent with a galactic halo environment which deflates the afterglow intensity.

We plan to fit the complete gamma- and X-ray light curves of the long duration GRB060218, including the prompt emission, in order to clarify the nature of the progenitors and the astrophysical scenario of the class of GRBs associated to SNe Ib/c. The initial total energy of the electron-positron plasma E_{e^\pm}^{tot}=2.32\times 10^{50} erg has a particularly low value similarly to the other GRBs associated with SNe. For the first time we observe a baryon loading B=10^{-2} which coincides with the upper limit for the dynamical stability of the fireshell. The effective CircumBurst Medium (CBM) density shows a radial dependence n_{cbm} \propto r^{-\alpha} with 1.0<\alpha<1.7 and monotonically decreases from 1 to 10^{-6} particles/cm^3. Such a behavior is interpreted as due to a fragmentation in the fireshell. Analogies with the fragmented density and filling factor characterizing Novae are outlined. The fit presented is particularly significant in view of the complete data set available for GRB060218 and of the fact that it fulfills the Amati relation. We fit GRB060218, usually considered as an X-Ray Flash (XRF), as a "canonical GRB" within our theoretical model. The smallest possible black hole, formed by the gravitational collapse of a neutron star in a binary system, is consistent with the especially low energetics of the class of GRBs associated with SNe Ib/c. We give the first evidence for a fragmentation in the fireshell. Such a fragmentation is crucial in explaining both the unusually large T_{90} and the consequently inferred abnormal low value of the CBM effective density.

We outline the main results of our GRB model, based on the three interpretation paradigms we proposed in July 2001, comparing and contrasting them with the ones in the current literature. Thanks to the observations by Swift and by VLT, this analysis points to a "canonical GRB" originating from markedly different astrophysical scenarios. The communality is that they are all emitted in the formation of a black hole with small or null angular momentum. The following sequence appears to be canonical: the vacuum polarization process creating an optically thick self accelerating electron-positron plasma; the engulfment of baryonic mass during the plasma expansion; the adiabatic expansion of the optically thick "fireshell" up to the transparency; the interaction of the remaining accelerated baryons with the interstellar medium (ISM). This leads to the canonical GRB composed of a proper GRB (P-GRB), emitted at the moment of transparency, followed by an extended afterglow. The parameters are the plasma total energy, the fireshell baryon loading and the ISM filamentary distribution around the source. In the limit of no baryon loading the total energy is radiated in the P-GRB. In this limit, the canonical GRBs explain as well the short GRBs.

GRB970228 is analyzed as a prototype to understand the relative role of short GRBs and their associated afterglows, recently observed by Swift and HETE-II. Detailed theoretical computation of the GRB970228 light curves in selected energy bands are presented and compared with observational BeppoSAX data.

The luminosity and the spectral distribution of the afterglow of GRB 031203 have been presented within our theoretical framework, which envisages the GRB structure as composed by a proper-GRB, emitted at the transparency of an electron-positron plasma with suitable baryon loading, and an afterglow comprising the "prompt emission" as due to external shocks. In addition to the GRB emission, there appears to be a prolonged soft X-Ray emission lasting for 10^6-10^7 seconds followed by an exponential decay. This additional source has been called by us URCA-3. It is urgent to establish if this component is related to the GRB or to the Supernova (SN). In this second case, there are two possibilities: either the interaction of the SN ejecta with the interstellar medium or, possibly, the cooling of a young neutron star formed in the SN 2003lw process. The analogies and the differences between this triptych GRB 031203 / SN 2003lw / URCA-3 and the corresponding ones GRB 980425 / SN 1998bw / URCA-1 and GRB 030329 / SN 2003dh / URCA-2, as well as GRB 060218 / SN 2006aj are discussed.