The long-lasting optical afterglow plateau of short burst GRB 130912A

The long-lasting optical afterglow plateau of short burst GRB 130912A

Key Words.:
Gamma-ray burst: general – Gamma-ray burst: individual: 130912A

Abstract

Context:The short burst GRB 130912A was detected by Swift, Fermi satellites, and several ground-based optical telescopes. Its X-ray light curve decayed with time normally. The optical emission, however, displayed a long-term plateau.

Aims:We examine the physical origin of the X-ray and optical emission of short GRB 130912A.

Methods:The afterglow emission was analysed and the light curve fitted numerically.

Results:The canonical forward-shock model of the afterglow emission accounts for the X-ray and optical data self-consistently, so the energy injection model that has been widely adopted to interpret the shallowly decaying afterglow emission is not needed.

Conclusions:The burst was born in a very-low density interstellar medium, which is consistent with the compact-object merger model. Significant amounts of the energy of the forward shock were given to accelerate the non-thermal electrons and amplify the magnetic fields (i.e., and , respectively), which are much more than those inferred in most short-burst afterglow modelling and can explain why the long-lasting optical afterglow plateau is rare in short GRBs.

1 Introduction

Gamma-ray bursts (GRBs) are short flashes of -rays from the outer/deep space. Based on the duration distribution of the prompt emission, GRBs can be divided into two sub-groups (Kouveliotou et al. 1993). One group has a typical duration of . The other has a much shorter duration s, centered at s. The long GRBs are usually related to the death of massive stars, and one smoking-gun signature of such a scenario is the bright supernova emission, as identified in most nearby long GRBs (Woosley 1993; MacFadyen & Woosley 1999; Hjorth et al. 2003). The short GRBs are rarer by a factor of about four according to the BATSE, or about ten from the Swift observation. The rate depends strongly on the energy bands and on sensitivities of the instruments (e.g., Qin et al. 2013, Zhang et al. 2012). The understanding of such a group of events had not been revolutionized until 2005 when Swift and HETE-II localized such events, and the long-lasting multi-wavelength afterglow emission had been detected (Gehrels et al. 2005; Fox et al. 2005).

Although not as abundant as that of long GRBs, the afterglow emission data of short GRBs by then are valuable for revealing the physical processes taking place in the central engine. For example, the peculiar X-ray emission, such as the X-ray flares and the X-ray plateau followed by abrupt quick decline (Villasenor et al. 2005; Barthelmy et al. 2005; Rowlinson et al. 2010; Margutti et al. 2011), was found to be inconsistent with the so-called standard forward-shock afterglow model. This emission instead implies the prolonged activity of the central engines (Fan et al. 2005; Dai et al. 2006; Perna et al. 2006; Gao & Fan 2006; Proga & Zhang 2006; Metzger et al. 2008; Rowlinson et al. 2013), which is possibly associated with non-ignorable gravitational wave radiation (Fan et al. 2013a; Zhang 2013). If these X-ray plateaus are indeed powered the internal energy dissipation of supra-massive neutron star (SMNS) wind and the sharp declines mark the collapse of the SMNSs (Gao & Fan 2006; Rowlinson et al. 2013), then the maximum gravitational mass of non-rotating neutron star can be estimated to be (Fan et al. 2013a; Li et al. 2014). Moreover, the detection of a weak infrared bump in GRB 130603B strongly favors the physical origin of the merger of a compact object binary (Tanvir et al. 2013; Berger et al. 2013; Jin et al. 2013). To account for the shallowly decaying X-ray emission, a millisecond magnetar central engine, which was still active in s after the short burst, is needed, and the progenitor stars are likely to be double neutron stars (Fan et al. 2013b; de Ugarte Postigo et al. 2014; Lü & Zhang 2014; Lü et al. 2015).

Motivated by this progress, in this work we examine the physical origin of the X-ray and optical emission of short GRB 130912A, which is characterized by a long-lasting optical plateau. In section 2 we introduce the observations of GRB 130912A. In section 3 we interpret the data. We summarize our results with some discussion in section 4.

2 Observations

2.1 Swift observations

The Swift Burst Alert Telescope (BAT) detected GRB 130912A at 08:34:57 UT on September 12, 2013 (D’Elia et al. 2013). The duration is second. The light curve shows two overlapping peaks. The time-averaged spectrum of the first 0.32 second is best fitted by a simple power law, with photon index . The total fluence is erg cm , and the peak photon flux is ph cm s (Krimm et al. 2013). All values are in the 15 – 150 keV energy band. There is no evidence of extended emission detected in the BAT energy range (Norris et al. 2010), which makes it an unambiguous short GRB.

The Swift X-ray Telescope (XRT) began to observe the field at 08:36:31.7 UT, 93.9 seconds after the BAT trigger. Kennea et al. (2013) analyzed the initial XRT data and report that the light curve can be modelled with a power-law decay with a decay index of . The spectrum formed from the PC mode data can be fitted with an absorbed power law with a photon spectral index of and cm, consistent with the Galactic value of cm (Kalberla et al. 2005).

The Swift UVOT took 150 seconds to find the chart exposure 98 seconds after the BAT trigger (D’Elia et al. 2013). No optical afterglow within the enhanced XRT position (Beardmore et al. 2013) has been detected in this initial and the subsequence exposures (Chester & D’Elia 2013).

2.2 Ground-based optical observations

The field of GRB 130912A was observed at early times by several ground-based optical telescopes, but only two detected the afterglow. GROND started observing at 08:50 UT on September 12, 2013 (16 min after the GRB trigger), and found the afterglow at coordinates RA (J2000) = 03:10:22.23 and Dec. (J2000) = 13:59:48.7 within the Swift XRT error box. Over the first hour-long observation in the filter, the afterglow seemed to be constant, with a magnitude (here and throughout this paper, magnitudes are in AB system) of r’=22.00.2 (Tanga et al. 2013). Another telescope, P60, started observing at 8:48 UT on September 12, 2013 (13 min after the Swift trigger) also in band, confirmed the GROND observation that the source was almost unchanged in the first hour. It led to the suspicion that the detected source is the host galaxy of GRB 130912A (Cenko et al. 2013); however, when RATIR observed the field again started on 2013 September 13.30, the source had faded down to a limited magnitude 23.89 (3-sigma) in SDSS band, so the previously detected source should be the afterglow (Butler et al. 2013). According to the RATIR observation, the host galaxy is fainter than 23.89, 23.79, 22.86, 22.38, 22.30, and 21.78 magnitudes in , , , , , and bands, respectively (Butler et al. 2013). None of these ground-based observations is publicly available in the form of raw data. We collected the public GCN data as shown in Table 1, and loose limits had been discarded. The magnitudes used in the following analysis were corrected for the Galactic extinction, assuming (Schlegel et al. 1998) and a ratio of total-to-selective extinction , the Galactic extinction in band is .

3 Interpreting the optical and X-ray afterglow

In the Swift-era, long-lasting plateau-like X-ray emission (i.e., the so-called shallow decay phase) was detected well in a good fraction of GRB afterglows (e.g., Zhang et al. 2006; Nousek et al. 2006). The leading interpretation of the long-lasting plateau-like X-ray emission is the energy injection model, which is valid if the central engine works continually, or alternatively, the bulk Lorentz factor of the outflow material has a wide distribution (e.g., Dai & Lu 1998; Zhang & Mészáros 2001; Fan & Xu 2006). A general prediction of the energy injection model is that the temporal behaviours of multi-wavelength afterglow emission will be shaped simultaneously (Fan & Piran 2006). However, the X-ray and optical observations of the GRB afterglow usually do not track each other. The X-ray afterglows usually displayed an early shallow decline, which is often not observed in the optical (e.g., Fan & Piran 2006; Panaitescu et al. 2006). The physical reason is still not clear for such a puzzle, for example, the widely-adopted energy injection model for the X-ray decline is found to be unable to account for the optical data.

Recently, there is another unusual situation. As shown in Fig.1, the X-ray afterglow emission of GRB 130912A can be fitted by a single power law (Beardmore et al. 2013), as found in most GRB X-ray emission. However, the optical emission is plateau-like on a very long timescale s, which is a very strange behaviour that is rarely observed in the optical afterglow. A long-lasting optical plateau was also detected in GRB 090510 (De Pasquale et al. 2010; Gao et al. 2009). The duration of the optical plateau of GRB 090510, however, is just about half of that of GRB 130912A. The optical plateau of GRB 130912A is thus very likely the longest one people have detected so far in short GRB afterglows. In reality, the lack of optical variability in the first hour after the trigger of GRB 130912A motivated the idea that these emission were from the host galaxy (Tanga et al. 2013; Cenko et al. 2013). The afterglow emission nature of the optical data had not been established until significant fading was identified about one day after the burst (Butler et al. 2013). The main purpose of this work is to interpret the “unusual” optical plateau and the X-ray emission self-consistently.

Because nothing is unusual in the X-ray band for GRB 130912A, we conclude that the energy injection model does not work for the current data, as discussed above. It is widely known that the forward-shock emission is governed by some physical parameters that can be parameterized as (e.g., Sari et al. 1998; Yost et al. 2003, Fan & Piran 2006)

 Fν,max=6.6 mJy (1+z2)D−2L,28.34ϵ1/2B,−2Ek,53n1/20, (1)
 νm=2.4×1016 Hz E1/2k,53ϵ1/2B,−2ϵ2e,−1C2p(1+z2)1/2t−3/2d,−3, (2)
 νc=4.4×1016 Hz E−1/2k,53ϵ−3/2B,−2n−10(1+z2)−1/2t−1/2d,−31(1+Y)2, (3)

where , () is the fraction of shock energy given to the electrons (magnetic field), is the time in days, the Compton parameter , , and . Here and throughout this text, the convention has been adopted.

At , if we have , the temporal behaviour of the optical emission would be , and the X-ray emission light curve should be in the case of or, alternatively, in the case of . While the temporal behaviour of optical emission agrees nicely with the data, the X-ray emission does not. The case of is clearly at odds with the data. The case of on the timescale of is also inconsistent with the X-ray spectrum spectra/).

If we, instead, have at , the temporal behaviour of the optical emission would be , and the X-ray emission light curve should be in the case of or, alternatively, in the case of . Now the temporal behaviours of both optical and X-ray emission are consistent with the data for , as are the spectral behaviours. Since an optical plateau lasting a few thousand seconds is rare in the short GRB afterglow data, it is highly necessary to examine whether the required forward-shock physical parameters are reasonable or not.

To self-consistently interpret the optical and X-ray data, Hz, Hz and mJy are needed. Following Zhang et al. (2015), we have

 ϵ1/2B,−2Ek,53n1/20≈a, (4)
 E1/2k,53ϵ1/2B,−2ϵ2e,−1≈b, (5)
 E−1/2k,53ϵ−3/2B,−2n−10(1+Y)−2≈c, (6)

where , , and .

Now we have three relations but four free parameters (i.e., ), which implies that these parameters cannot be uniquely determined. However, can be expressed as the functions of , and it is possible to reasonably constrain the range of , as found in Zhang et al. (2015).

In the case of (i.e., the synchrotron-self Compton cooling is unimportant), can be ignored, and we have (see also Zhang et al. 2015)

 ϵB,−2=a−25c−45n−350, (7)
 ϵe,−1=a−15b12c110n150, (8)
 Ek,53=a65c25n−150. (9)

These three variables weakly depend on .

In the case of (i.e., the synchrotron-self Compton radiation is important), we have and then

 ϵB,−2=(10cd)−85p−9a−2(p−1)5p−9b−4(p−1)5p−9n−3p+15p−90, (10)
 ϵe,−1=(10cd)15p−9a−p−25p−9b3p−55p−9np−15p−90, (11)
 Ek,53=(10cd)45p−9a6p−105p−9b−2(p−1)5p−9n−p−55p−90, (12)

where . Compared with the case of , and depend strongly on .

The redshift is unknown, and we assume that it is in the reasonable range of . However, for in solving eqs.(4-6), we do not find any reasonable values of the parameters with any given , while solutions are obtainable for larger . In Fig. 2 we present the physical parameters as functions of for , respectively. As already shown in the analytical approaches, both and evolve with quickly, while the dependence of on is rather weak. One can also find from the figure that there is a very interesting constraint that for a reasonable ( i.e., expected to be not much larger than the equipartition value ), though the afterglow physical parameters cannot be uniquely determined. GRB 130912A was therefore born in a very low-density medium, consistent with the compact-object merger model.

To better show whether the afterglow model can indeed reasonably account for the data or not, in Fig.1 we numerically fit the optical and X-ray data of GRB 130912A. The numerical calculation code was developed by Fan & Piran (2006) and Zhang et al. (2006). In it, (i) the dynamical evolution of the outflow formulated by Huang et al. (2000) that can describe the dynamical evolution of the outflow for both the relativistic and non-relativistic phases has been adopted. (ii) The energy distribution of the shock-accelerated electrons is calculated by solving the continuity equation with the power-law source function , normalized by a local injection rate (Moderski et al. 2000). (iii) The cooling of the electrons due to both synchrotron and synchrotron-self Compton has been taken into account. Assuming (following Rowlinson et al. (2013) we adopt the average value of redshift of short GRBs), the fit parameters are , where is the half opening angle of the GRB ejecta. An isotropic fireball is found to be unable to reproduce the data. The inferred and are at the high end of the distribution of the shock parameters of short GRB afterglows (e.g., Soderberg et al. 2006; De Pasquale et al. 2010), which is unexpected since the optical afterglow plateau of GRB 130912A has the longest duration people have ever detected in short events. We would like to point out that the above fit parameters are for , which is very high. If is intrinsically smaller, and hence are lowered accordingly. As shown in eqs.(7-9), or alternatively eqs.(10-12), and would increase, while would decrease. The contrary holds in the case of larger .

4 Summary and conclusions

The most remarkable feature of the short burst GRB 130912A is an optical plateau lasting about s, which is the longest one in current short GRB observations, and it is about twice longer than that of GRB 090510. In this work we examined whether any ”unusual” information can be extracted from the afterglow data of GRB 130912A. Though the energy injection model has been widely adopted to interpret the shallowly decaying afterglow emission of long and short GRBs (see Zhang et al. 2006; Nousek et al. 2006 and the references therein), it was found to be unable to account for the X-ray and optical data of GRB 130912A self-consistently. Instead the canonical afterglow emission of an ejecta with an opening angle can reasonably reproduce the data. The circum-burst medium is found to be ISM-like and has a very low density , consistent with the model of merger of binary compact objects (either double neutron stars or a neutron starblack hole). Significant amounts of the energy of the forward shock were given to accelerate the non-thermal electrons and amplify the magnetic fields (i.e., and , respectively), which are much more than those inferred in most short burst afterglow modelling and which can explain why the long-lasting optical afterglow plateau is rare in short GRBs.

Acknowledgements. We thank the anonymous referee for helpful comments and Dr. Y. Z. Fan for stimulating discussion. This work made use of data supplied by the UK Swift Science Data Centre at the University of Leicester, and is supported in part by 973 Program of China under grant 2014CB845800, National Natural Science Foundation of China under grants 11163003, 11273063, 11103084, 11303098, U1331101, and 11361140349, and the Chinese Academy of Sciences via the Strategic Priority Research Programme (Grant No. XDB09000000). F.-W.Z. also acknowledges the support by the Guangxi Natural Science Foundation (No. 2013GXNSFAA019002), the Open Research Programme of Key Laboratory for the Structure and Evolution of Celestial Objects (OP201207) and the project of outstanding young teachers’ training in higher education institutions of Guangxi.

Footnotes

1. email: fwzhang@hotmail.com (FWZ) and jin@pmo.ac.cn (ZPJ)
2. email: fwzhang@hotmail.com (FWZ) and jin@pmo.ac.cn (ZPJ)
3. email: fwzhang@hotmail.com (FWZ) and jin@pmo.ac.cn (ZPJ)
4. email: fwzhang@hotmail.com (FWZ) and jin@pmo.ac.cn (ZPJ)
5. email: fwzhang@hotmail.com (FWZ) and jin@pmo.ac.cn (ZPJ)

References

1. Barthelmy, S. D., Chincarini, G., Burrows, D. N., et al. 2005, Nature, 438, 994
2. Beardmore, A. P., Evans, P. A., Goad, M. R., & Osborne, J. P. 2013, GRB Coordinates Network, 15217, 1
3. Berger, E., Fong, W., & Chornock, R. 2013, ApJ, 774, L23
4. Butler, N., Watson, A. M., Kutyrev, A., et al. 2013, GRB Coordinates Network, 15226, 1
5. Cenko, S. B., Sudilovsky, V., Tanga, M., & Greiner, J. 2013, GRB Coordinates Network, 15222, 1
6. Chester, M. M., & D’Elia, V. 2013, GRB Coordinates Network, 15229, 1
7. Dai, Z. G., & Lu, T. 1998, A&A, 333, L87
8. Dai, Z. G., Wang, X. Y., Wu, X. F., & Zhang, B. 2006, Science, 311, 1127
9. D’Elia, V., Chester, M. M., Cummings, J. R., et al. 2013, GRB Coordinates Network, 15212, 1
10. De Pasquale, M., Schady, P., Kuin, N. P. M., et al. 2010, ApJ, 709, L146
11. de Ugarte Postigo, A., Thöne, C. C., Rowlinson, A., et al. 2014, A&A, 563, A62
12. Evans, P. A., Beardmore, A. P., Page, K. L., et al. 2009, MNRAS, 397, 1177
13. Fan, Y.-Z., & Piran, T. 2006, MNRAS, 369, 197
14. Fan, Y.-Z., Wu, X. F., & Wei, D. M. 2013a, Phys. Rev. D., 88, 067304
15. Fan, Y.-Z., & Xu, D. 2006, MNRAS, 372, L19
16. Fan, Y. Z., Zhang, B., & Proga, D. 2005, ApJ, 635, L129
17. Fan, Y.-Z., Yu, Y.-W., Xu, D., et al. 2013b, ApJ, 779, L25
18. Fox, D. B., Frail, D. A., Price, P. A., et al. 2005, Nature, 437, 845
19. Gao, W.-H., & Fan, Y.-Z. 2006, Chinese J. Astron. Astrophys., 6, 513
20. Gao, W.-H., Mao, J., Xu, D., & Fan, Y.-Z. 2009, ApJ, 706, L33
21. Gehrels, N., Sarazin, C. L., O’Brien, P. T., et al. 2005, Nature, 437, 851
22. Hjorth, J., Sollerman, J., Møller, P., et al. 2003, Nature, 423, 847
23. Huang, Y. F., Gou, L. J., Dai, Z. G., & Lu, T. 2000, ApJ, 543, 90
24. Jin, Z.-P., Xu, D., Fan, Y.-Z., Wu, X.-F., & Wei, D.-M. 2013, ApJ, 775, L19
25. Kalberla, P. M. W., Burton, W. B., Hartmann, D., et al. 2005, A&A, 440, 775
26. Kennea, J. A., Gompertz, B. P., Osborne, J. P., et al. 2013, GRB Coordinates Network, 15221, 1
27. Kouveliotou, C., Meegan, C. A., Fishman, G. J., et al. 1993, ApJ, 413, L101
28. Krimm, H. A., Barthelmy, S. D., Baumgartner, W. H., et al. 2013, GRB Coordinates Network, 15216, 1
29. Li, X., Zhou, B., He, H. N., Fan, Y. Z., & Wei, D. M.  2014, ApJ, 797, 33
30. Lü, H.-J., & Zhang, B. 2014, ApJ, 785, 74
31. Lü, H.-J., Zhang, B., Lei, W.-H., Li, Y., & Lasky, P. D 2015, arXiv:1501.02589
32. MacFadyen, A. I., & Woosley, S. E. 1999, ApJ, 524, 262
33. Margutti, R., Chincarini, G., Granot, J., et al. 2011, MNRAS, 417, 2144
34. Metzger, B. D., Quataert, E., & Thompson, T. A. 2008, MNRAS, 385, 1455
35. Moderski, R., Sikora, M., & Bulik, T. 2000, ApJ, 529, 151
36. Nousek, J. A., Kouveliotou, C., Grupe, D., et al. 2006, ApJ, 642, 389
37. Norris, J. P., Gehrels, N., & Scargle, J. D. 2010, ApJ, 717, 411
38. Panaitescu, A., Mészáros, P., Burrows, D., et al. 2006, MNRAS, 369, 2059
39. Perna, R., Bozzo, E., & Stella, L. 2006, ApJ, 639, 363
40. Proga, D., & Zhang, B. 2006, MNRAS, 370, L61
41. Qin, Y., Liang, E.-W., Liang, Y.-F., et al. 2013, ApJ, 763, 15
42. Rowlinson, A., O’Brien, P. T., Tanvir, N. R., et al. 2010, MNRAS, 409, 531
43. Rowlinson, A., O’Brien, P. T., Metzger, B. D., Tanvir, N. R., & Levan, A. J. 2013, MNRAS, 430, 1061
44. Sari, R., Piran, T., & Narayan, R. 1998, ApJ, 497, L17
45. Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525
46. Soderberg, A. M., Berger, E., Kasliwal, M., et al. 2006, ApJ, 650, 261
47. Tanga, M., Klose, S., Sudilovsky, V., Filgas, R., & Greiner, J. 2013, GRB Coordinates Network, 15214, 1
48. Tanvir, N. R., Levan, A. J., Fruchter, A. S., et al. 2013, Nature, 500, 547
49. Villasenor, J. S., Lamb, D. Q., Ricker, G. R., et al. 2005, Nature, 437, 855
50. Woosley, S. E. 1993, ApJ, 405, 273
51. Yost, S. A., Harrison, F. A., Sari, R., & Frail, D. A. 2003, ApJ, 597, 459
52. Zhang, B. 2013, ApJ, 763, L22
53. Zhang, B., Fan, Y. Z., Dyks, J., et al. 2006, ApJ, 642, 354
54. Zhang, B., & Mészáros, P. 2001, ApJ, 552, L35
55. Zhang, F.-W., Shao, L., Yan, J.-Z., & Wei, D.-M., 2012, ApJ, 750, 88
56. Zhang, S., Jin, Z.-P., & Wei, D.-M. 2015, ApJ, 798, 3
100685