Short gamma-ray bursts result from mergers of two neutron stars or from collapsars, but probably at a smaller rate. In the first case, a kilonova occurs while in the second case a Type Ic supernova is expected. } {Even if future observations of kilonovae in association with gravitational wave events provide better data, detecting a kilonova during an afterglow follow-up would remain useful for exploring the diversity of the kilonova phenomenon. As supernovae produce a weaker gravitational signal, afterglow follow-up will be the only possible method to find one. In this work, we identify the conditions of the burst energy, external density, kilonova mass, supernova luminosity, that are necessary for the detection of a kilonova or supernova in the follow-up of short GRB afterglows.} {We have used a simple kilonova model to obtain the peak luminosities and times as a function of mass, expansion velocity and ejected matter opacity. Afterglow light curves are computed for a uniform medium and a stellar wind, in the kilonova and supernova cases, respectively.} {We represent, using diagrams of the burst kinetic energy vs. density of the external medium, the domains where the kilonova or supernova at maximum is brighter than the afterglow. %\LEt{Please check I have retained your intended meaning -> It's OK

Ever since the insightful analysis of the durations of GRBs by [Kouveliotou:93], GRBs have most often been classified into two populations: "short bursts" (shorter than 2.0 seconds) and "long bursts" (longer than 2.0 seconds). However, recent works have suggested the existence of an intermediate population in the bursts observed by the Swift satellite. Moreover, some researchers have questioned the universality of the 2.0-second dividing line between short and long bursts: some bursts may be short but actually result from collapsars, the physical mechanism behind normally long bursts, and some long ones may originate from mergers, the usual progenitors of short GRBs. In this work, we focus on GRBs detected by the Fermi satellite and study the distribution of their durations measured in the observer's reference frame and, for those with known redshifts, in the bursts' reference frames. However, there are relatively few bursts with measured redshifts, and this makes an accurate study difficult. To overcome this problem, we follow [Zhang:2018] and determine a "pseudo-redshift" from the correlation relation between the luminosity $L_p$ and the energy $E_p$, both of which are calculated at the peak of the flux. Interestingly, we find that the uncertainties in the quantities observed and used in the determination of pseudo-redshifts, do affect the precision of the individual results significantly, but they keep the distribution of pseudo-redshifts very similar to that of the actual ones and thus allow us to use pseudo-redshifts for our statistical study. We use the reduced chi-square and the maximization of the log-likelihood to statistically analyze the distribution of Fermi GRB durations. Both methods show that the distribution of the observed and the intrinsic bursts durations are better represented by two groups, rather than three.

We use the correlation relation between the energy emitted by the GRBs in their prompt phases and the X-ray afterglow fluxes, in an effort to constrain cosmological parameters and aiming to construct a Hubble diagram at high redshifts, i.e. beyond those found with Type Ia supernovae. We use a sample of 126 \textit{Swift} GRBs, that we have selected among more than 800 long bursts observed until April 2015. The selection is based on a few observational constraints: GRB flux higher than 0.4 $photons/cm^2/s$ in the band 15-150 keV, spectrum fitted with simple power law, redshift accurately known and given, and X-ray afterglow observed and flux measured. The statistical method of maximum likelihood is then used to determine the best cosmological parameters ($\Omega_M$, $\Omega_{\Lambda}$) that give the best correlation for two relations: a) the Amati relation (between intrinsic spectral peak energy $E_{p,i}$ and the equivalent isotropic energy), b) the Dainotti relation, namely between the X-ray afterglow luminosity $L_X$ and the break time $T_a$, which is observed in the X-ray flux FX. Although the number of GRBs with high redshifts is rather small, and despite the notable dispersion found in the data, the results we have obtained are quite encouraging and promising. The results obtained using the Amati relation are close to those obtained using the Type Ia supernovae, and they appear to indicate a universe dominated by dark energy. However, those obtained with the correlation between the break time and the X-ray afterglow luminosity is consistent with the findings of the WMAP study of the cosmic microwave background radiation, and they seem to indicate a de Sitter-Einstein universe dominated by matter.

Studies of \textit{BATSE} bursts \citep{kouveliotou:93} have resulted in the widespread adoption of a two-group categorization: long bursts (those with durations $\geq 2$ seconds) and short bursts (those with durations $\leq 2$ seconds). This categorization, one must recall, used the observed $T_{90}$ time durations for bursts (during which 90\% of a burst's fluence is measured). In this work, we have explored two ideas: 1) a statistical search for a possible third, intermediate category of bursts (between the "short" and the "long" ones) among 2041 \textit{BATSE} GRBs and 757 \textit{Swift/BAT} ones, 2) a study of bursts' intrinsic durations, where durations in the bursts' reference frames (instead of the observed durations) are considered, for this, 248 \textit{Swift/BAT} bursts that have redshift measurements were statistically analyzed for the same categorization goal. We first use a Monte Carlo method to determine the proper binning of each GRB, considering that bursts come with different uncertainties on their durations. Then, using the method of minimization of chi-square $\chi^2$, we search for the best fit of the normalized frequency distributions $\frac{1}{N_0}\frac{dN}{d\ln{T}}$ of durations, this allows us to compare fits with two groups ("short" and "long") with fits with three groups ("short", "long", and "intermediate"). Our results indicate that the distributions of observed durations are better fitted by three groups than two groups for \textit{Swift/BAT} data, interestingly, the "intermediate" group appears rather clearly for both observed and intrinsic durations. For BATSE data, the statistical test does not prefer three groups over two. We discuss the results, their possible underlying causes, and reasonable interpretations.

We use a sample of \textit{Swift} gamma-ray bursts (GRBs) to analyze the Amati and Yonetoku correlations. The first relation is between $E_{p,i}$, the intrinsic peak energy of the prompt GRB emission, and $E_{iso}$, the equivalent isotropic energy. The second relation is between $E_{p,i}$ and $L_{iso}$, the isotropic peak luminosity. We select a sample of 71 \textit{Swift} GRBs that have a measured redshift and whose observed $E^{obs}_p$ is within the interval of energy 15-150 keV with a relative uncertainty of less than 70\%. We seek to find correlation relations for long-duration GRBs (LGRBs) with a peak photon flux $P_{ph}\geq 2.6~ \mathrm{ph/cm^{2}/s}$. Uncertainties (error bars) on the values of the calculated energy flux \textit{P}, the energy $E_{iso}$, and the peak isotropic luminosity $L_{iso}$ are estimated using a Monte Carlo approach. We find 27 \textit{Swift} LGRBs that satisfy all our constraints. Results of our analyses of the sample of 71 GRBs and the selected subsample (27 GRBs) are in good agreement with published results. The plots of the two relations for all bursts show a large dispersion around the best straight lines in the sample of 71 LGRBs but not so much in the subsample of 27 GRBs.

We compute the expected luminosity function of GRBs in the context of the internal shock model. We assume that GRB central engines generate relativistic outflows characterized by the respective distributions of injected kinetic power Edot and contrast in Lorentz factor Kappa = Gamma_max/Gamma_min. We find that if the distribution of contrast extends down to values close to unity (i.e. if both highly variable and smooth outflows can exist) the luminosity function has two branches. At high luminosity it follows the distribution of Edot while at low luminosity it is close to a power law of slope -0.5. We then examine if existing data can constrain the luminosity function. Using the log N - log P curve, the Ep distribution of bright BATSE bursts and the XRF/GRB ratio obtained by HETE2 we show that single and broken power-laws can provide equally good fits of these data. Present observations are therefore unable to favor one form of the other. However when a broken power-law is adopted they clearly indicate a low luminosity slope ~ -0.6 +- 0.2, compatible with the prediction of the internal shock model.