The possible ubiquity of energy injection in Gamma-Ray Burst afterglows
Since its launch in 2004, the Swift satellite has monitored the X-ray afterglows of several hundred
Gamma-Ray Bursts, and revealed that their X-ray light-curves are more complex than previously thought,
exhibiting up to three power-law segments. Energy injection into the relativistic blast-wave energizing
the burst ambient medium has been proposed most often to be the reason for the X-ray afterglow complexity.
We examine 117 light-curve breaks of 98 Swift X-ray afterglows, selected for their high-quality monitoring
and well-constrained flux decay rates. Thirty percent of afterglows have a break that can be an adiabatic
jet-break, in the sense that there is one variant of the forward-shock emission from a collimated outflow
model that can account for both the pre- and post-break flux power-law decay indices, given the measured
X-ray spectral slope. If allowance is made for a steady energy injection into the forward-shock, then
another 56 percent of X-ray afterglows have a light-curve break that can be explained with a jet-break.
The remaining 12 percent that are not jet-breaks, as well as the existence of two breaks in 19 afterglows
(out of which only one can be a jet-break), suggest that some X-ray breaks arise from a sudden change
in the rate at which energy is added to the blast-wave, and it may well be that a larger fraction of
X-ray light-curve breaks are generated by that mechanism.
The fractional increase in the shock energy, inferred from the energy injection required to account
for the observed X-ray flux decays, may be anticorrelated with the GRB prompt output, whether
the X-ray break is a jet-break or an energy-injection break. That anticorrelation can also be seen
as bursts with a higher energy output being followed by faster fading X-ray afterglows.
To test the above two mechanisms for afterglow light-curve breaks, we derive comprehensive analytical
results for the dynamics of outflows undergoing energy injection and for their light-curves, including
closure relations for inverse-Compton afterglows and for the emission from spreading jets interacting
with an wind-like ambient medium.
keywords:radiation mechanisms: non-thermal, relativistic processes, shock waves, gamma-ray bursts, ISM: jets and outflows
The power-law temporal decays displayed by the radio, optical, and X-ray fluxes of GRB afterglows () have provided strong support for the prediction (Meszaros & Rees 1997) that afterglows arise from the relativistic blast-wave generated as the GRB outflow ejecta interacts with the ambient medium. If the GRB outflow is a collimated jet then a break in the smooth power-law decay of the afterglow is expected for both either a conical jet (Panaitescu, Meszaros & Rees 1998) or a sideways-spreading jet (by Rhoads 1999). The breaks found in the flux decay light curves of most of the well-sampled afterglows that were localized by the BeppoSAX satellite seemed to provide solid observational support for the existence of these predicted jet-breaks (e.g. Frail et al 2001).
It was also realized early on that there should be a substantial diversity in the steepness of the afterglow flux decay, even when the spectral slope is constant, if the fundamental forward-shock parameters vary as the shock traverses the ambient medium. Rees & Meszaros (1998) derived the expected power-law flux decay for the forward-shock synchrotron flux in the case of an outflow whose kinetic energy-per-solid angle is not constant in time, either because the shock’s energy increases due to fresh ejecta arriving at the blast-wave or because is anisotropic. Two important microphysical parameters ( and ), which quantify the fraction of the post-shock energy acquired by electrons and magnetic field, also determine the intensity the afterglow synchrotron flux, and their variability could also yield a variety of decay indices .
Swift measurements of GRB X-ray afterglows provide support for that shock parameters vary. The standard forward-shock model with constant energy and shock micro-parameters predicts a closure relation between the temporal decay index and the slope of the X-ray continuum () of the form . The early X-ray flux of some Swift afterglows decays slower than expected for the measured spectral slope , which was attributed to a variable kinetic energy of the forward shock (Nousek et al 2006, Panaitescu et al 2006, Zhang et al 2006). Moreover, the flux decay indices displayed by Swift X-ray afterglows after the slow decay phase are not correlated with the X-ray spectral slope (Panaitescu 2007), as would be expected from the closure relation for the forward-shock emission.
When optical light-curves are taken into account, the standard jet model fares even worse, being inconsistent with practically all Swift afterglows (Liang et al 2008). While more than half of the well-sampled optical afterglows display coupled optical and X-ray light-curves (Panaitescu & Vestrand 2011), exhibiting similar flux decay indices and/or simultaneous light-curve breaks, many other afterglows show a decoupling of their optical and X-ray light-curves, with a light-curve break occurring only in the X-ray. At least four possibilities for this light-curve decoupling, which is not predicted by the standard jet model, have been proposed: ) the X-ray afterglow emission arises occasionally from the same not-yet-understood mechanism which produces the prompt burst emission (Ghisellini et al 2007), ) there is a dominant contribution of the reverse-shock to the afterglow emission (Uhm & Beloborodov 2007), ) the forward-shock emission is reprocessed by bulk-scattering and/or inverse-Compton scattering outside the forward-shock (Panaitescu 2008), or ) the optical and X-ray afterglow flux arise from different parts of the forward-shock, from different jets (e.g. Racusin et al 2008).
In this work, we assume that Swift X-ray afterglows arise from the forward-shock, and attribute the failure of the standard jet model in accounting for the X-ray flux decay indices to the assumption that the forward shock energy is constant.
Evolving micro-physical parameters (e.g. Ioka et al 2006) could, in principle, explain the measured X-ray flux decay indices, but the existence of light-curve breaks requires an unnatural, sudden change in their evolution. For that reason, we maintain in this work the basic assumption of the standard jet model that micro-physical parameters do not evolve.
We also assume that the ambient medium with which the post-GRB outflow interacts is uniform, without discontinuities that could yield breaks in the afterglow light-curve at a frequency. A notable model thus excluded is that where the afterglow-producing outflow interacts with a prior outflow (Ioka et al 2006) whose density increases radially outward up to its termination shock, yielding a slower decay of the afterglow flux than expected for a homogeneous ambient medium, followed by a break to a steeper decay.
One reason for relaxing the constant shock energy constraint is the observational evidence for X-ray flares in Swift afterglows (e.g. Chincarini et al 2007). The short timescale of those flares indicates that the GRB progenitor produces relativistic ejecta after the burst phase, whose arrival at the decelerating forward-shock could increase substantially the blast-wave’s energy. Our hypothesis is that such an injection of energy in the forward-shock is, in general, a more persistent and continuous process. However, we note that the fluence of X-ray flares is about 10 percent that of the burst (Falcone et al 2007), hence the late outflow powering the flares could carry a dynamically-important energy only if that outflow is much less efficient than the GRB outflow in producing X-ray emission.
We derive in §2 the closure relations for a continuous energy injection in the forward-shock. For the ease of extracting from observations the required injected power history, we assume that the forward-shock energy increases as a power with observer-time, , which is motivated by that the resulting afterglow flux should be power-law in observer time. This assumption circumvents a calculation of the kinematics for the ejecta – forward-shock interaction.
In §3, the closure relations for the forward-shock emission with energy injection are applied to a set of Swift afterglows with well sampled light-curves. There, we examine two plausible scenarios for the origin of afterglow X-ray light-curve breaks (break in energy injection, jet-break), three dynamical regimes of the forward-shock (spherical expansion, conical jet, spreading jet), and two radiation processes that may yield the afterglow X-ray flux (synchrotron, inverse-Compton).
2 Analytical flux decays for forward-shock models with energy injection
There are three mechanisms that could produce a change in the flux power-law decay index , i.e. an afterglow light-curve break.
(1) Passage of a spectral break. X-ray light-curve breaks are, generally, not accompanied by a spectral evolution (Nousek et al 2006; Liang, Zhang & Zhang 2007; Willingale et al 2007), hence their origin in the passage of the spectral feature is largely ruled out, but there are other reasons why such an X-ray light-curve break origin is unlikely. The slower evolution of the cooling frequency (the synchrotron characteristic frequency at which radiate the electrons that cool radiatively on the hydrodynamical timescale) leads to extended X-ray light-curve breaks when crosses the 0.3-10 keV band. The peak frequency of the synchrotron spectrum decreases faster, but it is likely to be lower than the X-rays even at early times, hence it cannot cross the X-ray band during afterglow observations.
(2) Jet-break. A light-curve break occurs at the ‘jet-break time’, when the outflow Lorentz factor drops below the inverse of the jet half-angle . For a ‘conical jet’ that does not spread laterally (as in the case of an outflow core embedded in an envelope that prevents the spreading of the core but yields a negligible emission), the steepening of the X-ray flux decay is solely due to that, after the jet-break time, the average kinetic energy per solid angle over the visible cone of angle (where relativistic Doppler boost is effective) is decreasing with observer-frame time. For a ‘spreading jet’ (as expected for a top-hat distribution of ), in addition to the above-described geometrical-relativistic effect, there is a second contribution to the X-ray flux decay steepening from the enhanced jet deceleration produced by the jet lateral spreading.
(3) Energy injection break. If energy injection occurs routinely in GRB afterglows, then a light-curve break could be produced by a change in the injected power. The dynamical behaviours of the forward-shock model to be considered are: ) a ‘spherical’ outflow, i.e. a wide jet whose edge becomes visible to the observer after the last X-ray measurement (at ), which occurs if jet half-angle satisfies , and ) a narrow jet whose edge is already visible at the first X-ray observation (at ), occurring when . For the latter model, we will consider both a conical jet and a spreading jet. These dynamical models are applicable also to X-ray afterglows without light-curve breaks, in which case the energy injection law must be unevolving (a single index ).
We will investigate the ability of the energy-injection break and jet-break models above to accommodate the observed X-ray flux decays before and after the break, with allowance for all three dynamical models (spherical outflow, conical jet, spreading jet), for two types of ambient medium (homogeneous and wind-like), two radiation mechanisms dominating the X-ray flux (synchrotron and inverse-Compton), and two possible locations of the X-ray relative to the cooling frequency but, first, a word about their rationale.
Wind-like media, having a particle density, are expected at 0.1-1 pc from the source of long GRBs, i.e. where the afterglow emission is produced, because the progenitors of such bursts are massive Wolf-Rayet stars that drive powerful winds. The termination-shock of such winds is of several pcs (Garcia-Segura et al 1996), thus a homogeneous medium at the afterglow location seems unjustified. However, if the GRB progenitor resides in a highly-pressurized bubble, the wind termination-shock radius could be as low as 0.5 pc (Chevalier et al 2004), with a further reduction of that radius if the WR wind is weak, as expected along the rotation axis (Ramirez-Ruiz et al 2005), or if the WR star moves sufficiently fast relative to the ISM (van Marle et al 2006), thus the afterglow could occur in shocked wind, which is quasi-homogeneous.
Synchrotron emission is generally believed to produce all the afterglow emission, but the likelihood of inverse-Compton emission being dominant is higher at higher photon frequency, like X-rays. Admittedly, that requires a high environment density, perhaps ten times higher (Panaitescu 2011) than implied by the mass-loss rates and terminal wind velocities measured for Galactic Wolf-Rayet winds (Nugis & Lamers 2000).
2.1 Spectral characteristics
The closure relations expected for a given variant of the forward-shock model are derived from the evolution of the spectral characteristics (peak flux), (frequency of the spectrum peak), and (cooling frequency), and from the shape of the spectrum at frequencies higher than its peak: for , for and for , where is the power-law index of the electron distribution with energy produced by the shock (). Therefore, the flux received at observing frequency satisfies for and for ( implies an X-ray spectrum with that is harder than measured for nearly all Swift X-ray afterglows).
For a spherical outflow (or a jet before the jet-break time ) and for synchrotron emission, , , . Here is the mass of the ambient medium swept-up by the forward-shock (i.e. the number of radiating electrons); is the post-shock magnetic field, whose energy density is assumed to be a fraction of the post-shock energy density (from the shock jump conditions); is the external density at the shock’s location; is the typical electron random Lorentz factor, parametrized as a fraction of the proton random Lorentz factor ; and is the random Lorentz factor of the electrons whose radiative cooling timescale equals the comoving-frame shock’s age , being the outflow radius and the Compton parameter. We assume the more likely situations of either electrons cooling mostly through synchrotron emission () or, if inverse-Compton is dominant (), that , in which case is constant and does not affect the evolution of .
For a jet (conical or spreading) after the jet-break time and for synchrotron emission, we have , with and as above. Thus
The characteristics of the inverse-Compton spectrum are immediately related with those of the synchrotron spectrum: (where is the optical thickness to electron scattering of the shocked ambient medium), , and , thus
2.2 Outflow dynamics
The above equations for the spectral characteristics show that their evolution and, thus, the resulting flux decay index , is set by the forward-shock dynamics: , , and . The energy-per-particle behind that shock is , hence the lab-frame energy of the forward-shock is and its dynamics is given by , where is the kinetic energy of the ejecta (all transferred to the forward-shock) at the time when the assumed power-law energy injection () begins. The other equations for the forward-shock dynamics are , with the density of the ambient medium at radius ( for a homogeneous medium, for a wind) and the increase of the jet opening at the comoving-frame sound speed (), and for the observer-frame arrival time of photons emitted along the direction toward the observer.
Defining , with (the radius where the lateral spreading almost doubles the initial opening of an adiabatic jet interacting with a homogeneous medium), and , the above equations for the outflow dynamics become:
Equation (10) allows the conversion of the jet dynamics [, , ] from lab frame to observer frame [, , ].
For a spherical outflow, a conical jet, or a spreading jet before the jet-break time, is constant and the above equations can be solved easily, leading to
which can be replaced in the equations for the synchrotron and inverse-Compton spectral characteristics (equations 1–3 and 4–6) to obtain their evolutions and, then, the power-law index of the afterglow flux decays given in lines 1, 2, 4, and 5 of Table 1 and Table 2, respectively.
|Model||Homogeneous medium||Wind-like medium|
|syn – spherical||-||-||-||-|
|syn – conical jet||-||-||-|
|syn – spreading jet||-|
|iC – spherical||-||-1|
|iC – conical jet||-||-|
|iC – spreading jet|
2.3 Spreading jets
2.3.1 Adiabatic jet, homogeneous medium
For a jet spreading laterally, the dynamics equations can be solved analytically in the case of no energy injection and for a homogeneous medium. That has been done by Rhoads (1999), who obtained an exponential jet deceleration for a lateral spreading given by : , , at , where . Because the differential lateral spreading is “on top” of a conical jet and not “over” a cylindrical jet (as in Rhoads’ formalism), the correct lateral expansion is . For and , equations (7)–(10) (which use the latter prescription for jet spreading) have the solution , , for , which can be shown to lead to the same flux power-law decay indices as for Rhoads’ prescription because the increase of the jet radius is only logarithmic in observer time (note that, for a nearly constant jet radius , implies ).
The above exponential solutions for the dynamics of an adiabatic jet interacting with a homogeneous medium can be expressed in observer time by integrating , which leads to for . Thus
from where it can be shown that
where is the jet energy and .
|Model||Homogeneous medium||Wind-like medium|
|syn – spherical|
|syn – conical jet|
|syn – spreading jet|
|iC – spherical|
|iC – conical jet|
|iC – spreading jet|
2.3.2 Energized jets
and integrated numerically, which showed that the left-hand sides are (asymptotically) constant, i.e. the dynamics [, , ] of a spreading jet undergoing energy injection according to is a power-law in observer time. That dynamics is between the extreme cases of an adiabatic jet interacting with a homogeneous medium (which yields the fastest deceleration) and a spherical outflow with energy injection (which yields the slowest deceleration). Given that, in both of those two extreme cases, the dynamical quantities have a power-law evolution with time, it is not surprising that an intermediate case manifests that feature too.
Replacing power-law solutions for , , and in equations (19)–(21), allows the calculation of the temporal exponents by requiring that the right-hand sides are also constants, i.e. without a dependence. The result is
It is worth mentioning that the above dynamics of spreading jets undergoing energy injection can be recovered from the dynamics of a spreading adiabatic jet interacting with a homogeneous medium (equations 16, 17) if one replaces (the prescription for energy injection) and . The latter relation is trivial for an outflow decelerating as a power-law with radius (as for a spherical outflow or a conical jet) and is approximately satisfied for an exponential deceleration (as for an adiabatic spreading jet interacting with a homogeneous medium), which suggests that it should also stand in any other case of an intermediate-strength deceleration (that is verified with the numerical integration of the jet dynamics). Either from equations (22) and (24)or by replacing and in equations (15) and (18) with the above prescriptions, it can be shown that
for a spreading jet and with energy injection.
By substituting equations (22)–(26) in the relations for the synchrotron and inverse-Compton spectral characteristics (equations 1–3 and 4–6), one can derive the evolutions from lines 3 and 6 of Table 1 and the power-law indices of the afterglow flux decay from lines 3 and 6 of Table 2.
We acknowledge that the dynamics of spreading jets undergoing energy injection was derived using a simplistic toy model for the jet lateral spreading, where the jet was assumed to have a uniform kinetic energy per solid angle , although both lateral spreading and energy injection could/should introduce a substantial angular gradient of , decreasing with angle measured from the jet’s spine.
3 Energy injection in Swift X-ray afterglows
3.1 Sample selection
Out of the 634 afterglow light-curves listed at the Swift-XRT light-curve repository (Evans et al 2007) posted until the end of year 2011, we have retained only those for which a light-curve break is obvious or seems possible, a power-law light-curve segment is monitored for longer than a factor 5 in time before and after the break, and the spectral slope listed at the Swift-XRT spectra repository (Evans et al 2007) has a 90 percent confidence level (CL) less than 0.16.
The light-curves have been piece-wise fit with power-laws, using measurements sufficiently far from the break time, to capture the asymptotic power-law flux decay. Shallow breaks for which the data could be fit satisfactorily (i.e. reduced ) with a single power-law were excluded. The 90 percent confidence level () was determined from the variation of around the best power-law fit.
As shown in Table 2, for an adiabatic outflow (), the power-law index is + const = /2 + const, with the coefficient of being as large as 9/8 for other models. Thus, an equal importance of measurement uncertainties and requires that we retain only afterglows whose decay indices can be determined with a 90 percent CL smaller than .
The above two basic selection criteria led to the sample of 98 afterglows and 117 X-ray light-curve breaks listed in Table 3. Forty-four afterglows have temporal coverage starting from the burst end (for them, we determine the fractional increase of the outflow energy during the afterglow phase) and a measured redshift, which allows the calculation of the GRB output and of the afterglow jet energy.
The light-curves of the afterglows listed in Table 3 are morphologically diverse. A minority display a continuous, slowly-decaying power-law flux from the burst phase to the afterglow, the afterglow beginning being chosen when the flux fluctuations disappear. The majority of GRBs display a sharp decay after the prompt phase (the GRB tail) followed by a much slower decay and spectral hardening; the afterglow beginning is chosen at the epoch when the slower flux decay starts.
The closure relations listed in Table 2 are used to calculate the required energy injection indices from the measured pre- and post-break X-ray flux decay indices and the X-ray slope (which gives the electron index ). For a model to be validated, the resulting energy injection index should be positive (i.e. , with the uncertainty calculated from the propagation of and ), and for jet-break models, should lead to a decelerating outflow, so that the jet boundary can become visible to the observer at some time. From equation 11, this condition is .
3.2 Adiabatic jets
For adiabatic jets (), the observables and yield two constraints. As there are no free model parameters, aside the various possibilities related to the jet dynamics, location of the cooling frequency, ambient medium type, and dominant X-ray emission process, the model is over-constrained. Thus, it should not be surprising that only a small fraction of Swift X-ray breaks satisfy both before and after the break, i.e. can be explained by any single variant of the adiabatic jet break model (upper part of Table 4).
|model||Homog medium||Wind-like medium|
|jet-break with or without steady EI|
Evidently, Swift X-ray afterglows may not be all of the same forward-shock model type. Even when we take into account all these possible variants, we find that adiabatic jet-break model can explain only 30 afterglows (indicated in boldface in Table 3), i.e. about 30 percent of the sample considered here. Six of those 30 breaks occur in afterglows with two light-curve breaks, the break compatible with an adiabatic jet-break being always the second. For all six, there is always a spherical energy-injection break model with after the break that accounts for the first break, and compatible with the jet-break model, i.e. the two break models have the same location of cooling frequency relative to the X-rays, the same type of external medium, and the same X-ray emission process.
3.3 Energized jets
For the jet-break model, we require that the X-ray light-curve break arises solely from the outflow collimation, without any contribution from a break in the injected power, i.e. we require a constant energy injection index across the X-ray light-curve break: . Thus, the model is still over-constrained, with two observables ( and ) that must be accommodated by a single model parameter (). Considering that the jet-break model with steady energy injection injection is less over-constrained than the adiabatic jet model, it may not be surprising that the former model can account for a larger number of Swift breaks (mid part of Table 4).
With allowance for all jet-break model variants listed in Table 4 (i.e. for either location of the cooling frequency relative to the X-ray, for either radiation process, for either type of medium, and either post jet-break dynamics, conical or spreading), the jet-break model with steady energy injection (including adiabatic jet-breaks, for which the energy injection parameter is consistent with zero) can account for 103 of the 117 X-ray light-curve breaks considered here (the unexplained breaks are shown in italics in Tables 3). 37 of those potential jet-breaks occur in afterglows with two light-curve breaks; with the exception of the first break of afterglow 081222, all other 36 jet-breaks can be paired consistently with an energy-injection break for the remaining break, i.e. the two break models have the same medium stratification, the same X-ray radiative process, the same location of X-ray relative to the cooling frequency, and compatible ejecta dynamics.
Because only one break can arise from the outflow collimation and 19 afterglows in our sample have two light-curve breaks, a more correct way of assessing the success of the jet-break model with steady energy injection is to say that it can explain at least one of the breaks observed for 85 afterglows out of 98 or, conversely, that it cannot account for any of the X-ray light-curve breaks of 12 percent of afterglows. Excluding the 30 breaks (of 30 afterglows) that can be accounted for with an adiabatic jet-break, the jet-break model with non-zero energy injection can explain one X-ray light-curve break for 55 of the 98 afterglows considered here, which argues strongly that energy injection may be often present in Swift afterglows.
3.3.2 Energy injection breaks
Given that 12 percent of afterglows have an X-ray light-curve break that cannot be explained with a jet-break even when a steady energy injection is allowed, and that some afterglows display two breaks, out of which only one can be a jet-break, it is natural to speculate that some X-ray light-curve breaks must arise from a change in the power at which energy is added to the forward-shock. In the extreme, one could speculate that all X-ray breaks could be such energy injection breaks, but the large fraction of breaks that can be explained as jet-breaks argues against that extreme speculation.
The broken power-law energy injection model is not over-constrained, as there are two observables ( and ) for two model parameters ( and ), the only restriction being that . Thus, it should not be surprising that there are variants of the energy-injection break model that can account for most of the Swift X-ray breaks (lower part of Table 4). There is only one X-ray break that cannot be explained with an energy injection break: GRB afterglow 060607. From Table 2, it can be seen that, for , the steepest possible decay for an adiabatic outflow is that for the inverse-Compton emission from a conical jet interacting with a wind, and for : 111 For , the fastest decay , exceeds the limit for the large-angle emission decay, resulting when the emission from a spherical outflow stops suddenly; however, for a jet seen after the jet-break time, the large-angle emission does not exist because there is no radiating fluid at angles larger than . For GRB afterglow 060607, , hence the steepest possible decay for the forward-shock emission has , which is slower than the measured post-break .
4 Implications of Energy Injection
4.1 Afterglow decay rate – GRB output correlation
After having identified the power-law energy injection that accommodates each of the X-ray light-curve breaks, one can calculate the increase in the forward-shock energy from first to last X-ray measurement. The interesting result is that the fractional increase of the forward-shock energy over a source-frame time-interval common to all afterglows (from 100 s to 100 ks), , is anticorrelated with the isotropic-equivalent GRB output , for any model variant that can account for a significant fraction of the afterglow breaks listed in Table 3.
The anticorrelation is illustrated in Figure 1 for two jet-break models and in Figure 2 for two energy-injection break models. For the synchrotron jet-break model (black points in Figure 1), the linear correlation coefficient in log-log space is , yielding a probability for a stronger anticorrelation in the null hypothesis. For the inverse-Compton jet-break model (red points in Figure 1), and . For the synchrotron energy-injection break model (black points in Figure 2), , and , while for the inverse-Compton energy-injection break model (red points in Figure 2), , and . Similarly high chance probabilities are obtained for other models; thus the anticorrelation is only “tentative”.
The above linear correlation coefficients (Pearson’s) have been calculated assuming a 50 percent uncertainty (in natural scale) for both quantities plotted, the chance probabilities of those anticorrelations have been estimated conservatively for a linear correlation coefficient . The best-fits shown in Figures 1 and 2 were obtained by minimizing using the uncertainties of both quantities.
That ) the fractional increase of the forward-shock energy over a fixed time interval is determined by the exponent of the power-law energy injection which accommodates the measured X-ray flux decay index , and that ) (a stronger energy injection yields a slower flux decay) suggest that the correlation may originate from a correlation. Indeed, as illustrated in Figure 3, more energetic bursts are followed by faster decaying X-ray afterglows. The correlation between these two observables is slightly stronger than the correlation for afterglows with jet-breaks at the X-ray light-curve epoch, but significantly weaker if the X-ray light-curve breaks are energy-injection breaks.
The correlation implies that either the efficiency of the GRB mechanism or the energy of the ejecta producing the GRB emission are anticorrelated with the kinetic energy of the ejecta interacting with the forward-shock after the burst phase.
If the dissipation mechanism leading to the production of -rays is internal shocks (Rees & Meszaros 1994) in a fluctuating relativistic wind, then the burst dissipation efficiency is determined by the ratio between the Lorentz factors of interacting shells. For an impulsive ejecta release (short-lived engine), all the ejecta are produced quasi-instantaneously and the energy injection occurs as the inner parts of the outflow catch-up with the ever-decelerating forward-shock (Rees & Meszaros 1998). In this model, the correlation implies an anticorrelation between the average shells’ Lorentz factor ratio with the energy of the shells behind the leading edge of the outflow, smoother outflows (i.e. with a smaller Lorentz factor contrast) having more energy localized toward the trailing edge. For an extended ejecta release (long-lived engine), where the ejecta are released over a source-frame duration comparable to the observer-frame time during which energy injection takes place, the correlation would imply an anticorrelation among bursts between the Lorentz factor contrast in the leading, GRB-producing part of the outflow and the kinetic energy of the outflow expelled by the central engine at later times. This is a puzzling implication, which makes the long-lived engine scenario less likely.
Irrespective of the dissipation mechanism, if the GRB output is correlated with the energy of the GRB-producing ejecta, then the correlation implies a relationship between the GRB ejecta energy and the radial distribution of the afterglow outflow energy. For a short-lived central engine, where the GRB and afterglow outflows are the same, the correlation would require that more energetic outflows have less energy behind the leading front. For a long-lived central engine, where the GRB outflow is the early afterglow blast-wave, the correlation implies that engines that are initially more energetic release less energy after the burst. Perhaps this scenario appears more plausible because it opens the possibility of a constant total ejecta energy. Unfortunately, we do not find that is universal; instead, its distribution extends over at least two decades, for the calculated in any light-curve break model.
We note that the detection of X-ray flares in some Swift afterglows indicates an extended ejecta release, but the energetic output of those flares is only about 10 percent of the prompt GRB output (Falcone et al 2007). In contrast, the energy injection required during the afterglow phase corresponds to an increase of the outflow energy after the burst by a factor up to (Figures 1 and 2). To reconcile the flare and afterglow energetics, the internal shocks that produce flares must be very inefficient X-ray radiators compared to the forward-shock (see also Maxham & Zhang 2009) or, else, the X-ray flares could be produced by the central engine, bearing little correlation with the energy of the ejecta released after the burst.
4.2 Afterglow energetics
The large increase of the forward-shock energy, by a factor (Fig 1) or (Fig 2), required to account for the light-curve breaks of Swift X-ray afterglows raises the question if the total outflow energy is not too large and incompatible with the expected energy reservoir for a few solar-masses black-hole and its sub solar-mass torus formed in the collapse of a WR massive star.
The jet energy can be estimated by assuming that, at the first X-ray measurement (), the isotropic-equivalent forward-shock energy is the GRB output . This assumption introduces an underestimation of the jet energy by a factor up to a few because energy could be injected into the forward-shock even before (100 – 1000 s), maybe starting as early as the burst end (10 – 100 s), at a rate (this is the range of energy injection laws that accommodate the early X-ray flux decays of Table 3), and with the estimated jet energy satisfying for a homogeneous medium and for a wind. For a short-lived engine, where all the outflow produces the burst emission and only its leading edge yields the early afterglow, we have ; in this case, the total jet energies calculated here may be substantially overestimated. For a long-lived engine, where the leading part of the outflow produces both the burst and early afterglow emissions, is equivalent to assuming a (reasonable) 50 percent GRB efficiency.
The jet initial aperture , needed to calculate the jet energy from isotropic equivalents, can be constrained as following. For the model of a narrow jet whose edge becomes visible before the first X-ray measurement and the X-ray light-curve break is an energy injection, . For the model where the X-ray light-curve break (at ) is a jet-break undergoing a steady energy injection, . For the model of a wide jet whose edge becomes visible after the last X-ray measurement (), . Prior to the jet-break, the jet Lorentz factor is given by the dynamics of a spherical blast-wave:
for a homogeneous medium of particle density per and
for a wind-like medium of density normalized to that corresponding to progenitor mass-loss rate-to-terminal speed ratio of , respectively, with being the blast-wave isotropic-equivalent energy at observer time . The latter is calculated assuming that and for the energy injection law required for each model to account for the measured X-ray flux decay. The observer time in the above equations corresponds to the arrival-time of photons emitted by the fluid moving at an angle relative to the direction toward the observer, from where arises most of the high-energy emission received by the observer at time .
Putting together all the elements described above, we have calculated the jet final energy (at ) for the X-ray afterglows whose light-curve breaks can be explained by the a jet-break model with steady energy injection or an energy-injection break model (the models listed in the mid and lower parts of Table 4). The jet energies for and , averaged over all the afterglows that can be accounted for by each model (fewer than indicated in Table 4, because the jet energy calculation requires the GRB output/burst redshift), are listed in Table 5. For jet-breaks models, the jet opening is determined by the light-curve break epoch. For energy injection-break models, the epoch of the first (last) X-ray measurement provide an upper (lower) limit on the jet-break epoch, which yields an upper (lower) limit on the jet opening and, consequently, on the jet energy. Numbers in round brackets give the dispersion of individual jet energies (or their lower/upper limits).
|model||Homog medium||Wind-like medium|
A potentially important result shown in Table 5 is that models where X-ray is inverse-Compton, the ambient medium is a wind, and X-ray is lower than the cooling frequency, can be ruled out on energetic grounds because, for them, the required jet total energy exceeds erg222 Magneto-hydrodynamic simulations of accretion disks (e.g. McKinney 2005, Hawley & Krolik 2006) show that the energy in relativistic matter (electromagnetic outflow) extracted from an accreting black-hole increases from 0.2 (0.03) percent of the accreted mass, for a non-rotating black-hole, to 15-40 (5-20) percent for a maximally rotating one. The maximal energy budget quoted here corresponds to a 20 percent MHD efficiency and a accretion disk. . Those models lead to a high final jet energy because a substantial increase of the jet energy is required to account for the measured X-ray flux decay rate, which stems from that the corresponding adiabatic jet model yields a steep X-ray flux decay (Table 2: , for a conical jet), from that energy injection has a weak effect in mitigating the flux decay (, for the spreading jet model), or a combination of those two factors (for the spherical outflow model).
Assuming the same jet energy at the beginning of the afterglow, one can also calculate the ratios of jet energies among the three types of X-ray break models and outflow dynamics: a narrow-jet model with an energy-injection break (models E-cjet and E-sjet in Table 5), a jet-break model with steady energy injection (models J-cjet and J-sjet), and a wide-jet model with energy-injection break (model E-sph). The interesting result is that, going from type to type models above increases on average (i.e over afterglows) the total jet energy by a factor 2-3, with another increase by a factor 2-3 required when going from type to models. These results hold for either type of ambient medium, either location of the cooling frequency relative to the X-ray, and for either X-ray radiative process, if the jet does not spread sideways, and only for inverse-Compton emission, if the jet is spreading laterally. Thus, the most economical model for a light-curve break is that of a narrow jet with an early jet-break and a broken power-law energy injection origin for the X-ray light-curve break. However, these conclusions do not hold for spreading jets if the X-ray radiation process is synchrotron. In this case, we find that the required jet energy is about the same for all model types above. Note that these relative energy requirements cannot be deduced from the average jet energies listed in Table 5. Here, we are referring to averages over afterglows of jet energy ratios between two models that explain the same afterglow break, while Table 5 gives jet energies for a given model averaged over all the afterglows that can be explained by that model.
4.3 Loose ends
The first caveat/weakness of the model testing presented above was already discussed: the energy-injection break origin for X-ray light-curve breaks is only weakly falsifiable by a comparison between the model closure relation and the observed flux power-law decay index and spectral slope . The only way for this model to fail is an X-ray flux that decays steeper than expected for an adiabatic shock, as energy injection can only mitigate the afterglow flux decay.
Furthermore, a comparison of the model light-curve breaks against observations is needed to validate possible X-ray light-curve break models, as the sharpness of X-ray light-curve breaks varies among afterglows. The majority of the well-sampled light-curve breaks used here are sharp, lasting less than a factor two in time. Quite likely, the sharpest breaks cannot be accounted by a wind-like medium, whatever is the break mechanism, owing to the slower jet deceleration produced by a decreasing external density (Kumar & Panaitescu 2000). Without taking into account such limitations, testing only the ability of some models to account for the observed flux power-law decay indices will overestimate the fraction of afterglows that can be explained by a given model.
We emphasize the importance of using afterglows whose light-curves have been sampled over a sufficiently long time before and after the break, so that the asymptotic flux decay indices can be measured accurately. Otherwise, shortly-monitored light-curves with breaks can lead to an underestimation of the break magnitude , which will favour a steady energy injection interpretation, because its effect is to reduce the break magnitude . We note that eighty percent of the power-law segments used in this work last a decade or more in time.
4.3.1 Optical afterglows
Inclusion of the optical light-curves could provide a further test to the origin of X-ray light-curve breaks (e.g. Liang et al 2007; Curran et al 2009; Liang et al 2008) in the two models discussed above (jet-break and energy-injection break). Both models yield achromatic light-curve breaks, appearing at the same time in all light-curves. For those afterglows where it can be established that the optical and X-ray emissions arise from the same source, either because the light-curve breaks are achromatic or because the optical and X-ray spectra are on the same power-law continuum, testing the light-curve break origin and determining the model details (location of cooling frequency, ambient medium stratification, radiation process, outflow dynamics) can be done by comparing the optical flux decay indices with the model expectations, for the energy injection parameter required by the X-ray flux decay indices .
In GRBlog database of optical light-curves (Quimby, McMahon & Jeremy 2004), we find five afterglows (050401, 061121, 070420, 090424, 090618) listed in Table 3 with chromatic X-ray breaks, hence their X-ray and optical emissions arise from different parts of the outflow, and eight with achromatic breaks, thus a comprehensive test of the X-ray break models cannot be performed at this time. We note that, of those 8 afterglows with achromatic breaks, two (060729, 081029) have , for five (050922C, 051109, 060614, 080710, 090424) , and for one (090510) .
In general, (as for 060729) is expected if the cooling frequency is not between the optical and X-ray and if is in between. According to Table 1, a decreasing (which leads to ), is obtained only for a homogeneous medium, thus all six afterglows with require such an ambient medium. Energy injection accelerates the evolution of and increases the index difference above the expectation for an adiabatic outflow: for a spherical/conical outflow and for a spreading jet. Thus, the energy injection indices in the range required by the X-ray decays of those six afterglows imply that , consistent with what is measured for five of them, while the large measured for 090510 indicates that the optical and X-ray are from different radiative processes (synchrotron and inverse-Compton, respectively).
4.3.2 Reverse shock
In this work, we have considered only the forward-shock emission and ignored the emission from the reverse-shock, which energizes the GRB/post-GRB ejecta, the agent of energy injection into the blast-wave. The two shocks dissipate equally the kinetic energy of the ejecta, which indicates that the reverse-shock could make a comparable contribution to the afterglow emission.
For a short-lived ejecta source, it can be shown from the kinematics of the catch-up between the GRB ejecta and the decelerating forward-shock that the ratio of the Lorentz factor of the ejecta arriving at the forward-shock at time to that of the swept-up ambient medium is (evidently, the ratio of those two Lorentz factors must be above unity, but cannot be much above unity, else the catch-up would have occurred earlier). In this case, the reverse shock is only mildly relativistic and the ejecta electrons accelerated by the reverse-shock could be much less energetic than the electrons behind the forward-shock, making the emission of the latter dominant at X-ray photon energies.
For a long-lived central source, the contrast between the Lorentz factor of the unshocked outflow and that of the forward-shock can be arbitrarily large, the reverse-shock could be highly relativistic, provided that the ejecta are not too dense (e.g. Sari & Piran 1995). The decay of the reverse-shock flux can be calculated in a way similar to that shown here for the forward-shock; however, the parametrization of the reverse-shock emission is not that simple because ) its dynamics depends on the Lorentz factor of the ejecta and their density, and ) its peak flux depends on the total ejecta mass having entered the reverse-shock. Thus, the calculation of the reverse-shock emission requires at least two parameters that quantify the energy and mass mass injection into the blast-wave, leading to a model even less constrainable (observationally) than the forward-shock (which requires only one parameter, , for the energy injection law).
By comparing the flux decay index () with the spectral slope () measured for a set of 98 well-sampled Swift X-ray afterglows with the expectations for the forward-shock model, we find that about a third of those afterglows display a break that could be the traditional ‘jet-break’ claimed to have been observed in a dozen pre-Swift optical afterglows (§2.3.1). For this test of adiabatic jet-breaks, we have considered various details of the forward-shock model: homogeneous or wind-like ambient media, spectral cooling frequency higher or lower than the X-rays, synchrotron or inverse-Compton as the dominant X-ray process, a collimated conical or spreading outflow.
If a steady energy injection into the forward-shock occurs, then the jet-break model may account for 88 percent of the 117 light-curve breaks or, else put, for at least one break in 87 percent of the 98 afterglows selected here (§3.3.1). This vast improvement over the adiabatic jet model suggests that energy injection could be a prevalent process in GRB afterglows.
The remaining 12 percent of breaks that cannot be explained by a jet-break undergoing a steady energy injection, and the 19 afterglows with two breaks (only one being, at most, a jet-break), indicate that some X-ray breaks may originate from a change in the rate at which energy is added to the forward-shock (an energy-injection break), and it is possible that some (many ?) of the identified jet-breaks are also energy-injection breaks (§3.3.2). It should be recognized that the complete success of the energy-injection break model over the traditional jet-break model could be due mostly to that the latter model is over-constrained (observations provide two constraints: flux pre- and post-break power-law decay indices – for one model parameter: the energy injection power-law exponent), while the former model is not.
We find a tentative correlation between the GRB output and the afterglow X-ray flux decay rate, more energetic bursts being followed by faster decaying afterglows (§4.1). If energy injection is a common process in GRB afterglows, then the slower decaying afterglows should be identified with a stronger total energy injection. This means that the GRB energy – afterglow (X-ray flux) decay rate correlation should lead to an anticorrelation of the GRB output with the increase of the outflow kinetic energy during the afterglow phase, i.e. more energetic bursts should be followed by afterglows with less energy being injected. This anticorrelation is, indeed, observed for both the jet-break and the energy injection break models for X-ray light-curve breaks, and for various model details (dynamics, X-ray afterglow emission process, radial stratification of the ambient medium). If the GRB progenitor is a long-lived source of relativistic ejecta, the above anticorrelation implies an anticorrelation of the kinetic energy of the leading outflow (producing the burst) with that of the following ejecta (which inject energy in the blast-wave during the afterglow phase), but the total ejecta energy (burst plus afterglow) is far from being universal.
The energy budget required by most of the possible variants of the forward-shock model ranges from to ergs, and only variants involving inverse-Compton emission in the X-rays and a wind-like medium can be excluded on energetic grounds because, in those models, the required jet-energy exceeds erg (§4.2). Considering both models for light-curve breaks (jet-break with steady energy-injection and energy-injection break), we find that the most energetically economical model is that of a narrow jet whose jet-break occurs early on (before first X-ray measurement) with the light-curve break being an energy-injection break, followed by a model where the X-ray light-curve break is a jet-break, with the most ”wasteful” model being that of a wide jet whose break occurs late (after last X-ray measurement), the X-ray break being an energy-injection break. However, the average (over afterglows) energy ratio among pairs of these three break models is only 2–3.
This work was supported by an award from the Laboratory Directed Research and Development program at the Los Alamos National Laboratory and made use of data supplied by the UK Science Data Center at the University of Leicester, by the “GRB log” site for optical light-curves (http://grblog.org/grblog.php), and by the GRB afterglow repository of GCN circulars (http://www.mpe.mpg.de/jcg/grbgen.html).
-  Chevalier R., Li Z., Fransson C., 2004, ApJ, 606, 369
-  Chincarini G. et al, 2007, ApJ, 671, 1903
-  Curran P. et al, 2009, MNRAS, 395, 580
-  Evans E. et al, 2007, A&A, 469, 379
-  Evans E. et al, 2009, MNRAS, 397, 1177
-  Falcone A. et al, 2007, ApJ, 671, 1921
-  Ghisellini G., Ghirlanda G., Nava L., 2007, ApJ, 658, L75
-  Hawley J. & Krolik J., 2006, ApJ, 641, 103
-  Ioka K. et al, 2006, A&A 458, 7
-  Kumar P., Panaitescu A., 2000, ApJ, 541, L9
-  Liang E., Zhang B.-B., Zhang B., 2007, ApJ, 670, 565
-  Liang E. et al, 2008, ApJ, 675, 528
-  van Marle A. et al, 2006, A&A, 460, 105
-  Maxham A. & Zhang B., 2009, ApJ, 707, 1623
-  McKinney J., 2005, ApJ, 630, L5
-  Mészáros P. & Rees M., 1997, ApJ, 476, 232
-  Nousek J. et al, 2006, ApJ, 642, 389
-  Nugis T. & Lamers H., 2000, A&A, 360, 227
-  Panaitescu A., Mészáros P., Rees M., 1998, ApJ, 503, 314
-  Panaitescu A. et al, 2006, MNRAS, 366, 1357
-  Panaitescu A., 2007, MNRAS, 379, 331
-  Panaitescu A., 2008, MNRAS, 383, 1143
-  Panaitescu A., 2011, MNRAS, 414, 1379
-  Panaitescu A. & Vestrand T., 2011, MNRAS, 414, 3537
-  Quimby R., McMahon E., Jeremy M., 2004, “GRBs: 30 years of Discovery”, AIP Conf Proc, 727, 529
-  Racusin J. et al, 2008, Nature, 455, 183
-  Ramirez-Ruiz E. et al, 2005, ApJ, 631, 435
-  Rees M. & Mészáros P., 1994, ApJ, 430, L93
-  Rees M. & Mészáros P., 1998, ApJ, 496, L1
-  Rhoads J., 1999, ApJ, 525, 737
-  Sari R. & Piran T., 1995, ApJ, 455, L143
-  Uhm Z. & Beloborodov A., 2007, ApJ, 665, L93
-  Willingale R. et al, 2007, ApJ, 662, 1093
-  Zhang B. et al, 2006, ApJ, 642, 354