Fermi satellite discovered that cosmological gamma-ray bursts (GRBs) are accompanied by long GeV flashes. In two GRBs, an optical counterpart of the GeV flash has been detected. Recent work suggests that the GeV+optical flash is emitted by the external blast wave from the explosion in a medium loaded with copious $e^\pm$ pairs. The full light curve of the flash is predicted by a first-principle radiative transfer simulation and can be tested against observations. Here we examine a sample of 7 bursts with best GeV+optical data and test the model. We find that the observed light curves are in agreement with the theoretical predictions and allow us to measure three parameters for each burst: the Lorentz factor of the explosion, its isotropic kinetic energy, and the external density. With one possible exception of GRB 090510 (which is the only short burst in the sample) the ambient medium is consistent with a wind from a Wolf-Rayet progenitor. The wind density parameter $A=\rho r^2$ varies in the sample around $10^{11}$g/cm. The initial Lorentz factor of the blast wave varies from 200 to 540 and correlates with the burst luminosity. Radiative efficiency of the prompt emission in the sample is between 0.1 and 0.8. For the two bursts with detected optical flash, GRB 120711A and GRB 130427A, we also estimate the magnetization of the external blast wave. Remarkably, the model reproduces the entire optical light curve of GRB 120711A (with its sharp peak, fast decay, plateau, and break) as well as the GeV data. The spectrum of GeV flashes is predicted to extend above 0.1 TeV, where they can be detected by ground-based Cherenkov telescopes.

The short GRB 120323A had the highest flux ever detected with the Fermi/GBM. Here we study its remarkable spectral properties and their evolution using two spectral models: (i) a single emission component scenario, where the spectrum is modeled by the empirical Band function, and (ii) a two component scenario, where thermal (Planck-like) emission is observed simultaneously with a non-thermal component (a Band function). We find that the latter model fits the integrated burst spectrum significantly better than the former, and that their respective spectral parameters are dramatically different: when fit with a Band function only, the Epeak of the event is unusually soft for a short GRB, while adding a thermal component leads to more typical short GRB values. Our time-resolved spectral analysis produces similar results. We argue here that the two-component model is the preferred interpretation for GRB 120323A, based on: (i) the values and evolution of the Band function parameters of the two component scenario, which are more typical for a short GRB, and (ii) the appearance in the data of a significant hardness-intensity correlation, commonly found in GRBs, when we employee two-component model fits; the correlation is non-existent in the Band-only fits. GRB 110721A, a long burst with an intense photospheric emission, exhibits the exact same behavior. We conclude that GRB 120323A has a strong photospheric emission contribution, first time observed in a short GRB. Magnetic dissipation models are difficult to reconcile with these results, which instead favor photospheric thermal emission and fast cooling synchrotron radiation from internal shocks. Finally, we derive a possibly universal hardness-luminosity relation in the source frame using a larger set of GRBs L,i=(1.59+/-0.84).10^50 (Epeak,i)^(1.33+/-0.07) erg/s), which could be used as a possible redshift estimator for cosmology.

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 outflow 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 dissipating 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 imply huge bulk Lorentz factors to avoid a large gamma gamma optical depth at high energy. Estimates can be as high as Gamma ~ 1000 in the most extreme cases. This puts severe constraints on models of the central engine and the jet acceleration in GRBs. These estimates are however obtained from a simplified single zone model. We present here a more realistic calculation which takes into account the time, space and direction dependent photon field existing in an outflow with several relativistically moving emitting zones. The formalism is general and can be applied to many models of the prompt GRB emission. We present results obtained for a numerical implementation in the framework of the internal shock model. We show that (i) the minimum Lorentz factor Gamma_min in bright LAT GRBs is reduced by a factor ~ 2-3 compared to previous estimates if the GeV and MeV emission are produced in the same region, and by an additional factor ~ 2-8 if the GeV emission is produced at larger radii. We provide an improved approximate formula for Gamma_min which is in good agreement with our numerical results and can be directly applied to LAT GRB data; (ii) a delayed GeV onset can be due to the time evolution of the opacity. As an illustration of these two first results, we present a synthetic GRB that reproduces most features of GRB 080916C with a mean Lorentz factor of ~ 340, an optically thin regime for gamma gamma opacity at 3GeV in bin 'b', a variability timescale of 0.5s in the MeV lightcurve and a delayed onset of ~ 5s of the GeV emission; (iii) the gamma gamma opacity can smooth the short timescale variability in the GeV lightcurve. This last result implies that the observed variability at high energy is not necessarily a good test to distinguish between an internal and an external origin for the GeV emission in GRBs. [abridged]

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.

A few long gamma-ray bursts such as GRB 050421 show no afterglow emission beyond the usual initial steep decay phase. It has been suggested that these events correspond to "naked" bursts that occur in a very low density environment. We reconsider this possibility in the context of various scenarios for the origin of the afterglow. In the standard model where the afterglow results from the forward shock as well as in the alternative model where the afterglow comes from the reverse shock, we aim to obtain constraints on the density of the environment, the microphysics parameters, or the Lorentz factor of the ejecta, which are imposed by the absence of a detected afterglow. For the two models we compute the afterglow evolution for different values of the external density (uniform or wind medium) and various burst parameters. We then compare our results to the Swift data of GRB 050421, which is the best example of a long burst without afterglow. In the standard model we show that consistency with the data imposes that the external density does not exceed 1E-5 cm-3 or that the microphysics parameters are very small with epsilon_e <~ 1E-2 and epsilon_B <~ 1E-4. If the afterglow is caused by the reverse shock, we find that its contribution can be strongly reduced if the central source has mainly emitted fast-moving material (with less than 10 - 30 % of the kinetic energy at Gamma<100 and was located in a dense environment. The two considered scenarios therefore lead to opposite constraints on the circumburst medium. The high-density environment, favored by the reverse shock model, better corresponds to what is expected if the burst progenitor was a massive star.

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

The huge optical brightness of GRB 080319B (the "Naked Eye Burst) makes this event really challenging for models of the prompt GRB emission. In the framework of the internal shock model, we investigate a scenario where the dominant radiative process is synchrotron emission and the high optical flux is due to the dynamical properties of the relativistic outflow: if the initial Lorentz factor distribution in the jet is highly variable, many internal shocks will form within the outflow at various radii. The most violent shocks will produce the main gamma-ray component while the less violent ones will contribute at lower energy, including the optical range.