The initial Lorentz factors of fireballs inferred from the early X-ray data of SWIFT GRBs
R.-R. Xue1,2 - Y.-Z. Fan3,1 - D.-M. Wei1,4
1 - Purple Mountain Observatory, Chinese Academy of Sciences, Nanjing 210008, PR China
2 - Graduate School, Chinese Academy of Sciences, Beijing, 100012, PR China
3 - Niels Bohr International Academy, Niels Bohr Institute, University of Copenhagen, Blegdamsvej 17, 2100 Copenhagen, Denmark
4 - Joint Center for Particle Nuclear Physics and Cosmology of Purple Mountain Observatory - Nanjing University, Nanjing 210008, PR China
Received 17 November 2008 / Accepted 10 February 2009
Aims. We intend to determine the type of circumburst medium and directly measure the initial Lorentz factor of GRB outflows.
Methods. If the early X-ray afterglow light curve has a peak and the whole profile across the peak is consistent with the standard external shock model, the early rise profile of light curves can be used to differentiate whether the burst was born in the interstellar medium (ISM) or in a stellar wind. In the thin shell case related to a subrelativistic reverse shock, the peak time occurring after the end of the prompt emission can be used to derive an accurate , especially for the ISM case. The afterglow light curves for a flat electron spectrum 1<p<2 have been derived analytically.
Results. In our GRB sample, we obtained for the bursts born in ISM. We did not find any good case for bursts born in stellar wind and behaving as a thin shell that can be used to constrain reliably.
Key words: gamma rays: bursts - ISM: jets and outflows - radiation mechanisms: non-thermal
Gamma-ray bursters are among the most mysterious celestial objects that have interested people since its first detection in 1967 (Klebesadel et al. 1973). The time variability of pulses, as short as a millisecond, limits this event to an object on a stellar scale. The random occurrence and the short time-duration of this kind of event has led to difficulties in detection. The dark era of research on gamma-ray bursts (GRBs) lasted until the release of X-ray afterglow data of GRB 970228, confirming GRBs at the cosmological distances (Costa et al. 1997).
The power-law decay of multi-waveband afterglows of many GRBs are consistent with the standard external shock model (Wijers et al. 1997; Waxman 1997). However, the multi-waveband afterglows were usually monitored several hours after the burst trigger before Swift era. The late afterglow, independent of the initial values of the fireball, cannot provide information about the fireball characteristics. X-ray Telescope (XRT) and Ultra-Violet Telescope (UVOT) onboard the Swift satellite (Gehrels et al. 2004) and other ground-based telescopes can slew to GRB within ten seconds and then begin observations. The early afterglow data released in Swift era provide us an opportunity to study properties of fireballs, e.g., the initial Lorentz factor of the fireball.
The fireball is expected to be a highly relativistic ejection from the central engine to avoid the ``compact problem'' (Shemi & Piran 1990; Lithwick & Sari 2001). After the radiation-dominated acceleration phase, the fireball goes into a matter-dominated phase when the fireball is no long accelerated. The fireball keeps an approximately invariable velocity until it sweeps up a considerable mass of ambient medium. We call this episode the coasting phase (Piran et al. 1993). Although the profile of the early afterglow appears quite different from burst to burst and also from the X-ray to infrared bands, the peaks in the early afterglow light curves may indicate the arrival of the deceleration radius in some GRBs. For example, Molinari et al. (2007) attributed the peaks in near-infrared afterglows of GRB 060418 and GRB 060607A to the end of the coasting phase and determined the initial Lorentz factors of the fireballs (see also Jin & Fan 2007). Recently, Oates et al. (2009) have analyzed the early afterglows of Swift-UVOT data and measured the initial Lorentz factor of GRBs for those showing an early power-law increase in flux.
In contrast to these works, we now use the early X-ray data of Swift GRBs to constrain the initial Lorentz factors (). As a probe of , the X-ray data is better than the optical data for the following reasons. (1) In the standard fireball model, the X-ray emission quickly decays with time after the outflow has decelerated, independent of the profile of the medium surrounding the progenitor (Fan & Wei 2005). This occurs because both the typical frequencies of the forward shock (FS) and the reverse shock (RS) emission are usually below the X-ray band (see Table 1 for the light curves). The optical emission, however, will increase until the typical synchrotron frequency of the FS drops below the optical band (Sari et al. 1998). (2) In the thin shell case of interest, usually the RS X-ray emission is not strong enough to outshine the FS emission component. The origin of the X-ray peak can thus be reliably established. (3) X-ray afterglows are hardly influenced by the self-absorption effect and dust extinction, different from the emission at lower frequencies. One disadvantage of our method is that the early X-ray emission of most GRBs have been polluted by the delayed flares (Falcone et al. 2007) powered by the prolonged activity of the central engine. Fortunately, the X-ray flares usually have a decline as steep as , which is significantly sharper than what the fireball model predicts, so one can distinguish the peak in FS emission from the peak of the flare in the X-ray band convincingly.
A very bright optical flash has been detected in GRB 990123 (Akerlof et al. 1999). The most widely discussed interpretation is the external RS model (Sari & Piran 1999; Mészáros & Rees 1999 however, it can also be produced by the internal shock model; e.g. Wei 2007). Since then, the RS emission in the optical band has been extensively investigated (see Zhang 2007, for a review). However, the RS X-ray emission has just been calculated by a few authors (Fan & Wei 2005; Zou et al. 2005). In this work we focus on the profile of early X-ray afterglow light curves in different cases.
In the Fermi acceleration process, the power-law index of shocked electrons p>2 is the result (Gallant 2002). It has been taken as the standard scenario and widely used in the afterglow calculation (Sari et al. 1998; Chevalier & Li 2000). Some afterglow modeling (Dai & Cheng 2001; Bhattacharya 2001), however, favor a flat electron spectrum 1<p<2, for which a reliable estimate of the afterglow emission is still unavailable. We discuss such a scenario in some detail in Sect. 2.2.
First, we discuss the ``thin shell case'' for which , where is the crossing time of RS and T90 is the duration of the burst, referring to a subrelativistic RS. Assuming typical parameters (e.g., the fraction of the shock energy given to the electrons , the fraction given to the magnetic field , the total energy of the fireball and p=2.3), we have the typical synchrotron radiation frequency and the cooling frequency of the FS emission , for bursts born in ISM (Sari et al. 1998) and , born in the wind medium (Chevalier & Li 2000), where n and A* are two medium parameters: n indicates the density of ISM, in units of , while A* is the dimensionless parameter of stellar wind environment. t indicates the observed time postburst. The convention Qx=Q/10x has been adopted in this paper, in units of cgs.
we find that both
of FS and RS emission (marked by the subscripts f and r, respectively) are below the X-ray band
Hz, assuming typical parameters. If the shock parameters are similar for the FS and the RS, the flux contrast between the RS and FS X-ray emission is thus (Fan & Wei 2005)
where the subscript represents the parameters measured at , and and represent the initial Lorentz factor and the instant Lorentz factor at , respectively. The corresponding relative Lorentz factor is .
In the so-called ``thin shell case'',
For a typical ,
Eq. (1) gives
for , which suggests that the RS X-ray emission can be ignored. This conclusion is unchanged if is actually above , because in such a case the FS X-ray emission would be stronger.
In the ``thick shell case'' (for which
), particularly for a relativistic RS satisfying
As a result, the RS X-ray emission may be able to outshine the FS component.
Kobayashi 2000) when . The simultaneous FS X-ray emission is . In the wind case, both the FS and the RS emission decrease slowly with time as t(2-p)/2 when (see also Fan & Wei 2005).
For , it is well known that , independent of the type of circumburst medium for both thin shell and thick shell cases.
Please note that we assume in our above analysis that both and of FS and RS are well below the XRT band. More general results have been summarized in Table 1. One can see that is also possible but only for a fast-cooling forward shock. Its spectrum should be , which can be distinguished from the shallow decline predicted for the FS and RS emission in the wind case at time .
Table 1: The temporal behavior of the X-ray afterglow light curves in the case of p>2.
With the shock jump conditions, Bhattacharya (2001) derived a minimum Lorentz factor ( ) of the electrons that depends on the maximal one ( ), i.e., , where is the Lorentz factor of the shock, and are the rest mass of protons and electrons, respectively. However, in reality, particle acceleration proceeds from low to high energy. Both and should be determined by the first shock-crossing and by radiative losses or escape from the acceleration region, respectively; hence, should have no ``knowledge'' of (Panaitescu 2006, private communication). Motivated by the above arguments, we assume in this work that , where f(p) is a function of p. Such treatment requires that only a small fraction of the upstream material has been accelerated; otherwise, the energy momentum conservation law will be violated.
Assuming that the shock-accelerated electrons have a power-law energy distribution
with the shock jump conditions that
We also assume that the maximal Lorentz factor is limited by the synchrotron losses and is given by (Cheng & Wei 1996)
where is the magnetic field strength of the shock region, the number density of the upstream medium measured in its rest frame, and c the speed velocity of light. The afterglow light curves in the case of 1<p<2 are different from those presented in Table 1 by a factor of given below.
Assuming the X-band is above
the flux contrast between the RS and the FS emission can be estimated (for a thin shell) by
As a result, in the thin shell case, the RS X-ray emission is usually outshone by the FS X-ray radiation for both a flat electron spectrum (1<p<2) and a standard electron spectrum (p>2).
For the FS before deceleration, we have , , and . For the GRBs born in ISM in the case of a thin shell, we have (Fan et al. 2002), , and for the RS. In the case of a thick shell, for the RS we have (Kobayashi 2000), , and . For the GRBs born in the wind and for the RS we have , , and .
For the FS taking a Blandford-McKee profile, we have , , and . As long as the edge of the ejecta is visible, the sideways expansion may be so strong that it cannot be ignored, and the dynamics is governed by (Rhoads 1999). We then have , independent of k. The detailed light curves are summarized in Table 2.
Table 2: Temporal behavior of the X-ray afterglow light curves in the case of 1<p<2.
In the thick shell case, the RS disappears almost simultaneously with the prompt X-ray and ray emission. As a result, the RS X-ray emission may be outshone by the prompt emission component. Moreover, for a thick shell, we need self-consistent modeling of the FS and RS emission to get and then . In such a process, quite a few free parameters are introduced and the constraint on is uncertain. In this work, we only focus on the thin shell case with .
As summarized in Tables 1 and 2, the outflow expanding into the ISM in the thin shell case will give rise to an increase no shallower than t2, while expanding into the wind cannot account for an increase steeper than t1/2. Therefore we can judge whether the GRB was born in the ISM or wind medium according to the sharpness of flux increase. One can also speculate that the X-ray data for the GRBs born in the wind are not very good probes of the initial Lorentz factor of the ejecta due to the lack of a distinguished peak for typical parameters that give rise to .
The number of GRBs recorded by XRT exceeds 160 till August 1, 2008. Most of them, however, play no role in constraining . For our purpose, the light curves have to be characterized by (i) a distinguished/single peak; (ii) across the peak, a smooth transition to a single power-law decay is followed and the whole profile must be consistent with the standard afterglow model (Sari et al. 1998; Chevalier & Li 2000).
We checked all the Swift X-ray afterglows and find 7 GRBs whose early peaks may be attributed to the arrival of , and the post-peak temporal behaviors of these bursts are consistent with the standard afterglow model. All these bursts except GRB 051111 have the X-ray afterglow peaks much later than the end of the prompt emission and thus fall in the thin shell case. A good example is that of GRB 080319C (see Fig. 1). We divide our GRB sample, in total 6 events, into two groups: 4 were born in ISM and 2 that might be born in wind, based on the rise behavior of the early X-ray afterglow light curve.
|Open with DEXTER|
Cucchiara et al. (2006). The XRT began taking data 63 s after the BAT trigger (Racusin et al. 2006a). The flux increased with time in the form of and peaked 110 s after the BAT trigger, followed by a decay of with time to 400 s from the trigger (Racusin et al. 2006b). The XRT spectrum is . The decline and the spectrum are consistent with the fireball afterglow model, supposing that and (see Table 1). It is straightforward to show that, for , the FS X-ray emission still dominates the RS X-ray component.
Holland et al. (2006) considered this peak as a flare, the post-peak flux decayed temporally in the same way as the late time afterglows. D'Elia et al. (2006) claim a redshift of z=3.20 based on a damped Ly line and various absorption features. The XRT spectrum is . The decline and the spectrum are consistent with the fireball afterglow model supposing that and (see Table 1).
Pagani et al. 2008). The redshift z=1.95 was found from absorption lines of the afterglow spectrum (Wiersema et al. 2008). As shown in Fig. 1, a quick increase in flux is obvious up to s after the GRB trigger. After the peak time, the late afterglow exhibited a broken power law - a relatively shallow decline followed by a profile. We attribute the late steep decline to the jet effect of this GRB. With a flat electron spectrum , both the temporal behavior of light curves in the two late episodes and the X-ray spectrum can be reproduced provided that (see Table 2).
Stamatikos et al. 2008b). The early XRT data recorded a steep increase in flux, peaking around 167 s from the trigger. The X-ray light curve is well-fitted by a simple power law, with a decay slope of 0.88 0.06 from 167 s to 104 s after the burst trigger (Troja & Stamatikos 2008). Based on the detection of numerous absorption features, including Fe II, Mg II, and Mg I lines, Vreeswijk et al. (2008) infer a redshift z = 1.10. The XRT spectrum is . The decline and the spectrum are roughly consistent with the fireball afterglow model, supposing that and (see Table 1).
2008). The X-ray spectrum is . As shown in Table 1 in the thin shell case for , the FS emission of an ejecta expanding into the stellar wind can give rise to an initial rise t0.5, roughly consistent with the data. However, this model gives a decline t-1/4 after the X-ray peak, deviating from the data significantly, so this event won't be included in our constraint.
Holland et al. 2008). The XRT spectrum is . Similar to GRB 080307, an initial rise t0.6 may be the FS emission of an ejecta born in wind. However, this model gives a decline t-1/4 after the X-ray peak unless we see the edge of the ejecta for which the decline can be as steep as t-1 (see Table 1). Again, we do not include it in the following constraint.
Table 3: GRB sample.
In our GRB sample, the peaks appearing in the early X-ray afterglow light curves are related to the position of the outflow
from the central engine. At this time, the swept medium has an energy comparable to that of the GRB outflow, and the instant Lorentz factor at
is about half of the initial one. Based on the assumption above, we can measure the initial Lorentz factor of GRBs in ISM case (Blandford & Mckee 1976):
where is the isotropic energy of prompt gamma-ray emission. Here, we take the radiation efficiency = 0.2 in the calculation according to Guetta et al. (2001) and Molinari et al. (2007). The densities in ISM case are assumed. Since is weakly dependent on the unknown density of the cicumburst medium in the ISM case, it can be measured relatively accurate. The results have been presented in Table 3. Usually we have in the case of ISM with a number density , consistent with some other independent probes (Molinari et al. 2007; Jin & Fan 2007; Zou et al. 2009; Pe'er et al. 2007; Zhang et al. 2006; Lithwick & Sari 2001; Sari & Piran 1999; Zou & Piran 2009), some of which are independent of the profile of the circumburst medium.
We used the current Swift-XRT data to constrain the profile of the circumburst medium and then measured the initial Lorentz factor of the fireballs. As a reliable probe, the X-ray light curves should have the following characters for our purpose. (i) There is a distinguished peak. (ii) Across the peak, a smooth transition to a single power-law flux is followed, and the whole profile must be consistent with the standard afterglow model. The early peaks accompanied by steep decay, usually steeper than t-3, are abandoned. Among the Swift bursts we checked, only 4 events meet such requests, because in most events the early X-ray emission are polluted by the emission powered by the prolonged activity of the central engine. In all 4 bursts, the initial increase in the X-ray flux is quicker than t2, strongly suggesting a constant low-density medium. In general, we find , consistent with the constraints obtained in other analysis. In this work all the events in the sample are bright GRBs. For the nearby sub-luminous events, such as GRB 980425, GRB 031203, GRB 060218, and GRB 060614, a reliable estimate of their initial Lorentz factors is still infeasible at present.
In our analysis we did not find a good case for the burst born in a stellar wind and behaving as a thin shell (see Sect. 3.2). One possible reason is that usually the wind medium is so dense that the outflow has decelerated significantly on a timescale t<T90, so the X-ray data is likely to be unsuitable as a probe of the initial Lorentz factor of GRB outflows for most of events occurring in the wind medium.
Note added in proof After the acceptance of the paper, the details of the X-ray afterglow data of GRB 090113 became available (Krimm et al. 2009): The 0.3-10 keV light curve shows an initial period of roughly constant emission. For t>530 s, the light curve can be modeled with a power-law decay as t-1.3. The X-ray spectrum is . One can see that, with , both the temporal and the spectral behaviors of this X-ray afterglow are consistent with the forward shock emission model in the case of a wind medium and . The initial Lorentz factor can be estimated as , where A is related to the regular wind parameter A* by A=3 1035A*. Assuming , we have 1051 erg and .
We thank an anonymous referee for helpful comments and Nathaniel R. Butler for providing some of the data used in this work. This work made use of data supplied by the UK Swift Science Data Centre at the University of Leicester. This work is supported by the National Science Foundation (grants 10673034 and 10621303) and National Basic Research Program (973 programs 2007CB815404 and 2009CB824800) of China. Y.Z.F. is also supported by Danish National Research Foundation and by Chinese Academy of Sciences.
- Akerlof, C., Balsano, R., Barthelmy, S., et al. 1999, Nature, 398, 400 [NASA ADS] [CrossRef] (In the text)
- Bhattacharya, D. 2001, Bull. Astr. Soc. India, 29, 107 [NASA ADS]
- Barthelmy, S. D., Baumgartner, W., Cummings, J. R., et al. 2008, GCN Circ., 7606
- Blandford, R. D., & McKee, C. F. 1976, Phys. Fluids, 19, 1130 [NASA ADS] [CrossRef] (In the text)
- Butler, N. R., Kocevski, D., Bloom, J. S., & Curtis, J. L. 2007, ApJ, 671, 656 [NASA ADS] [CrossRef]
- Cheng, K. S., & Wei, D. M. 1996, MNRAS, 283, L133 [NASA ADS] (In the text)
- Chevalier, R. A., & Li, Z. Y. 2000, ApJ, 536, 195 [NASA ADS] [CrossRef]
- Costa, E., Frontera, F., Heise, J., et al. 1997, Nature, 387, 783 [NASA ADS] [CrossRef] (In the text)
- Cucchiara, A., Fox, D. B., Berger, E., & Price, P. A. 2006, GCN Circ., 5470 (In the text)
- Cummings, J. R., Barbier, L., Barthelmy, S. D., et al. 2006, GCN Circ., 5621
- Dai, Z. G., & Cheng, K. S. 2001, ApJ, 558, L109 [NASA ADS] [CrossRef]
- D'Elia, V., Piranomonte, S., Covino, S., et al. 2006, GCN Circ., 5637 (In the text)
- Evans, P. A., Beardmore, A. P., Page, K. L., et al. 2007, A&A, 469, 379 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Falcone, A. D., Morris, D., Racusin, J., et al. 2007, ApJ, 671, 1921 [NASA ADS] [CrossRef] (In the text)
- Fan, Y. Z., & Wei, D. M. 2005, MNRAS, 364, L42 [NASA ADS] (In the text)
- Fan, Y. Z., Dai, Z. G., Huang, Y. F., & Lu, T. 2002, ChJAA, 2, 449 [NASA ADS] (In the text)
- Gallant, Y. A. 2002, Lecture Notes in Physics, 589, 24 [NASA ADS] (In the text)
- Gehrels, N., Chincarini, G., Giommi, P., et al. 2004, ApJ, 611, 1005 [NASA ADS] [CrossRef] (In the text)
- Guetta, D., Spada, M., & Waxman, E. 2001, ApJ, 557, 399 [NASA ADS] [CrossRef] (In the text)
- Holland, S. T., Barthelmy, S. D., Starling, R., et al. 2006, GCN Rep., 1.2 (In the text)
- Holland, S. T., Palmer, D. M., Starling, R. L. C., & Schady, P. 2008, GCN Rep., 128.1 (In the text)
- Jin, Z. P., & Fan, Y. Z. 2007, MNRAS, 378, 1043 [NASA ADS] [CrossRef] (In the text)
- Klebesadel, R. W., Strong, I. B., & Olson, R. A. 1973, ApJ, 182, L85 [NASA ADS] [CrossRef] (In the text)
- Kobayashi, S. 2000, ApJ, 545, 807 [NASA ADS] [CrossRef] (In the text)
- Krimm, H. A., Evans, P. A., Oates, S. R., et al. 2009, GCN Rep., 193.1 (In the text)
- Lithwick, Y., & Sari, R. 2001, ApJ, 555, 540 [NASA ADS] [CrossRef]
- Mészáros, P., & Rees, M. J. 1999, MNRAS, 306, L39 [NASA ADS] [CrossRef] (In the text)
- Molinari, E., Vergani, S. D., Malesani, D., et al. 2007, A&A, 469, L13 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Oates, S. R., et al. 2009, MNRAS, in press [arXiv:0901.3597] (In the text)
- Pagani, C., Racusin, J. L., Kennea, J. A., et al. 2008, GCN Circ, 7460 (In the text)
- Page, K. L., Osborne, J. P., & Holland, S. T. 2008, GCN Circ, 7376 (In the text)
- Pe'er, A., Ryde, F., Wijers, R. A. M. J., Mészáros, P., & Rees, M. J. 2007, ApJ, 664, L1 [NASA ADS] [CrossRef]
- Piran, T., Shemi, A., & Narayan, R. 1993, MNRAS, 263, 861 [NASA ADS] (In the text)
- Racusin, J. L., Barbier, L. M., Brown, P. J., et al. 2006a, GCN Circ., 5378 (In the text)
- Racusin, J. L., Grupe, D., Morris, D., & Stroh, M. 2006b, GCN Circ., 5382 (In the text)
- Rhoads, J. E. 1999, ApJ, 525, 737 [NASA ADS] [CrossRef] (In the text)
- Sari, R., & Piran, T. 1999, ApJ, 520, 641 [NASA ADS] [CrossRef] (In the text)
- Sari, R., Piran, T., & Narayan, R. 1998, ApJ, 497, L17 [NASA ADS] [CrossRef] (In the text)
- Sato, G., Barbier, L., Barthelmy, S. D., et al. 2006, GCN Circ., 5381
- Shemi, A., & Piran, T. 1990, ApJ, 365, L55 [NASA ADS] [CrossRef]
- Stamatikos, M., Barthelmy, S. D., Cummings, J., et al. 2008a, GCN Circ., 7483
- Stamatikos, M., Beardmore, A. P., Gehrels, N., et al. 2008b, GCN Circ., 7598 (In the text)
- Troja, E., & Stamatikos, M. 2008, GCN Circ., 7608 (In the text)
- Vreeswijk, P. M., Thoene, C. C., Malesani, D., et al. 2008, GCN Circ., 7601 (In the text)
- Waxman, E. 1997, ApJ, 485, L5 [NASA ADS] [CrossRef]
- Wei, D. M. 2007, MNRAS, 374, 525 [NASA ADS] [CrossRef] (In the text)
- Wiersema, K., Tanvir, N., Vreeswijk, P., et al. 2008, GCN Circ., 7517 (In the text)
- Wijers, R. A., Rees, M. J., & Mészáros, P. 1997, MNRAS, 288, L51 [NASA ADS]
- Zhang, B. 2007, ChJAA, 7, 1 [NASA ADS] (In the text)
- Zhang, B., Fan, Y. Z., Dyks, J., et al. 2006, ApJ, 642, 354 [NASA ADS] [CrossRef]
- Zou, Y. C., & Piran, T. 2009, in preparation
- Zou, Y. C., Wu, X. F., & Dai, Z. G. 2005, MNRAS, 363, 93 [NASA ADS] [CrossRef]
- Zou, Y. C., Fan, Y. Z., & Piran, T. 2009, MNRAS, submitted [arXiv:0811.2997]
- ... thus
- This expression is valid in the ``thin shell'' case. In the ``thick shell'' case, such an expression is an upper limit since the total number of the RS electrons is less than times that of the FS electrons (see Eq. (3)).
Table 1: The temporal behavior of the X-ray afterglow light curves in the case of p>2.
Table 2: Temporal behavior of the X-ray afterglow light curves in the case of 1<p<2.
Table 3: GRB sample.
|Open with DEXTER|
|In the text|
Copyright ESO 2009