Articles

SYNCHROTRON SELF-COMPTON EMISSION AS THE ORIGIN OF THE GAMMA-RAY AFTERGLOW OBSERVED IN GRB 980923

, , and

Published 2012 May 2 © 2012. The American Astronomical Society. All rights reserved.
, , Citation N. Fraija et al 2012 ApJ 751 33 DOI 10.1088/0004-637X/751/1/33

0004-637X/751/1/33

ABSTRACT

GRB 980923 was one of the brightest bursts observed by the Burst and Transient Source Experiment. Previous studies have detected two distinct components in addition to the main prompt episode, which is well described by a Band function. The first of these is a tail with a duration of ≃ 400 s, while the second is a high-energy component lasting ≃ 2 s. We summarize the observations and argue for a unified model in which the tail can be understood as the early γ-ray afterglow from forward shock synchrotron emission, while the high-energy component arises from synchrotron self-Compton from the reverse shock. Consistency between the main assumption of thick shell emission and agreement between the observed and computed values for fluxes, break energies, starting times, and spectral indices leads to a requirement that the ejecta must be highly magnetized.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

The most successful theory in terms of explaining gamma-ray bursts (GRBs) and their afterglows is the fireball model (see Mészáros 2006; Zhang & Mészáros 2004 for recent reviews). This model predicts an expanding ultrarelativistic shell that moves into the external surrounding medium. The collision of the expanding shell with another shell (internal shocks) or the interstellar medium (ISM; external shocks) gives rise to radiation emission through the synchrotron and synchrotron self-Compton (SSC) processes. In addition, when the expanding relativistic shell encounters the external medium two shocks are involved: an outgoing, or forward, shock (Rees & Mészáros 1994; Paczyński & Rhoads 1993) and another one that propagates back into the ejecta, the reverse shock (Mészáros & Rees 1994, 1997a).

According to the standard relativistic fireball model, the forward shock accelerates electrons (through first-order Fermi mechanism or electric fields associated with the Weibel instability) up to relativistic energies and generates magnetic fields (Medvedev & Loeb 1999). The afterglow emission is more likely to be synchrotron (Sari et al. 1998); however, inverse Compton (IC) scattering may affect the electron radiative cooling. The progressive, power-law deceleration of the forward shock leads to a continuous softening of the afterglow spectrum (Panaitescu 2007). This has been observed in the spectra and light curves of GRB 970228 (Sari et al. 1998), GRB 970508 (Panaitescu et al. 1998; Sari et al. 1998), GRB 980923 (Giblin et al. 1999), and among a dozen other afterglows (Zeh et al. 2006).

On the other hand, the reverse shock is predicted to produce a strong optical flash (Mészáros & Rees 1997a; Sari & Piran 1999a, 1999b). When it crosses the shell, the reverse shock heats it up and accelerates electrons, but it operates only once. Thus, unlike the forward shock emission that continues later at lower energies, the reverse shock emits a single burst. After the peak of the reverse shock no new electrons are injected and the shell material cools adiabatically (Vedrenne & Atteia 2009). However, this picture can be modified allowing a long-lived reverse shock if the central engine emits slowly moving material in which the reverse shock could propagate and survive for hours to days (Genet et al. 2007; Uhm & Beloborodov 2007). Now, although the contribution of the reverse shock synchrotron emission to the X-ray band is small, electrons in the reverse shock region can upscatter the synchrotron photons (SSC process) up to the X-ray band or even higher energies (Mészáros & Rees 1993; Wang et al. 2001a, 2001b, 2005; Granot & Guetta 2003; Kobayashi et al. 2007; Zhang et al. 2006).

The soft tail component seen in GRB 980923 has been described reasonably well by synchrotron emission from a decelerating relativistic shell that collides with an external medium (forward shock model), and has also been described as such in more than a dozen independent events (see, e.g., Chevalier & Li 1999; Rhoads 1999; Sari et al. 1998; Panaitescu 2007; Giblin et al. 1999; Shao & Dai 2005). The high-energy component, on the other hand, has been explored through IC emission models (Granot & Guetta 2003; Pe'er & Waxman 2004) and SSC emission from the reverse shock (Wang et al. 2005, 2001b; Kobayashi et al. 2007). In all cases a full description of the high-energy component features has not been achieved.

In this paper, we present a unified description of the tail and the high-energy component in GRB 980923 through forward and reverse shocks, accounting for their energy, spectral indices, fluxes, and duration. In order to do so consistently we find that a magnetized outflow is required.

2. GRB 980923

GRB 980923 was observed by the Burst and Transient Source Experiment (BATSE) on 1998 September 23 at 20:10:52 UT for 32.02 s. It was localized to 234° with respect to the pointing-axis direction of CGRO. In agreement with the light curve given by González et al. (2009, 2012; Figure 1), the event consists of three components. The first is related to the typical prompt emission, the second one to a smooth tail which lasts ∼400 s, and the last to a hard component extending up to ∼150 MeV. The smooth tail was well described by Giblin et al. (1999) as the evolution of a synchrotron cooling break in the slow-cooling regime starting at t0 = 32 s. This conclusion was based on two different aspects. First, the time-dependent break energy was modeled with a power law of the form E0(tt0)δ and for t0 = 32.109 s, they obtained δ = −0.52 ± 0.12. Second, in agreement with the computed value from the spectral index p = 2.4 ± 0.11 and the value of t0, they obtained the relationship between the temporal and spectral indices and observed that the index values were indeed closer to the slow cooling regime. Thus, Giblin et al. (1999) identified the evolution of the spectrum in the tail of the burst as the evolution of a synchrotron cooling break in the slow-cooling regime, implying also that the transition from fast to slow cooling could take place on short timescales, comparable to the duration of the burst. On the other hand, González et al. (2009, 2012) described the high-energy component peaking at t = 20 s as a power law with spectral index, γ of 1.439 ± 0.0687 and flux F = 49.2 ± 3.8 × 10−6 erg cm2 s−1 and also pointed out that the tail could begin before or at least about 14 s after the burst trigger. The last argument is mainly based on the monotonic time evolution of the low-energy spectral index (αBAND) during the main prompt episode, except around ≃14 s, when the value of αBAND is consistent with the value during the tail. We take the results of González et al. (2009, 2012) to develop a unified model where both components are related.

Figure 1.

Figure 1. BATSE SD7 count rates for GRB 980923 are plotted in four energy ranges. A peak at 20 s is apparent and stands out particularly at the low and high energy ranges.

Standard image High-resolution image

3. DYNAMICS OF THE FORWARD AND REVERSE SHOCKS

In this section, we first extend the work done by Giblin et al. (1999) on the smooth tail, and then proceed to compute the energy range for SSC emission from a thick shell of the reverse shock fireball to account for the hard component. We use the fact that the tail could have begun before or at least about 20 s (before the high-energy component) after the burst trigger, suggesting that before this time there already existed forward and reverse shocks, and find that the characteristics of the observed hard emission noted by González et al. (2009, 2012) could be accounted for consistently by a scenario where the reverse shock becomes relativistic during its propagation, with magnetization parameter close to unity. The subscripts f and r refer throughout this paper to the forward and reverse shocks, respectively.

3.1. Implications of t0

In general t = 0 is defined by the trigger time of the burst. Because of the connection between the very early afterglow (smooth tail) and this prompt emission, we plotted the power index δ as function of t0 (Figure 2) for the data and the power law of the form E0(tt0)δ (Giblin et al. 1999). It is clear that δ goes from a fast cooling regime to a slow cooling regime depending on the choice of t0. Also, we observe as a particular case that for t0 = 32 s, δ ∼ 0.52 is obtained (i.e., slow cooling). However, González et al. (2009, 2012) point out that the tail could begin before or at least about 14 s after the burst trigger. From Figure 3 we observe that the data used to define t0 and δ are not very restrictive, so values of t0 = 14 s and δ ∼ 0.9 are consistent with the data as well. In other words, the data cannot uniquely define t0, and a value smaller than 32 s is possible. Giblin found that after 32 s the spectrum of the tail appears to mimic the evolution of a synchrotron cooling break in the slow-cooling regime. This then implies that the transition from fast to slow cooling must take place in the first 32 s of the burst. The importance of the suggested value by Gonzalez for t0 = 14 s lies in the fact that it allows causality between the high-energy component and the tail. Our model assumes that both components have a common origin. Thus the exact value of t0 is not relevant as long as it is less than 20 s, when the high-energy component is apparent. A transition time from fast to slow cooling regime comparable to the duration of the burst would support our assumption (we calculate it in the following section). Hence, choosing t0 = 14 s gives a δ value in the transition from fast to slow cooling that can be explained if the transition time is sufficiently short so that the spectra would never be observed in the fast cooling regime (this is the case we calculate later). We can assume that about this time the GRB ejecta collides with the ISM generating forward and reverse shocks.

Figure 2.

Figure 2. Power-law decay index δ is plotted as a function of t0 for a power law of the form E0 (tt0)δ, where t0 is the time at which the tail begins. The error bars are obtained for the different fits, as shown in Figure 3. The index is sensitive to the value t0, and thus the choice of t0 is crucial and illustrates the effect on the light curves and therefore on the apparent evolution of the spectrum. The expected values for the slow and fast cooling regimes are indicated.

Standard image High-resolution image
Figure 3.

Figure 3. Break energy as a function of time t on logarithmic axes. The continuous line represents the fitted power-law slope 0.5186 ± 0.1113 for t0 = 32 s (Giblin et al. 1999), while the dashed line is the fitted power-law slope 0.9039 ± 0.1999 for t0 = 14 s. t = 0 corresponds to the burst trigger.

Standard image High-resolution image

3.2. Smooth Tail from Synchrotron Radiation Forward Shock

For the forward shock, we assume that electrons are accelerated in the shock to a power-law distribution of Lorentz factor γe with a minimum Lorentz factor γm: Ne)dγe∝γpedγe, γe ⩾ γm and that constant fractions epsilone, f and epsilonB, f of the shock energy go into the electrons and the magnetic field, respectively. Then

Equation (1)

where we have used the value of p = 2.4 ± 0.11 as obtained by Giblin et al. (1999). Adopting the notation of Sari et al. (1998) and ignoring self-absorption, the observed spectral flux in the fast-cooling regime is given by

Equation (2)

Similarly, the flux in the slow-cooling regime can be written as

Equation (3)

Using the typical parameters given by Björnsson (2001), we compute the typical and cooling frequencies of the forward shock synchrotron emission (Sari et al. 1998) which are given by

Equation (4)

where the convention Qx = Q/10x has been adopted in cgs units throughout this document unless otherwise specified. ttr, f is the transition time, when the spectrum changes from fast cooling to slow cooling, D is the luminosity distance, nf is the ISM density, t is the time of the evolution of the tail, E is the energy, and the term (1 + xf) was introduced because a once-scattered synchrotron photon generally has energy larger than the electron mass in the rest frame of the second-scattering electrons. Multiple scattering of synchrotron photons can be ignored. xf is given by Sari & Esin (2001) as

Equation (5)

where η = (γc, fm, f)2 − p for slow cooling and η = 1 for fast cooling.

From Equation (4), we observe directly that νm, f ⩽ νc, f, the break energy Ec, f ∼ 124.1 keV is consistent with the values given by Giblin et al. (1999) and ttr, f ∼ 8.7 s, implying also that our assumption that the tail and the high-energy component have a common origin is plausible, and that the transition from fast to slow cooling could take place on very short timescales, comparable to the duration of the burst, as expected.

3.3. X-Ray Flare from Thick Shell Reverse Shock

For the reverse shock, it is possible to obtain a simple analytic solution in two limiting cases, thin and thick shell (Sari & Piran 1995) by using a critical Lorentz factor Γc:

Equation (6)

where T90 is the time of the GRB, which is much larger that the peak time of the reverse shock emission and nr is the thick shell density. We consider the thick shell case in which the reverse shock becomes relativistic during the propagation and the shell is significantly decelerated by the reverse shock. Hence, the Lorentz factor at the shock crossing time tc is given by γd ∼ Γc (Kobayashi et al. 2007; Kobayashi & Zhang 2007), and for σ ∼ 1, where σ = Lpf/Lkn = B2r/4πnrmpc2Γ2r is the magnetization parameter, defined as the ratio of Poynting flux to matter energy flux, the crossing time tc is much shorter than T90, tcT90/6 (Fan et al. 2004; Fan 2008; Zhang & Kobayashi 2005; Drenkhahn 2002). Now, if the constant fractions, epsilone, r and epsilonB, r, of the reverse shock energy go into the electrons and magnetic field, respectively, we have

Equation (7)

where γr is the Lorentz factor of the thick shell. The spectral characteristics of the synchroton emission of the forward and reverse shocks are related (Zhang et al. 2003; Kobayashi et al. 2007; Fan & Wei 2005; Fan et al. 2004; Jin & Fan 2007; Shao & Dai 2005) by

Equation (8)

where $\mathcal {R}_B=\epsilon _{B,f}/\epsilon _{B,r}$, $\mathcal {R}_e=\epsilon _{e,r}/\epsilon _{e,f}$, $\mathcal {R}_x=(1+x_f)/(1+x_r+x_r^2),$ and $\mathcal {R}_M=\Gamma ^2_c/\gamma$. The previous relations tell us that including the re-scaling there is a unified description between both shocks (forward and reverse), and the distinction between forward and reverse magnetic fields considers that in some central engine models (Usov 1992; Mészáros & Rees 1997b; Wheeler et al. 2000) the fireball wind may be endowed with "primordial" magnetic fields. Also as the cooling Lorentz factor must be corrected, then $\mathcal {R}_x$ is introduced as a correction factor for the IC cooling, where xr is obtained by Kobayashi et al. (2007) as

Equation (9)

For fast cooling, we take η = 1, and hence with the standard values for epsilonB, r and epsilone, r, ηepsilonB, r/epsilone, r ∼ 4.8. Using Equations (4) and (8), the typical and cooling frequencies of the reverse shock synchrotron emission are

Equation (10)

From Equation (10) we see that νm, r and νc, r correspond to optical and IR frequencies, respectively, and that νm, r, which characterizes the frequency band, does not depend on xr. However these energies were not recorded. Instead, as higher energy photons were observed we compute the upscattering emission of the synchrotron radiations (Equations (2) and (3)) (rf) by relativistic electrons (fast cooling and slow cooling; Sari & Esin 2001). So the SSC spectrum in the fast cooling regime is

Equation (11)

and in the slow cooling regime we find

Equation (12)

ν(IC), ν(IC)c, and F(IC)max are given by

Equation (13)

where k = 4(p − 1)/(p − 2) and τ = (σTNe/4πRd) = c/3/(1 + z/2)−1 σTn Γ4c γ−1 is the optical depth of the shell. In agreement with Equations (10) and (13), in the self-Synchrotron Compton we have

Equation (14)

From Equation (14) we observe that the break energies and (νF)max = 21.2 × 10−6 erg cm−2 s−1 are within the range pointed out by González et al. (2009, 2012).

4. DISCUSSION AND CONCLUSIONS

As shown in Figure 2, the power-law decay index is sensitive to the chosen value of t0, so depending on our choice, we will apparently be closer to a fast or slow cooling regime. Choosing t0 ⩽ 20 s (González et al. 2009, 2012; Sacahui et al. 2011) and assuming that from this time until t ∼ 32 s the synchrotron emission was eclipsed by the prompt phase, we have calculated the transition time from fast to slow cooling as ∼8.7 s. This value is consistent with a common origin of the tail and the high-energy component, as well as with Giblin's results. If t0 = 14 s, we could also explain the diminishing flux at ≃ 14 s as due to pair production (γγ → e+e) between prompt emission and forward shock photons (Equation (4)) at the beginning of the afterglow. We have calculated the energy of the forward shock photons as 1.31 MeV and the time for this energy to decrease under the pair production threshold as 1.5 s, which is consistent with duration of the diminishing flux.

In the reverse shock, the synchrotron process emitted photons with νc, r ∼ 1.0 × 1010 Hz and νm, r ∼ 4.6 × 1016 Hz, which were not recorded but were upscattered by electrons up to break energies EICc ∼ 4.2 × 10−5 eV, EICm ∼ 427.9 MeV with a (νF)max = 21.2 × 10−6 erg cm−2 s−1, which were pointed out by González et al. (2009, 2012) and Sacahui et al. (2011). Now, in accordance with the observed value for the spectral index γ ∼ 1.44 ± 0.07, ν(IC)m, r > ν(IC)c, r we conclude the SSC spectrum corresponds to fast-cooling regime, very similar to GRB 941017 (González et al. 2003; Granot & Guetta 2003). For our case (thick case), the flare occurs during the prompt gamma-ray phase.

From the value of $\mathcal {R}_B$, we obtained that the forward and reverse magnetic fields are related by Bf = 0.9 × 10−3Br. The previous result indicates that there is a stronger magnetic field in the reverse-shock region than in the forward-shock region, which may suggest that the obtained results are given when the ejecta is highly magnetized, as in the interpretation of the early afterglow of GRB 990123 and GRB 021211 provided by Zhang et al. (2003).

Finally, because the Large Area Telescope covers the energy range from about 20 MeV to more than 300 GeV, we hope to detect other hard components in GRBs and so further constrain this model.

The current model accounts for the main characteristics of the burst: energies, spectral indices, fluxes, duration of the main components in a unified manner. The main requirements are that the ejecta be magnetized, leading to the formation of a reverse shock and that both the tail and the high-energy component have a common origin. The model has eight free parameters (equipartition magnetic field, equipartition electron energy, Lorentz factor, and densities all of them in the reverse and forward shocks), with standard values. It is important to point out that our model is not dependent on the value of the t0 as long as it is less than 20 s. The main difference between our model and previous interpretations (Wang et al. 2005, 2001b; Kobayashi et al. 2007) are the assumption of different equipartition values in the forward and reverse shocks, leading to the conclusion of magnetization of the jet.

This burst has similar characteristics to GRB 090926A (Ackermann et al. 2011), but the high energy extended emission component may require SSC from the forward shock J. R. Sacahui et al. (2012, in preparation).

We thank Enrico Ramírez-Ruiz and Charles Dermer for useful discussions as well as the referee for a critical reading of the paper and valuable suggestions. This work is partially supported by UNAM-DGAPA PAPIIT IN105211 (M.G.) and CONACyT 103520 (N.F.) and 83254 (W.L.).

Please wait… references are loading.
10.1088/0004-637X/751/1/33