• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Probing the Progenitor of High-z Short-duration GRB 201221D and its Possible Bulk Acceleration in Prompt Emission

    2022-08-01 01:47:58HaoYuYuanHouJunYeLiBinBinZhangHuiSunJaredRiceJunYangandEnWeiLiang

    Hao-Yu Yuan, Hou-Jun Lü, Ye Li, Bin-Bin Zhang, Hui Sun, Jared Rice, Jun Yang, and En-Wei Liang

    1 Guangxi Key Laboratory for Relativistic Astrophysics, Department of Physics, Guangxi University, Nanning 530004, China; lhj@gxu.edu.cn

    2 Purple Mountain Observatory, Chinese Academy of Sciences, Nanjing 210008, China

    3 Key Laboratory of Modern Astronomy and Astrophysics (Nanjing University), Ministry of Education, Nanjing 210093, China

    4 Key Laboratory of Space Astronomy and Technology, National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100101, China

    5 Department of Mathematics and Physical Science, Southwestern Adventist University, Keene, TX 76059, USA

    Received 2022 February 22; revised 2022 May 9; accepted 2022 May 16; published 2022 June 14

    Abstract The growing observed evidence shows that the long-and short-duration gamma-ray bursts(GRBs)originate from massive star core-collapse and the merger of compact stars, respectively. GRB 201221D is a short-duration GRB lasting ~0.1 s without extended emission at high redshift z=1.046. By analyzing data observed with the Swift/BAT and Fermi/GBM, we find that a cutoff power-law model can adequately fit the spectrum with a soft Ep = keV,and isotropic energyEγ,iso=1.36×10 51erg.In order to reveal the possible physical origin of GRB 201221D, we adopted multi-wavelength criteria (e.g., Amati relation, ε-parameter, amplitude parameter,local event rate density, luminosity function, and properties of the host galaxy), and find that most of the observations of GRB 201221D favor a compact star merger origin.Moreover,we find that α?is larger than2+β? in the prompt emission phase which suggests that the emission region is possibly undergoing acceleration during the prompt emission phase with a Poynting-flux-dominated jet.

    Key words: (stars:) gamma-ray burst: individual (GRB 201221D) – stars: massive – acceleration of particles

    1. Introduction

    Observationally, supernovae (SNe) are associated with some long GRBs (or without short GRBs; Galama et al.1998; Hjorth et al. 2003; Soderberg et al. 2004; Campana et al. 2006) and the host galaxies of long (or short) GRBs are typically associated with irregular galaxies with intense (or little) star formation (Tanvir et al. 2005; Fruchter et al. 2006).These lines of observational evidence, as well as the joint detection of the gravitational-wave event GW170817 and the short GRB 170817A (Abbott et al. 2017; Goldstein et al.2017; Savchenko et al. 2017; Zhang et al. 2018b), encourage people to believe that long and short GRBs are likely related to the deaths of massive stars (Type II) and the merger of two compact stellar objects (Type I), respectively (Eichler et al. 1989; Woosley 1993; Zhang 2006). However, some apparently long-duration (or short-duration high-z) GRBs have been suggested to originate in compact stars mergers (or massive star core-collapse). Two counterexamples are GRB 060614 (Gehrels et al. 2006), and GRB 090426 (Levesque et al. 2010;Xin et al. 2011). In particular, Zhang et al.(2021)recently discovered another peculiarly short-duration GRB 200826A which seems to originate in massive star corecollapse.

    In either the death of massive star or merger of compact stars, the catastrophic event leaves behind a hyper-accreting black hole or a rapidly rotating highly magnetized neutron star(called a magnetar), which serves as the central engine of a collimated outflow (or jet) with a relativistic speed toward Earth(Usov 1992;Dai&Lu 1998b;Zhang&Mészáros 2001;Zhang 2011; Lü & Zhang 2014; Kumar & Zhang 2015; Lü et al.2015;Chen et al.2017).One basic question is what is the composition of the relativistic jets? There are two models widely discussed in the literature (Lei et al. 2013). One is the fireball model with a matter-dominated outflow which dissipates its kinetic energy in internal shocks or external shocks to produce the observed GRB emission (Rees &Meszaros 1992, 1994; Meszaros et al. 1993; Kobayashi et al.1997). Within this model, the fireball has a rapid acceleration early on and can only reduce its kinetic energy at large radii from the central engine. The other one is Poynting-fluxdominated outflow. Within this scenario, the Poynting flux energy can be converted to kinetic energy(Drenkhahn&Spruit 2002; Komissarov et al. 2009), and then converted to particle energy and radiation via magnetic dissipation, such as reconnection, current instabilities, and internal collisions(Zhang & Yan 2011). In comparison to the fireball model, a Poynting-flux-dominated jet can undergo gradual acceleration in a large range of emission region (Gao & Zhang 2015; Uhm& Zhang 2015).

    Traditionally, it is assumed that the curvature effect6The curvature effect is due to the observer receiving progressively delayed emission from higher latitudes (Liang et al. 2006; Zhang et al. 2007).can be used to interpret the pulse decay(including both prompt emission and X-ray flare) if the emission region moves with a constant Lorentz factor. One can measure the decay index (α?) and the spectral index () during the pulse decay, and they satisfy a simple relationship7Throughout the paper, the notationfv (t )∝t- α?v -β?is adopted.in the Lab frame (Kumar &Panaitescu 2000).However,it is difficult to interpret the observed spectral lags of pulses with the curvature effect(Uhm&Zhang 2016b).Indeed,Uhm&Zhang(2016b)found that the dissipation of magnetic field energy in the shell via reconnection of magnetic field lines would result in a faster decrease than indicated by flux conservation if the magnetic field strength in the emitting region decreases with radius.This means that the emission region does not have a constant Lorentz factor but is accelerated, and the decay slopeα?should be steeper than2+β?(Uhm & Zhang 2015). Afterwards, evidence for rapid bulk acceleration was discovered in observations of both X-ray flares(Jia et al. 2016;Uhm&Zhang 2016a)and GRB prompt emission(Uhm&Zhang 2016b;Li&Zhang 2021).In fact,Jia et al.(2016)found that a large fraction of X-ray flares in GRBs are inconsistent with the predicted relation from the curvature effect.Li&Zhang(2021)invoked the same method to analyze the prompt emission lightcurves of single-pulse GRBs, and suggested that the emission region of at least some GRBs is undergoing acceleration during the prompt emission phase.

    A bright short-duration GRB 201221D, triggered the Swift Burst Alert Telescope at 23:06:34 UT on 2020 December 21(BAT; Krimm et al. 2020) and located the source at R.A.=11h24m12sand=+ ° ′ ″decl.42 08 39 (J2000).This GRB was also detected by the Fermi Gamma-ray Burst Monitor (GBM;Hamburg et al. 2020) and Konus-Wind (Frederiks et al. 2020).Based on spectroscopy of the optical counterpart, it was measured at a redshift of 1.046 (de Ugarte Postigo et al. 2020)which is larger than 95% of short-duration GRBs without extended emission (Dichiara et al. 2020). In order to test the physical origin of the high-z short-duration GRB 201221D, we perform a comprehensive analysis of Fermi and Swift data on this burst as shown in Section 2.In Section 3,we compare some statistical relations of this burst with those of other long and short GRBs,and discuss its physical origin. In Section 4, we discuss the possible bulk acceleration in the prompt emission of this burst. A summary and conclusions are presented in Section 5.Throughout the paper, a concordance cosmology with parameters H0=71 km s?1Mpc?1, ΩM=0.30, and ΩΛ=0.70 is adopted.

    2. Data Reduction and Analysis

    2.1. Swift Data Reduction

    GRB 201221D first triggered the Swift/BAT at 23:06:34 UT on 2020 December 21 (Page et al. 2020). The BAT data were processed using the HEASOFT package (v6.28). The light curves in different energy bands and spectra were extracted by running batbinevt (Sakamoto et al. 2008). The time bin size is fixed to 64 ms in this case due to the short duration, and the light curve shows a short-pulse with duration T90=0.15±0.04 s in the 15–350 keV (see Figure 1). The time-averaged spectrum from T0?0.06 to T0+0.17 s is best fit by a simple power-law model with spectral index 1.56±0.13 due to the narrow energy band.Moreover,we do not find any signature of extended emission even up to 100 s following the burst. The X-ray Telescope (XRT) began observing the field at 87 s after the BAT trigger, but the source is too faint to be detected with photon counting (Evans et al. 2020).

    2.2. Fermi Data Reduction

    We downloaded the corresponding Time-Tagged-Event data from the public data site for Fermi/GBM.8https://heasarc.gsfc.nasa.gov/FTP/fermi/data/gbm/daily/For more details on the light curve and spectra data reduction procedure see the discussion in Zhang et al. (2016). The light curves of the n8 and b1 detectors are shown in Figure 1,and the background is modeled via applying the “baseline” method (Zhang et al.2011) to a wide time interval before and after the signal and subtracting the GBM light curve. The lightcurves show a single-pulse emission with a duration of T90=0.13±0.01 in 50?300 keV. No significant signatures of precursor emission before the burst and extended emission(EE)after the burst are found in the GBM temporal analysis.

    Figure 1. Swift/BAT and Fermi/GBM light curves of GRB 201221D in different energy bands with a 64 ms time bin.

    We also extract both time-integrated and time-dependent spectral analyses of GRB 201221D between T0?0.03 and T0+0.1. This time interval is divided into four slices (see Table 1) based on brightness and the count statistical significance of the spectral fitting (Zhang et al. 2016, 2018a).The background spectra are extracted from the time intervals before and after the prompt emission phase and modeled with an empirical function (Zhang et al. 2011), and the spectral fitting is performed by using our automatic code “McSpecfit”in Zhang et al.(2018a).Several spectral models can be selected to test the spectral fitting of the burst,such as power-law (PL),cutoff power-law(CPL),Band function(Band),and Blackbody(BB).In order to test which model is the best fit of the data,weinvoke the Bayesian Information Criteria (BIC)9BIC is a criterion for model selection among a finite set of models,and it is defined asBIC=χ2+ k· ln(n ),where k is the number of model parameters,and n is the number of data points.The model with the lowest BIC is preferred.to judge the best model among different models. The comparison of the goodness of the fits for different models is shown in Table 2.We find that the CPL model is the best one for adequately describing the observed data. The CPL model fit of the timeintegrated spectrum is shown in Figure 2 for parameter constraints of the fit. It gives peak energyEp=113keV,and a lower energy spectral index ofΓ = -0.26.The bestfit parameters of the CPL fits are listed in Table 3. The CPL model is expressed as

    Table 1 Properties of GRB 201221D

    Table 2 BIC Values for Different Models We Adopted to Fit within Time-dependent Spectral Fitting

    where Γ and N0are the photon index and the CPL spectral fitting normalization, respectively. To extract the time-dependent spectrum, we use a similar method to the one mentioned above. We find that the CPL model is also the best fit,and the fitting results are shown in Table 3. One can see that the tracking spectral evolution is observed during the burst (see Figure 3).

    2.3. Host Galaxy and Other High-z Short GRBs

    GRB 201221D was initially localized to R.A.=11h24m14 19, decl. = +42d08′35 5 with 3 9 uncertainty by Swift/XRT(Evans et al.2020).Later,with ther′ band afterglow from GTC/OSIRIS, it was localized to R.A.=11h24m14 09,decl.=+42d08′40 0 with 1″ uncertainty (Agüí Fernández et al.2021).A faint galaxy around it was identified as the host galaxy.Kilpatrick et al. (2020)analyzed stacked images of the Pan-STARRS data release (Flewelling & Alatalo 2016), and found the galaxy to have a g-band magnitude g=23.2±0.2 mag. This source was also observed by the Nordic Optical Telescope with r=23.1±0.3 mag (Malesani & Knudstrup 2020) and the Lowell Discovery Telescope with r=23.9 mag(Dichiara et al. 2020) which is consistent with emission identified in Kilpatrick et al.(2020).Recently,Agüí Fernández et al.(2021)observed the host and afterglow of GRB 201221D with GTC, and measured the redshift z=1.045±0.0008 which is consistent with the GCN report. In addition, it is detected in the Dark Energy Spectroscopic Instrument Legacy Survey (DESI/LS), and listed in the LS DR8 catalog. The DESI/LS g-band image with the localization from Swift/XRT(white dashed circle), the optical afterglow (green circle and dot), and the host galaxy (red cross) are presented in Figure 4.We analyzed the images and the host galaxy properties as follows.

    Figure 2.The parameter constraints of the spectral fit with the CPL model for GRB 201221D.Histograms and contours in the corner plots show the likelihood map of constrained parameters by using our McSpecFit package. Red crosses are the best-fitting values, and pink, yellow, and green circles are the 1σ, 2σ, and 3σ uncertainties, respectively.

    de Ugarte Postigo et al.(2020)obtained spectroscopy of the optical counterpart of GRB 201221D with the 10.4 m Telescope, and measured the redshift at z=1.046 based on a prominent emission feature. The high-z short GRBs play an important role in understanding the age of stellar progenitors,the cosmic chemical evolution, and formation channels of binary systems if we believe that short GRBs originate in the mergers of compact stars (Zheng & Ramirez-Ruiz 2007;Dominik et al. 2012; Anand et al. 2018). To date, more than 130 short GRBs have been detected by Swift/BAT,but less than 5% are found at z>1 (Dichiara et al. 2021).Dichiara et al. (2021) studied the properties of high-z short GRBs with z>1. They find that there are eight short GRBs with redshift z>1, with five short GRBs having EE and three short GRBs (GRBs 090426, 111117A, and 121226A) without EE. However, the short-duration GRB 090426 with z=2.609 seems to originate in a massive star core-collapse based on the properties of host galaxy and afterglow (Antonelli et al. 2009;Th?ne et al. 2011; Xin et al. 2011).

    2.4. Burst Energy

    Figure 3.Zooming in to the GBM lightcurve of the pulse with photon and energy flux in time interval[T0 ?0.4,T0+0.5],FRED fit(red solid line),and power-law fit (red dashed line), respectively (Top panel). The evolution of photon index (middle panel) and peak energy (bottom panel) with CPL model is presented.

    Table 3 Time-dependent Spectral Fitting Results of GRB 201221D with CPL Model

    3. What is the Physical Origin of GRB 201221D?

    The classification of GRBs remains an open question(Zhang 2011).Our purpose is to investigate the physical origin of GRB 201221D. In this section, we will discuss the origin of GRB 201221D by comparing the properties of GRB 201221D with those of other long- and short-duration GRBs, e.g., the Amati relation (Amati et al. 2002), luminosity function, properties of host galaxy, “tip of the iceberg” effect (Lü et al. 2014), and εclassification method (Lü et al. 2010).

    3.1.Comparisons of the Empirical Relationships of GRB 201221D with other Type I/II GRBs

    Figure 4. The DESI/LS g-band image with the localization of Swift/XRT (white dashed circle), optical afterglow (green circle and dot) and the host galaxy (red cross) of GRB 201221D.

    3.1.1. Ep–Eγ,isoRelation at the outlier of 3σ uncertainty of fits for both Type I and Type II GRBs, but seems to be closer to the distribution of Type I GRBs in comparison to Type II GRBs. It is difficult to judge the progenitor based only on this empirical correlation.

    3.1.2. Local Event Rate Density

    The event rate density describes how many events happen per volume per unit of time.The observed event rate density is redshift-dependent,luminosity-dependent,and beaming-dependent.Most short-duration GRBs are believed to originate in the merger of compact stars, and the observed local event rate density has a large uncertainty. Sun et al. (2015) estimated the local event rate density of short GRBs with ~(0.5–3)Gpc?3yr?1above 1050erg s?1. We estimate the local event rate density of the source fromρ0=1(VmaxT), whereVmaxis the maximum volume from which the source can be detected weighted by the redshift evolution (see Equation (3) in Sun et al.2015)and T is the total exposure time of the telescope or survey. Assuming a flux limit of Fth=10?7erg cm?2s?1, and total operation time of 16 yr for Swift, the local event rate density of GRB 201221D is ~3.0×10-3Gpc-3yr-1for the peak bolometric luminosity of ~2.09×1052erg s?1in 1–104keV. The 1σ errors are derived from Gehrels (1986)based on one detection. This value is lower than that of other Type II GRBs within the same luminosity range.A comparison of the event rate density of GRB 201221D with other Type I/II GRBs is shown in Figure 6.Moreover,we do not consider the beaming factor effect in our calculations due to uncertainty of jet opening-angle in this case.12By considering the effect of beaming factor,one can roughly estimate local event rate densities which are 25 and 500 times greater than those ignoring the beaming factor,respectively when we adopt the typical beaming factors ~0.04 and ~0.002 for short- and long-duration GRBs cases (Frail et al. 2001; Fong et al. 2015). If this is the case, the local event rate density is nearly consistent with other Type I and II GRBs within the same luminosity range.

    3.1.3. Luminosity Function

    The luminosity function of high-luminosity (HL) longduration GRBs can be characterized as a broken power law function via a large enough sample of redshift measurements,but the luminosity function of low-luminosity (LL) longduration GRBs is not well constrained with a small sample size due to the detectors’ sensitivity limit (Liang et al. 2007;Virgili et al. 2009; Sun et al. 2015). In comparison with longduration GRBs, the luminosity function of short-duration GRBs is less well constrained with a small fraction of the data. Sun et al. (2015) found that the luminosity function of short-duration GRBs can be roughly fitted with a simple power law in the luminosity range of [7×1049, 1050] erg s?1by assuming that all short-duration GRBs have a compact star merger origin.GRB 201221D broadly follow the distributions of both the long and short GRB populations, as shown in the right panel of Figure 6. Moreover, Paul (2018) analyzes the luminosity function with a large sample of short-duration GRBs observed by CGRO/BATSE, Swift/BAT, and Fermi/GBM. They found that the luminosity function can be described with the exponential cutoff power law and broken power law models. However, the redshifts of most short GRBs in their samples are not measured, and they adopt pseudo-redshifts which are derived by the empirical relationship. Therefore, we still used the results from Sun et al. (2015).

    Figure 5. Top panel: Ep and Eiso correlation diagram. Black points and gray diamonds correspond to Type I and Type II GRBs,respectively.The red star is GRB 201221D,and other data are taken from Zhang et al.(2009).The best-fit Ep–Eiso correlations for both Type II(gray diamonds)and Type I(black points)GRBs are plotted (solid lines) with the 3σ boundary (dashed line) marked.Bottom panel: Ep distribution of GRB 201221D and other Type I/II GRBs observed by Fermi/GBM. The data of Ep values of other Type I/II GRBs are taken from Lu et al. (2017) and von Kienlin et al. (2020), respectively.

    3.1.4. ε-parameter

    3.1.5. Amplitude Parameter

    Figure 6.Local event rate density(top)and normalized luminosity function(bottom)distribution for GRB 201221D(red),LL-LGRBs(blue),HL-LGRBs(olive),and short-duration GRBs (black) inferred Sun et al. (2015).

    Figure 7. (a)1D and 2D distributions of GRBs samples in T90–ε space.The dotted line is ε=0.03.(b) f(feff)–T90 for type I and Type II GRBs which are from Lü et al. (2014). The vertical solid line is T90=2 s. Blue (black) solid circles represent the Type I(II) GRB candidates, and green triangles denote the nearby lowluminosity long GRBs. The red star is GRB 201221D.

    3.1.6. Host Galaxy Properties—the Probability to Be a Type I GRB—log Odds

    3.2. Merger of Compact Star Scenario

    A leading progenitor for producing short-duration GRBs is binary neutron star mergers which were confirmed via the detection of the associated of gravitational-wave (GW)170817 and GRB 170817A (Abbott et al. 2017; Zhang et al. 2018b). Neutron stars likely receive a kick at birth from SNe explosions, and mergers of binary neutron stars are usually found to have a large offset from the galaxy center or even be outside of the host galaxy (Bloom et al.1999; Li & Zhang 2020). As a caution, Mandhai et al. (2021)point out that the observational offset is the minimum value due to projection effects, namely the real-offset is much larger. Most observations (amplitude parameter, properties of host galaxy, local event rate) support a compact star merger origin.

    Moreover,black hole and white dwarf merger systems were also proposed as possible progenitors of short-duration GRBs(Fryer et al.1999),but these systems may not be favorable for launching a GRB jet with a small fraction of accreted mass(Narayan et al. 2001). MacFadyen et al. (2005) proposed that accretion-induced collapse of a neutron star can power a GRB,but this model can produce GRBs in non-star-forming galaxies with a small offset of GRB location. It is also inconsistent with the observations of the host galaxy of GRB 201221D.White dwarf-white dwarf mergers producing a short GRB were discussed by Lyutikov & Toonen (2017), but this model predicts the powering of long-lasting EE after the prompt emission, which is disfavored for GRB 201221D because there is no observed EE.

    3.3. Massive Star Core-collapse Scenario

    Figure 9.One-and two-dimensional projections of the posterior distributions of the Equation(2)free parameters by MCMC fit with pulse.Vertical dashed lines mark the the median and 1σ range. The contours are drawn at 68%, 95%, and 99% credible levels.

    Figure 10. Similar with Figure 9, but adopting a power-law model to fit the decay segment of pulse.

    Such a puzzle might be solved in three ways if we believe that GRB 201221D originated from massive star core collapse.One is that GRB 201221D is produced by the collapse of a supramassive NS (or magnetar) into a black hole (Rezzolla &Kumar 2015),and the supramassive NS is the initial remnant of massive star collapse due to the stiff equation of state of NS,but does not power the GRB during this process. If this is the case, it should leave clues in the X-ray band when it loses rotational energy via dipole radiation. Unfortunately, there is not enough X-ray observational data after the prompt emission of GRB 201221D. The second possibility is that the total central engine timescale(teng)is much longer than that of what we observed, but the majority of its time (tjet) is used in breaking the envelope out of the stellar surface,so the observed duration time of GRB (tGRB=teng?tjet) could be as short as 1 s(Bromberg et al.2012).However,this timescale is still one order of magnitude longer than the observed duration of GRB 201221D. The third possibility is that GRB 201221D is not intrinsically short but is actually long duration due to the tip of the iceberg effect (Lü et al. 2014). Within this scenario, the observed fluence is underestimated, which means that the intrinsic isotropic energy Eγ,isoshould be larger than the current value what we calculated in GRB 201221D.If this is the case,the location of GRB 201221D in both Ep–Eγ,isoand εparameter diagrams will move into the long GRB population.In that case, it is natural to explain the observed empirical correlations. However, the f-parameter of GRB 201221D should be consistent with that of other Type II GRBs, but the observed f-parameter of GRB 201221D does not support this hypothesis.

    4. Possible Bulk Acceleration in Prompt Emission?

    Uhm & Zhang (2016b) proposed that the relationship betweenα?and2+β?is not consistent with the predicted of curvature effect when the emission region itself is undergoing acceleration or deceleration in the prompt emission phase. In this section, we describe how to test whether the prompt emission of GRB 201221D is undergoing acceleration or deceleration to diagnose the composition of the jet.

    As the light curve of GRB 201221D shows a single-peaked structure,we adopt a fast-rising and exponential decay(FRED)empirical function to fit in order to describe the shape of a pulse(Kocevski et al. 2003). It reads as,

    Figure 11.α? –β?relation of the curvature effect during the initial decay phase of prompt emission. The solid line is = 2+β?.

    5. Conclusions

    Combined with the above analysis, the origin of GRB 201221D is favored to be a binary neutron star merger. This model can be used to interpret some observations of this case,except for the ε-parameter and the inconclusiveness of the Amati relation.If this is the case,the remnant of such a merger should leave either a magnetar or black hole which produces a short-duration GRB with a Poynting-flux-dominated jet.It also naturally explains undergoing bulk acceleration in the prompt emission phase. Moreover, another possibility is that GRB 201221D originates in a black hole–neutron star (BH–NS)merger at high redshift. Mandhai et al. (2021) pointed out that the probabilities of BH–NS systems for producing GRBs peak at higher redshifts (e.g., z>1), although production of shortduration GRBs is rarer with BH–NS mergers (Shibata &Taniguchi 2011).

    On the other hand, the high-z short-duration GRBs without extended emission are very rare,with only five GRBs(060801,090426, 111117A, 121226A, and 190627A) observed so far.Within these five GRBs,GRB 111117A has been discussed as having a high probability of originating in a compact binary merger(Dichiara et al.2021).Moreover,Wiggins et al.(2018)estimated the rate for short GRBs peaks in the redshift range z=0.6–1 by studying the population synthesis, and this is similar to observations of short GRBs with a known redshift distribution(Fong et al.2015).Also,we expect more and more short GRBs with high-redshift to be observed in the future,and they can be used to study the star formation delay time(Wanderman&Piran 2015).If the binary neutron star mergers indeed occurred in high-z region of early universe, there are important implications for understanding binary stellar evolution, heavy element nucleosynthesis, and chemical evolution,or the possibility of the high-z short-duration GRBs belonging to a different population of bursts (e.g., NS–BH or WD–WD merger systems). However, owing to the relative lack of expected mergers in the high-z region, to understand those implications remains an open question. We hope that there are more neutron star mergers at high redshifts than that expected in the future. The third generation of GW detectors, such as Cosmic Explorer(Reitze et al.2019)or the Einstein Telescope(Punturo et al. 2010), may play a crucial role in understanding the mergers of these objects out to redshift z ~2–3.

    Acknowledgments

    We acknowledge the use of the public data from the Swift and Fermi data archive.This work is supported by the National Natural Science Foundation of China (grant Nos. 11922301 and 12133003), the Guangxi Science Foundation (grant Nos.2017GXNSFFA198008 and AD17129006). LHJ thanks support by the Program of Bagui Young Scholars Program, and special funding for Guangxi distinguished professors (Bagui Yingcai & Bagui Xuezhe). B.B.Z. acknowledges support by the National Key Research and Development Programs of China (2018YFA0404204), the National Natural Science Foundation of China (grant Nos. 11833003 and U2038105),and the Program for Innovative Talents, Entrepreneur in Jiangsu. H.S. acknowledges support by the Strategic Priority Research Program of the Chinese Academy of Sciences (grant Nos. XDA15310300, XDA15052100 and XDB23040000).Y.L. acknowledges support by the National Natural Science Foundation of China (grant Nos. 12041306 and 12103089),and the Natural Science Foundation of Jiangsu Province(grant No. BK20211000).

    在线观看美女被高潮喷水网站| 久久久久精品国产欧美久久久| 欧美xxxx黑人xx丫x性爽| 国产一区二区三区在线臀色熟女| 22中文网久久字幕| 欧美极品一区二区三区四区| 欧美xxxx黑人xx丫x性爽| 国产精品爽爽va在线观看网站| 日本色播在线视频| 日本黄色视频三级网站网址| 最新在线观看一区二区三区| 99久久精品国产国产毛片| 热99在线观看视频| 精品一区二区三区人妻视频| 人人妻人人看人人澡| 精品午夜福利视频在线观看一区| 免费人成视频x8x8入口观看| 国产一区亚洲一区在线观看| 日韩人妻高清精品专区| 91狼人影院| 色尼玛亚洲综合影院| 99久国产av精品国产电影| 国产一区二区三区在线臀色熟女| 最新中文字幕久久久久| 99视频精品全部免费 在线| 联通29元200g的流量卡| aaaaa片日本免费| 男女视频在线观看网站免费| 国产欧美日韩精品亚洲av| 国产午夜精品久久久久久一区二区三区 | av免费在线看不卡| 97碰自拍视频| 三级经典国产精品| 99久国产av精品| 69人妻影院| 中文字幕熟女人妻在线| 国产精品精品国产色婷婷| or卡值多少钱| ponron亚洲| 伦理电影大哥的女人| 午夜福利成人在线免费观看| 欧美丝袜亚洲另类| 色播亚洲综合网| 日本成人三级电影网站| 淫妇啪啪啪对白视频| 免费看av在线观看网站| 精品久久久久久成人av| 亚洲精品乱码久久久v下载方式| 一卡2卡三卡四卡精品乱码亚洲| 男女之事视频高清在线观看| 美女免费视频网站| 精品人妻视频免费看| 波多野结衣高清无吗| 亚洲经典国产精华液单| 在线播放无遮挡| 亚洲av中文字字幕乱码综合| 久久久国产成人免费| 天美传媒精品一区二区| 我要搜黄色片| 噜噜噜噜噜久久久久久91| 久久中文看片网| 我的女老师完整版在线观看| 在线免费观看的www视频| 久久精品影院6| 最近视频中文字幕2019在线8| 亚洲精品在线观看二区| 亚洲丝袜综合中文字幕| 亚洲欧美日韩高清在线视频| 啦啦啦韩国在线观看视频| 亚洲欧美日韩高清在线视频| 亚洲欧美日韩高清专用| 美女xxoo啪啪120秒动态图| 国产精品三级大全| 在线免费观看的www视频| 久久这里只有精品中国| 少妇被粗大猛烈的视频| 久久欧美精品欧美久久欧美| 不卡一级毛片| 欧美区成人在线视频| 悠悠久久av| 日本黄色视频三级网站网址| 亚洲欧美日韩东京热| 国产精品电影一区二区三区| 欧美另类亚洲清纯唯美| 99热6这里只有精品| 午夜福利高清视频| 日韩精品中文字幕看吧| 国产精品人妻久久久久久| 久久久久久久亚洲中文字幕| 欧美成人免费av一区二区三区| 国产精品美女特级片免费视频播放器| 一个人观看的视频www高清免费观看| 久久久精品94久久精品| 精品乱码久久久久久99久播| 97人妻精品一区二区三区麻豆| 99久久久亚洲精品蜜臀av| 人人妻人人澡人人爽人人夜夜 | 国产精品爽爽va在线观看网站| 女生性感内裤真人,穿戴方法视频| 色哟哟·www| 人人妻人人澡欧美一区二区| .国产精品久久| 看十八女毛片水多多多| 国产精品亚洲一级av第二区| 中文字幕av成人在线电影| 日产精品乱码卡一卡2卡三| 国产亚洲精品久久久久久毛片| 精品一区二区三区视频在线观看免费| 久久热精品热| 少妇丰满av| 黑人高潮一二区| 国产精品一区二区三区四区久久| 又粗又爽又猛毛片免费看| 国产亚洲精品久久久久久毛片| 看十八女毛片水多多多| 国产午夜福利久久久久久| 日韩精品青青久久久久久| 亚洲成人中文字幕在线播放| 成人高潮视频无遮挡免费网站| 亚洲欧美成人精品一区二区| 日韩欧美在线乱码| 日韩在线高清观看一区二区三区| 人妻丰满熟妇av一区二区三区| 欧美激情国产日韩精品一区| 三级男女做爰猛烈吃奶摸视频| 久久精品人妻少妇| 在线观看av片永久免费下载| 色在线成人网| 久久久成人免费电影| 丰满人妻一区二区三区视频av| 色吧在线观看| 色综合站精品国产| 国产黄色小视频在线观看| 亚洲av熟女| 中国国产av一级| 麻豆国产97在线/欧美| 久久精品综合一区二区三区| 91在线精品国自产拍蜜月| 性色avwww在线观看| or卡值多少钱| 国产一区亚洲一区在线观看| 免费高清视频大片| 国产国拍精品亚洲av在线观看| 亚洲国产欧美人成| 日韩强制内射视频| 卡戴珊不雅视频在线播放| 91久久精品国产一区二区三区| 国产av不卡久久| 国产精品一及| 嫩草影院精品99| 午夜亚洲福利在线播放| 两个人视频免费观看高清| 综合色丁香网| 日韩成人av中文字幕在线观看 | 国产精品久久电影中文字幕| 乱码一卡2卡4卡精品| 国内久久婷婷六月综合欲色啪| 亚洲美女搞黄在线观看 | 神马国产精品三级电影在线观看| 日本-黄色视频高清免费观看| 变态另类成人亚洲欧美熟女| .国产精品久久| 免费观看在线日韩| 亚洲精品成人久久久久久| 麻豆国产97在线/欧美| 国国产精品蜜臀av免费| 露出奶头的视频| 99riav亚洲国产免费| 国产视频内射| 中文在线观看免费www的网站| 婷婷亚洲欧美| 欧美日韩综合久久久久久| 国产精品亚洲美女久久久| 色5月婷婷丁香| 欧美日韩国产亚洲二区| 亚洲中文字幕一区二区三区有码在线看| 最好的美女福利视频网| 波多野结衣高清作品| 美女内射精品一级片tv| 在线免费观看不下载黄p国产| 国产欧美日韩一区二区精品| 久久精品国产亚洲网站| 2021天堂中文幕一二区在线观| 亚洲经典国产精华液单| 99视频精品全部免费 在线| 午夜影院日韩av| 有码 亚洲区| 国产亚洲精品久久久com| 中文字幕免费在线视频6| 精品久久久久久久久久免费视频| 天堂av国产一区二区熟女人妻| 国产毛片a区久久久久| 亚洲电影在线观看av| 久久欧美精品欧美久久欧美| 久久精品国产清高在天天线| 欧美日韩乱码在线| 欧美色视频一区免费| 久久久久久久久久久丰满| 黄色欧美视频在线观看| 久久人妻av系列| 国产色爽女视频免费观看| 天堂av国产一区二区熟女人妻| 99久久九九国产精品国产免费| 人人妻,人人澡人人爽秒播| 成人毛片a级毛片在线播放| 国产高清不卡午夜福利| 成人av一区二区三区在线看| 亚洲美女黄片视频| 免费av毛片视频| 国产一区二区亚洲精品在线观看| av在线播放精品| 俺也久久电影网| 男人狂女人下面高潮的视频| 成人综合一区亚洲| 久久久久免费精品人妻一区二区| 国产午夜精品论理片| 黄色一级大片看看| 人人妻人人澡人人爽人人夜夜 | 国产精品美女特级片免费视频播放器| 欧美另类亚洲清纯唯美| 色5月婷婷丁香| 天天躁夜夜躁狠狠久久av| 老熟妇仑乱视频hdxx| 日本-黄色视频高清免费观看| 亚洲自拍偷在线| 有码 亚洲区| 亚洲欧美中文字幕日韩二区| 乱码一卡2卡4卡精品| 九九在线视频观看精品| 在线观看一区二区三区| 女生性感内裤真人,穿戴方法视频| 一边摸一边抽搐一进一小说| 校园春色视频在线观看| 成人永久免费在线观看视频| 22中文网久久字幕| 长腿黑丝高跟| 国产真实乱freesex| 精品久久久久久成人av| 国产 一区 欧美 日韩| 国产精品一区二区免费欧美| 国产欧美日韩精品一区二区| av天堂中文字幕网| 国产高清激情床上av| 一边摸一边抽搐一进一小说| 欧美高清性xxxxhd video| 久久久久免费精品人妻一区二区| 欧美最新免费一区二区三区| 免费观看在线日韩| 国产午夜精品久久久久久一区二区三区 | 夜夜看夜夜爽夜夜摸| 成人三级黄色视频| 九色成人免费人妻av| 久久天躁狠狠躁夜夜2o2o| 亚洲人成网站在线播| 精品少妇黑人巨大在线播放 | 亚洲国产欧美人成| 国产精品久久视频播放| 亚洲一区高清亚洲精品| 午夜精品一区二区三区免费看| 人妻制服诱惑在线中文字幕| 国产男靠女视频免费网站| 久久久精品94久久精品| 国模一区二区三区四区视频| 国产成人aa在线观看| 成人二区视频| 一本精品99久久精品77| 国产精品1区2区在线观看.| 精品久久久久久久久亚洲| 免费观看人在逋| 少妇猛男粗大的猛烈进出视频 | 身体一侧抽搐| 日韩精品有码人妻一区| 国产精品精品国产色婷婷| 在线免费十八禁| 亚洲人成网站高清观看| 欧美三级亚洲精品| 久久精品国产自在天天线| 亚洲熟妇熟女久久| 日韩av在线大香蕉| 国产黄色视频一区二区在线观看 | av中文乱码字幕在线| 香蕉av资源在线| 免费观看人在逋| 简卡轻食公司| 尤物成人国产欧美一区二区三区| 久久久久免费精品人妻一区二区| 国产极品精品免费视频能看的| 亚洲av一区综合| 日本一二三区视频观看| 免费看日本二区| av专区在线播放| 欧美区成人在线视频| 成人毛片a级毛片在线播放| 免费不卡的大黄色大毛片视频在线观看 | 亚洲国产精品sss在线观看| 在线免费观看的www视频| 听说在线观看完整版免费高清| 一区二区三区免费毛片| 成熟少妇高潮喷水视频| 国产精品一区二区三区四区免费观看 | 九色成人免费人妻av| 最近的中文字幕免费完整| av在线老鸭窝| av天堂中文字幕网| 丰满的人妻完整版| 午夜免费激情av| a级一级毛片免费在线观看| 国产三级在线视频| 狠狠狠狠99中文字幕| 精品午夜福利视频在线观看一区| 91狼人影院| 国产精品乱码一区二三区的特点| 男女啪啪激烈高潮av片| 亚洲18禁久久av| 麻豆国产av国片精品| 国产精品野战在线观看| 蜜臀久久99精品久久宅男| 99九九线精品视频在线观看视频| 成人二区视频| 久久久久久久久久成人| 亚洲成人av在线免费| 日韩欧美精品v在线| 国产男人的电影天堂91| 麻豆成人午夜福利视频| 一本精品99久久精品77| 国产熟女欧美一区二区| 黑人高潮一二区| 免费在线观看成人毛片| 日韩,欧美,国产一区二区三区 | 夜夜爽天天搞| 国产久久久一区二区三区| 日本-黄色视频高清免费观看| 亚洲精品色激情综合| 亚洲va在线va天堂va国产| 欧美国产日韩亚洲一区| 国产白丝娇喘喷水9色精品| 天堂√8在线中文| 99热6这里只有精品| 国产精品亚洲美女久久久| 婷婷色综合大香蕉| 狠狠狠狠99中文字幕| 一进一出抽搐动态| 日韩欧美国产在线观看| 国产高清不卡午夜福利| 亚洲国产欧洲综合997久久,| 别揉我奶头~嗯~啊~动态视频| 天堂av国产一区二区熟女人妻| 精品一区二区三区人妻视频| 美女大奶头视频| 狂野欧美激情性xxxx在线观看| 国产黄色视频一区二区在线观看 | 国产成人91sexporn| 国模一区二区三区四区视频| 色尼玛亚洲综合影院| 国产精品1区2区在线观看.| 深夜精品福利| 国产日本99.免费观看| 91狼人影院| 亚洲激情五月婷婷啪啪| 国产一区二区激情短视频| 在线看三级毛片| 看十八女毛片水多多多| 亚州av有码| 日韩欧美精品免费久久| 别揉我奶头 嗯啊视频| 亚洲av二区三区四区| 少妇丰满av| 69av精品久久久久久| 欧美极品一区二区三区四区| 国产亚洲精品综合一区在线观看| 大型黄色视频在线免费观看| 国国产精品蜜臀av免费| 最近2019中文字幕mv第一页| 免费一级毛片在线播放高清视频| 色哟哟哟哟哟哟| 亚洲熟妇熟女久久| 乱系列少妇在线播放| 国产激情偷乱视频一区二区| 日韩一本色道免费dvd| 色哟哟·www| 51国产日韩欧美| 国产精品人妻久久久久久| 国产乱人视频| 久久久久久久午夜电影| 热99re8久久精品国产| 内射极品少妇av片p| 波多野结衣高清无吗| 观看美女的网站| 桃色一区二区三区在线观看| 大香蕉久久网| 国产高清视频在线观看网站| 麻豆乱淫一区二区| 97超碰精品成人国产| 亚洲成人久久性| 亚洲av成人av| 麻豆乱淫一区二区| 国语自产精品视频在线第100页| 亚洲熟妇熟女久久| 晚上一个人看的免费电影| 午夜福利在线观看免费完整高清在 | 久久久久久久久中文| 成人三级黄色视频| 久久久色成人| 伊人久久精品亚洲午夜| 日本精品一区二区三区蜜桃| 又爽又黄无遮挡网站| 亚洲国产高清在线一区二区三| 九色成人免费人妻av| 亚洲精品456在线播放app| 伦精品一区二区三区| 18禁黄网站禁片免费观看直播| 一夜夜www| 国产 一区 欧美 日韩| 国产精品一区二区性色av| 观看美女的网站| 亚洲在线观看片| 国产不卡一卡二| 午夜激情福利司机影院| 白带黄色成豆腐渣| 亚洲综合色惰| 在线观看美女被高潮喷水网站| 久久久久久久久久久丰满| 国产高潮美女av| 校园人妻丝袜中文字幕| 国产v大片淫在线免费观看| 网址你懂的国产日韩在线| 97在线视频观看| 直男gayav资源| 99国产极品粉嫩在线观看| 国产亚洲av嫩草精品影院| 欧美绝顶高潮抽搐喷水| 色综合亚洲欧美另类图片| 国产片特级美女逼逼视频| 欧美激情国产日韩精品一区| 亚州av有码| 99热全是精品| 精品福利观看| 变态另类丝袜制服| 欧美高清性xxxxhd video| 97碰自拍视频| 午夜老司机福利剧场| 亚洲欧美成人精品一区二区| 免费电影在线观看免费观看| 12—13女人毛片做爰片一| 亚洲美女视频黄频| 简卡轻食公司| 亚洲人成网站在线播放欧美日韩| 亚洲性夜色夜夜综合| 亚洲国产精品国产精品| 成人特级黄色片久久久久久久| 亚洲精品粉嫩美女一区| 成年女人毛片免费观看观看9| 精品不卡国产一区二区三区| 亚洲国产欧美人成| 精品日产1卡2卡| 国产一级毛片七仙女欲春2| 欧美性猛交╳xxx乱大交人| 欧美高清性xxxxhd video| 国产免费一级a男人的天堂| 99九九线精品视频在线观看视频| 久久久久久久久久久丰满| 麻豆久久精品国产亚洲av| 一个人免费在线观看电影| 国产aⅴ精品一区二区三区波| 嫩草影视91久久| 毛片一级片免费看久久久久| 欧美潮喷喷水| 久久精品国产自在天天线| 老司机影院成人| 免费看光身美女| 国产男人的电影天堂91| 成人二区视频| 国产成人精品久久久久久| 观看美女的网站| 久久久精品大字幕| 日本色播在线视频| 一本一本综合久久| 国产精品乱码一区二三区的特点| 色吧在线观看| a级毛色黄片| 亚洲电影在线观看av| 女人十人毛片免费观看3o分钟| 日韩一区二区视频免费看| 丰满人妻一区二区三区视频av| 国产精品久久久久久久久免| 国产熟女欧美一区二区| 国产精品爽爽va在线观看网站| 精品99又大又爽又粗少妇毛片| 悠悠久久av| 嫩草影院入口| 国产精品亚洲美女久久久| 国产精品日韩av在线免费观看| 久久人妻av系列| 一进一出好大好爽视频| 国产精品一区www在线观看| h日本视频在线播放| 欧美性猛交黑人性爽| 99热这里只有是精品50| 亚洲欧美日韩无卡精品| 免费观看的影片在线观看| 久久鲁丝午夜福利片| 舔av片在线| av专区在线播放| 国产成年人精品一区二区| 亚洲国产精品国产精品| 97人妻精品一区二区三区麻豆| 国产亚洲av嫩草精品影院| 人妻制服诱惑在线中文字幕| 日本精品一区二区三区蜜桃| 国产69精品久久久久777片| 色哟哟哟哟哟哟| av在线播放精品| 97超碰精品成人国产| 久久人人爽人人爽人人片va| 老女人水多毛片| 久久99热这里只有精品18| 亚洲av电影不卡..在线观看| 国产精品福利在线免费观看| 能在线免费观看的黄片| 伦理电影大哥的女人| 听说在线观看完整版免费高清| 国产探花在线观看一区二区| 国产精品,欧美在线| 老司机午夜福利在线观看视频| 久久久久久九九精品二区国产| 在现免费观看毛片| 欧美最新免费一区二区三区| 少妇裸体淫交视频免费看高清| 日韩成人伦理影院| 久久久久久大精品| 亚洲av中文字字幕乱码综合| 久久人人爽人人爽人人片va| 精品久久久久久久久久免费视频| 国产一区二区在线观看日韩| 国产亚洲欧美98| 精品一区二区三区视频在线观看免费| 久久天躁狠狠躁夜夜2o2o| 老熟妇乱子伦视频在线观看| 国产黄色小视频在线观看| 18+在线观看网站| 国产男人的电影天堂91| 午夜久久久久精精品| 久久精品国产亚洲av香蕉五月| 尤物成人国产欧美一区二区三区| 香蕉av资源在线| 天堂影院成人在线观看| 精品一区二区免费观看| 亚洲精品乱码久久久v下载方式| 一区福利在线观看| 啦啦啦啦在线视频资源| 中文字幕av成人在线电影| 99久久精品国产国产毛片| 99久久九九国产精品国产免费| 舔av片在线| 午夜视频国产福利| 成人特级av手机在线观看| 少妇的逼水好多| 露出奶头的视频| 中文字幕av在线有码专区| 成年女人永久免费观看视频| 黄色视频,在线免费观看| 免费一级毛片在线播放高清视频| 精品乱码久久久久久99久播| 中文字幕免费在线视频6| av国产免费在线观看| 国产精品久久久久久久久免| 亚洲欧美日韩无卡精品| 免费看av在线观看网站| 亚洲成人精品中文字幕电影| 日本黄色片子视频| 国产三级在线视频| 国产69精品久久久久777片| 人人妻人人澡欧美一区二区| 少妇熟女欧美另类| 97碰自拍视频| 一进一出抽搐动态| 国产精品久久久久久亚洲av鲁大| 狂野欧美白嫩少妇大欣赏| 99国产极品粉嫩在线观看| 麻豆av噜噜一区二区三区| 国产一区二区三区av在线 | 成人午夜高清在线视频| 精品少妇黑人巨大在线播放 | 日韩一区二区视频免费看| 国产亚洲91精品色在线| 亚洲在线观看片| 国产男人的电影天堂91| 欧美一级a爱片免费观看看| 国产一区亚洲一区在线观看| 又黄又爽又刺激的免费视频.| 免费电影在线观看免费观看| 中文在线观看免费www的网站| 高清日韩中文字幕在线| 精品人妻一区二区三区麻豆 | 亚洲七黄色美女视频| 亚洲最大成人中文| 日日干狠狠操夜夜爽| 又黄又爽又免费观看的视频| 亚洲aⅴ乱码一区二区在线播放| 看免费成人av毛片| 亚洲aⅴ乱码一区二区在线播放| 大香蕉久久网| 亚洲五月天丁香| 欧美区成人在线视频| 男插女下体视频免费在线播放| av.在线天堂| 一本一本综合久久| 91久久精品电影网| 热99re8久久精品国产| 人人妻,人人澡人人爽秒播| 国产亚洲欧美98| 搡老岳熟女国产| 免费黄网站久久成人精品| 亚洲综合色惰| 国产色爽女视频免费观看|