NDAFs for LGRBs

Testing black hole neutrino-dominated accretion discs for long-duration gamma-ray bursts

Cui-Ying Song, Tong Liu, Wei-Min Gu, and Jian-Xiang Tian
Shandong Provincial Key Laboratory of Laser Polarization and Information Technology, Department of Physics,
Qufu Normal University, Jining, Shandong 273165, China
Department of Astronomy, Xiamen University, Xiamen, Fujian 361005, China
Department of Physics and Astronomy, University of Nevada, Las Vegas, NV 89154, USA
E-mail: tongliu@xmu.edu.cn
Accepted XXX. Received YYY; in original form ZZZ
Abstract

Long-duration gamma-ray bursts (LGRBs) are generally considered to originate from the massive collapsars. It is believed that the central engine of gamma-ray bursts (GRBs) is a neutrino-dominated accretion flow (NDAF) around a rotating stellar-mass black hole (BH). The neutrino annihilation above the NDAF is a feasible mechanism to power GRB. In this work, we analyse the distributions of the isotropic gamma-ray radiated energy and jet kinetic energy of 48 LGRBs. According to the NDAF and fireball models, we estimate the mean accreted masses of LGRBs in our sample to investigate whether the NDAFs can power LGRBs with the reasonable BH parameters and conversion efficiency of neutrino annihilation. The results indicate that most of the values of the accreted masses are less than for the extreme Kerr BHs and high conversion efficiency. It suggests that the NDAFs may be suitable for most of LGRBs except for some extremely high energy sources.

keywords:
accretion, accretion discs - black hole physics - gamma-ray burst: general - neutrinos
pubyear: 2016pagerange: Testing black hole neutrino-dominated accretion discs for long-duration gamma-ray burstsLABEL:lastpage

1 Introduction

Gamma-ray bursts (GRBs) are extremely energetic transient events and isotropically distribute over the sky. Depending on a separation at 2 s of their duration , they are generally divided into two classes, i.e., long- and short-duration GRBs (LGRBs and SGRBs, see Kouveliotou et al., 1993). The percentage rates for LGRBs and SGRBs in the BATSE sample are about and (Sakamoto et al., 2008, 2011), respectively, and the LGRB fraction is further larger for other space detectors, such as Swift and Fermi (Zhang et al., 2012). The progenitors of LGRBs are usually considered as the collapse of massive stars with rapid spin, low metallicity, and stripped of their helium and hydrogen envelope (e.g., Woosley, 1993; Kumar & Zhang, 2015). In this circumstance, the stellar-mass Kerr black hole (BH) surrounded by a debris torus or disc forms quickly in the central parts of the stellar core. Moreover, SGRBs are related to the mergers of two neutron stars (NSs) or NS-BH binaries (e.g., Eichler et al., 1989; Narayan et al., 1992; Nakar, 2007). A similar accretion system may form in SGRBs due to the loss of the angular momentum and mass of the remnants.

Neutrino annihilation and Blandford-Znajek (BZ, see Blandford & Znajek, 1977) mechanisms are proposed to power GRBs if there is a hyperaccretion disc in the centre of GRBs. A typical cosmological GRB requires the value of the accretion rate ranging from a fraction of solar mass per second to several solar masses per second. If the neutrino annihilation is the dominant cooling mechanism, a geometrically and optically thick disc is formed and a large number of free baryons appear in the inner region of the disc. The temperature and density are so high that photons are trapped and neutrino cooling becomes effective. This disc is named the neutrino-dominated accretion flow, which has been widely studied (NDAF, see, e.g., Popham et al., 1999; Narayan et al., 2001; Di Matteo et al., 2002; Kohri & Mineshige, 2002; Kohri et al., 2005; Gu et al., 2006; Chen & Beloborodov, 2007; Kawanaka & Mineshige, 2007; Liu et al., 2007, 2008, 2012a, 2012b, 2013, 2015d; Li & Liu, 2013; Janiuk et al., 2013; Kawanaka et al., 2013; Xue et al., 2013; Hou et al., 2014a, b).

Another generally accepted model to power GRBs is the magnetar model (e.g., Usov, 1992; Thompson et al., 2004; Metzger et al., 2008, 2011; Mao et al., 2010; Lü & Zhang, 2014). The released energy via its spinning down can effortlessly provide the requirements of GRBs (e.g., Usov, 1992; Wheeler et al., 2000), and this model may offer a plausible explanation for the shallow decay phases in some GRB afterglow light curves (e.g., Dai, 2004; Zhang et al., 2006). No matter what kind of models, the puzzled central engines of GRBs are hidden by the electromagnetic radiation (e.g., Liu et al., 2015d), and the BH-NDAF systems may be probed directly by gravitational waves (e.g., Suwa & Murase, 2009; Romero et al., 2010; Sun et al., 2012) and MeV neutrino emission (e.g., Liu et al., 2015e).

Fan & Wei (2011) and Liu et al. (2015c) investigated the capabilities of the NDAFs to power SGRBs. As a result, the disc mass of a certain SGRB mainly depends on the output energy of the NDAF, jet opening angle, and characteristics of central BH. Almost all of SGRBs are satisfied with the reasonable disc masses, only a few approach or exceed the limits in simulations, (e.g., Kluźniak & Lee, 1998; Lee & Kluźniak, 1999; Popham et al., 1999; Liu et al., 2012b), even with the extreme BH parameters. Compared with SGRBs, LGRBs need more massive accreted mass and more explosive energy, so we want to figure out whether the NDAFs could power LGRBs.

This paper is organized as follows. In Section 2, we introduce the NDAF and fireball models to connect the observational values and model parameters. Based on the collected LGRB data, we analyse the distributions of the isotropic radiation energy and jet kinetic energy and show the distribution of the accreted masses with different parameters of the BH-NDAF system in Section 3. The last section is devoted to the conclusions and discussion.

2 Model

The NDAF and fireball models for GRBs have been widely discussed for several decades. Fan & Wei (2011) proposed a method to estimate the mass accreted by the central BH of progenitors of SGRBs where it is assumed that the outflow is driven from the neutrino-antineutrino annihilation steaming from an underlying NDAF model. Liu et al. (2015c) developed this method further and enlarged the number of SGRB events to which the method is applied. Here we reuse the same methodology as in the former papers but applied to LGRBs. Generally, the mean output power from the central engine can be expressed as

 ˙E≈(1+z)(Eγ,iso+Ek,iso)θ2j2T90, (1)

where is the redshift, can roughly be considered as the duration of the activity of the central engine, is the isotropic kinetic energy estimated by the X-ray luminosity during the afterglow phase with the standard afterglow model, is the isotropic equivalent energy radiated in gamma-ray band, and is the opening angle of the ejecta.

Meanwhile, is a fraction of the total neutrino annihilation luminosity averaged in the whole accretion process, i.e.,

 ˙E=ηLν¯ν, (2)

where is a dimensionless number that represents the conversion efficiency the neutrino annihilation. (e.g., Eichler et al., 1989; Aloy et al., 2005; Liu et al., 2012b, 2015c). In this paper, we take = 0.2 and 0.5. The approximate analytical formula of neutrino annihilation luminosity we adopt here is proposed by Zalamea & Beloborodov (2011),

 Lν¯ν (erg s−1)≈5.7×1052x−4.8ms(MBH/M⊙)−3/2 ×{0˙M<˙Mign(˙M/M⊙ s−1)9/4˙Mign<˙M<˙Mtrap(˙Mtrap/M⊙ s−1)9/4˙M>˙Mtrap, (3)

where is the mass accretion rate, , is radius of the marginally stable orbit, and for . is the Schwarzschild radius, and and are the critical ignition accretion rate and the accretion rate when neutrino trapping event occurs in the inner region (e.g., Chen & Beloborodov, 2007; Zalamea & Beloborodov, 2011).

Hence, the mean accretion rate for the case can be defined as (e.g., Fan & Wei, 2011; Liu et al., 2015c)

 ˙M≈0.12 [(1+z)(Eγ,iso,51+Ek,iso,51)θ2jηT90,s]4/9 ×  x2.1ms (MBHM⊙)2/3 M⊙ s−1, (4)

where and . Of course, we must ensure that the NDAF can be ignited, i.e., . Then we can estimate the accreted mass , i.e., (e.g., Liu et al., 2015c)

 Macc≈0.12 [(Eγ,iso,51+Ek,iso,51)θ2jη]4/9(T90,s1+z)5/9 ×  x2.1ms (MBHM⊙)2/3 M⊙. (5)

Although the accreted masses can be estimated with the above equations and the observational data, it is still worth noting that there are some uncertainties. Firstly, The hyperaccretion may result in the violent evolution of the BH’s characteristics, which further leads to the evolution of the neutrino annihilation luminosity during the long duration. Fortunately, the extreme Kerr BHs and low mean accretion rates are considered in the central engines of LGRBs (e.g., Popham et al., 1999; Woosley & Heger, 2006). Recently, we found that this evolution has little effect on the total energy of the neutrino annihilation for LGRBs (Song et al., 2015). Secondly, we replace the interval of the central engine with in this work. In fact, the duration of a GRB should be longer than the time interval of the activity of central engine if the radial expansion of the fireball is considered, and this effect is potentially more important in SGRBs than in LGRBs (Aloy et al., 2005; Janka et al., 2006). Thirdly, we do not consider the outflow from the disc, which may seriously influences the measurement of the disc mass (Liu et al., 2012b; Janiuk et al., 2013). Therefore the accreted mass we calculated can be considered as a lower limit of mass of the original accretion matter.

3 Results

Since no LGRB has a positively identified origin of accretion or magnetar models, in order to test the ability of the NDAF model, all LGRBs with the data of , , , , and should be collected. There are some selection criterions should be stressed. Firstly, in the collapse events, we can consider that all Type Ib/c supernovae (SNe) have an engine-driven BH/NS-NDAF (e.g., Liu et al., 2015e). If the ejections from NDAFs direct to the observers, various kinds of GRB light curves should be detected. So the LGRBs associated with SNe should be included, which are marked in Table 1. Secondly, the ultra-long GRBs (ULGRBs) may derive from the blue supergiants (e.g., Nakauchi et al., 2013), which is different from the origin of “normal” LGRBs. For the disc model, they cannot be explained by neutrino annihilation mechanism but BZ mechanism (e.g., Nathanail & Contopoulos, 2015). So these sources are not included. Thirdly, GRB 060614 is a typical short-long GRBs, which may originate from a event of compact objects merger (e.g., Gehrels et al., 2006; Yang et al., 2015). It should not be considered.

Our sample contains 48 LGRBs as shown in Table 1. All data are taken from the literatures and Swift data archives. They are discovered and measured by the different telescopes including BeppoSAX, BATSE, HETE, Swift and Fermi. The isotropic radiated energy in the prompt emission phase , redshift , and the duration are directly available from measurements. and can be estimated by the X-ray afterglow phase with the standard afterglow model. Moreover, due to the restricted observations of LGRB afterglows, it is difficult to obtain the accurate jet opening angles. In Table 1, the lower limits of the jet opening angles in some LGRB data are shown, which cause the calculated accreted masses lower. In this sample, we collect the data of , , , and with errors because the results may depend on them.

Figure 1 shows the distributions of the accreted masses for the different typical mean BH masses and spins and conversion efficiencies, which are respectively set to = 3, 5, and 10, = 0.9, 0.95, and 0.998, and =0.2 and 0.5 corresponding to Figures 1 (a-i). It needs to be emphasized that the BHs in the centre of the collapsars should be rotating very rapidly, thus the spin parameters are set larger than 0.9 (e.g., Popham et al., 1999; Woosley & Heger, 2006). It is obvious that larger central BHs will require more accreted mass. Conversely, the accreted mass is negatively correlated with the mean BH spin and conversion efficiency. In these cases, the mean spin parameters are definitely more effective on the values of the accreted masses than the mean BH masses. For the cases with = 0.9, 0.95, and 0.998, , and = 0.2, about , , and LGRBs’ accreted masses are less than , respectively. If we consider = 0.998, , and =0.5, all the values of the accreted masses are less than . Obviously, if the accreted mass is larger than , the mean BH mass must be much larger than .

According to these results, the accreted mass of most LGRBs are less than the simulation results of the disc mass, (e.g., MacFadyen & Woosley, 1999; Popham et al., 1999; Zhang et al., 2003). Of course, there may still exist some LGRBs without the extreme Kerr low-mass BH in the centre or high efficiency, which need invoking alternative energy extraction mechanisms, such as BZ mechanism (e.g., Blandford & Znajek, 1977; Lee et al., 2000a, b; Yuan & Zhang, 2012; Kawanaka et al., 2013; Hou et al., 2014a; Liu et al., 2015b).

4 Conclusions and discussion

The goal of this work is to verify whether the NDAFs can account for the total energy of the prompt emission and afterglow phases in LGRBs. We collected data and calculated the distributions of the mean accreted masses for the reasonable states of the BH hyperaccretion system. As a result, most of LGRBs can be satisfied with the NDAFs except for some extremely high energy cases. Additionally, one should stress that the accretion rate and the BH mass and spin we adopted are the mean quantity throughout the accretion process. Recently, we studied the effects of the evolution of the BH-NDAF system on the neutrino annihilation energy (Song et al., 2015), which has slight influences on the current results.

As mentioned above, increasing the neutrino emission rate or enhancing annihilation efficiency are the effective ways to improve the ability of NDAF model. For example, Liu et al. (2015a) presented that the vertical convection can suppress the radial advection, which effectively increase the neutrino emission rate. Moreover, Liu et al. (2010) investigated the vertical structure of the NDAF model. We noticed that the half-opening angle of the disc is actually very large, . As mentioned in Birkl et al. (2007), the BH gravity can affect the traces of the neutrinos, especially for the neutrinos launched from the inner region of the disc, so they considered that the thin disc geometries is more efficient to raise the energy of the neutrino annihilation. But once the disc is extremely geometrically thick, the annihilable efficiency could be greatly enhanced due to the neutrinos trapped in a narrow space, even though some neutrinos falling into the BH.

BZ mechanism is another popular candidate for the energy sources of GRBs (e.g., Blandford & Znajek, 1977; Lee et al., 2000a, b). Moreover, Yuan & Zhang (2012) proposed that the episodic magnetic reconnection from the disc, which is very similar to the events at the solar surface, can produce GRBs and subsequent flares. Additionally, the closed magnetic field lines connecting a BH with its surrounding disc can transfer the angular momentum and the energy, which may significantly enhance the neutrino luminosity (e.g., Lei et al., 2009; Luo et al., 2013). It should be emphasized that magnetar model is more dynamic than ever since ULGRBs and super-luminous SNe were observed (e.g., Gao et al., 2015; Metzger et al., 2015; Wang & Dai, 2013; Wang et al., 2015a, b).

Furthermore, X-ray flares as the common features in GRB afterglows (e.g., Gehrels et al., 2004; Burrows et al., 2005; Chincarini et al., 2007) are not considered in this framework. They may have the similar origin with GRBs and ask for more extreme dynamic conditions and more massive accreted mass of the systems. Besides, the long ‘plateau’ phases (shallow decay phases) with more than thousands of seconds are generally observed in the X-ray afterglows (e.g., Zhang et al., 2006), which may require the energy injection from the central engine (e.g., Dai & Lu, 1998; Zhang et al., 2006). Thus the NDAF model has to confront more rigorous challenges.

Acknowledgements

We thank Bing Zhang, Ye Li, Shu-Jin Hou, and Hui-Jun Mu for beneficial discussion, and the anonymous referee for very useful suggestions and comments. This work was supported by the National Basic Research Program of China (973 Program) under grant 2014CB845800, the National Natural Science Foundation of China under grants 11274200, 11333004, 11373002, 11473022, 11573023, and U1331101, the Fundamental Research Funds for the Central Universities under grant 20720160024. TL acknowledges financial support from China Scholarship Council to work at UNLV.

References

• Aloy et al. (2005) Aloy, M. A., Janka, H.-T., & Müller, E. 2005, A&A, 436, 273
• Amati (2006) Amati, L. 2006, MNRAS, 372, 233
• Birkl et al. (2007) Birkl, R., Aloy, M. A., Janka, H.-T., & Müller, E. 2007, A&A, 463, 51
• Blandford & Znajek (1977) Blandford, R. D., & Znajek, R. L. 1977, MNRAS, 179, 433
• Bloom et al. (2003) Bloom, J. S., Frail, D. A., & Kulkarni, S. R. 2003, ApJ, 594, 674
• Burrows et al. (2005) Burrows, D. N., Romano, P., Falcone, A., et al. 2005, Science, 309, 1833
• Cenko et al. (2010) Cenko, S. B., Frail, D. A., Harrison, F. A., et al. 2010, ApJ, 711, 641
• Cenko et al. (2011) Cenko, S. B., Frail, D. A., Harrison, F. A., et al. 2011, ApJ, 732, 29
• Chen & Beloborodov (2007) Chen, W. X., & Beloborodov, A. M. 2007, ApJ, 657, 383
• Chincarini et al. (2007) Chincarini, G., Moretti, A., Romano, P., et al. 2007, ApJ, 671, 1903
• Dai (2004) Dai, Z. G. 2004, ApJ, 606, 1000
• Dai & Lu (1998) Dai, Z. G., & Lu, T. 1998, A&A, 333, L87
• Di Matteo et al. (2002) Di Matteo, T., Perna, R., & Narayan, R. 2002, ApJ, 579, 706
• Eichler et al. (1989) Eichler, D., Livio, M., Piran, T., & Schramm, D. N. 1989, Nature, 340, 126
• Fan & Wei (2011) Fan, Y.-Z., & Wei, D.-M. 2011, ApJ, 739, 47
• Gao et al. (2015) Gao, H., Ding, X., Wu, X.-F., Dai, Z.-G., & Zhang, B. 2015, ApJ, 807, 163
• Gehrels et al. (2004) Gehrels, N., Chincarini, G., Giommi, P., et al. 2004, ApJ, 611, 1005
• Gehrels et al. (2006) Gehrels, N., Norris, J. P., Barthelmy, S. D., et al. 2006, Nature, 444, 1044
• Ghirlanda et al. (2012) Ghirlanda, G., Nava, L., Ghisellini, G., et al. 2012, MNRAS, 420, 483
• Gu et al. (2006) Gu, W.-M., Liu, T., & Lu, J.-F. 2006, ApJ, 643, L87
• Hjorth& Bloom (2012) Hjorth, J., & Bloom, J. S. 2012, Chapter 9 in “Gamma-Ray Bursts”, Cambridge Astrophysics Series 51, eds. Kouveliotou, C., Wijers, R. A. M. J., and Woosley, S., Cambridge University Press (Cambridge), 169-190
• Hou et al. (2014a) Hou, S.-J., Gao, H., Liu, T., et al. 2014a, MNRAS, 441, 2375
• Hou et al. (2014b) Hou, S.-J., Liu, T., Gu, W.-M., et al. 2014b, ApJ, 781, L19
• Hurkett et al. (2006) Hurkett, C. P., Osborne, J. P., Page, K. L., et al. 2006, MNRAS, 368, 1101
• Janiuk et al. (2013) Janiuk, A., Mioduszewski, P., & Moscibrodzka, M. 2013, ApJ, 776, 105
• Janka et al. (2006) Janka, H.-T., Aloy, M.-A., Mazzali, P. A., & Pian, E. 2006, ApJ, 645, 1305
• Japelj et al. (2014) Japelj, J., Kopač, D., Kobayashi, S., et al. 2014, ApJ, 785, 84
• Kawanaka & Mineshige (2007) Kawanaka, N., & Mineshige, S. 2007, ApJ, 662, 1156
• Kawanaka et al. (2013) Kawanaka, N., Piran, T., & Krolik, J. H. 2013, ApJ, 766, 31
• Kluźniak & Lee (1998) Kluźniak, W., & Lee, W. H. 1998, ApJ, 494, L53
• Kohri & Mineshige (2002) Kohri, K., & Mineshige, S. 2002, ApJ, 577, 311
• Kohri et al. (2005) Kohri, K., Narayan, R., & Piran, T. 2005, ApJ, 629, 341
• Kouveliotou et al. (1993) Kouveliotou, C., Meegan, C. A., Fishman, G. J., et al. 1993, ApJ, 413, L101
• Kumar & Zhang (2015) Kumar, P., & Zhang, B. 2015, Phys. Rep., 561, 1
• Laskar et al. (2014) Laskar, T., Berger, E., Tanvir, N., et al. 2014, ApJ, 781, 1
• Laskar et al. (2015) Laskar, T., Berger, E., Margutti, R., et al. 2015, ApJ, 814, 1
• Lee et al. (2000a) Lee, H. K., Brown, G. E., & Wijers, R. A. M. J. 2000a, ApJ, 536, 416
• Lee et al. (2000b) Lee, H. K., Wijers, R. A. M. J., & Brown, G. E. 2000b, Phys. Rep., 325, 83
• Lee & Kluźniak (1999) Lee, W. H., & Kluźniak, W. 1999, ApJ, 526, 178
• Lei et al. (2009) Lei, W.-H., Wang, D.-X., Zhang, L., et al. 2009, ApJ, 700, 1970
• Li & Liu (2013) Li, A., & Liu, T. 2013, A&A, 555, A129
• Liang et al. (2008) Liang, E.-W., Racusin, J. L., Zhang, B., Zhang, B.-B., & Burrows, D. N. 2008, ApJ, 675, 528
• Liu et al. (2010) Liu, T., Gu, W.-M., Dai, Z.-G., & Lu, J.-F. 2010, ApJ, 709, 851
• Liu et al. (2015a) Liu, T., Gu, W.-M., Kawanaka, N., & Li, A. 2015a, ApJ, 805, 37
• Liu et al. (2007) Liu, T., Gu, W.-M., Xue, L., & Lu, J.-F. 2007, ApJ, 661, 1025
• Liu et al. (2012a) Liu, T., Gu, W.-M., Xue, L., & Lu, J.-F. 2012a, Ap&SS, 337, 711
• Liu et al. (2008) Liu, T., Gu, W.-M., Xue, L., Weng, S.-S., & Lu, J.-F. 2008, ApJ, 676, 545
• Liu et al. (2015b) Liu, T., Hou, S.-J., Xue, L., & Gu, W.-M. 2015b, ApJS, 218, 12
• Liu et al. (2015c) Liu, T., Lin, Y.-Q., Hou, S.-J., & Gu, W.-M. 2015c, ApJ, 806, 58
• Liu et al. (2012b) Liu, T., Liang, E.-W., Gu, W.-M., et al. 2012b, ApJ, 760, 63
• Liu et al. (2013) Liu, T., Xue, L., Gu, W.-M., & Lu, J.-F. 2013, ApJ, 762, 102
• Liu et al. (2015d) Liu, T., Xue, L., Zhao, X.-H., Zhang, F.-W., & Zhang, B. 2015d, arXiv:1511.02800
• Liu et al. (2015e) Liu, T., Zhang, B., Ma, R.-Y., & Xue, L. 2015e, arXiv:1512.07203
• Lloyd-Ronning & Zhang (2004) Lloyd-Ronning, N. M., & Zhang, B. 2004, ApJ, 613, 477
• Lu et al. (2012) Lu, R.-J., Wei, J.-J., Qin, S.-F., & Liang, E.-W. 2012, ApJ, 745, 168
• Luo et al. (2013) Luo, Y., Gu, W.-M., Liu, T., & Lu, J.-F. 2013, ApJ, 773, 142
• Lü & Zhang (2014) Lü, H.-J., & Zhang, B. 2014, ApJ, 785, 74
• MacFadyen & Woosley (1999) MacFadyen, A. I., & Woosley, S. E. 1999, ApJ, 524, 262
• Mao et al. (2010) Mao, Z., Yu, Y. W., Dai, Z. G., Pi, C. M., & Zheng, X. P. 2010, A&A, 518, A27
• Melandri et al. (2014) Melandri, A., Pian, E., D’Elia, V., et al. 2014, A&A, 567, A29
• Metzger et al. (2011) Metzger, B. D., Giannios, D., Thompson, T. A., Bucciantini, N., & Quataert, E. 2011, MNRAS, 413, 2031
• Metzger et al. (2015) Metzger, B. D.,Margalit, B., Kasen, D., & Quataert, E. 2015, MNRAS, 454, 3311
• Metzger et al. (2008) Metzger, B. D., Quataert, E., & Thompson, T. A. 2008, MNRAS, 385, 1455
• Nakar (2007) Nakar, E. 2007, Phys. Rep., 442, 166
• Nakauchi et al. (2013) Nakauchi, D., Kashiyama, K., Suwa, Y., & Nakamura, T. 2013, ApJ, 778, 67
• Narayan et al. (1992) Narayan, R., Paczyński, B., & Piran, T. 1992, ApJ, 395, L83
• Narayan et al. (2001) Narayan, R., Piran, T., & Kumar, P. 2001, ApJ, 557, 949
• Nathanail & Contopoulos (2015) Nathanail, A., & Contopoulos, I. 2015, MNRAS, 453, L1
• Nava et al. (2012) Nava, L., Salvaterra, R., Ghirlanda, G., et al. 2012, MNRAS, 421, 1256
• Popham et al. (1999) Popham, R., Woosley, S. E., & Fryer, C. 1999, ApJ, 518, 356
• Romano et al. (2006) Romano, P., Campana, S., Chincarini, G., et al. 2006, A&A, 456, 917
• Romero et al. (2010) Romero, G. E., Reynoso, M. M., & Christiansen, H. R. 2010, A&A, 524, A4
• Sakamoto et al. (2008) Sakamoto, T., Barthelmy, S. D., Barbier, L., et al. 2008, ApJS, 175, 179
• Sakamoto et al. (2011) Sakamoto, T., Barthelmy, S. D., Baumgartner, W. H., et al. 2011, ApJS, 195, 2
• Song et al. (2015) Song, C.-Y., Liu, T., Gu, W.-M., et al. 2015, ApJ, 815, 54
• Sun et al. (2012) Sun, M.-Y., Liu, T., Gu, W.-M., & Lu, J.-F. 2012, ApJ, 752, 31
• Suwa & Murase (2009) Suwa, Y., & Murase, K. 2009, Phys. Rev. D, 80, 123008
• Thompson et al. (2004) Thompson, T. A., Chang, P., & Quataert, E. 2004, ApJ, 611, 380
• Usov (1992) Usov, V. V. 1992, Nature, 357, 472
• Wang & Dai (2013) Wang, L.-J., & Dai, Z.-G. 2013, ApJ, 774, L33
• Wang et al. (2015a) Wang, S. Q., Liu, L. D., Dai, Z. G., Wang, L. J., & Wu, X. F. 2015a, arXiv:1509.05543
• Wang et al. (2015) Wang, X.-G., Zhang, B., Liang, E.-W., et al. 2015, ApJS, 219, 9
• Wang et al. (2015b) Wang, S. Q., Wang, L. J., Dai, Z. G., & Wu, X. F. 2015b, ApJ, 799, 107
• Wheeler et al. (2000) Wheeler, J. C., Yi, I., Höflich, P., & Wang, L. 2000, ApJ, 537, 810
• Woosley (1993) Woosley, S. E. 1993, ApJ, 405, 273
• Woosley & Heger (2006) Woosley, S. E., & Heger, A. 2006, ApJ, 637, 914
• Xue et al. (2013) Xue, L., Liu, T., Gu, W.-M., & Lu, J.-F. 2013, ApJS, 207, 23
• Yang et al. (2015) Yang, B., Jin, Z.-P., Li, X., et al. 2015, Nature Communications, 6, 7323
• Yuan & Zhang (2012) Yuan, F., & Zhang, B. 2012, ApJ, 757, 56
• Zalamea & Beloborodov (2011) Zalamea, I., & Beloborodov, A. M. 2011, MNRAS, 410, 2302
• Zhang et al. (2006) Zhang, B., Fan, Y. Z., Dyks, J., et al. 2006, ApJ, 642, 354
• Zhang et al. (2012) Zhang, F.-W., Shao, L., Yan, J.-Z., & Wei, D.-M. 2012, ApJ, 750, 88
• Zhang et al. (2003) Zhang, W., Woosley, S. E., & MacFadyen, A. I. 2003, ApJ, 586, 356
• Zhang et al. (2011) Zhang, Z., Zhao, Y., & Zhang, Y. 2011, Ap&SS, 333, 89
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