While the Band function or other phenomenological functions are commonly used to fit GRB prompt emission spectra, we propose a new parametric function that is inspired by an internal shock physical model. We use this function as a proxy of the model
to confront it easily to GRB observations. We built a parametric function that represents the spectral form of the synthetic bursts provided by our internal shock synchrotron model (ISSM). We simulated the response of the Fermi instruments to the synthetic bursts and fitted the obtained count spectra to validate the ISSM function. Then, we applied this function to a sample of 74 bright GRBs detected by the Fermi/GBM, and we computed the width of their spectral energy distributions around their peak energy. For comparison, we fitted also the phenomenological functions that are commonly used in the literature. Finally, we performed a time-resolved analysis of the broadband spectrum of GRB 090926A, which was jointly detected by the Fermi GBM and LAT. The ISSM function reproduces 81% of the spectra in the GBM bright GRB sample, versus 59% for the Band function, for the same number of parameters. It gives also relatively good fits to the GRB 090926A spectra. The width of the MeV spectral component that is obtained from the fits of the ISSM function is slightly larger than the width from the Band fits, but it is smaller when observed over a wider energy range. Moreover, all of the 74 analysed spectra are found to be significantly wider than the synthetic synchrotron spectra. We discuss possible solutions to reconcile the observations with the internal shock synchrotron model, such as an improved modeling of the shock micro-physics or more accurate spectral measurements at MeV energies.
Bright X-ray flares are routinely detected by the Swift satellite during the early afterglow of gamma-ray bursts, when the explosion ejecta drives a blast wave into the external medium. We suggest that the flares are produced as the reverse shock pro
pagates into the tail of the ejecta. The ejecta is expected to contain a few dense shells formed at an earlier stage of the explosion. We show an example of how such dense shells form and describe how the reverse shock interacts with them. A new reflected shock is generated in this interaction, which produces a short-lived X-ray flare. The model provides a natural explanation for the main observed features of the X-ray flares --- the fast rise, the steep power-law decline, and the characteristic peak duration Delta t /t= (0.1-0.3).
The cosmic evolution of the neutron star merger (NSM) rate can be deduced from the observed cosmic star formation rate. This allows to estimate the rate expected in the horizon of the gravitational wave detectors advanced Virgo and ad LIGO and to com
pare those rates with independent predictions. In this context, the rapid neutron capture process, or r process, can be used as a constraint assuming NSM is the main astrophysical site for this nucleosynthetic process. We compute the early cosmic evolution of a typical r process element, Europium. Eu yields from NSM are taken from recent nucleosynthesis calculations. The same approach allows to compute the cosmic rate of Core Collapse SuperNovae (CCSN) and the associated evolution of Eu. We find that the bulk of Eu observations at high iron abundance can be rather well fitted by either CCSN or NSM scenarios. However, at lower metallicity, the early Eu cosmic evolution favors NSM as the main astrophysical site for the r process. A comparison between our calculations and spectroscopic observations at very low metallicities allows to constrain the coalescence timescale in the NSM scenario to about 0.1 to 0.2 Gyr. These values are in agreement with the coalescence timescales of some observed binary pulsars. Finally, the cosmic evolution of Eu is used to put constraints on the NSM rate, the merger rate in the horizon of the gravitational wave detectors advanced Virgo/ad LIGO, as well as the expected rate of electromagnetic counterparts to mergers (kilonovae) in large near-infrared surveys.
The early X-ray afterglow of gamma-ray bursts revealed by Swift carried many surprises. We focus in this paper on the plateau phase whose origin remains highly debated. We confront several newly discovered correlations between prompt and afterglow qu
antities (isotropic emitted energy in gamma-rays, luminosity and duration of the plateau) to several models proposed for the origin of plateaus in order to check if they can account for these observed correlations. We first show that the scenario of plateau formation by energy injection into the forward shock leads to an efficiency crisis for the prompt phase and therefore study two possible alternatives: the first one still takes place within the framework of the standard forward shock model but allows for a variation of the microphysics parameters to reduce the radiative efficiency at early times; in the second scenario the early afterglow results from a long-lived reverse shock. Its shape then depends on the distribution of energy as a function of Lorentz factor in the ejecta. In both cases, we first present simple analytical estimates of the plateau luminosity and duration and then compute detailed light curves. In the two considered scenarios we find that plateaus following the observed correlations can be obtained under the condition that specific additional ingredients are included. In the forward shock scenario, the preferred model supposes a wind external medium and a microphysics parameter epsilon_e that first varies as n^{- u} (n being the external density), with u~1 to get a flat plateau, before staying constant below a critical density n_0. To produce a plateau in the reverse shock scenario the ejecta must contain a tail of low Lorentz factor with a peak of energy deposition at Gamma >~ 10.
We study the observed correlations between the duration and luminosity of the early afterglow plateau and the isotropic gamma-ray energy release during the prompt phase. We discuss these correlations in the context of two scenarios for the origin of
the plateaus. In the first one the afterglow is made by the forward shock and the plateau results from variations of the microphysics parameters while in the second one the early afterglow is made by a long-lived reverse shock propagating in a low Lorentz factor tail of the ejecta.
The peak time of optical afterglow may be used as a proxy to constrain the Lorentz factor Gamma of the gamma-ray burst (GRB) ejecta. We revisit this method by including bursts with optical observations that started when the afterglow flux was already
decaying; these bursts can provide useful lower limits on Gamma. Combining all analyzed bursts in our sample, we find that the previously reported correlation between Gamma and the burst luminosity L_gamma does not hold. However, the data clearly shows a lower bound Gamma_min which increases with L_gamma. We suggest an explanation for this feature: explosions with large jet luminosities and Gamma < Gamma_min suffer strong adiabatic cooling before their radiation is released at the photosphere; they produce weak bursts, barely detectable with present instruments. To test this explanation we examine the effect of adiabatic cooling on the GRB location in the L_gamma - Gamma plane using a Monte Carlo simulation of the GRB population. Our results predict detectable on-axis orphan afterglows. We also derive upper limits on the density of the ambient medium that decelerates the explosion ejecta. We find that the density in many cases is smaller than expected for stellar winds from normal Wolf-Rayet progenitors. The burst progenitors may be peculiar massive stars with weaker winds or there might exist a mechanism that reduces the stellar wind a few years before the explosion.
GRB spectra appear non-thermal, but recent observations of a few bursts with Fermi GBM have confirmed previous indications from BATSE of the presence of an underlying thermal component. Photospheric emission is indeed expected when the relativistic o
utflow emerging from the central engine becomes transparent to its own radiation, with a quasi-blackbody spectrum in absence of additional sub-photospheric dissipation. However, its intensity strongly depends on the acceleration mechanism - thermal or magnetic - of the flow. We aim to compute the thermal and non-thermal emissions produced by an outflow with a variable Lorentz factor, where the power injected at the origin is partially thermal (fraction epsilon_th) and partially magnetic (fraction 1-epsilon_th). The thermal emission is produced at the photosphere, and the non-thermal emission in the optically thin regime. Apart from the value of epsilon_th, we want to test how the other model parameters affect the observed ratio of the thermal to non-thermal emission. If the non-thermal emission is made by internal shocks, we self-consistently obtained the light curves and spectra of the thermal and non-thermal components for any distribution of the Lorentz factor in the flow. If the non-thermal emission results from magnetic reconnection we were unable to produce a light curve and could only compare the respective non-thermal and thermal spectra. In the different considered cases, we varied the model parameters to see when the thermal component in the light curve and/or spectrum is likely to show up or, on the contrary, to be hidden. We finally compared our results to the proposed evidence for the presence of a thermal component in GRB spectra. Focussing on GRB 090902B and GRB 10072B, we showed how these observations can be used to constrain the nature and acceleration mechanism of GRB outflows.
The Swift-XRT observations of the early X-ray afterglow of GRBs show that it usually begins with a steep decay phase. A possible origin of this early steep decay is the high latitude emission that subsists when the on-axis emission of the last dissip
ating regions in the relativistic outflow has been switched-off. We wish to establish which of various models of the prompt emission are compatible with this interpretation. We successively consider internal shocks, photospheric emission, and magnetic reconnection and obtain the typical decay timescale at the end of the prompt phase in each case. Only internal shocks naturally predict a decay timescale comparable to the burst duration, as required to explain XRT observations in terms of high latitude emission. The decay timescale of the high latitude emission is much too short in photospheric models and the observed decay must then correspond to an effective and generic behavior of the central engine. Reconnection models require some ad hoc assumptions to agree with the data, which will have to be validated when a better description of the reconnection process becomes available.
GRB 080503, detected by Swift, belongs to the class of bursts whose prompt phase consists of an initial short spike followed by a longer soft tail. It did not show any transition to a regular afterglow at the end of the prompt emission but exhibited
a surprising rebrightening after one day. We aim to explain this rebrightening with two different scenarios - refreshed shocks or a density clump in the circumburst medium - and two models for the origin of the afterglow, the standard one where it comes from the forward shock, and an alternative one where it results from a long-lived reverse shock. We computed afterglow light curves either using a single-zone approximation for the shocked region or a detailed multizone method that more accurately accounts for the compression of the material. We find that in several of the considered cases the detailed model must be used to obtain a reliable description of the shock dynamics. The density clump scenario is not favored. We confirm previous results that the presence of the clump has little effect on the forward shock emission, except if the microphysics parameters evolve when the shock enters the clump. Moreover, we find that the rebrightening from the reverse shock is also too weak when it is calculated with the multi-zone method. On the other hand, in the refreshed-shock scenario both the forward and reverse shock models provide satisfactory fits of the data under some additional conditions on the distribution of the Lorentz factor in the ejecta and the beaming angle of the relativistic outflow.
Recent detections of GeV photons in a few GRBs by Fermi-LAT have led to strong constraints on the bulk Lorentz factor in GRB outflows. To avoid a large gamma-gamma optical depth, minimum values of the Lorentz factor have been estimated to be as high
as 800-1200 in some bursts. Here we present a detailed calculation of the gamma-gamma optical depth taking into account both the geometry and the dynamics of the jet. In the framework of the internal shock model, we compute lightcurves in different energy bands and the corresponding spectrum and we show how the limits on the Lorentz factor can be significantly lowered compared to previous estimates. Our detailed model of the propagation of high energy photons in GRB outflows is also appropriate to study many other consequences of gamma-gamma annihilation in GRBs: (i) the gamma-gamma cutoff transition in a time-integrated spectrum is expected to be closer to a power-law steepening of the spectrum than to a sharp exponential decay; (ii) the temporal evolution of the gamma-gamma opacity during a burst favors a delay between the MeV and GeV light curves; (iii) for complex GRBs, the gamma-gamma opacity suppresses the shortest time-scale features in high energy light curves (above 100 MeV). Finally we also consider GRB scenarii where MeV and GeV photons are not produced at the same location, showing that the gamma-gamma opacity could be further lowered, reducing even more the constraint on the minimum Lorentz factor.
