EDP Sciences logo
Subscriber Authentication Point
Search
Menu
Home All issues Volume 576 (April 2015) A&A, 576 (2015) A71 Full HTML
Free Access
Issue
A&A
Volume 576, April 2015
Article Number A71
Number of page(s) 5
Section Astrophysical processes
DOI https://doi.org/10.1051/0004-6361/201424993
Published online 01 April 2015

© ESO, 2015

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 ~20 s. The other has a much shorter duration <2 s, centered at ~0.1 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 ~2.3 M⊙ (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 ~103 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).

Table 1

Optical observations of the field of GRB 130912A.

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 Sect. 2 we introduce the observations of GRB 130912A. In Sect. 3 we interpret the data. We summarize our results with some discussion in Sect. 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 T90 duration is 0.28 ± 0.03 second. The light curve shows two overlapping peaks. The time-averaged spectrum of the first 0.32 s is best fitted by a simple power law, with photon index Γγ = 1.20 ± 0.20. The total fluence is 1.7 ± 0.2 × 10-7 erg cm -2, and the peak photon flux is 2.2 ± 0.3 ph cm-2 s-1 (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 α = 1.20 ± 0.04. The spectrum formed from the PC mode data can be fitted with an absorbed power law with a photon spectral index of and cm-2, consistent with the Galactic value of NH = 1.2 × 1021 cm-2 (Kalberla et al. 2005).

The Swift UVOT took 150 seconds to find the chart exposure 98 s 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 r′ filter, the afterglow seemed to be constant, with a magnitude (here and throughout this paper, magnitudes are in AB system) of r′ = 22.0 ± 0.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 r′ 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 r 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 mag in r, i, Z, Y, J, and H 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 E(B − V) = 0.28 (Schlegel et al. 1998) and a ratio of total-to-selective extinction RV = 3.1, the Galactic extinction in r′ band is Ar′ = 0.78.

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 ~3.2 ×103 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.

thumbnail Fig. 1

X-ray (square) and optical (circle) light curves of GRB 130912A and the theoretical model curves. The X-ray data are from the UK Swift Science Data Center (Evans et al. 2009) and transformed to 1.732 keV. In order to reduce the influence of the error of the spectral index on the flux calculation, we used the geometric mean of the lower and upper boundaries of the corresponding X-ray energy band 0.3–10 keV. Proper corrections for extinction in the Milky Way Galaxy have been made. The solid and dashed curves are the theoretical optical and X-ray afterglow prediction with a forward shock.
Open with DEXTER

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) where Cp ≡ 13(p − 2)/ [ 3(p − 1) ], ϵe (ϵB) is the fraction of shock energy given to the electrons (magnetic field), td is the time in days, the Compton parameter , , and . Here and throughout this text, the convention Qx = Q/ 10x has been adopted.

At t ~ 3.2 × 103 s, if we have min { νc,νm } = νc ≳ 5 × 1014 Hz, the temporal behaviour of the optical emission would be , and the X-ray emission light curve should be in the case of νc<νx<νm or, alternatively, in the case of νc<νm<νx. While the temporal behaviour of optical emission agrees nicely with the data, the X-ray emission does not. The case of νc<νx<νm is clearly at odds with the data. The case of νc<νm<νx on the timescale of 100 s <t< 3.2 × 103 s is also inconsistent with the X-ray spectrum Fν ∝ ν− 0.50 ± 0.16 1.

If we, instead, have min { νc,νm } = νm ≳ 5 × 1014 Hz at t ~ 3.2 × 103 s, the temporal behaviour of the optical emission would be , and the X-ray emission light curve should be in the case of νm<νx<νc or, alternatively, in the case of νm<νc<νx. Now the temporal behaviours of both optical and X-ray emission are consistent with the data for p ~ 2.3, 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, νc(t ~ 1.0 × 104 s) ≈7.254 × 1016 Hz, νm(t ~ 3.2 × 103 s) ≈ 5 × 1014 Hz and Fν,max ~ 0.02 mJy are needed. Following Zhang et al. (2015), we have

where , , and .

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

In the case of Y ≤ 1 (i.e., the synchrotron-self Compton cooling is unimportant), (1 + Y)-2 can be ignored, and we have (see also Zhang et al. 2015) These three variables weakly depend on n0.

In the case of Y ≥ 1 (i.e., the synchrotron-self Compton radiation is important), we have (1 + Y)-2 ≈ Y-2 and then where . Compared with the case of Y ≤ 1, ϵB and Ek depend strongly on n0.

thumbnail Fig. 2

(Ek,ϵe,ϵB) as functions of n obtained in solving Eqs. (4)–(6). The solid, dotted, and dashed lines represent ϵB, ϵe, and Ek,53, respectively. The thin black line (i.e., Value = 0.4) is the reasonable upper limit of ϵB and ϵe, above which the solution is unphysical.
Open with DEXTER

The redshift z is unknown, and we assume that it is in the reasonable range of 0.1 ≤ z ≤ 1.4. However, for z< 0.4 in solving Eqs. (4)–(6), we do not find any reasonable values of the parameters with any given n0, while solutions are obtainable for larger z. In Fig. 2 we present the physical parameters (Ek,ϵe,ϵB) as functions of n0 for z = (0.4, 0.6, 0.8, 1.0, 1.2, 1.4), respectively. As already shown in the analytical approaches, both Ek and ϵB evolve with n0 quickly, while the dependence of ϵe on n0 is rather weak. One can also find from the figure that there is a very interesting constraint that n0 ≤ 3 × 10-3 cm-3 for a reasonable ϵe ≤ 0.4 (i.e., expected to be not much larger than the equipartition value ~1/3), 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 z = 0.72 (following Rowlinson et al. (2013) we adopt the average value of redshift of short GRBs), the fit parameters are (Ek,n0,ϵe,ϵB,p,θj) = (1.7 × 1051 erg, 0.002, 0.37, 0.16, 2.3, 0.03), where θj is the half opening angle of the GRB ejecta. An isotropic fireball is found to be unable to reproduce the data. The inferred ϵe and ϵB 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 Ar′ = 0.78, which is very high. If Ar′ is intrinsically smaller, Fνmax and hence a are lowered accordingly. As shown in Eqs. (7)–(9), or alternatively Eqs. (10)–(12), ϵB and ϵe would increase, while Ek would decrease. The contrary holds in the case of larger Ar′.

4. Summary and conclusions

The most remarkable feature of the short burst GRB 130912A is an optical plateau lasting about 4000 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 θj ~ 0.03 can reasonably reproduce the data. The circum-burst medium is found to be ISM-like and has a very low density ~10-3 cm-3, consistent with the model of merger of binary compact objects (either double neutron stars or a neutron star-black 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., ϵe ~ 0.37 and ϵB ~ 0.16, 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.

Acknowledgments

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.

References

All Tables

Table 1

Optical observations of the field of GRB 130912A.

In the text

All Figures

thumbnail Fig. 1

X-ray (square) and optical (circle) light curves of GRB 130912A and the theoretical model curves. The X-ray data are from the UK Swift Science Data Center (Evans et al. 2009) and transformed to 1.732 keV. In order to reduce the influence of the error of the spectral index on the flux calculation, we used the geometric mean of the lower and upper boundaries of the corresponding X-ray energy band 0.3–10 keV. Proper corrections for extinction in the Milky Way Galaxy have been made. The solid and dashed curves are the theoretical optical and X-ray afterglow prediction with a forward shock.
Open with DEXTER
In the text
thumbnail Fig. 2

(Ek,ϵe,ϵB) as functions of n obtained in solving Eqs. (4)–(6). The solid, dotted, and dashed lines represent ϵB, ϵe, and Ek,53, respectively. The thin black line (i.e., Value = 0.4) is the reasonable upper limit of ϵB and ϵe, above which the solution is unphysical.
Open with DEXTER
In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.