Late activity of the central engine is often invoked in order to explain the flares observed in the early X-ray afterglow of gamma-ray bursts, either in the form of an active neutron star remnant or (fall-back) accretion onto a black hole. However, these scenarios are not always plausible, in particular when flares are delayed to very late times after the burst. Recently, a new scenario was proposed that suggests X-ray flares can be the result of the passing of a long-lived reverse shock through a stratified ejecta, with the advantage that it does not require late-time engine activity. In this work, we numerically demonstrate this scenario to be physically plausible, by performing onedimensional simulations of ejecta dynamics and emission using our novel moving-mesh relativistic hydrodynamics code. Improved efficiency and precision over previous work enables the exploration of a broader range of setups. We can introduce a more physically realistic description of the circumburst medium mass density. We can also locally trace the cooling of electrons when computing the broadband emission from these setups. We show that the synchrotron cooling timescale can dominate the flare decay time if the stratification in the ejecta is constrained to a localised angular region inside the jet, with size corresponding to the relativistic causal connection angle, and that it corresponds to values reported in observations. We demonstrate that this scenario can produce a large range of observed flare times, suggesting a connection between flares and initial ejection dynamics rather than with late-time remnant activity.

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.

Using multiple observational arguments, recent work has shown that cosmological GRBs are typically viewed at angles within, or close to the cores of their relativistic jets. One of those arguments relied on the lack of tens-of-days-long periods of very shallow evolution that would be seen in the afterglow light-curves of GRBs viewed at large angles. Motivated by these results, we consider that GRBs efficiently produce $\gamma$-rays only within a narrow region around the core. We show that, on these near-core lines-of-sight, structured jets naturally produce shallow phases in the X-ray afterglow of GRBs. These plateaus would be seen by a large fraction of observers and would last between $10^2-10^5$ s. They naturally reproduce the observed distributions of time-scales and luminosities as well as the inter-correlations between plateau duration, plateau luminosity and prompt $\gamma$-ray energy. An advantage of this interpretation is that it involves no late time energy injection which would be both challenging from the point of view of the central engine and, as we show here, less natural given the observed correlations between plateau and prompt properties.

The up-to-now only binary neutron star merger gravitational wave event with detected electromagnetic counterparts, GRB170817A, occurred in a rarefied medium with a density smaller than $10^{-3}-10^{-2}~{\rm cm}^{-3}$. As neutron star binaries are imparted kicks upon formation, and due to their long delay-times before merger, such low-density circum-merger media are generally expected. However, there is some indirect evidence for a class of faster-merging binaries, which would coalesce in denser environments. Nonetheless, astronomical data is largely inconclusive on the possibility of these high-density mergers. We describe a method to directly probe this hypothetical population of high-density mergers through multi-messenger observations of binary neutron star merger afterglows, exploiting the sharp sensitivity of these to the circum-merger medium density. This method is based on a yet-to-be-collected sample of merger afterglows. Its constraining power is large even with a small sample of events. We discuss the method's limitations and applicability. In the upcoming era of 3$^{\rm rd}$-generation gravitational wave detectors, this method's potential will be fully realized as it will allow to probe mergers having occurred soon after the peak of cosmic star formation.

Following the historical observations of GW170817 and its multi-wavelength afterglow, more radio afterglows from neutron star mergers are expected in the future as counterparts to gravitational wave inspiral signals. We wish to describe these events using our current knowledge of the population of neutron star mergers coming from gamma-ray burst science, and taking into account the sensitivities of current and future gravitational wave and radio detectors. We combine analytical models for the merger gravitational wave and radio afterglow signals to a population model prescribing the energetics, circum-merger density and other relevant parameters of the mergers. We report the expected distributions of observables (distance, orientation, afterglow peak time/flux, etc.) from future events and study how these can be used to further probe the population of binary neutron stars, their mergers and related outflows during future observing campaigns. In the case of the O3 run of the LIGO-Virgo Collaboration, the radio afterglow of one third of gravitational-wave-detected mergers should be detectable (and detected if the source is localized thanks to the kilonova counterpart) by the Very Large Array, and these events should have viewing angles similar to that of GW170817. These findings confirm the radio afterglow as a powerful insight on these events, though some key afterglow-related techniques, such as very long baseline interferometry imaging of the merger remnant, may no longer be feasible as the gravitational wave horizon increases.

Following the historical observation of GW170817 and its electromagnetic follow-up, new neutron star merger afterglows are expected to be observed as counterparts to gravitational wave signals during the next science runs of the gravitational interferometer network. The diversity of the observed population of afterglows of these future events is subject to various factors, which are (i) intrinsic, such as the energy of the ejecta, (ii) environmental, such as the ambient medium density or (iii) observational, such as the viewing angle and distance of the source. Through prescribing a population of mergers and modelling their afterglows, we study the diversity of those events to be observed jointly in gravitational waves and electromagnetic bands. In the future, observables of detected events such as viewing angle, distance, afterglow peak flux or proper motion will form distributions which together with predictions from our study will provide insight on neutron star mergers and their environments.

The X-ray light curves of the early afterglow phase from gamma-ray bursts present a puzzling variability, including flares. The origin of these flares is still debated, and often associated with a late activity of the central engine. We discuss an alternative scenario where the central engine remains short-lived and flares are produced by the propagation of a long-lived reverse shock in a stratified ejecta. Here we focus on the hydrodynamics of the shock interactions. We perform one-dimensional ultrarelativistic hydrodynamic simulations with different initial internal structure in the gamma-ray burst ejecta. We use them to extract bolometric light curves and compare with a previous study based on a simplified ballistic model. We find a good agreement between both approaches, with similar slopes and variability in the light curves, but identify several weaknesses in the ballistic model: the density is underestimated in the shocked regions, and more importantly, late shock reflections are not captured. With accurate dynamics provided by our hydrodynamic simulations, we confirm that internal shocks in the ejecta lead to the formation of dense shells. The interaction of the long-lived reverse shock with a dense shell then produces a fast and intense increase of the dissipated power. Assuming that the emission is due to the synchrotron radiation from shock-accelerated electrons, and that the external forward shock is radiatively inefficient, we find that this results in a bright flare in the X-ray lightcurve, with arrival times, shapes, and duration in agreement with the observed properties of X-ray flares in GRB afterglows.

The most energetic gamma-ray bursts (GRBs) are remarkable sources releasing huge amounts of energy on short timescales. Their prompt emission, which usually lasts few seconds, is so bright that it is visible across the whole observable universe. Studying these extreme events may provide clues on the nature of GRB progenitors and on the physical processes at work in relativistic jets. In this paper, we study the bright end of the isotropic energy distribution of long GRBs. We use two samples of long GRBs with redshift detected by Fermi/GBM or Konus-Wind, two instruments which measure the spectral shape and the energetics of the prompt emission accurately. We focus on GRBs within a range of redshifts z = 1 -- 5, a volume that contains a large number of energetic GRBs, and we propose a simple method to reconstruct the bright end of the GRB energy distribution from the observed one. We find that the GRB energy distribution cannot be described by a simple power law but requires a strong cutoff above $1-3 \times 10^{54}$ erg. We attribute this feature to an intrinsic limit on the energy per unit of solid angle radiated by gamma-ray bursts.

ThunderKAT is the image-plane transients programme for MeerKAT. The goal as outlined in 2010, and still today, is to find, identify and understand high-energy astrophysical processes via their radio emission (often in concert with observations at other wavelengths). Through a comprehensive and complementary programme of surveying and monitoring Galactic synchrotron transients (across a range of compact accretors and a range of other explosive phenomena) and exploring distinct populations of extragalactic synchrotron transients (microquasars, supernovae and possibly yet unknown transient phenomena) - both from direct surveys and commensal observations - we will revolutionise our understanding of the dynamic and explosive transient radio sky. As well as performing targeted programmes of our own, we have made agreements with the other MeerKAT large survey projects (LSPs) that we will also search their data for transients. This commensal use of the other surveys, which remains one of our key programme goals in 2016, means that the combined MeerKAT LSPs will produce by far the largest GHz-frequency radio transient programme to date.

We report on a new time-resolved analysis of the GRB 090926A broad-band spectrum during its prompt phase, and on its interpretation in the framework of prompt emission models. We characterize the high-energy extra component of GRB 090926A using Pass 8 data from the Fermi Large Area Telescope (LAT), which offer a greater sensitivity than any data set used in previous studies of this burst, particularly in the 30-100 MeV energy band. Then, we combine the LAT data with the Fermi Gamma-ray Burst Monitor (GBM) in joint spectral fits in order to characterize the time evolution of the broad-band spectrum from keV to GeV energies. Careful attention is paid to the LAT systematic effects. Finally, we perform a temporal analysis of the light curves and we compute the variability time scales during and after the bright spike at ~10 s post-trigger. Our analysis confirms and better constrains the spectral break in the high-energy extra component, which has been previously reported during the bright spike. Furthermore, it reveals that the spectral attenuation persists at later times, with an increase of the break characteristic energy up to the GeV domain until the end of the prompt phase. We interpret the high-energy spectral break as caused by photon opacity to pair creation. Requiring that all emissions are produced above the photosphere of GRB 090926A, we compute the bulk Lorentz factor of the outflow, $\Gamma$. The latter decreases from 230 during the spike to 100 at the end of the prompt emission. Assuming, instead, that the spectral break reflects the natural curvature of the inverse Compton spectrum, lower limits corresponding to larger values of $\Gamma$ are also derived. Combined with the extreme temporal variability of GRB 090926A, these Lorentz factors lead to emission radii $R\sim10^{14}$ cm which are consistent with an internal origin of both the keV-MeV and GeV prompt emissions.

To take advantage of the astrophysical potential of Gamma-Ray Bursts (GRBs), Chinese and French astrophysicists have engaged the SVOM mission (Space-based multi-band astronomical Variable Objects Monitor). Major advances in GRB studies resulting from the synergy between space and ground observations, the SVOM mission implements space and ground instrumentation. The scientific objectives of the mission put a special emphasis on two categories of GRBs: very distant GRBs at z$>$5 which constitute exceptional cosmological probes, and faint/soft nearby GRBs which allow probing the nature of the progenitors and the physics at work in the explosion. These goals have a major impact on the design of the mission: the on-board hard X-ray imager is sensitive down to 4 keV and computes on line image and rate triggers, and the follow-up telescopes on the ground are sensitive in the NIR. At the beginning of the next decade, SVOM will be the main provider of GRB positions and spectral parameters on very short time scale. The SVOM instruments will operate simultaneously with a wide range of powerful astronomical devices. This rare instrumental conjunction, combined with the relevance of the scientific topics connected with GRB studies, warrants a remarkable scientific return for SVOM. In addition, the SVOM instrumentation, primarily designed for GRB studies, composes a unique multi-wavelength observatory with rapid slew capability that will find multiple applications for the whole astronomy community beyond the specific objectives linked to GRBs. This report lists the scientific themes that will benefit from observations made with SVOM, whether they are specific GRB topics, or more generally all the issues that can take advantage of the multi-wavelength capabilities of SVOM.

We briefly present the science capabilities, the instruments, the operations, and the expected performance of the SVOM mission. SVOM (Space-based multiband astronomical Variable Objects Monitor) is a Chinese-French space mission dedicated to the study of Gamma-Ray Bursts (GRBs) in the next decade. The SVOM mission encompasses a satellite carrying four instruments to detect and localize the prompt GRB emission and measure the evolution of the afterglow in the visible band and in X-rays, a VHF communication system enabling the fast transmission of SVOM alerts to the ground, and a ground segment including a wide angle camera and two follow-up telescopes. The pointing strategy of the satellite has been optimized to favor the detection of GRBs located in the night hemisphere. This strategy enables the study of the optical emission in the first minutes after the GRB with robotic observatories and the early spectroscopy of the optical afterglow with large telescopes to measure the redshifts. The study of GRBs in the next decade will benefit from a number of large facilities in all wavelengths that will contribute to increase the scientific return of the mission. Finally, SVOM will operate in the era of the next generation of gravitational wave detectors, greatly contributing to searches for the electromagnetic counterparts of gravitational wave triggers at Xray and gamma-ray energies.

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

Gamma-Ray Bursts (GRBs) are the strongest explosions in the Universe, which due to their extreme character likely involve some of the strongest magnetic fields in nature. This review discusses the possible roles of magnetic fields in GRBs, from their central engines, through the launching, acceleration and collimation of their ultra-relativistic jets, to the dissipation and particle acceleration that power their $\gamma$-ray emission, and the powerful blast wave they drive into the surrounding medium that generates their long-lived afterglow emission. An emphasis is put on particular areas in which there have been interesting developments in recent years.

Evidence has been accumulated on the existence of a thermal-like component during the prompt phase of GRBs. This component, often associated with the GRB jet's photosphere, is usually subdominant compared to a much stronger non-thermal one. The prompt emission of Fermi GRB 131014A provides a unique opportunity to study this thermal-like component. Indeed, the thermal emission in GRB 131014A is much more intense than in other GRBs and a pure thermal episode is observed during the initial 0.16 s. The thermal-like component cools monotonically during the first second while the non-thermal emission kicks off. The intensity of the non-thermal component progressively increases until being energetically dominant at late time. This is a perfect scenario to disentangle the thermal component from the non-thermal one. A low-energy spectral index of +0.6 better fit the thermal component than the typical index value +1 corresponding to a pure Planck function. The non-thermal component is adequately fitted with a Band function whose low and high energy power law indices are ~-0.7 and <~-3, respectively; this is also statistically equivalent to a cutoff power law with a ~-0.7 index. This is in agreement with our previous results. Finally, a strong correlation is observed between the time-resolved luminosity of the non-thermal component, L$_i^{nTh}$, and its corresponding rest frame spectral peak energy, E$_{peak,i}^{rest,nTh}$, with a slope similar to the one reported in our previous articles. Assuming this relation to be universal for all GRBs we estimate a redshift of ~1.55 for GRB 131014A that is a typical value for long GRBs. These observational results are consistent with the models in which the non-thermal emission is produced well above the GRB jet photosphere but they may also be compatible with other scenarios (e.g., dissipative photosphere) that are not discussed in this article.

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 propagates 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).

We reanalyze the prompt emission of two of the brightest Fermi GRBs (080916C and 090926A) with a new model composed of 3 components: (i) a thermal-like component--approximated with a black body (BB)--interpreted as the jet photosphere emission of a magnetized relativistic outflow, (ii) a non-thermal component--approximated with a Band function--interpreted as synchrotron radiation in an optically thin region above the photosphere either from internal shocks or magnetic field dissipation, and (iii) an extra power law (PL) extending from low to high energies likely of inverse Compton origin, even though it remains challenging. Through fine-time spectroscopy down to the 100 ms time scale, we follow the smooth evolution of the various components. From this analysis the Band function is globally the most intense component, although the additional PL can overpower the others in sharp time structures. The Band function and the BB component are the most intense at early times and globally fade across the burst duration. The additional PL is the most intense component at late time and may be correlated with the extended high-energy emission observed thousands of seconds after the burst with the Fermi/Large Area Telescope (LAT). Unexpectedly, this analysis also shows that the additional PL may be present from the very beginning of the burst. We investigate the effect of the three components on the new time-resolved luminosity-hardness relation in both the observer and rest frames and show that a strong correlation exists between the flux of the non-thermal component and its E$_{peak}$ only when the three components are fitted simultaneously to the data (i.e., F$_i^{NT}$-E$_{peak,i}^{NT}$ relation). In addition, this result points toward a universal relation between those two quantities for all GRBs when transposed to the central engine rest frame (i.e., L$_i^{NT}$-E$_{peak,i}^{rest,NT}$ relation).

Several trends have been identified in the prompt gamma-ray burst (GRB) emission: e.g. hard-to-soft evolution, pulse width evolution with energy, time lags, hardness-intensity/-fluence correlations. Recently Fermi has significantly extended the spectral coverage of GRB observations and improved the characterization of this spectral evolution. We study how internal shocks can reproduce these observations. In this model the emission comes from the synchrotron radiation of shock accelerated electrons, and the spectral evolution is governed by the evolution of the physical conditions in the shocked regions. We present a comprehensive set of simulations of a single pulse and investigate the impact of the model parameters, related to the shock microphysics and to the initial conditions in the ejecta. We find a general qualitative agreement between the model and the various observations used for the comparison. All these properties or relations are governed by the evolution of the peak energy and photon indices of the spectrum. In addition, we identify the conditions for a quantitative agreement. We find that the best agreement is obtained for (i) steep electron slopes (p>~2.7), (ii) microphysics parameters varying with shock conditions so that more electrons are accelerated in stronger shocks, (iii) steep variations of the initial Lorentz factor in the ejecta. When simulating short GRBs by contracting all timescales, all other parameters being unchanged, we show that the hardness-duration correlation is reproduced, as well as the evolution with duration of the pulse properties. Finally, we investigate the signature at high energy of these different scenarios and find distinct properties - delayed onset, longer emission, and flat spectrum in some cases - suggesting that internal shocks could have a significant contribution to the prompt LAT emission. [abridged]

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 quantities (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^{-\nu} (n being the external density), with \nu~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.

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.

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.

Gamma-Ray Bursts (GRBs) are the most powerful cosmic explosions since the Big Bang, and thus act as signposts throughout the distant Universe. Over the last 2 decades, these ultra-luminous cosmological explosions have been transformed from a mere curiosity to essential tools for the study of high-redshift stars and galaxies, early structure formation and the evolution of chemical elements. In the future, GRBs will likely provide a powerful probe of the epoch of reionisation of the Universe, constrain the properties of the first generation of stars, and play an important role in the revolution of multi-messenger astronomy by associating neutrinos or gravitational wave (GW) signals with GRBs. Here, we describe the next steps needed to advance the GRB field, as well as the potential of GRBs for studying the Early Universe and their role in the up-coming multi-messenger revolution.

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.

We present a mechanism based on internal self-annihilation of dark matter accreted from the galactic halo in the inner regions of neutron stars that may trigger the full or partial conversion into a quark star. We explain how this effect may induce a gamma ray burst that could be classified as short, according to the usual definition based on time duration of the prompt gamma-ray emission. This mechanism differs in many aspects from the most discussed scenario associating short gamma-ray bursts with compact object binary mergers. We list possible observational signatures that should help to distinguish between these two possible classes of progenitors.

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.

In this contribution we briefly introduce a mechanism for short gamma ray burst emission different from the usually assumed compact object binary merger progenitor model. It is based on the energy release in the central regions of neutron stars. This energy injection may be due to internal self-annihilation of dark matter gravitationally accreted from the galactic halo. We explain how this effect may trigger its full or partial conversion into a quark star and, in such a case, induce a gamma ray burst with isotropic equivalent energies in agreement with those measured experimentally. Additionally, we show how the ejection of the outer crust in such events may be accelerated enough to produce Lorentz factors over those required for gamma ray emission.

We perform a detailed study on the dynamics of a relativistic blast wave with the presence of a long-lived reverse shock (RS). Although a short-lived RS has been widely considered, the RS is believed to be long-lived as a consequence of a stratification expected on the ejecta Lorentz factors. The existence of a long-lived RS makes the forward shock (FS) dynamics to deviate from a self-similar Blandford-McKee solution. Employing the "mechanical model" that correctly incorporates the energy conservation, we present an accurate solution for both the FS and RS dynamics. We conduct a sophisticated calculation of the afterglow emission. Adopting a Lagrangian description of the blast wave, we keep track of an adiabatic evolution of numerous shells between the FS and RS. An evolution of the electron spectrum is also followed individually for every shell. We then find the FS and RS light curves by integrating over the entire FS and RS shocked regions, respectively. In particular, we make use of an analytic expression for observed spectral flux, which we derive here in terms of an observed frequency and observer time. Exploring a total of 20 different ejecta stratifications, we explain in detail how a stratified ejecta affects its blast wave dynamics and afterglow light curves. We show that, while the FS light curves are not sensitive to the ejecta stratifications, the RS light curves exhibit much richer features, including steep declines, plateaus, bumps, re-brightenings, and a variety of temporal decay indices. These distinctive RS features may be observable if the RS has higher values of the micophysics parameters than the FS. We discuss possible applications of our results in understanding the GRB afterglow data.

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.

GRB 041219A is one of the longest and brightest gamma-ray bursts (GRBs) ever observed. It was discovered by the INTEGRAL satellite, and thanks to a precursor happening about 300 s before the bulk of the burst, ground based telescopes were able to catch the rarely-observed prompt emission in the optical and in the near infrared bands. Here we present the detailed analysis of its prompt gamma-ray emission, as observed with IBIS on board INTEGRAL, and of the available X-ray afterglow data collected by XRT on board Swift. We then present the late-time multi-band near infrared imaging data, collected at the TNG, and the CFHT, that allowed us to identify the host galaxy of the GRB as an under-luminous, irregular galaxy of about 5x10^9 M_Sun at best fit redshift of z=0.31 -0.26 +0.54. We model the broad-band prompt optical to gamma-ray emission of GRB 041219A within the internal shock model. We were able to reproduce the spectra and light curve invoking the synchrotron emission of relativistic electrons accelerated by a series of propagating shock waves inside a relativistic outflow. On the other hand, it is less easy to simultaneously reproduce the temporal and spectral properties of the infrared data.

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.

From July 2008 to October 2009, the Gamma-ray Burst Monitor (GBM) on board the Fermi Gamma-ray Space Telescope (FGST) has detected 320 Gamma-Ray Bursts (GRBs). About 20% of these events are classified as short based on their T90 duration below 2 s. We present here for the first time time-resolved spectroscopy at timescales as short as 2 ms for the three brightest short GRBs observed with GBM. The time-integrated spectra of the events deviate from the Band function, indicating the existence of an additional spectral component, which can be fit by a power-law with index ~-1.5. The time-integrated Epeak values exceed 2 MeV for two of the bursts, and are well above the values observed in the brightest long GRBs. Their Epeak values and their low-energy power-law indices ({\alpha}) confirm that short GRBs are harder than long ones. We find that short GRBs are very similar to long ones, but with light curves contracted in time and with harder spectra stretched towards higher energies. In our time-resolved spectroscopy analysis, we find that the Epeak values range from a few tens of keV up to more than 6 MeV. In general, the hardness evolutions during the bursts follows their flux/intensity variations, similar to long bursts. However, we do not always see the Epeak leading the light-curve rises, and we confirm the zero/short average light-curve spectral lag below 1 MeV, already established for short GRBs. We also find that the time-resolved low-energy power-law indices of the Band function mostly violate the limits imposed by the synchrotron models for both slow and fast electron cooling and may require additional emission processes to explain the data. Finally, we interpreted these observations in the context of the current existing models and emission mechanisms for the prompt emission of GRBs.

Observations of GRB 100724B with the Fermi Gamma-Ray Burst Monitor (GBM) find that the spectrum is dominated by the typical Band functional form, which is usually taken to represent a non-thermal emission component, but also includes a statistically highly significant thermal spectral contribution. The simultaneous observation of the thermal and non-thermal components allows us to confidently identify the two emission components. The fact that these seem to vary independently favors the idea that the thermal component is of photospheric origin while the dominant non-thermal emission occurs at larger radii. Our results imply either a very high efficiency for the non-thermal process, or a very small size of the region at the base of the flow, both quite challenging for the standard fireball model. These problems are resolved if the jet is initially highly magnetized and has a substantial Poynting flux.

(abridged)Prompt GRB emission is often interpreted as synchrotron radiation from high-energy electrons accelerated in internal shocks. Fast synchrotron cooling predicts that the photon index below the spectral peak is alpha=-3/2. This differs significantly from the observed median value alpha \approx -1. We quantify the influence of inverse Compton and adiabatic cooling on alpha to understand whether these processes can reconcile the observations with a synchrotron origin. We use a time-dependent code that follows both the shock dynamics and electron energy losses. We investigate the dependence of alpha on the parameters of the model. Slopes between -3/2 and -1 are reached when electrons suffer IC losses in the Klein-Nishina regime. This does not necessarily imply a strong IC component in the Fermi/LAT range because scatterings are only moderately efficient. Steep slopes require that a large fraction (10-30%) of the dissipated energy is given to a small fraction (<~1%) of the electrons and that the magnetic energy density fraction remains low (<~ 0.1%). Values of alpha up to -2/3 can be obtained with relatively high radiative efficiencies (>50%) when adiabatic cooling is comparable with radiative cooling (marginally fast cooling). This requires collisions at small radii and/or with low magnetic fields. Amending the standard fast cooling scenario to account for IC cooling naturally leads to alpha up to -1. Marginally fast cooling may also account for alpha up to -2/3, although the conditions required are more difficult to reach. About 20% of GRBs show spectra with slopes alpha>-2/3. Other effects, not investigated here, such as a thermal component in the electron distribution or pair production by HE photons may further affect alpha. Still, the majority of observed GRB prompt spectra can be reconciled with a synchrotron origin, constraining the microphysics of mildly relativistic internal shocks.

The curvature of a relativistic blast wave implies that its emission arrives to observers with a spread in time. This effect is believed to wash out fast variability in the lightcurves of GRB afterglows. We note that the spreading effect is reduced if emission is anisotropic in the rest-frame of the blast wave (i.e. if emission is limb-brightened or limb-darkened). In particular, synchrotron emission is almost certainly anisotropic, and may be strongly anisotropic, depending on details of electron acceleration in the blast wave. Anisotropic afterglows can display fast and strong variability at high frequencies (above the 'fast-cooling' frequency). This may explain the existence of bizarre features in the X-ray afterglows of GRBs, such as sudden drops and flares. We also note that a moderate anisotropy can significantly delay the 'jet break' in the lightcurve, which makes it harder to detect.

We present the Space-based multi-band astronomical Variable Object Monitor (SVOM), a future satellite mission for Gamma-Ray Burst (GRB) studies, developed in cooperation between the Chinese National Space Agency (CNSA), the Chinese Academy of Science (CAS), the French Space Agency (CNES) and French research institutes. The scientific objectives of the SVOM GRB studies cover their classification (GRB diversity and unity of the model), their physics (particle acceleration and radiation mechanisms), their progenitors, cosmological studies (host galaxies, star formation history, re-ionization, cosmological parameters), and fundamental physics (origin of cosmic rays, Lorentz invariance, gravitational wave sources). From 2015 on, SVOM will provide fast and accurate localizations of all known types of GRB, and determine the temporal and spectral properties of the GRB emission, thanks to a set of four onboard instruments. The trigger system of the coded-mask telescope ECLAIRs onboard SVOM images the sky in the 4-120 keV energy range, in order to detect and localize GRB in its 2 sr-wide field of view. The low-energy threshold of ECLAIRs is well suited for the detection of highly red-shifted GRB. The high-energy coverage is extended up to 5 MeV thanks to a non-imaging gamma-ray spectrometer. GRB alerts are sent in real-time to the ground observers community, and a spacecraft slew is performed in order to place the GRB within the narrow fields of view of a soft X-ray telescope and a visible-band telescope, to refine the GRB position and study its early afterglow. Ground-based robotic telescopes and wide-angle cameras complement the onboard instruments. A large fraction of GRB will have redshift determinations, thanks to an observing strategy optimized to facilitate follow-up observations by large ground-based spectroscopic telescopes.

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.

Using a detailed model of the internal shock phase, we discuss the origin of the prompt emission in gamma-ray bursts. We focus on the identification of the dominant radiative process (Fermi-GBM range) and propose an explanation for some features observed by Fermi-LAT at high energy in some GRB lightcurves.

New observations from the Hubble ultra deep field suggest that the star formation rate at z>7 drops off faster than previously thought. Using a newly determined star formation rate for the normal mode of Population II/I stars (PopII/I), including this new constraint, we compute the Thomson scattering optical depth and find a result that is marginally consistent with WMAP5 results. We also reconsider the role of Population III stars (PopIII) in light of cosmological and stellar evolution constraints. While this input may be needed for reionization, we show that it is essential in order to account for cosmic chemical evolution in the early Universe. We investigate the consequences of PopIII stars on the local metallicity distribution function of the Galactic halo (from the recent Hamburg/ESO survey of metal-poor stars) and on the evolution of abundances with metallicity (based on the ESO large program on very metal-poor stars), with special emphasis on carbon-enhanced metal-poor stars. Our most important results show that the nucleosynthetic yields of PopIII stars lead to abundance patterns in agreement with those observed in extremely metal-poor stars. In this chemical approach to cosmic evolution, PopIII stars prove to be a compulsory ingredient, and extremely metal-poor stars are inevitably born at high redshift. (Abridged)

The prompt GRB emission is thought to arise from electrons accelerated in internal shocks propagating within a highly relativistic outflow. The launch of Fermi offers the prospect of observations with unprecedented sensitivity in high-energy (>100 MeV) gamma-rays. The aim is to explore the predictions for HE emission from internal shocks, taking into account both dynamical and radiative aspects, and to deduce how HE observations constrain the properties of the relativistic outflow. The emission is modeled by combining a time-dependent radiative code with a dynamical code giving the evolution of the physical conditions in the shocked regions.Synthetic lightcurves and spectra are compared to observations. The HE emission deviates significantly from analytical estimates, which tend to overpredict the IC component, when the time dependence and full cross-sections are included. The exploration of the parameter space favors the case where the dominant process in the BATSE range is synchrotron emission. The HE component becomes stronger for weaker magnetic fields. The HE lightcurve can display a prolonged pulse duration due to IC emission, or even a delayed peak compared to the BATSE range.Alternatively, having dominant IC emission in the BATSE range requires most electrons to be accelerated into a steep power-law distribution and implies strong 2nd order IC scattering. In this case, the BATSE and HE lightcurves are very similar. The combined dynamical and radiative approach allows a firm appraisal of GRB HE prompt emission. A diagnostic procedure is presented to identify from observations the dominant emission process and derive constrains on the bulk Lorentz factor, particle density and magnetic field of the outflow.

Polarization measurements provide direct insight into the nature of astrophysical processes. Unfortunately, only a few instruments are available for this kind of measurements at gamma-ray energies, and the sources need to be very bright. Gamma-Ray Bursts (GRBs) are ideal candidates due to their large flux over limited time intervals, maximizing the available signal-to-noise ratio. To date a few polarization measurements have been reported, claiming of a high degree of polarization in the prompt emission of GRBs but with low statistical evidence. We used the IBIS telescope on board the INTEGRAL satellite to measure the polarization of the prompt gamma-ray emission of the long and bright GRB 041219A in the 200-800 keV energy band. We find a variable degree of polarization ranging from less than 4% over the first peak to 43+/-25% for the whole second peak. Time resolved analysis of both peaks indicates a high degree of polarization, and the null average polarization in the first peak can be explained by the rapid variations observed in the polarization angle and degree. Our results are consistent with different models for the prompt emission of GRBs at these energies, but they favor synchrotron radiation from a relativistic outflow with a magnetic field which is coherent on an angular size comparable with the angular size of the emitting region (~1/Gamma) . Indeed this model has the best capabilities to maintain a high polarization level, and to produce the observed variability.

The paucity of reliable achromatic breaks in Gamma-Ray Burst afterglow light curves motivates independent measurements of the jet aperture. Serendipitous searches of afterglows, especially at radio wavelengths, have long been the classic alternative. These survey data have been interpreted assuming a uniformly emitting jet with sharp edges (``top-hat'' jet), in which case the ratio of weakly relativistically beamed afterglows to GRBs scales with the jet solid angle. In this paper, we consider, instead, a very wide outflow with a luminosity that decreases across the emitting surface. In particular, we adopt the universal structured jet (USJ) model, that is an alternative to the top-hat model for the structure of the jet. However, the interpretation of the survey data is very different: in the USJ model we only observe the emission within the jet aperture and the observed ratio of prompt emission rate to afterglow rate should solely depend on selection effects. We compute the number and rate of afterglows expected in all-sky snapshot observations as a function of the survey sensitivity. We find that the current (negative) results for OA searches are in agreement with our expectations. In radio and X-ray bands this was mainly due to the low sensitivity of the surveys, while in the optical band the sky-coverage was not sufficient. In general we find that X-ray surveys are poor tools for OA searches, if the jet is structured. On the other hand, the FIRST radio survey and future instruments like the Allen Telescope Array (in the radio band) and especially GAIA, Pan-Starrs and LSST (in the optical band) will have chances to detect afterglows.

The early optical afterglow of GRB 050820a recorded by the RAPTOR telescope shows both a contribution from the prompt emission and the initial rise of the afterglow. It is therefore well-suited for the study of the dynamical evolution of the GRB ejecta when it first undergoes the decelerating effect of the environment. This is a complex phase where the internal, reverse, and forward shocks can all be present simultaneously. We have developed a simplified model that can follow these different shocks in an approximate, but self-consistent way. It is applied to the case of GRB 050820a to obtain the prompt and afterglow light curves. We show that the rise of the afterglow during the course of the prompt emission has some important consequences. The reverse shock propagates back into the ejecta before internal shocks are completed, which affects the shape of the gamma-ray profile. We get the best results when the external medium has a uniform density, but obtaining a simultaneous fit of the prompt and afterglow emission is not easy. We discuss a few possibilities that could help to improve this situation.

(abridged) The association of GRB 980425 with SN 1998bw at z=0.0085 implies the existence of a population of GRBs with an isotropic-equivalent luminosity which is about 10^4 times smaller than in the standard cosmological case. We investigate two scenarios to explain a weak GRB : a normal (intrinsically bright) GRB seen off-axis or an intrinsically weak GRB seen on-axis. For each of these two scenarios, we first derive the conditions to produce a GRB 980425-like event and we then discuss the consequences for the event rate. If we exclude the possibility that GRB 980425 is an occurence of an extremely rare event observed by chance during the first eight years of the 'afterglow era', the first scenario implies that (i) the local rate of standard bright GRBs is much higher than what is usually expected; (ii) the typical opening angle is much narrower than what is derived from observations of a break in the afterglow lightcurve. In addition, we show that the afterglow of GRB 980425 in this scenario should have been very bright and easily detected. For these reasons the second scenario appears more realistic. We show that the parameter space of the internal shock model indeed allows GRB 980425-like events, in cases where the outflow is only mildly-relativistic and mildly-energetic. The rate of such weak events in the Universe has to be much higher than the rate of standard bright GRBs to allow the discovery of GRB 980425 during a short period of a few years. However it is still compatible with the observations as the intrinsic weakness of these GRB 980425-like bursts does not allow detection at cosmological redshift. We finally briefly discuss the consequences of such a high local rate of GRB 980425-like events for the future prospects of detecting non-electromagnetic radiation, especially gravitational waves.

We propose to explain the recent observations of GRB early X-ray afterglows with SWIFT by the dissipation of energy in the reverse shock which crosses the ejecta as it is decelerated by the burst environment. We compute the evolution of the dissipated power and discuss the possibility that a fraction of it can be radiated in the X-ray range. We show that this reverse shock contribution behaves in a way very similar to the observed X-ray afterglows if the following two conditions are satisfied: (i) the Lorentz factor of the material which is ejected during the late stages of source activity decreases to small values Gamma < 10 and (ii) a large part of the shock dissipated energy is transferred to a small fraction (zeta < 0.1) of the electron population. We also discuss how our results may help to solve some puzzling problems raised by multiwavelength early afterglow observations such as the presence of chromatic breaks.

We predict the redshift distribution of long Gamma-Ray Bursts (GRBs) with Monte Carlo simulations. Our improved analysis constrains free parameters with three kinds of observation: (i) the log(N)-log(P) diagram of BATSE bursts; (ii) the peak energy distribution of bright BATSE bursts; (iii) the HETE2 fraction of X-ray rich GRBs and X-ray flashes. The statistical analysis of the Monte Carlo simulation results allow us to carefully study the impact of the uncertainties in the GRB intrinsic properties on the redshift distribution. The comparison with SWIFT data then leads to the following conclusions. The Amati relation should be intrinsic, if observationally confirmed by SWIFT. The progenitor and/or the GRB properties have to evolve to reproduce the high mean redshift of SWIFT bursts. Our results favor an evolution of the efficiency of GRB production by massive stars, that would be ~6-7 times higher at z~7 than at z~2. We finally predict around 10 GRBs detected by SWIFT at redshift z>6 for a three year mission. These may be sufficient to open a new observational window over the high redshift Universe.

We compute the afterglow of gamma-ray bursts produced by purely electromagnetic outflows to see if it shows characteristic signatures differing from those obtained with the standard internal/external shock model. Using a simple approach for the injection of electromagnetic energy to the forward shock we obtain the afterglow evolution both during the period of activity of the central source and after. Our method equally applies to a variable source. Afterglow light curves in the visible and X-ray bands are computed both for a uniform medium and a stellar wind environment. They are brighter at early times than afterglows obtained with the internal/external shock model but relying only on these differences to discriminate between models is not sufficient.

We consider the stochastic background of gravitational waves produced by an early generation of Population III stars coupled with a normal mode of star formation at lower redshift. The computation is performed in the framework of hierarchical structure formation and is based on cosmic star formation histories constrained to reproduce the observed star formation rate at redshift $z \la 6$, the observed chemical abundances in damped Lyman alpha absorbers and in the intergalactic medium, and to allow for an early reionization of the Universe at $z\sim 10-20$ as indicated by the first year results released by WMAP. We find that the normal mode of star formation produces a gravitational wave background which peaks at 300-500 Hz and is within LIGO III sensitivity. The Population III component peaks at lower frequencies (30-100 Hz depending on the model), and could be detected by LIGO III as well as the planned BBO and DECIGO interferometers.

Gamma-ray bursts (GRB), at least those with a duration longer than a few seconds are the most energetic events in the Universe and occur at cosmological distances. The ECLAIRs micro-satellite, to be launched in 2009, will provide multi-wavelength observations of GRB, to study their astrophysics and to use them as cosmological probes. Furthermore in 2009 ECLAIRs is expected to be the only space borne instrument capable of providing a GRB trigger in near real-time with sufficient localization accuracy for GRB follow-up observations with the powerful ground based spectroscopic telescopes available by then. A "Phase A study" of the ECLAIRs project has recently been launched by the French Space Agency CNES, aiming at a detailed mission design and selection for flight in 2006. The ECLAIRs mission is based on a CNES micro-satellite of the "Myriade" family and dedicated ground-based optical telescopes. The satellite payload combines a 2 sr field-of-view coded aperture mask gamma-camera using 6400 CdTe pixels for GRB detection and localization with 10 arcmin precision in the 4 to 50 keV energy band, together with a soft X-ray camera for onboard position refinement to 1 arcmin. The ground-based optical robotic telescopes will detect the GRB prompt/early afterglow emission and localize the event to arcsec accuracy, for spectroscopic follow-up observations.

We study the evolution of the circumstellar medium of massive stars. We pay particular attention to Wolf-Rayet stars that are thought to be the progenitors of some long Gamma-Ray Bursts. We detail the mass-loss rates we use in our stellar evolution models and how we estimate the stellar wind speeds during different phases. With these details we simulate the interactions between the wind and the interstellar medium to predict the circumstellar environment around the stars at the time of core-collapse. We then investigate how the structure of the environment might affect the GRB afterglow. We find that when the afterglow jet encounters the free-wind to stalled-wind interface that rebrightening occurs and a bump is seen in the afterglow light curve. However our predicted positions of this interface are too distant from the site of the GRB to reach while the afterglow remains observable. The values of the final-wind density, A_{*}, from our stellar models are similar to the values inferred from observed afterglow lightcurves and those from observed Wolf-Rayet stars. However we do not reproduce the lowest observed A_{*} values below 0.3. For these cases we suggest that the progenitors could have been a WO type Wolf-Rayet star, be in a close binary or very low metallicity star. Finally we turn our attention to the matter of stellar wind material producing absorption lines in the afterglow spectra. We discuss the observational signatures of two Wolf-Rayet stellar types, WC and WO, in the afterglow lightcurve and spectra. We also indicate how it may be possible to constrain the initial mass and metallicity of a GRB progenitor by using the inferred wind density and wind velocity.

Evidence from the WMAP polarization data indicates that the Universe may have been reionized at very high redshift. It is often suggested that the ionizing UV flux originates from an early population of massive or very massive stars. Depending on their mass, such stars can explode either as type II supernovae or pair-instability supernovae, or may entirely collapse into a black hole. The resulting neutrino emission can be quite different in each case. We consider here the relic neutrino background produced by an early burst of Population III stars coupled with a normal mode of star formation at lower redshift. The computation is performed in the framework of hierarchical structure formation and is based on cosmic star formation histories constrained to reproduce the observed star formation rate at redshift z \la 6, the observed chemical abundances in damped Lyman alpha absorbers and in the intergalactic medium, and to allow for an early reionization of the Universe at z \sim 10-20. We find that although the high redshift burst of Population III stars does lead to an appreciable flux of neutrinos at relatively low energy (E_\nu \approx 1 MeV), the observable neutrino flux is dominated by the normal mode of star formation. We also find that predicted fluxes are at the present level of the SuperK limit. As a consequence, the supernova relic neutrino background has a direct impact on models of chemical evolution and/or supernova dynamics.

The cosmic star formation histories are evaluated for different minimum masses of the initial halo structures, with allowance for realistic gas outflows. With a minimum halo mass of 10^{7} - 10^{8} M_odot and a moderate outflow efficiency, we reproduce both the current baryon fraction and the early chemical enrichment of the IGM. The intensity of the formation rate of ``normal'' stars is also well constrained by the observations: it has to be dominated by star formation in elliptical galaxies, except perhaps at very low redshift. The fraction of baryons in stars is predicted as are also the type Ia and II supernova event rates. Comparison with SN observations in the redshift range z=0-2 allows us to set strong constraints on the time delay of type Ia supernovae (a total delay of \sim 4 Gyr is required to fit the data), the lower end of the mass range of the progenitors (2 - 8 M_odot) and the fraction of white dwarfs that reproduce the type Ia supernova (about 1 per cent). The intensity in the initial starburst of zero metallicity stars below 270 M_\odot must be limited in order to avoid premature overenrichment of the IGM. Only about 10 - 20 % of the metals present in the IGM at z = 0 have been produced by population III stars at very high z. The remaining 80 - 90 % are ejected later by galaxies forming normal stars, with a maximum outflow efficiency occurring at a redshift of about 5. We conclude that 10^{-3} of the mass in baryons must lie in first massive stars in order to produce enough ionizing photons to allow early reionization of the IGM by z \sim 15.

We discuss the origin of X-Ray Flashes (XRFs), a recently discovered class of Gamma-Ray Bursts (GRBs). Using a simplified model for internal shocks we check if XRFs can be intrinsically soft due to some specific values of the parameters describing the relativistic outflow emerging from the central engine. We generate a large number of synthetic events and find that XRFs are obtained when the contrast Gamma_max/Gamma_min of the Lorentz factor distribution is small while the average Lorentz factor Gamma is large. A few XRFs may be GRBs at large redshifts but we exclude this possibility for the bulk of the population. If outflows with a small contrast are commonly produced, even a large population of XRFs could be explained. If conversely the Lorentz factor distribution within the wind is broad, one should then rely on extrinsic causes, such as viewing angle effects or high redshift.

The cosmological revolution of 1997 has established that (at least long duration) gamma-ray bursts (GRB) are among the most energetic events in the Universe and occur at cosmological distances. The ECLAIRs micro-satellite, to be launched in 2009, will provide multi-wavelength observations for astrophysical studies of GRB and for their possible use as cosmological probes. It is expected to be the only space borne GRB trigger available for ground based robotic telescopes operational at that time. This paper presents the ECLAIRs project and its status. An X/gamma-ray camera onboard ECLAIRs with a wide field of view of 2 sr, will detect ~100 GRB/yr in the 4-50 keV energy range, localize the GRB with a precision of ~10 arcmin on the sky, and transmit this information to the ground in near real-time, as a GRB trigger for ground based optical telescopes. Inspired by the INTEGRAL imager IBIS, it is based on a CdTe detection plane covering 1000 cm^2, placed 35 cm below a coded mask. An optical camera, sensitive to magnitude-15 stars, covering up to 1/4th of the X/gamma-ray camera's field of view, will observe the prompt emission and a possible precursor of ~10 GRB/yr in the visible-band. Used in a continuous acquisition mode at a rate of ~5 images/s dumped into an on-board memory, a GRB event sent by the X/gamma-ray camera triggers a seek-back in memory for the GRB optical precursor. The full X/gamma-ray and visible-band data of a GRB are sent to ground when a high data-rate telemetry ground receiver is reachable.

Motivated by the WMAP results indicating an early epoch of reionization, we consider alternative cosmic star formation models which are capable of reionizing the early intergalactic medium. We develop models which include an early burst of massive stars (with several possible mass ranges) combined with standard star formation. We compute the stellar ionizing flux of photons and we track the nucleosynthetic yields for several elements: D, He4, C, N, O, Si, S, Fe, Zn. We compute the subsequent chemical evolution as a function of redshift, both in the intergalactic medium and in the interstellar medium of forming galaxies, starting with the primordial objects which are responsible for the reionization. We apply constraints from the observed abundances in the Lyman alpha forest and in Damped Lyman alpha clouds in conjunction with the ability of the models to produce the required degree of reionization. We also consider possible constraints associated with the observations of the two extremely metal-poor stars HE 0107-5240 and CS22949-037. We confirm that an early top-heavy stellar component is required, as a standard star formation model is unable to reionize the early Universe and reproduce the abundances of the very metal-poor halo stars. A bimodal (or top-heavy) IMF (40 - 100 M_\odot) is our preferred scenario compared to the extreme mass range (\ga 100 M_\odot) often assumed to be responsible for the early stages of reionization. A mode of even more extreme stellar masses in the range (\ge 270 M_\odot) has also been considered. All massive stars in this mode collapse entirely into black holes, and as a consequence, chemical evolution and reionization are de-correlated. [Abstract abbreviated.]

We present results from a comprehensive number of relativistic, time-dependent, axisymmetric simulations of the runaway instability of non-constant angular momentum thick discs around black holes. This second paper extends earlier results where only constant angular momentum discs were considered. All relevant aspects of the theory of stationary thick discs around rotating black holes, necessary to build the initial state in our simulations, are presented in great detail. The angular momentum of the discs is assumed to increase outwards with the radial distance according to a power law. The main simplifying assumptions of our approach are not to include magnetic fields and self-gravity in the discs. Furthermore, the dynamics of the spacetime is accounted for by computing the transfer of mass and angular momentum from the disc to the black hole through the event horizon : the evolution of the central black hole is assumed to follow a sequence of Kerr black holes of increasing mass and spin. In agreement with previous results based on stationary models we find that by allowing the mass and the spin of the black hole to grow, constant angular momentum discs rapidly become unstable on a dynamical timescale. The comparison with the results of paper I shows that the effect of the angular momentum transfer from the torus to the black hole is to make constant angular momentum discs less unstable, increasing the timescale of the instability. However, we find that non-constant angular momentum discs are dramatically stabilized for very small values of the angular momentum slope. Our time-dependent simulations confirm, thus, the predictions of stationary studies concerning the stabilizing effect of non-constant angular momentum distributions.

We report the results of XMM-Newton observations of Sgr A*, the radiative counterpart of the massive black hole at the nucleus of our Galaxy, performed in the frame of the guaranteed time survey program of the Galactic Center region. The discovery of bright X-ray flares from Sgr A* with Chandra in October 2000 have opened new perspectives to understand the processes at work in this object and in general in black holes accreting at low accretion rates. We report here the important results obtained with XMM-Newton on the Sgr A* high-energy flaring activity and we discuss the implications on the models and the future observational perspectives.

The compact radio source Sgr A*, believed to be the counterpart of the massive black hole at the Galactic nucleus, was observed to undergo rapid and intense flaring activity in X-rays with Chandra in October 2000. We report here the detection with XMM-Newton EPIC cameras of the early phase of a similar X-ray flare from this source, which occurred on 2001 September 4. The source 2-10 keV luminosity increased by a factor ~ 20 to reach a level of 4 10$^{34}$ erg s$^{-1}$ in a time interval of about 900 s, just before the end of the observation. The data indicate that the source spectrum was hard during the flare and can be described by simple power law of slope ~ 0.7. This XMM-Newton observation confirms the results obtained by Chandra, suggests that, in Sgr A*, rapid and intense X-ray flaring is not a rare event and therefore sets some constraints on the emission mechanism models proposed for this source.

We present a simple, semi-analytical model to explain GRB temporal and spectral properties in the context of the internal shock model. Each individual pulse in the temporal profiles is produced by the deceleration of fast moving material by a comparatively slower layer within a relativistic wind. The spectral evolution of synthetic pulses is first obtained with standard equipartition assumptions to estimate the post-shock magnetic field and electron Lorentz factor. We get Ep propto t^-delta with delta=7/2 which is much steeper than the observed slopes delta(obs) <= 1.5. We therefore consider the possibility that the equipartition parameters depend on the shock strength and post-shock density. We then get a much better agreement with the observations and our synthetic pulses satisfy both the hardness-intensity and hardness-fluence correlations. We also compute time lags between profiles in different energy channels and we find that they decrease with increasing hardness. We finally compare our predicted time lag - luminosity relation to the Norris et al. (2000) result obtained from 6 bursts with known redshift.

The prompt emission of gamma-ray bursts (hereafter GRBs) probably comes from a highly relativistic wind which converts its kinetic energy into radiation via the formation of shocks within the wind itself. Such "internal shocks" can occur if the wind is generated with a highly non uniform distribution of the Lorentz factor Gamma. Taking into account such a variable distribution of Gamma, we estimate the expected thermal emission of the relativistic wind when it becomes transparent. We compare this emission (temporal profile + spectrum) to the emission produced by the internal shocks. In most cases we predict a rather bright thermal emission that could easily be detected. This favors acceleration mechanisms for the wind where the main energy reservoir is under magnetic rather than thermal form. Such scenarios can produce thermal X-ray precursors comparable to those observed by GINGA and WATCH/GRANAT.

We have developed a toy model for internal shocks which has been used to generate a large number of synthetic GRBs in order to find in the parameter space the conditions which can lead to the formation of X-ray flashes. The key condition appears to be a small contrast of the Lorentz factor in the relativistic wind emitted by the central engine.

Heise et al. (2001) have reported the identification of a new class of sources by Beppo-SAX, the so-called X-ray flashes (XRFs), which have many common properties with gamma-ray bursts (GRBs) but are not detected in the gamma-ray range. In the framework of the internal shock model, we investigate the possibility that XRFs have the same origin than GRBs but are intrinsically softer due to different shock parameters.

Sgr A*, the compact radio source, believed to be the counterpart of the massive black hole at the galactic nucleus, was observed to undergo rapid and intense flaring activity in X-rays with Chandra in October 2000. We report here the detection with XMM-Newton EPIC cameras of the early phase of a similar X-ray flare from this source, which occurred on September 4, 2001. The source 2-10 keV luminosity increased by a factor about 20 to reach a level of 4 10^{34} erg s^{-1} in a time interval of about 900 s, just before the end of the observation. The data indicate that the source spectrum was hard during the flare. This XMM-Newton observation confirms the results obtained by Chandra and suggests that, in Sgr A*, rapid and intense X-ray flaring is not a rare event. This can constrain the emission mechanism models proposed for this source, and also implies that the crucial multiwavelength observation programs planned to explore the behaviour of the radio/sub-mm and hard X-ray/gamma-ray emissions during the X-ray flares, have a good chance of success.

Discerning the likelihood of the so-called runaway instability of thick accretion disks orbiting black holes is an important issue for most models of cosmic gamma-ray bursts. To this aim we investigate this phenomenon by means of time-dependent, hydrodynamical simulations of black hole plus torus systems in general relativity. The evolution of the central black hole is assumed to be that of a sequence of Kerr black holes of increasing mass and spin, whose growth rate is controlled by the transfer of mass and angular momentum from the material of the disk spiralling in through the event horizon of the black hole. The self-gravity of the disk is neglected. We find that when the black hole mass and spin are allowed to increase, constant angular momentum disks undergo a runaway instability on a dynamical timescale (a few orbital periods). However, our simulations show that a slight increase of the specific angular momentum of the disk outwards has a dramatic stabilizing effect. Our results, obtained in the framework of general relativity, are in broad agreement with earlier studies based both on stationary models and on time-dependent simulations with Newtonian and pseudo-Newtonian gravitational potentials.

The prompt emission of gamma-ray bursts probably comes from a highly relativistic wind which converts part of its kinetic energy into radiation via the formation of shocks within the wind itself. Such "internal shocks" can occur if the wind is generated with a highly non uniform distribution of the Lorentz factor. We estimate the expected photospheric emission of such a wind when it becomes transparent. We compare this thermal emission (temporal profile + spectrum) to the non-thermal emission produced by the internal shocks. In most cases, we predict a rather bright thermal emission that should already have been detected. This favors acceleration mechanisms for the wind where the initial energy input is under magnetic rather than thermal form. Such scenarios can produce thermal X-ray precursors comparable to those observed by GINGA and WATCH/GRANAT.

We present results from a numerical study of the runaway instability of thick discs around black holes. This instability is an important issue for most models of cosmic gamma-ray bursts, where the central engine responsible for the initial energy release is such a system consisting of a thick disc surrounding a black hole. We have carried out a comprehensive number of time-dependent simulations aimed at exploring the appearance of the instability. Our study has been performed using a fully relativistic hydrodynamics code. The general relativistic hydrodynamic equations are formulated as a hyperbolic flux-conservative system and solved using a suitable Godunov-type scheme. We build a series of constant angular momentum discs around a Schwarzschild black hole. Furthermore, the self-gravity of the disc is neglected and the evolution of the central black hole is assumed to be that of a sequence of exact Schwarzschild black holes of varying mass. The black hole mass increase is thus determined by the mass accretion rate across the event horizon. In agreement with previous studies based on stationary models, we find that by allowing the mass of the black hole to grow the disc becomes unstable. Our hydrodynamical simulations show that for all disc-to-hole mass ratios considered (between 1 and 0.05), the runaway instability appears very fast on a dynamical timescale of a few orbital periods, typically a few 10 ms and never exceeding 1 s for our particular choice of the mass of the black hole ($2.5 \mathrm{M_\odot}$) and a large range of mass fluxes ($\dot{m} \ga 10^{-3} \mathrm{M_{\odot}/s}$). The implications of our results in the context of gamma-ray bursts are briefly discussed.

Most models for the central engine of gamma-ray bursts involve a stellar mass black hole surrounded by a thick disk formed after the merging of a system of compact objects or the collapse of a massive star. Energy released from the accretion of disk material by the black hole or from the rotation of the hole itself extracted by the Blandford-Znajek mechanism powers a relativistic wind along the system axis. Lorentz factors of several hundreds are needed to solve the compactness problem in the wind which implies the injection of a tremendous power into a very small amount of matter. The Blandford-Znajek mechanism, where the outflow follows magnetic field lines anchored to the black hole is probably the best way to prevent baryonic pollution and can even initially produce a purely leptonic wind. In this paper we rather study the wind emitted from the inner part of the disk where the risk of baryonic pollution is much larger since the outflow originates from high density regions. We show that the baryonic load of this wind sensitively depends on the disk temperature and magnetic field geometry and that the outflow can become ultra-relativistic (Lorentz factor Gamma>100) under quite restrictive conditions only. Conversely, if Gamma remains of the order of unity the dense wind emitted from the inner disk could help to confine the central jet but may also represent a source of baryon contamination for the Blandford-Znajek mechanism.

The XMM-Newton X-ray observatory performed a pointed observation of the bursting pulsar GRO J1744-28 in April 2001 for about 10 ks during a program devoted to the scan of the Galactic center region. After the discovery of this source by BATSE in December 1995 during a very active bursting phase, it has been in quiescence since April 1997. We present here the first detection of GRO J1744-28 in its quiescent state at a position which is consistent with previous high-energy positions (but not consistent with a proposed IR counterpart). The observed luminosity of the source in quiescence is about 6 orders of magnitude weaker than the luminosity in outburst.

We derive a new formulation of the fully general relativistic equations describing a stationary equatorial MHD outflow from a rotating central object. The wind solution appears as a level contour of a `Bernoulli' function fixed by the requirements that it must pass through the slow and fast critical points. This approach is the general relativistic extension to the classical treatment of Sakurai (1985). We discuss in details how the efficiency of the magnetic to kinetic energy conversion depends mainly on the geometry of the flux tubes and show that the magnetic acceleration can work very well under some conditions. We show how this tool can be used for the study of several astrophysical phenomena, among which gamma-ray bursts.

We investigate a possible nucleosynthetic signature of highly energetic explosions of C-O cores ("hypernovae," HNe) which might be associated with gamma-ray bursts (GRBs). We note that the direct impact of C- and O-enriched hypernova ejecta on the ambient hydrogen and helium leads to spallation reactions which can produce large amounts of the light nuclides lithium, beryllium, and boron (LiBeB). Using analytic velocity spectra of the hypernova ejecta, we calculate the LiBeB yields of different exploding C-O cores associated with observed hypernovae. The deduced yields are $\sim 10^3$ times higher than those produced by similar (direct) means in normal Type II supernovae, and are higher than the commonly used ones arising from shock wave acceleration induced by Type II supernova (SN) explosions. To avoid overproduction of these elements in our Galaxy, hypernovae should be rare events, with $\la 10^{-3}$ hypernovae per supernova, assuming a constant HN/SN ratio over time. This rate is in good agreement with that of long duration GRBs if we assume that the gamma-ray emission is focussed with a beaming factor $\Omega/4\pi \la 10^{-2}$. This encouraging result supports the possible HN-GRB association. Thus, Galactic LiBeB abundance measurements offer a promising way to probe the HN rate history and the possible HN-GRB correlation. On the other hand, if hypernovae are associated to very massive pregalactic stars (Population III) they would produce a LiBeB pre-enrichment in proto-galactic gas, which could show up as a plateau in the lowest metallicities of the Be-Fe relation in halo stars.

The internal shock model for gamma-ray bursts involves shocks taking place in a relativistic wind with a very inhomogeneous initial distribution of the Lorentz factor. We have developed a 1D lagrangian hydrocode to follow the evolution of such a wind and the results we have obtained are compared to those of a simpler model presented in a recent paper (Daigne & Mochkovitch 1998) where all pressure waves are suppressed in the wind so that shells with different velocities only interact by direct collisions. The detailed hydrodynamical calculation essentially confirms the conclusion of the simple model: the main temporal and spectral properties of gamma-ray bursts can be reproduced by internal shocks in a relativistic wind.

We consider possible geometries of magnetic fields in GRB outflows, and their evolution with distance from the source. For magnetically driven outflows, with an assumed ratio of magnetic to kinetic energy density of order unity, the field strengths are sufficient for efficient production of gamma-rays by synchrotron emission in the standard internal shock scenario, without the need for local generation of small scale fields. In these conditions, the MHD approximation is valid to large distances (>10^19cm). In outflows driven by nonaxisymmetric magnetic fields, changes of direction of the field cause dissipation of magnetic energy by reconnection. This dissipation takes place outside the photosphere of the outflow, and can convert a significant fraction of the magnetic energy flux into radiation.

The source of cosmic gamma-ray bursts (hereafter GRBs) is usually believed to be a stellar mass black hole accreting material from a thick disk. The mechanism for the production of a relativistic wind by such a system is still unknown. We investigate here one of the proposal where the disk energy is extracted by a magnetic field amplified to very large values B \sim 10^15 G. Using some very simple assumptions we compute the mass loss rate along magnetic field lines and then estimate the Lorentz factor \Gamma at infinity. We find that \Gamma can reach high values only if severe constraints on the field geometry and the conditions of energy injection are satisfied. We discuss the results in the context of different scenarios for GRBs.

We compute the mass loss rate from a disk orbiting a stellar mass black hole assuming the flow is guided along magnetic field lines attached to the disk. We then estimate the Lorentz factor of the wind at infinity. We find that the Lorentz factor can reach high values only if severe constraints on the field geometry and the conditions of energy injection are satisfied. We discuss our results in the context of different scenarios for gamma-ray bursts. We mention a risk of 'contamination' of the Blandford-Znajek process by wind material emitted from the disk.

BeppoSAX observations of the single pulse burst GRB960720 have allowed a detailed study of its X-ray properties: pulse width in different energy bands, spectral evolution from 2 to 700 keV, etc. We show that the early (0-5s) X-ray emission and the gamma-rays are well explained by internal shocks in a relativistic wind while the late (5-20s) X-ray emission could come from the reverse shock generated in the wind when it interacts with the external medium. The results for a medium of uniform density are compared to the observations.

We construct models for gamma-ray bursts where the emission comes from internal shocks in a relativistic wind with a highly non uniform distribution of the Lorentz factor. We follow the evolution of the wind using a very simplified approach where a large number of layers interact by direct collisions but where all pressure waves have been suppressed. We suppose that the magnetic field and the electron Lorentz factor reach large equipartition values in the shocks. Synchrotron photons emitted by the relativistic electrons have a typical energy in the gamma-ray range in the observer frame. Synthetic bursts are constructed as the sum of the contributions from all the internal elementary shocks and their temporal and spectral properties are compared to the observations. We reproduce the diversity of burst profiles, the ``FRED'' shape of individual pulses and the short time scale variability. Synthetic bursts also satisfy the duration-hardness relation and individual pulses are found to be narrower at high energy, in agreement with the observations. These results suggest that internal shocks in a relativistic wind may indeed be at the origin of gamma-ray bursts. A potential problem is however the relatively low efficiency of the dissipation process. If the relativistic wind is powered by accretion from a disc to a stellar mass black hole it implies that a substantial fraction of the available energy is injected into the wind.