GRB 130831A central engine

# The central engine of GRB 130831A and the energy breakdown of a relativistic explosion.

M. De Pasquale , S. R. Oates, J. L. Racusin, D. A. Kann, B. Zhang, A. Pozanenko, A. A. Volnova, A. Trotter, N. Frank, A. Cucchiara, E. Troja, B. Sbarufatti, N. R. Butler, S. Schulze, Z. Cano, M. J. Page, A. J. Castro-Tirado, J. Gorosabel, A. Lien, O. Fox, O. Littlejohns, J. S. Bloom, J. X. Prochaska, J. A. de Diego, J. Gonzalez, M. G. Richer, C. Román-Zúñiga, A. M. Watson, N. Gehrels, H. Moseley, A. Kutyrev, S. Zane, V. Hoette, R. R. Russell, V. Rumyantsev, E. Klunko, O. Burkhonov, A. A. Breeveld, D. E. Reichart, J.B. Haislip
Mullard Space Science Laboratory, University College London, Dorking, United Kingdom E-mail: m.depasquale@ucl.ac.uk
Istituto Astrofisica Spaziale Fisica Cosmica, Palermo, Italy
Istituto Euro Mediterraneo di Scienza e Tecnologia, Palermo, Italy
Istituto de Astrofisica de Andalucia (CSIC), Granada, Spain
Center for Research and Exploration in Space Science And Technology (CRESST) and NASA Goddard Space Flight Center, Greenbelt
MD 20771, USA
Thüringer Landessternwarte Tautenburg, Sternwarte 5, 07778 Tautenburg, Germany
Department of Physics and Astronomy, University of Nevada Las Vegas, Las Vegas, USA
Space Research Institute (IKI), Moscow, Russia
University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, USA
Pennsylvania State University, University Park, PA, USA
Arizona State University, Tempe, Arizona, USA
Instituto de Astrofísica, Facultad de Física, Pontificia Universidad Católica de Chile, Vicuña Mackenna 4860, 7820436 Macul, Santiago, Chile
Millennium Institute of Astrophysics, Vicuña Mackenna 4860, 7820436 Macul, Santiago, Chile
Centre for Astrophysics and Cosmology, Science Institute, University of Iceland, 107 Reykjavik, Iceland.
Vasco UPV/EHU, Alameda de Urquijo s/n, E-48013 Bilbao, Spain
Ikerbasque, Basque Foundation for Science, Alameda de Urquijo 36-5, E-48008 Bilbao, Spain
Deceased
Department of Physics, University of Maryland, Baltimore County, Baltimore, MD 21250, USA
University of California Berkeley, USA
University of California Santa Cruz, USA
Instituto de Astronomía, Universidad Nacional Autónoma de México, Apartado Postal 70-264 04510 México, DF, Mexico
The University of Chicago, IL, USA
Crimean Astrophysical Observatory, 98409, pgt. Nauchny, Crimea
Institute of Solar-Terrestrial Physics, Russian Academy of Sciences, 664033, p/o box 291; Lermontov st., 126a, Irkutsk, Russia
Ulugh Beg Astronomical Institute, 100052, 33 Astronomicheskaya str., Tashkent, Uzbekistan
National Research Nuclear University MEPhI (Moscow Engineering Physics Institute), 115409 Moscow, Russia
###### Abstract

Gamma-ray bursts (GRBs) are the most luminous explosions in the universe, yet the nature and physical properties of their energy sources are far from understood.Very important clues, however, can be inferred by studying the afterglows of these events. We present optical and X-ray observations of GRB 130831A obtained by Swift, Chandra, Skynet, RATIR, Maidanak, ISON, NOT, LT and GTC. This burst shows a steep drop in the X-ray light-curve at s after the trigger, with a power-law decay index of . Such a rare behaviour cannot be explained by the standard forward shock (FS) model and indicates that the emission, up to the fast decay at  s, must be of “internal origin”, produced by a dissipation process within an ultrarelativistic outflow. We propose that the source of such an outflow, which must produce the X-ray flux for day in the cosmological rest frame, is a newly born magnetar or black hole. After the drop, the faint X-ray afterglow continues with a much shallower decay. The optical emission, on the other hand, shows no break across the X-ray steep decrease, and the late-time decays of both the X-ray and optical are consistent. Using both the X-ray and optical data, we show that the emission after s can be explained well by the FS model. We model our data to derive the kinetic energy of the ejecta and thus measure the efficiency of the central engine of a GRB with emission of internal origin visible for a long time. Furthermore, we break down the energy budget of this GRB into the prompt emission, the late internal dissipation, the kinetic energy of the relativistic ejecta, and compare it with the energy of the associated supernova, SN 2013fu.

###### keywords:
Gamma-Ray Burst, general – Gamma-Ray Burst, individual (GRB 130831A) - magnetar

## 1 Introduction

The study of gamma-ray bursts (GRBs) has received an exceptional boost thanks to the Swift mission (Gehrels et al. 2004), which has enabled rapid follow-up radio to X-ray observations of GRBs. However, despite a very large number of such fast follow-up observations performed by this spacecraft and by ground observatories, the characteristics of the “central engine” that produces the GRB are still unclear. The prevailing model (Woosley 1993; MacFadyen 1999, 2001; Thompson 2007) predicts the formation of a compact object, either a black hole or a magnetar, surrounded by an accretion disk. The compact object is the result of the core collapse of a very massive star or the final merger of two neutron stars (NS-NS) or NS - black hole binary. Under the correct conditions of angular momentum and magnetic field, such a system would launch a collimated, ultra-relativistic jet (e.g. via the Blanford-Znajek mechanism; Blandford & Znajek 1977). Within these jets, one or more “internal dissipation” processes take place (for a review, see Zhang 2011), converting part of the kinetic and/or magnetic energy into radiation. A burst of gamma-rays will then be visible if the observer is placed within the opening angle of the outflow.

Following the prompt gamma-ray emission, long-lived afterglow emission is detected in several energy bands (radio, optical, X-ray and probably high-energy gamma-rays). The consensus is that the afterglow radiation is emitted when ultra-relativistic ejecta interact with the circumburst medium, driving a forward shock (FS), which moves into the medium, and a reverse shock (RS), which propagates backwards through the ejecta. In particular, the emission due to the FS can in principle last indefinitely. A hallmark of RS and FS afterglow emission is that the flux density behaves as a power-law both in time and frequency, being described as , where is the time from the GRB trigger (Kobayashi & Zhang 2007) and is the frequency.
However, Swift observations have produced evidence of more complex phenomena during the afterglow phase, such as X-ray and optical flares (Falcone et al. 2006; Margutti et al 2011; Swenson et al. 2013), that cannot be attributed to the FS due to their fast temporal variability. For a few kiloseconds (ks) after the trigger, the afterglow flux often decays in a slow fashion (phase II of the canonical XRT light-curve, also known as “plateau”; O ’ Brien et al. 2006, Nousek et al. 2006; Racusin et al. 2009; Margutti et al. 2010) which cannot be understood if the fireball follows an adiabatic evolution; a process of energy injection (Zhang et al. 2006) into the fireball is usually invoked to explain such a feature. These observations strongly suggest that the GRB central engine is still active, and produces energy and/or a relativistic outflow.

In a small subset of Swift GRB afterglows, observations have shown slow decline phases of the X-ray flux which terminate with an abrupt fall in the emission, with slopes , sometimes approaching (Troja et al. 2007, Liang et al. 2007, Lyons et al. 2010, Rowlinson et al. 2013; Lü & Zhang 2014). Again, the FS model does not predict such behaviour. Instead, such steep decay can be expected when the central engine is a newly born magnetar, which emits a very high luminosity outflow due to the spin-down process (Usov 1992, Zhang & Meszaros 2001, Zhang & Meszaros 2002, Dall’Osso et al. 2011), which in turn produces emission we can directly observe in the X-ray. A very energetic outflow could also be produced by a newly formed stellar mass black hole surrounded by an accretion disk. The electromagnetic luminosity is expected to fall rapidly after the time-scale of the spin-down process, if the magnetar collapses into a black hole, or accretion onto the black hole stops.

The outflow produced by the spin-down process should generate emission via the synchrotron process. This kind of emission seems to be produced mostly in the X-ray band, while it is not usually detected at lower frequencies such as the optical (see Troja et al. 2007, Zhang et al. 2009, Rowlinson et al. 2013; however see Cano et al. 2014 for GRB130215A), perhaps because the optical band is at or below the synchrotron self-absorption frequency (Zhang 2009; see also Shen & Zhang 2009). At low frequencies, the dominant emission mechanism seems to be the standard FS, with its power-law decays and slow flux variations. FS emission, however, is still expected to be present even in the X-ray band, and it should emerge once the X-ray emission from the outflow produced by the magnetar or the black hole drops. So far, this has been seen clearly in the case of the long GRB 070110 (Troja et al. 2007). In this event, the X-ray light-curve showed a plateau lasting for  20 ks, after which the flux fell quickly with a slope of , and then resumed a power-law decay with a much shallower slope of . Such a late slope appears in most GRB afterglows and is likely produced by the FS.

In this paper, we present the well-sampled X-ray, UV/Optical and NIR observations of the afterglow of GRB 130831A, and show that its behaviour can be interpreted as a superposition of FS emission and “internal emission”, the latter of which suddenly ceases at  ks. This article is organized as follows. In Sect. 2, we present the observations taken by different instruments and observing facilities, and show how the data were reduced and analyzed. We recall the results of the observations taken by other groups on GRB130831A as well. In Sect. 3, we show and fit the resulting light-curves and spectral energy distributions of this GRB. In Sect. 4, we model the afterglow of GRB 130831A in the context of FS and internal emission models, we discuss the possible origins of the observed emission and the properties of the object that produced the explosion. Finally, we present our conclusions in Sect. 5. Throughout the paper, errors are expressed at the confidence level (CL) unless stated otherwise, and we assume a CDM Cosmology with km Mpc, and (Jarosik et al. 2011).

## 2 Observations and Analysis

### 2.1 X-ray data

#### 2.1.1 Swift-Bat

The Swift Burst Alert Telescope (BAT; Barthelmy et al. 2005) triggered on GRB 130831A at 13:04:16.54 UT, 2013 August 31. The mask-weighted light-curve in the 15-350 keV energy range (see Fig. 1) shows a main pulse with a fast-rise-and-exponential-decay (FRED) shape. It starts at T-2 s and peaks at T+3 s. The major pulse structure ends at T+12 s. However, there is some extended emission that lasts until T+41 s, with two additional peaks at T+20 s and T+32 s. T (15-350 keV), the interval during which from 5% to 95% of the total emission is recorded, is  s (error includes systematics).

The time-averaged spectrum between T1.9 s to T+41.4 s can be fitted by a simple power-law model with a spectral index ( for 57 degrees of freedom, ) in the energy band 15-150 keV. This model gives a fluence of erg cm in the same energy range. The 1-s peak spectrum can also be fitted by a simple power-law model with () in the energy band 15-150 keV, which gives a peak flux of erg cm s in the same band. The BAT analysis uses the event data from -240 s to +963 s.

#### 2.1.2 Swift-XRT Observations

After the BAT trigger, the Swift satellite promptly slewed to point the narrow field instruments at the source. The X-ray Telescope (XRT; Burrows et al. 2005) began observations of GRB 130831A 125.8 s after the trigger in the 0.3–10 keV energy band. The XRT monitored the source until 2013 September 14, collecting 84.9 s of data in Window Timing (WT) mode (11.6 s were taken while the spacecraft was settling and the remaining 73.3 s pointing mode) and 59.7 ks in Photon Counting (PC) mode.

XRT data were processed using the HEASoft software packages, version 6.15 and version 20130313 of the XRT Calibration DataBase222www.swift.ac.uk/analysis/xrt/files/SWIFT-XRT-CALDB-09_v17.pdf, applying calibrations and standard filtering and screening criteria (for more details, see Evans et al. 2009). WT data were extracted in an interval centred on the source, 20 pixels on each side, and the background estimated using intervals between 40 and 60 pixels from the source. The PC data were initially affected by pile-up, and corrected by excluding the central core of 6 pixel radius. The remaining PC data were extracted using a circular region with 30-pixel radius, and when the count rate dropped below  counts s, within a 10 pixel radius. During the extraction of the light-curve, a minimum signal-to-noise ratio (S/N) criterion was applied for the re-binning, such that we required a S/N greater than 3 for all data points, with the exception of the last data point. The 90% confidence level intervals were estimated following the Kraft et al. (1991) technique for low count rate sources. The X-ray light-curve (Fig. 2), which shows the temporal evolution of the flux, was constructed once the spectral information was obtained (see Sect.s 3.1 and 3.2, Tables 1 and 2).

#### 2.1.3 Chandra observations

After the very steep break, we requested and obtained two Chandra Director’s Discretionary Time (DDT) observations (PI: De Pasquale). Meanwhile, Swift continued observing for several more days detecting a dim X-ray afterglow, which suggested that the steep decay had broken to a gentler slope. The Chandra observations of GRB 130831A took place on 2013 September 17 (T+16.6 days, or 1430 ks) and 2013 October 3 (T+33.1 days, or 2860 ks) with an exposure of 15 ks each. The location of GRB 130831A was imaged on the ACIS S3 chip in both observations. The ACIS data were processed using ciao version 4.6 and version 4.5.9 of the Calibration Database (Fruscione et al. 2006) using standard tools. We filtered the events for grades 0, 2, 3, 4, and 6 and energy from 0.5–8 keV. We extracted the source counts within a radius region centered on the best known position of the GRB, yielding 8 and 1 counts in the two observations, respectively. The first epoch yielded a detection with a significance of determined using the method described in Kraft et al. (1991). The second epoch was only above background, thus it should not be regarded as a real detection.

### 2.2 Optical observations

Here we describe how we collected, reduced and calibrated the ultraviolet, optical and infrared photometry of GRB 130831A. All magnitudes have been calibrated to the Vega system. The full set of photometric measurements is provided in an on-line table; a sample is shown in Table 3.

#### 2.2.1 Uvot

Swift/UVOT (Roming et al. 2005) began observing the field of GRB 130831A 114 s after the trigger (Hagen et al., GCN Circular 15139) and started settled observations 191 s after the trigger, with a finding chart exposure in the band. The afterglow was detected in all 7 UVOT filters. Observations were taken in both imaging and event modes. Before extracting count rates from the event lists, the astrometry was refined following the method described in Oates et al. (2009). The source counts were extracted initially using a source region of 5” radius. When the count rate dropped to below 0.5 counts per second we then used a source region of 3” radius. In order to be consistent with the UVOT calibration, these count rates were then corrected to 5” using the curve of growth contained in the calibration files. Background counts were extracted using a circular region of radius 20” from a blank area of sky situated near to the source position. The count rates were obtained from the event and image lists using the Swift tools uvotevtlc and uvotsource, respectively. They were converted to magnitudes using the UVOT photometric zero points (Breeveld et al., 2011). The analysis pipeline used version 20130118 of the UVOT Calibration Database.

#### 2.2.2 Ratir

The Reionization And Transients Infra-Red camera (RATIR; Butler et al. 2012) observed GRB 130831A over a period of seven hours, beginning 15.8 hours after the Swift trigger, with follow-up observations on six nights over the next month. RATIR is mounted on the 1.5-meter Harold L. Johnson telescope of the Observatorio Astronï¿½omico Nacional on Sierra San Pedro Mï¿½artir in Baja California (Mexico). This facility, which became fully operational in December 2012, conducts autonomous observations (Watson et al. 2012; Klein et al. 2012) of its targets in the six photometric bands , , , , and simultaneously. RATIR captured 80 s exposure frames in and 67 s exposure frames in due to additional overhead. We applied their standard image reduction pipeline with twilight flat division and bias subtraction routines written in Python and using astrometry.net (Lang et al. 2010) for image alignment and swarp (Bertin 2010) for image co-addition. Aperture photometry was calculated using sextractor (Bertin & Arnouts 1996).

#### 2.2.3 NOT and LT

The 2.5 m Nordic Optical Telescope (NOT) and the 2 m Liverpool Telescope (LT), located at Roque de los Muchachos Observatory, La Palma in the Canary Islands (Spain), obtained three epochs of photometry presented in this paper: 2013 September ( and ), 2013 September 3 (), 2014 January 5 (). Several additional epochs of NOT and LT photometry were obtained as part of a campaign to observe SN 2013fu, which are presented in Cano et al. (2014). Image reduction of the and photometric data was performed using standard techniques in iraf333iraf is distributed by the National Optical Astronomy Observatory, which is operated by the Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation.: bias combine, bias-subtract, flat-field co-add, flat-field normalise, flat-field divide, align and co-add (Tody 1986, 1993). The optical data were calibrated using Sloan Digital Sky Survey (SDSS; Ahn et al. 2012) stars in the GRB field, and converted into Vega magnitudes.

#### 2.2.4 Skynet

Skynet (Reichart et al. 2005) obtained images of the field of GRB 130831A on 2013 August 31 with four 17” telescopes of the Panchromatic Robotic Optical Monitoring and Polarimetry Telescope (PROMPT) array at Siding Spring Observatory, New South Wales, Australia. Further observations were taken on September 1, this time observing also with the four 16” telescopes of the PROMPT array at the Cerro Tololo Inter-American Observatory, Chile and the 41” telescope at Yerkes Observatory, Wisconsin, USA. Beginning at 13:39 UT (T-T min), exposures ranging from 60 to 180 s were obtained in the (PROMPT), and (Yerkes) bands. Bias subtraction and flat-fielding were performed by Skynet’s automated pipeline. Post-processing occured in Skynet’s guided analysis pipeline, using both custom algorithms and ones based on iraf. Differential aperture photometry was performed on single and stacked images, with effective exposure times from 60 s to 2.7 h. Photometry was calibrated against the catalogued magnitudes of six APASS DR7 stars in the field (Henden & Munari 2014). The calibrations for magnitudes were derived using transformations obtained from prior observations of Landoldt stars (A. Henden, private communication).

#### 2.2.5 IKI Network for Transients

The field of GRB 130831A was observed by several facilities of the follow-up network organized by the Space Research Institute (IKI) of Moscow, where observations of GRBs are planned and data reduction is carried out. We detail the observations of GRB 130831A of the individual observatories below.

The International Scientific Optical-Observation Network (ISON; Molotov et al. 2008, Pozanenko et al. 2013c) started observing on 2013 August 31 at 13:14:32 UT, i.e. minutes after the trigger, with the 0.65-m telescope SANTEL-650 (Volnova et al., 2013a) of ISON-Ussuriysk observatory. Thirty unfiltered images, each with 120 seconds exposure, were taken. The 50-cm telescope VT-50 of ISON-Ussuriysk observatory started to observe at 13:26:10 UT, 22 min after the trigger, taking 384 unfiltered images with exposures of 30 seconds in two epochs with a gap of about 1.8 h (Volnova et al., 2013a). Starting at 19:12:30 UT during the first day, the 40-cm SANkovich-TELescope (SANTEL) -400AN of ISON-Kislovodsk observatory took 34 frames with exposures between 60 and 120 s (Volnova et al., 2013b).

The Astronomicheskii Zerkalnyi Telescope - 8 (AZT-8) 0.7-m telescope of Gissar (Tajikistan) observatory took 57 frames in band with 60 s exposure time each starting at 17:47:54 UT.

The AZT-8 0.7-m telescope of the Chuguev Observational Station (Institute of Astronomy, Kharkiv National University) observed the afterglow in filter starting on September 1 at 19:37:49 UT (Volnova et al. 2013c).

The optical afterglow was also observed by the 1.5-m telescope AZT-22 of Maidanak observatory (Uzbekistan) on 2013 September 1, 2, 4, 5, 7–11, 15, 16, 22, 27, and 29. Every observational night, frames were taken in the filter with an exposure time of 600 s each. On September 1 we also took several frames in the , and bands with the same exposure, and on September 2 we obtained additional observations in the filter.

Most of the data obtained by Maidanak are presented in Cano et al. (2014) and used to study SN2013fu; in this article we use the early data which were not contaminated by the SN emission (see Sect. 3.3).

Observations were also taken with the 1.6-m telescope AZT-33IK of Sayan observatory (Mondy, Russia) on September 3, 4, and October 10, 11 and 14, taking several frames in the -band with an exposure of 60 seconds each. The optical afterglow was also imaged with the 2.6-m Shajn telescope of the Crimean Astrophysical Observatory on September 5, taking several frames in the band with an exposure of 120 seconds each.

All data obtained by the facilities indicated above were processed with the same initial reduction including dark frame subtraction and flat-fielding, and using the iraf packages apphot and daophot. For the photometric calibrations we used four stars from the SDSS, indicated in Table 4. The magnitudes of the reference stars were transformed to the photometric system using the transformation equations attributed to Robert Lupton in the SDSS online documentation555www.sdss.org/DR7/algorithms/sdssUBVRITransform.html/
Lupton2005
.

#### 2.2.6 Gtc

The 10.4m Gran Telescopio Canarias (GTC; (Canary Island, Spain), equipped with the Optical System for Imaging and Low-Intermediate Resolution Imaging Spectroscopy (OSIRIS) instrument (Cepa et al. 2000), performed deep -band imaging of the GRB field 13330 ks after the trigger (2014 February 1), in order to obtain photometry for the host galaxy. Nine 60 s images were acquired in binning, providing a pixel scale of 0.25 arcsec pixel. The images were dark-subtracted and flat-fielded using custom iraf routines. Aperture photometry was performed using daophot tasks as implemented in iraf. Photometric calibration was based on SDSS standard stars present in the OSIRIS un-vignetted field of view (7.8 arcmin 7.8 arcmin).

#### 2.2.7 Results from other facilities

In this subsection, we describe results obtained by other teams, not involved in our work. Nonetheless, we will use their results in the next sections.

Konus-Wind Observations. GRB 130831A was observed by Konus-Wind onboard the WIND spacecraft. Golenetskii et al. (2013) found that this burst had a duration of 35 seconds; it was detected up to MeV. In addition, the spectrum is relatively soft; it can be fitted between 20 keV and 15 MeV with the Band model (Band et al. 1993), yielding a low-energy spectral index , a high-energy spectral index and peak energy ). The fluence between 20 keV-10 MeV is  erg cm (90% CL). We use the results of the Konus-Wind data analysis presented in Golenetskii et al. (2013), which spans a much wider spectral range than the BAT data, to derive the energetics of this burst. Thus, they are better suited to assess the energy emitted by the burst without using extrapolation. Using a cosmological k-correction (Bloom et al. 2001) and the redshift of this burst (Cucchiara et al. 2013), we derive 1–10000 keV rest-frame energetics of erg. This burst follows the Amati relation (Amati et al. 2006; Amati et al. 2009; see also Cano et al. 2014, their Fig. 12).

Radio. GRB 130831A field was also observed by the Jansky Very Large Array (EVLA) and Combined Array for Research in Millimeter Astronomy (CARMA), Zauderer et al. 2013. According to the analysis reported by Laskar et al. (2013), no significant radio emission was detected: the upper limits are 38 Jy and 71 Jy at 5.8 and 21.8 GHz, respectively, at 0.64 days after the trigger. Similarly, CARMA obtained a 3 upper limit of 780 Jy 0.76 days after the trigger (Zauderer et al. 2013).

#### 2.2.8 Building the UVOIR light-curves

In the subsections above we summarized the observations of the optical instruments and how the data produced by each one were reduced. We shall now describe how we combined these different datasets into a homogeneous set of flux light-curves. Magnitudes were translated to fluxes using the zero-magnitude flux densities listed in Table 5. Our overall approach has been that when data have been obtained from multiple observatories in the same (or almost the same) band, we scale all the data to match the dataset which provided the largest number of measurements in that band. The largest set of measurements in and come from RATIR. We normalized the NOT and LT and data points to RATIR using the calibration stars common to the two datasets (see Table 4). The comparison of the calibration stars led us to scale the flux from the NOT band data point by a factor of 0.99 and, based on the scatter between the NOT and RATIR measurements of the calibration stars we added a systematic error of 2 per cent to the errors on the NOT fluxes. The NOT and LT band flux was scaled by a factor of 1.01 and a systematic error of 1 per cent was added. The Skynet Yerkes and data were matched to RATIR using the common reference stars listed in Table 4. The and fluxes from Skynet were scaled by factors of 0.94 and 0.91 respectively and systematics of 3 per cent () and 2 per cent () were added to the uncertainties on the Skynet fluxes.

The largest set of measurements in , , and come from the Skynet PROMPT observations. To incorporate the UVOT and data into the and light-curves, we first transformed the UVOT magnitudes of the GRB afterglow and the Skynet calibration stars into Johnson magnitudes using the appropriate colour transformations in Poole et al. (2008), before translating to flux densities. The UVOT and Skynet photometry of the Skynet calibration stars (see Table 4) were then compared to determine the appropriate scaling factors and systematic errors for UVOT. The UVOT V and B fluxes were scaled by factors of 1.13 and 1.07 respectively, and systematic errors of 9 per cent and 6 per cent were added respectively to their flux uncertainties.

Based on the photometry of the calibration stars in common between the Maidanak and Skynet PROMPT observations (see Table 4) in the , and bands, we scaled the Maidanak , and fluxes by factors of 1.01, 1.01 and 1.06 respectively and added systematic errors of 10, 5 and 7 per cent. For the -band data from the IKI Gissar, Maidanak, Chuguev and Sayan observatories, the calibration stars in common with Skynet suggested a scaling factor of 1.01 for these data and a 6 per cent systematic error, but the resulting GRB -band light-curve from these observatories appeared to be systematically lower than that obtained from Skynet. A power-law fit to the -band flux light-curve derived from IKI observations between 15 and 100 ks, performed simultaneously with a power law fit to the Skynet band data in the same time interval, and with the power-law slope tied between the two datasets, gives a best-fitting normalisation for the IKI data which is 0.79 times that of the Skynet data. Therefore the IKI fluxes were scaled by this factor before combining them with the Skynet PROMPT data.

## 3 Results

### 3.1 Spectral analysis of X-ray data

The XRT spectral data (0.3–10 keV) were extracted in nine different temporal segments, with boundaries selected according to the different observing modes and breaks in the light-curve (see Table 2). Ancillary response files and exposure maps were created using the heasoft software for each segment, as well as appropriate response matrices from calibration database. All the spectra were fitted using xspec v12.8.1 (Arnaud et al. 1996) with an absorbed power-law model. Some indication of spectral evolution was found, with the spectral index trending from soft, during the first orbit, to hard from the second orbit onwards. Detailed results are shown in Table 1.

The X-ray spectrum between 9 ks and 132 ks, i.e. the end of the unusual late steep decay, can be modeled with a power-law with Galactic absorption and intrinsic absorption at the redshift of the burst, . The Galactic absorption has been fixed to cm (Kalberla et al. 2005). The best fit parameters are: spectral index , intrinsic cm which is consistent with 0 at CL.

To infer the late X-ray flux from Chandra measurements, we applied the following procedure. First, we corrected for the portion of the PSF excluded from the 1.5 arcsec radius aperture. Using the PC mode XRT spectral fit parameters (spectral index , intrinsic cm), we derived 0.3–10 keV fluxes for the two epochs of erg cm s and erg cm s, respectively. To derive the latter, we used the estimates on confidence limits for small numbers of events in astrophysical data (Gehrels 1986). The upper limit corresponding to the second Chandra observation is erg cm s. Both XRT and Chandra fluxes are corrected for absorption. The results of the spectral analysis were used to compute the rate to flux conversion factors employed to build the flux light-curve.

### 3.2 The X-ray light-curve of GRB 130831A.

The X-ray light-curve is shown in Fig. 2. After an initial fast decay with slope ending at  s, there is an X-ray flare starting at about 500 s after the trigger and lasting at least up to the end of the first orbit at 900 s. The highest flux recorded was about erg cm s at the beginning of XRT observations. When GRB observations resumed,  ks after the trigger, the decay slope was slower than that at the very beginning of observations. This decay terminated at about ks after the trigger, when the flux showed a surprisingly steep drop.

We fit the X-ray data (see Table 2) from the very beginning with the sum of an early power-law, a broken power-law decay, a Gaussian flare superimposed on the second segment, and another, final power-law. The first power-law basically represents the tail of the prompt emission, probably due the curvature effect, which then gives way to a shallower decay usually attributed to a different emission mechanism (Zhang et al. 2006). When this second process ends as well, the flux falls quickly again. The late power-law may represent emission powered by the FS mechanism, which can emerge once the emission from the previous process is over. The fit is acceptable, yielding . We find a decay slope between and  ks. Such a slope is intermediate between the typical shallow decay phase and the “normal” decay phase seen in a wide sample of GRB afterglows (see Evans et al. 2009). We can still find a value of , however, in the distributions of decay indices of both phases. After the break at  ks, the best fitting temporal slope is , much faster than is usually observed in late X-ray afterglows, even in the case of a jet break, when (Sari et al. 1999; Racusin et al. 2009). This steep drop suggests that the X-ray afterglow might have been, until then, produced by some internal dissipation mechanism rather than the typical FS emission. With this fit model, the late power-law component has a decay slope . We note that this decay slope of the X-ray flux after ks is steeper than the decay slope during the shallow decay phase. If the late X-ray flux were FS emission, it might have begun hundreds or even thousands of seconds after the trigger. However, the fitted model above does not include this possibility, and it may yield late power-law component slopes flatter than the real ones in order not to over-predict the flux at very early epochs. To better investigate the late emission, we use a different time interval. If we fit the X-ray data points from 100 ks onwards with a simple power-law model, we obtain a poor fit, . The best fit decay slope is . However, if we fit the same data points with a power-law + power-law model, we obtain a much better fit, with . This model yields a decay slope for the first power-law of , with a lower limit (with ) of . Such a value, though, is still too steep for the FS model, see above. The best fit value for the late power-law component is , while its flux at 2 days (173 ks) after the trigger is   erg cm s.

An analysis with the F-test suggests that the two-power-law model is not necessary, because the probability of an improvement by chance is , which is not negligible. However, if we adopted the simple power-law model with the steep decay slope after 100 ks, then the flux at the time of the first Chandra observation would be erg cm s. Such extremely low flux corresponds to less than counts with a ACIS-S 15 ks observation; thus our first Chandra observation would most likely yield 0 counts. We can place a 99.5% CL upper limit of 5.3 counts, according to Gehrels (1986). This prediction, however, is in disagreement with the fact that the observation actually produced 8 counts. According to Kraft et al. (1991), these 8 counts represent a detection. We therefore conclude that the X-ray afterglow decay has become much shallower at late epochs, and we will adopt the results of the fit of the two-power-law model after 100 ks.

### 3.3 Optical light-curve of GRB 130831A

The combined optical light-curves in Fig. 2 show an initial short plateau, which lasts until s, followed by a steep rise and a peak at  s. This optical flare is basically concurrent with the X-ray flare. Following the flare, there is another plateau that in turn gives way to a steeper decay at  s. In this phase, fitting the optical light-curves with simple models such as a broken power-law does not provide a statistically acceptable fit. For example, the best fit of band data between 3.5 and 15 ks with a broken power-law yields an early decay slope of , break time  ks, and post-break slope with . Fitting the other light-curves in this interval yields similar results. We nonetheless plot the best-fit curves in Fig. 2, as an indication for the behaviour of the optical afterglow, and report the results of fitting the light-curves in the interval between 3.5 and 15 ks in Table 6. Such high are due to some “wiggles” of the densely sampled light-curves in this phase, which has been seen in other GRBs (e.g. Swenson et al. 2013; Matheson et al. 2003). The post-plateau decay is moderately steep and does not seem to change its slope at the epoch of the X-ray drop. However, about days after the trigger, the optical emission starts to rise due to light coming from SN 2013fu, the supernova associated with GRB 130831A (Klose et al. 2013; Cano et al. 2014). Given the complication of considering the SN flux, we have excluded all data points that had a contribution from the SN, basically those after  ks and before  ks. The SN contribution at these epochs have been estimated using the SN 1998bw template program presented in Cano (2013).
In addition, we have observations at days in and filters, taken with Gran Telescopio Canarias (GTC) and LT (see Sect. 2.2.3 and 2.2.6). These late time data do not suffer significant contamination from the SN and the afterglow, which have faded away. They correspond to the magnitude of the host galaxy, and have been used to determine the optical afterglow behaviour (see below). Vega magnitudes of the host galaxy are , . Assuming the conversions from Jordi et al. 2006 (their Table 1), we find a Vega magnitude . Spectral observation taken at the Gemini North Observatory revealed that GRB 130831A occurred at redshift of (Cucchiara et al. 2013). For such a redshift, the magnitude of the host corresponds to a luminosity in the B band (Hjorth et al. 2012).

We know the host contribution in , and only, from our late-time GTC and LT images. We fitted the light-curves in these filters with the same model, namely a power-law plus constant, from 15 ks up to 100 days after the trigger. We ignored the data before 15 ks because they would lead to a very bad fit, as we previously noted (see above). The best fit decay slopes and are , ; , ; and , (see Table 6). These fits are statistically acceptable and consistent within . We find that the weighted mean is . The decay slopes of the flux in other filters are consistent with this value within as well. We note that Cano et al. (2014), in their analysis of the afterglow light-curves, find that the optical decay slope is , consistent with our analysis. We remark upon the fact that we can fit the optical light-curves with an uninterrupted power-law, even across the X-ray break. This feature strongly suggests that optical and X-ray emission (at least part of it) have different origins.

### 3.4 Spectral Energy Distributions

To test the hypothesis that the late emission is entirely due to FS, we built a spectral energy distribution (SED) with the available UVOIR + X-ray data at 2 days after the trigger (173 ks).

First, we calculated the count rates at 2 days. We used the data between 15 ks and ks, since no colour evolution was detected in the UV to the near IR, and the count rate in all light-curves could be fitted as a power-law with a common decay index between these two epochs. The UVOT data were translated to Xspec-compatible files using the standard ftool uvot2pha. Then, we adjusted the count rates of these files to the values determined by fitting the light-curves. Each of the ground based optical and near IR photometric data points were imported into Xspec using bespoke software, as follows. Each photometric data point was recorded as a single-channel spectral file containing a count rate and count-rate uncertainty, with a corresponding response file. To produce the response files, the responsivity of the filter/telescope combination as a function of wavelength was converted to a normalised effective area as a function of energy. As for the X-ray, we first determined the light-curve count rate at 2 days, . We then determined a new exposure time for which , where is the count rate of the source after background subtraction. We then imported the source and background XRT spectral files with the changed exposure times into Xspec. To build the SED, we used the XRT data after the steep decay slope, but did not use the Chandra data because XRT and Chandra fluxes would have to be renormalized to the same value and XRT has more counts.

We fitted the optical and X-ray data with Xspec (Arnaud et al. 1996). We adopted a simple power-law model with two absorbers and two zdust components, one at and another one at , i.e. the redshift of the burst (Cucchiara et al. 2013). The values of the absorbers are the same given in the X-ray data analysis section 3.1. The Galactic reddening was fixed at mag according to the map of Schlegel et al. (1998). As for the extinction at the redshift of the burst, we tried the Milky Way, Large and Small Magelanic Clouds (MW, LMC and SMC) extinction laws as in Pei (1992). We found that all of them yield acceptable and similar fiting results. However, the MW extinction law provides the best fit, so we have adopted the results of the fit with this law. The SED of GRB 130831A at 2 days is shown in Fig. 3. The best fit parameters are an index and a small or absent amount of extragalactic reddening mag (Cano et al. 2014 find a negligible rest-frame extinction as well); this fit yields . A broken power-law model does not significantly improve the fit, since it yields and the break energy is unconstrained. Moreover, fits with LMC and SMC extinction law result in a break energy above the X-ray band. All in all, we believe that a simple power-law model is adequate to describe the SED at this late epoch. Results of the fit of the 173 ks SED are shown in Table 7. There are only counts collected by XRT after the fast drop, and to build the SED and obtain the quoted results we constructed a single X-ray data bin spanning from 0.3 to 10 keV. We were concerned that the use of such a wide bin might be not the optimum in the SED fitting process. Therefore, we repeated the fit using standard statistics with optical data and Cash statistics for the X-ray data, where the bins were constituted of single counts. We obtained very similar results. The fit with a simple power-law model yielded , mag with total statistics of and 26 dof. A fit with a broken power-law model was marginally better, yielding statistics of 22.2 with 25 dof, but the F-test indicates that the probability of an improvement by chance was high, with a probability of .

We also show that the same model does not apply to earlier data. Following the same procedure outlined above, we built a SED with UVOIR and X-ray data at 80 ks (Fig. 3), before the steep X-ray drop. Then, we changed the normalization of the 173 ks power-law fit, following a decay slope of . We then plotted such a re-normalized power-law model with onto the 80 ks SED. We see that, while the optical emission is easily matched, the observed X-ray flux lies well above the model prediction. This finding is confirmed by the fact that if we fit the 80 ks SED with the same power-law model used for the 173 ks SED, we obtain a best-fit spectral slope . This is harder than the slope found at 2 days and inconsistent with it. This result confirms that the spectrum at 80 ks across the X-ray and optical bands is not consistent with the spectrum at 173 ks and lends credence to the idea that there is an additional component in the X-ray band at early epochs.

## 4 Discussion

### 4.1 Modeling of GRB 130831A and the Efficiency of its “Central Engine”.

We found that the X-ray light-curve, after the steep drop at 100 ks, resumes a slower decay slope of ; such a decay slope is statistically consistent with the optical decay slope from 15 ks onwards, , at CL. The SED at 2 days (173 ks) after the trigger is adequately fit with a simple power-law model, in which the X-ray and the optical band lie on the same spectral segment. Values of the spectral and temporal indices seem to be typical of GRB late afterglow emission, and are easily explainable by the FS model.
The FS model predicts basic relations between spectral and decay indices (for a review, see Zhang et al. 2006), which depend on the type of expansion (spherical or jet) and the spectral regime, i.e. where the observing bands are located relative to the synchrotron peak frequency and cooling frequency, and respectively, and the density profile of the circumburst medium. We can use these closure relations to find the conditions that apply to the case at hand. We adopted as decay slope, since it is better constrained than the late X-ray decay index, and as the spectral slope. If the observing bands (both the X-ray and optical) were below and the medium had a constant density profile, we would have ; this is in agreement with observations at . If the medium had a stellar wind profile, with density decreasing with radius as , then ; this is ruled out at level. If the observing frequency is above the cooling frequency, for both the constant medium and stellar wind profile cases the relation should be satisfied; but this is rejected at . Finally, if the outflow is collimated and has decelerated enough so that the observer detects emission from the edges, the decay should become steeper (“jet break”). After the jet break, the relations to satisfy are and if the observing frequency is below or above . These relations are ruled out at and respectively. The only relation of those above fulfilled within is that for the observing frequency below , constant density medium, and pre-jet break expansion.

In the context of the FS model, the flux produced by the afterglow depends on several factors, such as the kinetic energy of the outflow, the fraction of energy given to radiating electrons and to the magnetic field , the density of the environment , the slope of the power-law energy distribution of electrons , the type of expansion and where the observing bands are located. It is possible to derive the kinetic energy through simple relations once the flux is known, but one has to make assumptions of the values of the other parameters.

Following Zhang et al. (2007), the flux density of the FS emission is given by

 Fν(νm<ν<νc)=1600D−228(1+z)3+p4ϵp+14B,−2ϵp−1e,−1× E3+p4K,52n1/2t34(p−1)dfp−1p(ν3.3×1012ν)p−12 μJy (1)

Where is a parameter depending on , and is the luminosity distance in cm. Subindices indicate normalized quantities, in cgs units. A spectral index for an observing frequency between and implies . Such a large value of is not very common; typically . However, is still within two standard deviations from the centre of the distribution of this parameter, as implied from the analysis of more than Swift GRBs (Curran et al. 2010), and we will adopt this value. As for the observing frequency, we take  Hz, i.e. 1 keV. We convert the X-ray flux between 0.3 and 10 keV into a flux density of  Jy at this frequency.

From the equation above, we obtain

 ϵ1.02B,−2ϵ2.06e,−1E1.52K,52n1/2≃2.0 (2)

We infer if we assume , , and . These low values of and are required to have well above the X-ray band 2 days after the trigger. A lower value of may imply significant Inverse Compton flux, which is usually not detected in GRBs. The low density is not totally unprecedented in GRB afterglow modeling, since it has already been found for events in constant density media (Panaitescu & Kumar 2002, Cenko et al. 2011; with radio observations) and stellar wind media (Cenko et al. 2011, Perley et al 2014). On the other hand, a density lower than is not usually expected in long GRBs, which occur next to, or within, dense star forming regions. We also note that has a weak dependance on , so our inferred would not be dramatically different if the value of this parameter were within one order of magnitude from what we use. As for , the value we have chosen is close to that of equipartion, and is obtained by modeling GRB afterglows (see aforementioned references). However, we now show that in our modeling we cannot have . From Eq. 2, we have roughly. FS theory predicts that and , where is the peak synchrotron flux reached at . For , ; the flux below will therefore be . For the values of the parameters quoted above, the radio flux at 22 and 5.8 GHz predicted 0.67 days after the trigger would be 56 Jy and 36 Jy respectively, below the upper limits determined by the VLA (Laskar et al. 2013). However, if were less than the predicted radio flux would instead be larger than the above limits. We note, though, that interstellar scintillation could suppress the observed radio flux as well. For the values of parameters we’ve constrained, the FS peak flux will be  Jy. Such a flux is quite typical for GRB FS peak; see Chandra et al. (2012) and De Ugarte-Postigo (2012). We can assume that our estimate on is robust at least to an order of magnitude.

The modeling above enables us to determine the efficiency of the conversion of the energy of the outflow into energy emitted in high energy photons during the prompt emission. The efficiency is defined as . Given  erg and the value of constrained above, we find . By means of both optical and X-ray data, we have thus measured the efficiency of the central engine of a GRB with a prolonged internal emission episode. The optical light-curves were well sampled, and late Chandra observations constrained the X-ray flux light-curve after the steep drop. Thus, we could establish with little doubt that the late emission was entirely consistent with the FS model. Zhang et al (2007) and, more recently, Lü & Zhang (2014) carried out a similar task by using XRT light-curves; such a study with X-ray data alone, especially without precise late measurements, may be more ambiguous.

An efficiency of is lower than that of those few GRBs that unambiguously show a plateau of “internal origin” (see figure 12 of Lü & Zhang 2014). The efficiency of GRB 130831A is more characteristic of those GRBs that have their afterglow emission entirely explained by the FS mechanism, but with the presence of energy injection into the ejecta perhaps powered by a magnetar. Typically, GRBs with a plateau of internal origin have , while those explained by FS have a range of with most clustering around a few  .

### 4.2 Origins of the X-ray radiation between the end of the prompt emission and up to the 100 ks drop.

#### 4.2.1 Observations

The sudden drop in the X-ray flux at ks, with a decay index , cannot be interpreted as FS emission. The steepest decay in this model is , where is the index of the power-law energy distribution of radiating electrons, which occurs during a jet-break expansion phase. However, is not predicted at all on theoretical grounds (e.g. see Rieger et al. 2007); thus it is very difficult to explain such a steep decay index at late times. Theoretically, an index as steep as can be achieved by taking into account relativistic effects in simulations (e.g. Granot 2006, van Eerten et al. 2010) even if , but the value of the decay index reached in the case of GRB 130831A is greater than this prediction. Duffell &MacFadyen (2014) explored the possibility that the plateaux we see in GRB afterglows are produced by a jetted outflow before deceleration, followed by a steeper decay, which flags the Blandford & McKee deceleration of the ejecta. However, such a decay in their model does not reach a value as steep as detected in GRB 130831A between 100 and 200 ks.

By assuming that the steep decay is due to the curvature effect (Kumar & Panaitescu 2000 ; see also Uhm & Zhang 2014), we followed Liang et al. (2006) to test the internal origin of X-ray emission up to the steep break. We found that to satisfy the condition of the curvature effect, the zero time of this emission needs to be ks after the trigger, slightly before the beginning of the flux drop. The very steep decay component therefore strongly suggests that the X-ray emission up to 100 ks is of internal origin, since is allowed to occur at large time intervals from the initial prompt emission in this case (Zhang et al. 2006; Liang et al. 2006).

In a more general study of the curvature effect emission (Uhm & Zhang 2014), the emitting ejecta are assumed to accelerate or decelerate while producing the radiation. If the ejecta are accelerating, which may be the case for the magnetically dominated jet of the ICMART model (Zhang & Yan 2011), the decay index may temporarily reach a value of we observe even assuming as the trigger of the prompt emission and . A magnetically dominated jet thus appears to be a reasonable solution for the GRB at hand.

However, whichever solution applies, one would always have to assume that the emission is produced inside the ejecta and not in the medium surrounding the explosion as in the case of the FS scenario. Finally, in the previous section we made it clear that the final part of the X-ray light-curve, after the end of the steep decay, appears to be FS emission.

We conclude that the X-ray emission up to the 100 ks break is not produced by the external, circumburst medium energized by the FS, but it is instead of “internal origin”, generated directly within the explosion outflow. This component stops abruptly and the X-ray flux drops rapidly, until the FS emission in the X-ray band prevails. The optical emission is basically dominated by the FS mechanism (with the possible exception of the early flare; see Sect. 4.5). In the next sections, we shall discuss the possible origins of the high-energy emission between the end of the prompt emission and the fall of the X-ray flux at 100 ks.

#### 4.2.2 A Magnetar Central Engine

The nature of the X-ray afterglow and the nature of the central engine are two of the many open questions in the contemporary field of GRBs (see Zhang 2011 for a review), even more so because this GRB shows that we may be misinterpreting the behaviour of other X-ray light-curves solely or primarily attributed to FS. For example, if we had not observed the steep drop at 100 ks for 130831A, we might have easily mistaken the relatively ordinary decay and spectral slopes as being produced by the standard FS-emission.

The core of the stellar progenitor of GRB 130831A may have collapsed into a magnetar (see for example Thompson et al. 2010). The rotational kinetic energy of such objects is

 Erot=3×1052(M1.4M⊙)(R12km)2P−2ms erg (3)

where , and are the mass, radius and period in ms of the object, respectively. Assuming unitary values for the parameters above, erg, which is enough to power the supernova explosion and ultra-relativistic outflow. If the mass M is closer to the Tolman - Oppenheimer - Volkoff limit, which is thought to be slightly larger than (Antoniadis et al. 2013), the parameter might be slightly different (Metzger et al. 2015).

In a simple scenario, the magnetar would initially tap into rotational energy and produce a very energetic outflow, likely roughly collimated into bi-polar jets. Such a wind would be produced through the process of dipole spin-down. The energy of the outflow is initially imparted to the stellar envelope, causing, or at least contributing to, the supernova explosion. Moreover, the magnetar outflow may be long-lived, and produce radiation that we observe (see below). In this scenario, the luminosity of the magnetar could be roughly constant, even for a relatively long timescale (Zhang & Mészáros 2002), depending on physical parameters. After , or if the magnetar collapses into a black hole, the light-curve would show a flux drop. This magnetar model is relevant for a few GRBs that show a plateau with approximately constant flux, with , followed by the steep slope segment (see cases studied by Liang et al. 2007, Lü & Zhang 2014, Cano et al. 2014, Bernardini et al. 2014, Lü et al. 2015). However, the early decay slope of 130831A is , which is in contrast to the aforementioned cases, and requires a more complicated model to explain our observations.

#### 4.2.3 A Magnetar with decaying magnetic field

In the simple spin-down calculation above, one assumes that the magnetic field is constant and independent of the period of the magnetar . Metzger et al. (2011), however, assume that is linked to , taking into account that the magnetic field is generated from the energy available in differential rotation. They thus estimate

 B=1016R−1/26P−1−3 G (4)

Under this assumption, decreases with increasing , i.e., with increasing time since the explosion. In such conditions, the jet luminosity is not constant. As shown in Figures 2 and 5 of Metzger et al. (2011), decreases with time, scaling approximately as from 100–1000 s to tens of ks after the collapse for reasonable values of parameters at the beginning of the spin-down, such as G and  ms. The predicted luminosity of the jet seems to be in the right range to explain the X-ray emission of GRB 130831A during the slow decay phase, too. For the values quoted above, Metzger et al. (2011) predict a few erg s at 10 ks; at the same epoch, GRB 130831A had a luminosity of erg s.

We note that a jet luminosity does not convert immediately into X-ray luminosity ; one has to take into account that the radiation mechanism will have a certain efficiency. Moreover, we have not yet considered that GRB emission is beamed. We can write

 LX=ηXf−1b,xLjet (5)

where and are the efficiency in converting the jet luminosity into X-ray radiation and the correction for the beaming, respectively.

According to Metzger et al. (2011), at late epochs ( s after the collapse), the magnetar outflow is highly relativistic and Poynting-flux dominated; in such conditions, internal shocks and reconnections within the jet itself are not possible. Forced reconnection, however, can occur at large radii, when the outflow collides with the circumburst medium and/or the previous ejecta, and convert the jet energy into X-ray emission. Assuming that the efficiency of this process is similar to the one which generates the prompt emission, and the correction for beaming is 10 times lower than that during the prompt emission, Metzger et al. (2011) find that they can explain the observed and plateau durations of several GRBs similar to GRB 130831A.

The newly born magnetar may also provide a large energy input in the exploding progenitor, powering an energetic and luminous supernova explosion. This is in agreement with observations of Cano et al. (2014), who find a kinetic energy for the (non-relativistic) ejecta of SN 2013fu of erg and a peak absolute magnitude . We note that Greiner et al. 2015 found erg for SN2011 kl, a very bright supernova associated with GRB 111209A, whose properties can be explained by the energy injection of a newly born magnetar. In addition, the energetics of SN2013fu is quite typical of other SNe associated with GRBs (Cano et al. 2015).

As for the abrupt end of the X-ray emission of “internal origin”, with a very steep slope, we may attribute it to the delayed collapse of the magnetar into a black hole (Vietri & Stella 1998, Lyons et al. 2010, Rowlinson et al. 2013). Once the magnetar has lost much of its rotational energy to power the jet, the weakened centrifugal forces may not be able to avoid the collapse. The timescale of such event, for a neutron star mass of , initial period of s and magnetic field of Gauss would be  s (Vietri & Stella 1998, their Eq. 1), comparable to the epoch of the steep drop in the X-ray light curve of GRB 130831A. Such a collapse should be relatively quick, and rapidly stop energy emission from the central object. If a magnetar collapsing into a black-hole is the right model, GRB 130831A might be a candidate for the production of Fast-Radio Bursts (FRBs), as described in Zhang (2014), as thus a target for observational campaigns in radio aimed at understanding the origin of FRBs.

#### 4.2.4 A black hole with a fall-back accretion disk

Another possibility we consider is that the central engine of this GRB could be a stellar black hole with a fall-back accretion disk (Kumar et al. 2008). Basically, in the supernova explosion associated with the GRB, the innermost part of the star collapses into a black hole. A continued fall-back of matter - directly from the progenitor envelope or from supernova ejecta that failed to reach escape velocity - occurs at the centre. Some of this material does not accrete directly onto the black hole, but creates an accretion disk around it. Depending on the fall-back rate and the accretion time onto the black hole, the material of this disk can power relativistic ejecta from the black hole for a long time, which may produce both the prompt emission and a long-lived, slowly decaying X-ray flux.

For the latter, Kumar et al. (2008) envisage a few possibilities. One is that the accreting disk has a low viscosity. It will thus take a long time, of the order of  s, for the all the disk material to accrete onto the black hole and power the jet. This model explains a long-lived plateau, but it predicts a flux decay with slope at the end of this phase, which is not observed in GRB 130831A and other GRBs with similar features.

Another possibility is that the disk has high viscosity. Then is much less than the fall-back time, and the emission basically traces the rate of the matter falling back onto the accretion disk. However, such a scenario cannot explain why the fall-back rate is less steep than expected: on theoretical grounds, we expect a fall-back rate to vary as (Chevalier 1989). Secondly, the plateau slope and the sudden cut-off can be explained only by a particular density and angular momentum profile of the matter of the progenitor. A luminosity that evolves close to , as in the case of GRB 130831A, might be explained if the density profile of the stellar envelope has a profile of approximately , where is the radius from the centre of the star. A steep drop at the end of the plateau might be achieved if the material of the stellar envelope, which is the last to be accreted, has a relatively small angular momentum: this will cause the matter to fall rapidly onto the black hole and shut off the emission. Such a peculiar configuration, however, seems somewhat contrived and at odds with models of stellar progenitors (e.g. Woosley et al. 2011). Moreover, to support accretion for  s, the disk should be unusually large and massive, leaving little mass for ejecta. Cano et al. (2014) find instead that the SN 2013fu, associated with GRB 130831A, has an ejecta mass of , which is typical of other GRB-SNe (e.g. Cano 2013).

For the model discussed, the jet luminosity is expected to be in the range of erg s at the end of the plateau, which is similar to the X-ray luminosity of the GRB 130831A afterglow at the end of the shallow decline phase.

#### 4.2.5 A Binary Origin

Barkov & Komissarov (2010) conceive another scenario in which a black hole might power a GRB and long-lived outflow (see also Komissarov & Barkov, 2009). If a compact object and a Wolf-Rayet (WR) star form a very close binary, such a system can go through a common envelope phase in which the compact object spirals in and can accrete the matter of the companion. The common envelope matter will have very high angular momentum and can take a long time to accrete. According to Barkov & Kommisarov, the accretion timescale of such a system is:

 td≃8000(α0.01)−1(RsR⊙)3/2(Mc2M⊙)2(Ms10M⊙)−7/2s (6)

where represents the viscosity, and the radius and mass of the WR star, and is the mass of the compact object. For massive compact objects, an accretion timescale of several tens of ks is not impossible666Note, though, that Eq. 6 is valid if is considerably smaller than . However, if they become comparable, e.g. when we have a black hole of 10 and Wolf-Rayet star of similar mass, one has to find a different method to estimate .. According to Barkov & Kommisarov (2010), during accretion the compact object will produce jets via the Blandford-Znajek mechanism, and the jet luminosity will be in the range of erg s, more than enough to explain the luminosity during the slow decline phase. Nonetheless, this scenario might suffer from the same problems as the previous one. If the viscosity of the accretion disc is low, which is required to keep the material around the black hole for  day, then the flux should not decrease quickly at the end of the plateau. The binary origin scenario could not reproduce the steep X-ray flux drop and/or it would require a peculiar structure of the WR star to explain the temporal dependence of the luminosity. We note, however, that in their papers, Barkov and Kommissarov do not discuss what happens at the end of the plateau, so we can only postulate.

Finally, we remark that this binary origin model can predict a supernova, like the one associated with GRB 130831A.

#### 4.2.6 Concluding remarks on the origin of the X-ray emission

We have shown that the FS scenario (or any refreshed shock scenario) cannot explain the X-ray emission of GRB 130831A between the end of prompt emission and  ks, since such a model cannot entail the steep decay of the flux at that epoch.

We have investigated whether such early X-ray emission can be attributed to dissipation processes occurring in the outflow of a newly born magnetar, produced via spin-down energy extraction of the compact object. We have found that the simplest model cannot explain the observations, because it predicts a flat X-ray light-curve followed by a steep drop. In the case of GRB 130831A, the decay slope before the break is , which is much steeper than what we expect in this model. However, a more elaborate model of the magnetar spin-down, in which the magnetic field is expected to decay as the rotation time increases, predicts a luminosity decay more consistent with observations and a duration that can extend up to tens of ks. Moreover, the anticipated X-ray luminosity is in the right range. The assumed initial parameters - initial period and magnetic field - for the newly born magnetar that would produce such a X-ray light-curve are  ms and G, which are expected on theoretical grounds for such an object.

We have also discussed whether the compact object could be a black hole rather than a magnetar. In order to power emission for such a long time, the black hole should be surrounded by an extended disk, perhaps produced by fall-back material of the supernova associated with the GRB. While this model still may predict the right luminosity and duration for the “internal” X-ray emission, it would require an anomalous distribution of angular momentum in the progenitor star, and peculiar fall-back rates. Similar advantages and disadvantages may be present in a model in which the compact object forms a binary with a WR star and spirals in, blowing up the star into a massive disk. All in all, we deem the magnetar model with a decaying magnetic field the most plausible of those presented so far to explain the properties of GRB 130831A.

### 4.3 Energy budget of the X-ray radiation between the end of the prompt emission and the 100 ks drop.

In the cosmological rest frame, the X-ray internal emission begins no later than  s and lasts until  s. Taking into account cosmological corrections, the 0.3–10 keV luminosity at 135 s is  erg s; we assume that from this epoch the luminosity decreases as up to 66.4 ks. We subtract the amount of X-ray emission produced by the FS (see section 4.5), and we obtain a total energy erg. Such a value represents of the energy emitted during the prompt phase, and only of the kinetic energy of the relativistic ejecta producing the FS emission. We note, however, that the internal emission may extend below 0.3 keV and above 10 keV. Thus, and the percentage we have determined represent a lower limit.

### 4.4 Energy breakdown of GRB 130831A and the associated supernova

We know that the kinetic energy of the ejecta of SN 2013fu, the supernova associated with GRB 130831A, is erg (Cano et al. 2014). This is a factor lower than the kinetic energy of the relativistic ejecta and comparable to the energy emitted during the prompt emission, and times higher than the energy associated with the X-ray emission of “internal origin”. We do caution, however, that the estimated SN kinetic energy is isotropic, while the energy of the relativistic ejecta, the prompt and the early X-ray energetics must be corrected for an unknown beaming factor. GRB 130831A does not show the signature of a jet break (Sari, Piran & Halpern 1999) within our observations, so we cannot derive the beaming angle of the ejecta and thus the beaming factor. However, we can set limits. The detection of the X-ray afterglow by Chandra at 1430 ks after the trigger indicates that the FS X-ray afterglow had a decay slope of up to that epoch and no jet break had yet occurred.
Following Zhang et al. 2009, the beaming angle of the ejecta in a medium with constant density can be estimated as

where is the jet break time in days. Our first Chandra observation took place 16.6 days after the trigger and we adopt the values and determined from our modeling; we note that the exact value of depends only weakly on the value of these parameters anyway. We infer that  rad; thus the lower limit on the beaming factor is . By definition, the maximum value for the beaming factor is when the source is isotropic.
If the beaming factor of GRB 130831A emission were , the total energy budget of this event and its supernova would be erg. The prompt energy , the energy emitted in X-rays up to 100 ks , and the kinetic energy of the relativistic ejecta would be , and of the total energy budget. If, as an extreme and unlikely case, GRB 130831A were isotropic, its total energy budget would be erg; the above percentages would become , and . The kinetic energy of the relativistic ejecta is at least of the total energy produced by the GRB and the SN. Moreover, the fraction of energy going into the “internal emission” X-rays is always rather small, being substantially less than in both cases.
If GRB 130831A and its SN are powered by a magnetar, the total energy budget cannot be erg (see section 4.2). To not exceed this limit, the beaming factor of the GRB must be . If , , , and represent , and of the total energy. In reality, we should expect these percentages to be between those of the and  cases if GRB 130831A is actually powered by a magnetar. The breakdown is presented again in Table 8.

### 4.5 Early Afterglow

In our analysis, we have focused on the afterglow emission between 15 and 230 ks. It is worth exploring whether our model can explain the interesting features of the early afterglow, especially the optical band.

The initial flare peaks at 730 s. It takes place both in the X-ray and optical bands, but it is very pronounced in the latter. Its rapid temporal evolution (the optical flux increased by a factor of between 400 and 800 s after the trigger) suggests that it could be explained in the context of internal dissipation processes which occurred in the outflow, when the Lorentz factor is very high and relativistic effects cause rapid variations of the observed flux. Thus, GRB 130831A may show a clear example of internal dissipation that produces strong emission in the optical other than in the X-ray. Alternatively, the optical flare might flag the onset of FS emission. However, if this were the case, the decay slope after the flare peak would be consistent with the decay slope of the late optical afterglow . Instead, the decay rate after the flare peak is , which is inconsistent with found later. After the flare, the early optical emission shows a plateau up to a few ks (see Section 3). Typically, an early slow optical decay is interpreted as energy injection, which ends at the time of the break. However, this interpretation might be difficult in our scenario, because the energy injection would have to stop at ks while, according to our analysis, the GRB outflow is still active at ks.

A possibility is that the optical plateau basically results from the combination of the decaying optical flare and the rising of the FS peak, which has been shown for a number of Swift bursts for which early optical light-curves are available (Oates et al. 2009). The peak Lorentz factor of the ejecta then can be calculated (Molinari et al. 2007 and references therein) as

 Γ=160(Eγ,53(1+z)3η0.2 n t3dec,2)1/8 (8)

For a deceleration time  s and adopting the values of energy and density we have determined above, the resulting peak Lorentz factor would be , which is within the overall distribution of Lorentz factors for GRB afterglows (Oates et al. 2009). The reason for such a late deceleration, in our model, comes naturally given the low density of the circumburst medium.

The late emission, which we attribute to FS in our modeling, seems to have a relatively steep decay slope. So one may wonder whether it could give some important contribution to the X-ray flux as well at an earlier epoch. If we extrapolate the late X-ray flux to earlier epochs using a decay slope of , it would become comparable to or even higher than the observed X-ray flux at the end of the first orbit (at  s) and the shape of the X-ray light-curve would differ from what we see. However, if the FS onset occurs at  s, this problem is avoided. As for the observations from  ks (i.e. the beginning of the second orbit) onwards, at 10 ks the flux by FS emission is erg cm s erg cm s, which is times weaker than the observed X-ray flux. Afterwards, the FS flux decreases faster than that of “internal origin” that we see, which decays with a slope of . Similarly, for a and flux density in the R band ( Hz) of Jy at 10 ks, the expected X-ray 0.3-10 keV flux is erg cm s. This is again times lower than the X-ray flux observed at 10 ks. Thus, the X-ray flux from 9 ks up to the steep drop is not predominantly produced by the FS.

## 5 Conclusions

We have discussed the case of the long Swift GRB 130831A. The X-ray afterglow of this burst initially shows a shallow decay. However, at ks the X-ray light-curve breaks to an unusually steep decay slope , which cannot be explained by the standard FS model. Late XRT and especially Chandra observations show that the X-ray afterglow has a successive break to a more sedate decay with a slope .

The well-sampled optical afterglow shows no change of slope concurrent with the steep break in the X-ray band, which we interpret as arising from a different mechanism. By means of data taken by Swift, Chandra and several ground-based observatories, we have shown that both the optical and late X-ray emissions, after  ks after the trigger, have the typical decay and spectral slopes of GRB afterglows explained by the FS model.

We have interpreted the X-ray and optical afterglow as the superposition of two emission components. One component is of “internal” origin, generated within a relativistic outflow and responsible for the early X-ray emission up to  ks; the outflow is produced by either the spin-down of a newly formed magnetar or a black hole feeding on fall-back matter. A second component, responsible for the late X-ray and optical from a few ks after the trigger, is the typical FS emission. When the magnetar has lost much of its rotational energy or the black hole does not accrete and does not power the outflow any longer, the first component dies off and we see the steep X-ray decay that lasts until the standard FS emission emerges. We believe that the magnetar model is favoured, since the other scenario would require a rather peculiar stellar progenitor structure and fall-back process.

Modeling the late optical and X-ray afterglow, we have inferred the kinetic energy of the relativistic ejecta and thus an efficiency of the “central engine” of this GRB to produce -ray emission. This efficiency is smaller than that of other bursts that show emission of internal origin (Lü & Zhang 2014; although these were examined in the X-ray band only), and more typical of those GRBs in which no internal emission is clearly visible. Thus, GRB 130831A may represent a “trait d’union” between GRBs with different dominant emission processes.

More importantly, gathering the information on the kinetic energy of the supernova associated with 130831A, we have provided a breakdown of the energetics of the GRB and its associated supernova. We have found that, regardless of the nature of the central engine and unknown collimation of the ejecta, at least of the total energy of the event is coupled with relativistic ejecta; and less (probably significantly less) than of the energy goes into X-ray emission of “internal origin” lasting up to 100 ks in our case; this component produces a factor less energy than that released in -ray during the prompt emission.

Showing several emission processes at work, GRB 130831A has offered us the opportunity to investigate the complete phenomenon of a supernova with a central engine that produces the explosion and drives an energetic relativistic outflow, where dissipation processes take place for a long time.

## Acknowledgments

We thank H. Tananbaum for granting us DDT observations of GRB 130831A with Chandra. This research has made use of data obtained from the Chandra Data Archive and the Chandra Source Catalog, and software provided by the Chandra X-ray Center (CXC) in the application packages CIAO, ChIPS, and Sherpa.

MDP, MJP, SRO and AAB thank UK Space Agency for financial support. This work made use of data supplied by the UK Swift Science Data Centre at the University of Leicester.

AP, AV acknowledge partial support by RFBR grants 12-02-01336, 13-01-92204, 14-02-10015 and 15-02-10203.

AJCT and SRO acknowledges support from the Spanish Ministry Grant AYA 2012-39727-C03-01

SS acknowledges support from CONICYT-Chile FONDECYT 3140534, Basal-CATA PFB-06/2007, and Project IC120009 ”Millennium Institute of Astrophysics (MAS)” of Iniciativa Científica Milenio del Ministerio de Economía, Fomento y Turismo.

DAK acknowledges financial support by the Thüringer Landessternwarte Tautenburg, and the Max-Planck Institut für Extraterrestrische Physik.

ZC gratefully acknowledges support by a Project Grant from the Icelandic Research Fund.

We thank the RATIR project team and the staff of the Observatorio Astronï¿½mico Nacional on Sierra San Pedro Mártir. RATIR is a collaboration between the University of California, the Universidad Nacional Autonóma de México, NASA Goddard Space Flight Center, and Arizona State University, benefiting from the loan of an H2RG detector and hardware and software support from Teledyne Scientific and Imaging. RATIR, the automation of the Harold L. Johnson Telescope of the Observatorio Astronï¿½mico Nacional on Sierra San Pedro Mï¿½rtir, and the operation of both are funded through NASA grants NNX09AH71G, NNX09AT02G, NNX10AI27G, and NNX12AE66G, CONACyT grants INFR-2009-01-122785 and CB-2008-101958 , UNAM PAPIIT grant IN113810, and UC MEXUS-CONACyT grant CN 09-283.

Partly based on observations carried out with the 10.4m Gran Telescopio Canarias (GTC) installed in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofísica de Canarias in the island of La Palma.

## References

• Ahn (2012) Ahn C.P, Alendandroff, R., Allende-Prieto, C., 2012, ApJS, 203, 21
• Amati (2006) Amati L., 2006, MNRAS, 372, 233
• Amati et al. (2009) Amati L., Frontera F., Guidorzi C. 2009, A&A, 508, 173
• Antoniadis et al. (2013) Antoniadis J., Freire P., Wex C. 2013, Science, 508, 173
• Arnaud et al. (1996) Arnaud K. A., 1996, Astronomical Data Analysis Software and Systems V, eds. Jacoby G. and Barnes J., p17, ASP Conf. Series volume 101.
• Band et al. (1993) Band D., Matteson J., Ford L., et al., 1993, ApJ, 413, 281
• Barkov & Komissarov et al. (2010) Barkov M., & Komissarov S., 2010, MNRAS, 401, 1644
• Barthelmy et al. (2005) Barthelmy S. D., Barbier L. M., Cummings J. R., et al., 2005, SSR, 120, 143
• Bernardini et al. (2014) Bernardini M. G, Campana S., Ghisellini G., et al., 2014, MNRAS, 439L, 80
• Bertin & Arnouts (1996) Bertin E., Arnouts S., 1996, A&A, 117, 393
• Bertin (2010) Bertin E., et al., 2010, Astrophysics Source Code Library
• Blandford & Znajek (1977) Blandford R. D. & Znajek R. L., 1977, MNRAS, 179, 433
• Bloom et al. (2001) Bloom J. S., Frail D. A., Sari R. 2001, AJ, 121, 2879
• Burrows et al. (2005) Burrows D. N., Hill J. E., Nousek J. A., et al., 2005, SSR, 120, 165
• Butler et al. (2012) Butler N., Klein C., Fox O., et al., 2012, SPIE, 8446, 10
• Butler et al. (2013) Butler N., Watson A., Kutyrev A., et al., 2013, GCN Circ. 15165
• Breeveld et al. (2011) Breeveld A. A., Landsman W., Holland S. T., et al., 2011, AIPC, 1358, 373
• Cano et al. (2014) Cano Z., De Ugarte Postigo A., Pozanenko A., et al., 2014, A&A, 568, 19
• Cano (2013) Cano Z., 2013, MNRAS, 434, 1098
• Cenko et al. (2011) Cenko S. B., Bloom J. S., Kulkarni S. R., et al., 2011, MNRAS, 420, 2684
• Chandra (2012) Chandra P. & Frail, D. A. ApJ, 746, 156
• Cepa et al. (2000) Cepa J., Aguiar, M., Escalera, V. G., et al., 2000, SPIE 4008, 623
• Chevalier (1989) Chevalier R. A., 1989, ApJ, 346, 847
• Curran et al. (2010) Curran P. A., Evans P. A., De Pasquale M., et al., ApJL, 716, 135
• Cucchiara et al. (2013) Cucchiara A., Perley D., 2013, GCN Circ. 15144
• Dall’Osso et al. (2011) Dall’Osso S., Stratta G., Guetta D., et al., 2011, A&A, 526 121
• De Ugarte-Postigo (2012) De Ugarte Postigo A., Lundgren A., Martin S., et al. A&A, 538, 2012
• Duffel & MacFadyen (2014) Duffel P. C. & MacFadyen A. 2014, ApJL, 791, 1
• Evans et al. (2009) Evans P. A., Beardmore A. P, Osborne J. P., et al., 2009, MNRAS, 399, 1167
• Gao et al. (2013) Gao, H., Lei, W.-H., Zou, Y-.C. et al., 2013, NewAR, 57, 141
• Gehrels (1986) Gehrels N., 1986, ApJ, 303, 336
• Gehrels et al. (2004) Gehrels N., Chincarini G., Giommi P., et al., 2004, ApJ, 611, 1005
• Golenetskii et al. (2013) Golenetskii S., Aptekar R., Fredericks D., et al., 2013, GCN Circ., 15145
• Granot (2006) Granot J., 2006, RMXAC, 270, 140.
• Greiner et al. (2015) Greiner J., Mazzali, P. A., Kann, D. A., et al., 2015, Nature, 523, 189
• Falcone et al. (2007) Falcone A. D., Morris D., Racusin J. L., et al., 2007, ApJ, 671, 1921
• Fruscione et al. (2006) Fruscione A., et al. 2006, SPIE proc., 6270
• Henden & Munari (2014) Henden A. & Munari U., 2014, Contributions of the Astronomical Observatory Skalnate Pleso, 43, 518
• Jarosik et al. (2011) Jarosik N., Bennet, C. L., Dunkley J., et al., 2011, ApJS, 192, 14.
• Hjorth et al. (2012) Hjorth J., Malesani D., Jakobsson P., et al., 2012, ApJ, 756, 187
• Jordi et al. (2016) Jordi K., Grebel E. K., Ammon K., 2006, 460, 339
• Kalberla et al. (2005) Kalberla, P. M. W., Burton W. B, Hartmann D., et al., 2005, ApJ, 440, 775
• Klein et al. (2012) Klein, C. R., Kubanek P., Butler N. R., et al., 2012, 8453, 2
• Klose et al. (2013) Klose S., Nicuesa Guelbenzu A., Krühler T., et al., 2013, CBET 3677
• Kobayashi & Zhang (2007) Kobayashi S. & Zhang B., 2007, ApJ, 655, 973
• Komissarov & Barkov (2009) Komissarov S., Barkov M., 2009, MNRAS, 397, 1153
• Kraft et al. (1991) Kraft, R. P., Burrows D. N., Nousek J. A. 1991, ApJ, 374, 344
• Kumar et al. (2008) Kumar P., Narayan R., Johnson J. L., et al., 2008, MNRAS, 388, 1729
• Kumar & Panaitescu (2000) Kumar P., Panaitescu A., 2000, ApJL, 541, 51
• Lang et al. (2010) Lang D., Hogg D. J., Mierle, K. et al., 2010, AJ, 139, 1782
• Laskar et al. (2013) Laskar T., Zauderer A., Berger E., 2013, GCN Circ. 15162
• Liang et al. (2006) Liang E.-W., Zhang B., O’Brien, P., et al., 2006, ApJ, 646, 351
• Liang et al. (2007) Liang E.-W., Zhang B.-B., Zhang B., 2007, ApJ, 670, 565
• Lü & Zhang (2014) Lü, H. J., Zhang, B. 2014, ApJ, 785, 74
• Lü et al. (2015) Lü, H. J., Zhang, B., Lei, W.-H., ApJ, 805, 89
• Lyons et al. (2010) Lyons N., O’Brien P. T., Zhang B., et al., 2010, MNRAS, 402, 705
• Margutti et al. (2011) Margutti R., Bernardini G., Barniol-Duran R., et al., 2011, MNRAS, 410, 1064
• Matheson et al. (2003) Matheson T., Garnavich P. M., Stanek K. Z., et al., 2003, ApJ, 599, 394
• Metzger et al. (2011) Metzger B. D., Giannios D., Thompson T. A., et al., 2011, MNRAS, 413, 2031
• Metzger et al. (2015) Metzger B. D., Margalit, B., Kasen, D., et al., 2015, MNRAS Letters submitted, arXiv:1508.02712
• Molinari et al. (2011) Molinari E., Vergani S., Malesani D., et al., 2007, A&A, 469, 13
• Molotov et al. (2008) Molotov I., Agapov V., Titenko V. D., et al., 2008, AdSpR, 41, 1022
• Nousek et al. (2006) Nousek J. A., Kouveliotou C., Grupe D., et al., 2006 ApJ, 1, 389
• Oates et al. (2009) Oates S. R., Page M. J., Schady P., et al., 2009, MNRAS, 395, 490
• Panaitescu & Kumar (2002) Panaitescu A. & Kumar P., 2002, ApJ, 571, 779
• Pei (1992) Pei Y.-C., 1992, ApJ, 395, 130
• Perley et al. (2014) Perley D. A., Cenko S. B., Corsi A., et al., 2014, ApJ, 781, 37
• Poole et al. (2008) Poole T. S., Breeveld A. A., Page M. J., et al., 2008, MNRAS, 383, 627
• Pozanenko et al. (2013a) Pozanenko A., Volnova A., Hafizov B., et al., 2013a, GCN Circ. 15190
• Pozanenko et al. (2013b) Pozanenko A., Volnova A., Hafizov B., et al., 2013b, GCN Circ. 15237
• Pozanenko et al. (2013c) Pozanenko A., Elenin L., Litvinenko E., et al., 2013c, EAS Pub. Series, 61, 259
• Racusin et al. (2009) Racusin J. L., Liang E. W., Burrows D. N., et al., 2009, ApJ, 698, 43
• Reichart et al. (2005) Reichart D., Nysewander M., Moran J., et al., 2005, NCimC, 28, 767
• Rieger et al. (2007) Rieger F., Bosch-Ramon F., Duffy P., 2007, Ap&SS, 309, 119
• Rowlinson et al. (2013) Rowlinson A., O’Brien P. T., Metzger B. D., et al., 2013, MNRAS, 430, 1061
• Roming et al. (2005) Roming P. W. A., Kennedy T. E., Mason K. O., et al., 2005, SSRV, 120, 95
• Sari et al. (1999) Sari R., Piran T., Halpern, J. P., 1999, ApJL, 519, 17
• Swenson et al. (2013) Swenson C. A., Roming P. W. A. De Pasquale, M., et al., 2013 ApJ, 774, 2
• Schady et al. (2010) Schady P., Page M. J., Oates S. R., et al., 2010, MNRAS, 401, 2773
• Schlegel et al. (1998) Schlegel D. J., Finkbeiner D. P., Davis, M., 1998, ApJ, 500, 525
• Shen & Zhang (2009) Shen R.-F. & Zhang B. Mészáros P., 2009, MNRAS, 398, 1936
• Thompson (2007) Thomson T. A., 2007, RMxAC, 27, 80
• Thompson et al. (2010) Thomson T. A., Metzger B. D., Bucciantini, N., 2010, AIPC, 1279, 81
• Troja et al. (2007) Troja E., Cusumano, G., O’Brien P. T., et al., 2007 ApJ, 665, 599
• Tody (1986) Tody D., 1986, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 627, 733
• Tody (1993) Tody D., 1993 Astronomical Data Analysis Software and Systems II, A.S.P. Conference Ser., Vol 52, eds. R.J. Hanisch, R.J.V. Brissenden, & J. Barnes, 173
• Trotter et al. (2013) Trotter A., Haislip J., Lacluyze A., et al., 2013, GCN Circ. 15148
• Usov (1992) Usov V. V., Nature, 357, 472
• Van Eerten et al. (2010) van Eerten H., Zhang, W., MacFadyen, A., 2010, ApJ, 722, 235
• Vietri & Stella (1998) Vietri M., Stella, L. 1998, ApJL, 507, 45
• Volnova et al. (2013a) Volnova A., Matkin A., Stepura A., et al., 2013a, GCN Circ. 15185
• Volnova et al. (2013b) Volnova A., Linkov V., Polyakov K., et. al., 2013b, GCN Circ. 15188
• Volnova et al. (2013c) Volnova A., Krugly Y., Slyusarev I., et. al., 2013c, GCN Circ. 15189
• Watson et al. (2012) Watson A. M., Richer M. G., Bloom J. S., et al., 2012, SPIE, 8444, 5
• Woosley (2011) Woosley S.E., 2011, in ”Gamma-Ray Bursts”, Cambridge University Press, http://arxiv.org/abs/1105.4193
• Zauderer (2013) Zauderer B. A., Laskar T., & Berger E., 2013, GCN Circ. 15159
• Zhang et al. (2006) Zhang B., Fan Y.-Z., Dyks J., et al., 2006, ApJ, 642, 345
• Zhang (2011) Zhang B., 2011, CRPhy, 12, 206
• Zhang et al. (2007) Zhang B., Liang E.-W., Page K., et al., 2007, ApJ, 655, 989
• Zhang & MacFadyen (2009) Zhang W.-Q., MacFadyen, A., 2009, ApJ, 698, 1261
• Zhang & Mészáros (2001) Zhang B. & Mészáros P., 2001, ApJL, 552, 35
• Zhang & Mészáros (2002) Zhang B. & Mészáros P., 2002, ApJ, 566, 712
• Zhang (2009) Zhang X.-H., 2009, Research in Astronomy and Astrophysics, 9, 213
You are adding the first comment!
How to quickly get a good reply:
• Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
• Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
• Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters