PRIOR EMISSION MODEL FOR X-RAY PLATEAU PHASE OF GAMMA-RAY BURST AFTERGLOWS

Published 2008 December 19 © 2009. The American Astronomical Society. All rights reserved.
, , Citation Ryo Yamazaki 2009 ApJ 690 L118 DOI 10.1088/0004-637X/690/2/L118

1538-4357/690/2/L118

ABSTRACT

The two-component emission model to explain the plateau phase of the X-ray afterglows of gamma-ray bursts (GRBs) is proposed. One component, which is responsible for the plateau and subsequent normal decay phase of the X-ray afterglow, is the prior emission via outflow ejected from the central engine before the main burst. The other is the main outflow, which causes the prompt GRB emission and the initial steep decay phase of the X-ray afterglow. In this model, the transition from the plateau to the subsequent normal decay phase is an artifact of the choice of the zero of time. For events with distinct plateau phase, the central engine is active 103–104 s before the launch of the main outflow. According to this model, a prior emission in the X-ray and/or optical bands 103–104 s before the prompt GRB emission is possibly seen, which will be tested by near-future instruments such as Monitor of All-sky X-ray Image (MAXI), WIDe-field telescope for GRB Early Timing (WIDGET), and so on.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

The X-Ray Telescope (XRT) onboard Swift has revealed complex temporal behavior of the X-ray afterglows of gamma-ray bursts (GRBs) in the first few hours (Burrows et al. 2005; Tagliaferri et al. 2005; Nousek et al. 2006; O'Brien et al. 2006; Willingale et al. 2007; Zhang et al. 2007a; Liang et al. 2007, 2008). This time window had been largely unexplored before the Swift era, and studies of early afterglows have revealed many questions concerning GRBs, such as the emission mechanism, nature of the central engine, and burst environment (e.g., Zhang et al. 2006; Zhang 2007).

Early X-ray afterglows have three phases,1 which were not predicted by the standard model from the pre-Swift era (Nousek et al. 2006; Zhang 2007). Phase I: Initial Steep Decay Phase. Initially, the X-ray afterglow decays very steeply; the most popular interpretation is that this is the tail emission of the prompt GRB (Kumar & Panaitescu 2000; Zhang et al. 2006; Yamazaki et al. 2006; Zhang et al. 2009), although other possibilities have been proposed (e.g., Zhang et al. 2007a). Phase II: Plateau Phase. At several hundreds of seconds after the burst trigger, this phase begins until ∼103–104 s; its origin is quite uncertain. This is the main topic of this Letter. Phase III: Normal Decay Phase. After the plateau phase ends, the X-rays subsequently decay with the decay index usually steeper than unity, as expected in the pre-Swift era. This decay behavior is well explained by the classical external shock model (Sari et al. 1998), in which neither the delayed energy injection nor the time dependency of shock microphysics is considered.

Phase II is the most enigmatic in early X-ray afterglows. So far, various kinds of models have been proposed such as the energy injection model (Nousek et al. 2006; Zhang et al. 2006; Granot & Kumar 2006), the inhomogeneous or two-component jet model (Toma et al. 2006; Eichler & Granot 2006; Granot et al. 2006), the time-dependent microphysics model (Ioka et al. 2006; Granot et al. 2006; Fan & Piran 2006), the reverse shock-dominated afterglow model (Genet et al. 2007), the prior activity model (Ioka et al. 2006), the internal engine model (Ghisellini et al. 2007), the cannonball model (Dado et al. 2006), and the dust scattering model (Shao & Dai 2007). In this Letter, another model of phase II is proposed.

2. TWO-COMPONENT EMISSION MODEL

As described in the following, we consider a two-component emission model in which a prior and the main outflows emit X-rays independently (see Figure 1). One component is the prior emission via outflow ejected from the central engine before the main burst. This is responsible for phases II and III of the X-ray afterglow. This emission component arrives at the observer before the main burst triggering prompt GRB detectors (i.e., the onset of GRB) such as BAT onboard Swift. It decays with time simply in a single power-law form

Equation (1)

where A0 is a constant and the time coordinate, t, is measured in the rest frame of the observer. The epoch t = 0 is taken around the time of arrival at the observer of the (unseen) information of the launch of the outflow at the central engine. This kind of a choice of the time zero is seen in many references (e.g., Panaitescu & Kumar 2000; Kobayashi & Zhang 2007). The origin of the prior emission, f0, is not discussed in detail here. It can be either internal engine activity or external shock emission of the outflow.

Figure 1.

Figure 1. Schematic diagram of the model presented in this paper. In terms of the coordinate t, the prior emission component f0(t) takes a single power-law form throughout the burst, but in terms of T(= tT0), the function f0(T) has an artificial plateau phase (see the text for details).

Standard image High-resolution image

We set an observer time, T, where T = 0 corresponds to the onset of the prompt GRB. The interval between the time t = 0 and T = 0 is assumed to be T0 seconds, that is, t = T + T0. Then, one obtains

Equation (2)

which becomes constant if TT0, while $f_0\propto T^{-\alpha _0}$ when TT0. In order to explain phases II and III of the X-ray afterglow, T0 should be 103–104 s, and α0 should be the temporal decay index of phase III. It is noted that the onset time of the prior emission is unknown. The first detectable X-rays from the prior component arrive at the observer in the time range 0 ≲ tT0 (−T0T ≲ 0). This fact will be further discussed in Section 3.5. Another remark is that from Equation (2) alone, the introduction of T0 shifts the origin of time. It has been known that, for phase I, the choice of time zero affects the decay slope (Zhang et al. 2006; Yamazaki et al. 2006). The same argument for phase II is used for the first time in this Letter.

The other component is the main outflow, which causes the prompt GRB emission and the subsequent phase I of the X-ray afterglow. With time coordinate T, phase I of the X-ray afterglow is well approximated by a single power-law model given by

Equation (3)

where A1 is a constant and α1≈ 3–6.

The whole light curve of the X-ray afterglow from phase I to III is described by the sum of the two components introduced above, that is,

Equation (4)

This is our formula for observed X-ray afterglows, where f0(T) describes phases II and III, while f1(T) fits phase I.

We find that by choosing appropriate values of the parameters α0, α1, T0, A0, and A1, the observed light curves of X-ray afterglows are well described with Equation (4). Figure 2 shows some examples of the fit.

Figure 2.

Figure 2. Comparison of the observed light curves of X-ray afterglows with Equation (4), where all phases I, II, and III are well described. The adopted parameters are (α0, α1, T0) = (1.2, 7.0, 8000 s), (1.0, 4.0, 5000 s), (1.4, 4.0, 7000 s), and (1.6, 5.5, 4000 s), for GRB 051016B, 060428A, 060814, and 061121, respectively. Data of X-ray afterglows are taken from the Swift online repository (Evans et al. 2007).

Standard image High-resolution image

It is not surprising that our formula, Equation (4), well explains the observational results of X-ray afterglows. The functional form of Equation (2) is a good approximation of that introduced by Willingale et al. (2007):

Equation (5)

They have shown that the observed light curves of phases II and III are well fitted with Equation (5). The function fa(T) becomes constant with T if taTTa, while $f_a\propto T^{-\alpha _a}$ if TaT. This behavior is quite similar to that of f0(T). Indeed, one can find that over a wide parameter range, fa(T) ≈ f0(T) for taT if we take taTaT0 and αa ≈ α0.

3. DISCUSSION

3.1. Overall Shape of the X-ray Afterglow

In this Letter, we have seen that for most events, phases II and III of the observed light curves of the X-ray afterglow can be well fitted with a very simple formula, Equation (2). In particular, the observational facts are that, for most events, the transition from phase II to III is slow and fairly smooth, and that the X-ray spectrum remains unchanged across the transition (Nousek et al. 2006); this can be naturally explained by our model because the transition from phase II to III is an artifact of the choice of the time zero. A typical exception is GRB 070110; at the end of the plateau phase II, the light curve shows an abrupt drop (Troja et al. 2007). Such events may need other explanations (e.g., Kumar et al. 2008).

After phase III, a subsequent fourth phase is occasionally observed: the so called post jet break phase (Phase IV; Zhang 2007). Its typical decay index is ∼−2, satisfying the predictions of the jetted afterglow model (Sari et al. 1999). This phase IV can also be clearly seen in our sample (e.g., GRB 060428A; upper right of Figure 2). If the prior emission is from an external shock of the prior outflow, the jetted afterglow model is applied to the present case. As will be discussed in Section 3.4, the origin of the optical afterglow is different from the X-ray one, which explains the fact that the epoch of the jet break in the X-ray band is not generally the same as the optical one (Sato et al. 2007; Liang et al. 2008).

3.2. Observed f0(T0)–T0 Correlation

Sato et al. (2008) investigated the characteristics of the transition from phase II to III for 11 events with known redshifts. They derived the transition time, T0brk, and the isotropic luminosity at that time LX,end, and found that LX,end is well correlated with T0brk as LX,end ∝ (T0brk)−1.4. They adopt a broken power-law form to fit the light curve in phases II and III which is different from that considered in this paper. However, we expect that their T0brk and LX,end roughly correspond to T0 and f0(T0), respectively, that are considered in Section 2. Indeed, one finds $f_0(T_0)\propto T_0{}^{-\alpha _0}$ from Equation (2). Hence, if α0 ≈ 1.4, which is a typical number for the decay index of phase III, we can reproduce the observed result of Sato et al. (2008).

A similar analysis has been done but with Equation (5) as the fitting formula (Dainotti et al. 2008). For 32 events with measured redshifts, they found a correlation between Ta and the X-ray luminosity at the time Ta, LX(Ta), as LX(Ta) ∝ T−βa with β = 0.6–0.74. Since TaT0 and LX(Ta) ≈ f0(T0), their correlation indicates f0(T0) ∝ T0−β. The index β is smaller than the typical value of α0 (≈1.0–1.5). However, the claimed correlation has large scatter (see Figures 1 and 2 of Dainotti et al. (2008)), which may be explained by the scatter of α0 in our model.

There is a difference between the results of Sato et al. (2008) and Dainotti et al. (2008). At present, the number of events with known redshifts may be small, so this discrepancy may be resolved if the number of events increases.

3.3. Link Between T0 and the Prompt Emission Properties

Nava et al. (2007) studied the properties of prompt emission and X-ray afterglows of 23 GRBs with known redshifts. They adopted Equation (5) in fitting phases II and III of the X-ray afterglow, and found that for events with measured spectral peak energy Ep, the time Ta weakly correlates with the isotropic equivalent energy Eγ,iso of the prompt GRB emission. One can find from Figure 6 of Nava et al. (2007) that Ta seems to be roughly proportional to Eγ,iso. At present, this correlation is not firmly established because, as noted by Nava et al. (2007), there is no correlation between Ta and the isotropic equivalent energy of prompt GRBs in the 15–150 keV band for a larger sample of GRBs with known redshift but unknown Ep (hence without k-correction).

The quantity Eγ,iso correlates with Ep (Amati et al. 2002). Hence, if the TaEγ,iso correlation exists, the bright GRBs with large Eγ,iso and Ep have large Ta, which is responsible for the distinct plateau phase. On the other hand, the X-ray flashes or X-ray-rich GRBs (e.g., Heise et al. 2001; Barraud et al. 2003; Sakamoto et al. 2008) have small Ta (≈T0), and have X-ray afterglow without phase II. This tendency could have been seen in Sakamoto et al. (2008).

The above arguments may lead a link between long GRBs and X-ray flashes/X-ray-rich GRBs. Suppose that the outflow ejection is not continuous but intermittent, i.e., the central engine ejects two distinct outflows with a time interval of ∼T0. Just after the launch of the prior outflow, the central engine does not have enough energy for another outflow, so that it needs to store an additional one. During the time interval ∼T0, matter surrounding the central engine is accreted, increasing the gravitational binding energy. This energy is released as the main outflow causing the prompt GRB. It is expected that the larger is T0, the larger is the stored gravitational energy, resulting in a brighter burst with large Eγ,iso. On the other hand, if T0 is small the energy of the main outflow becomes small; this is responsible for the X-ray flash or X-ray-rich GRBs. Further details will be discussed in the near future.

3.4. External Shock Emission from the Main Outflow

The main outflow that is responsible for the prompt GRB might cause external shock X-ray emission, fX,ext(T). In the present two-component emission model, however, fX,ext(T) must be dimmer than the prior X-ray emission f0(T) throughout phases II and III. Let us consider the simplest model of external shock emission of the main outflow. The relativistically expanding shell with energy EK interacts with the surrounding medium with uniform density n0,2 and emits synchrotron radiation with microphysics parameters at the shock, p, εe, and εB (Sari et al. 1998). Then, in the case of slow cooling and νc < νX, the X-ray light curves are analytically calculated as fX,ext(T) ∝ E(p+2)/4Kεp−1eε(p−2)/4BT(2−3p)/4νp/2X, which is independent of n0 (e.g., Panaitescu & Kumar 2000). If p ≈ 2, then fX,ext(T) hardly depends on εB. Before the Swift era, typical values had been EK ∼ 1052–1053 erg, εe ∼ 10−1, and εB ∼ 10−2 so that the external shock emission reproduced the observed late-time X-ray afterglow. In the present case, fX,ext(T) must be dim, which implies small EK and/or εe. A similar discussion was made by Ghisellini et al. (2007, 2008). In some models, such as the energy injection model and the inhomogeneous jet model, prompt GRB emission needs high radiation efficiency, which is defined by εγ = Eγ,iso/(Eγ,iso + EK), because EK is small at the epoch of the prompt GRB emission (Fan & Piran 2006; Granot et al. 2006; Ioka et al. 2006; Zhang et al. 2007b). As discussed here, fX,ext(T) can be dim if εe is small while EK remains large, ∼1052–1053 erg. Hence, the present model could avoid a serious efficiency problem. On the other hand, if EK of the outflow is small, the efficiency εγ should be high. Then the mechanism of prompt GRB emission is unlike a classical internal shock model (e.g., Thompson et al. 2007; Ioka et al. 2007).

The observed optical afterglow comes mainly from the prior outflow. For some events, the rising part of the early optical afterglow proceeds until T ∼ 102 s (Molinari et al. 2007), which is difficult to explain with prior emission. In the case of prior emission, the time zero would be shifted T0 ∼ 103–104 s before the burst trigger, making the light curve extraordinarily spiky. Furthermore, in most cases, the transition from phase II to III is chromatic, i.e., the optical light curves do not show any break at that epoch (Panaitescu et al. 2006), although there exist a few exceptions (Liang et al. 2007; Grupe et al. 2007; Mangano et al. 2007). Hence, at least in the early epoch, the observed optical afterglow arises from the main outflow component or others, most likely an external shock emission. Indeed, there have been some observational facts that indicate different origins of X-ray and optical afterglows (e.g., Oates et al. 2007; Sato et al. 2007; Urata et al. 2007).3

3.5. Predicted Precursor Emission?

A possible prediction of the present model is a bright X-ray precursor before the prompt GRB emission. Let us assume that the prior X-ray emission starts at t ∼ 102 s, although its onset time is fairly uncertain (see Section 2). Then, from Equation (1), the X-ray flux in the 2–10 keV band is estimated as f0(t) ∼ 4 × 10−9(t/102 s)−1.2 erg s−1 cm−2, where α0 ≈ 1.2 is taken as a typical value and the flux normalization constant A0 is determined so that f0 ∼ 1 × 10−12 erg s−1cm−2 at t ∼ 105 s (Gendre et al. 2008). Such emission will be detected by Monitor of All-sky X-ray Image (MAXI).4 The expected event rate should be a few events per year (Suzuki et al. 2008).

However, if the emission starts with t ∼ 102 s, the predicted flux might be large enough to be detected by current instruments like BAT onboard Swift. In order to avoid this problem, the starting time of the emission should be comparable to or later than t ∼ 103 s so that the peak flux is smaller than the detection limit of the prompt GRB emission monitors. One possibility is off-axis jet emission in the context of the external shock model (Granot & Kumar 2003; Granot 2005). Due to the relativistic beaming effect, the observed X-ray emission is dim as long as the bulk Lorentz factor of the emitting outflow is larger than the inverse of the angle between the emitting matter and the observer's line of sight. This effect shifts the peak time of the X-ray emission toward the later epoch. For an appropriate choice of parameters, we may adjust the starting time (the peak time) of the X-ray emission.

A signal of the onset of the prior emission might also be seen in the optical band. If the prior emission is extremal shock origin, the reverse shock emission might cause a bright optical flash (Sari & Piran 1999). So far, for some events, WIDGET5 has given an observational upper limit on the prior optical emission, V > 10 mag, about 750 s before the prompt emission (e.g., Abe et al. 2006). Further observations will constrain the model parameters. In summary, in order to test the model presented in this paper, the search for a signal in the data of sky monitors in the optical and X-ray bands is crucial.

R.Y. would like to thank Takashi Nakamura, Kunihito Ioka, Takanori Sakamoto, Atsumasa Yoshida, Motoko Suzuki, Takeshi Uehara, and the anonymous referee for useful comments and discussions. This work was supported in part by grant-in-aid from the Ministry of Education, Culture, Sports, Science, and Technology (MEXT) of Japan, No. 18740153, No. 19047004.

Footnotes

  • For simplicity, the X-ray flares are not considered in this Letter.

  • Since the prior outflow may modify the circumburst medium density profile, the external shock emission from the main outflow deviates from the case of a uniform density profile. Nevertheless, we adopt the uniform density model here for simplicity.

  • One may expect that the optical emission also arises from the prior outflow. However, it may be outshone by the main outflow component at least in the early epoch. This condition will constrain the mechanisms of the prior X-ray emission as well as the optical one.

Please wait… references are loading.
10.1088/0004-637X/690/2/L118