GRB 200219A is a short gamma-ray burst (GRB) with an extended emission (EE) lasting $\sim 90$s. By analyzing data observed with the {\em Swift}/BAT and {\em Fermi}/GBM, we find that a cutoff power law (CPL) model can adequately well fit the spectra of the initial short pulse with $\rm E_{p}=1387^{+232}_{-134}$ keV. More interestingly, together with the EE component and early X-ray data, it exhibits a plateau emission smoothly connected with a $\sim t^{-1}$ segment, and followed by an extremely steep decay. The short GRB composed of those three segments is unique in the {\em Swift} era, and is very difficult to explain with the standard internal/external shock model of a black hole central engine, but could be consistent with the prediction of a magnetar central engine from the merger of an NS binary. We suggest that the plateau emission followed by a $\sim t^{-1}$ decay phase is powered by the spin-down of millisecond magnetar, which loses its rotation energy via GW quadrupole radiation. Then, the abrupt drop decay is caused by magnetar collapsing into black hole before switching to the EM-dominated emission. This is the first short GRB for which the X-ray emission has such an intriguing feature powered by magnetar via GW-dominated radiation. If this is the case, one can estimate the physical parameters of a magnetar, the GW signal powered by magnetar and the merger-nova emission are also discussed.

The LIGO-Virgo and Fermi collaborations recently reported a possible joint detection of a sub-threshold gravitational wave (GW) event and a sub-threshold gamma-ray burst (GRB), GBM-190816, that occurred 1.57 s after the merger. We perform an independent analysis of the publicly available data and investigate the physical implications of this potential association. By carefully studying the following properties of GBM-190816 using Fermi/GBM data, including signal-to-noise ratio, duration, f-parameter, spectral properties, energetic properties, and its compliance with some GRB statistical correlations, we confirm that this event is likely a typical short GRB. Assuming its association with the sub-threshold GW event, the inferred luminosity is $1.47_{-1.04}^{+3.40} \times 10^{49}$ erg s$^{-1}$. Based on the available information of the sub-threshold GW event, we infer the mass ratio q of the compact binary as $q=2.26_{-1.43}^{+2.75}$ according to the reported range of luminosity distance. If the heavier compact object has a mass > 3 solar masses, q can be further constrained to $q=2.26_{-0.12}^{+2.75}$. The leading physical scenario invokes an NS-BH merger system with the NS tidally disrupted. Within this scenario, we constrain the physical properties of such a system to produce a GRB. The GW data may also allow an NS-BH system with no tidal disruption of the NS or a BH-BH merger. We apply the charged compact binary coalescence (cCBC) theory (for both a constant charge and an increasing charge for the merging members) to derive the model parameters to account for GBM-190816 and found that the required parameters are extreme. Finally, we argue that the fact that the observed GW-GRB delay time scale is comparable to that of GW170817/GRB 170817A suggests that the GW-GRB time delay of these two cases is mainly defined by the time scale for the jet to propagate to the GRB emission site.

The observed internal plateau of X-ray emission in some short GRBs suggests the formation of a remnant supra-massive magnetar following a double neutron star (NS) merger. In this paper, we assume that the rotational energy is lost mainly via gravitational wave (GW) radiation instead of magnetic dipole (MD) radiation, and present further constraints on the NS nuclear equation of state (EoS) via mass quadrupole deformation and r-mode fluid oscillations of the magnetar. We present two short GRBs with measured redshifts, 101219A and 160821B, whose X-ray light curves exhibit an internal plateau. This suggests that a supra-massive NS may survive as the central engine. By considering twelve NS EoSs, within the mass quadrupole deformation scenario we find that the GM1, DD2, and DDME2 models give an $M_{\rm p}$ band falling within the 2$\sigma$ region of the proto-magnetar mass distribution for $\varepsilon=0.01$. This is consistent with the constraints from the MD radiation dominated model of rotational energy loss. However, for an r-mode fluid oscillation model with $\alpha=0.1$ the data suggest that the NS EOS is close to the Shen and APR models, which is obviously different from the MD radiation dominated and mass quadrupole deformation cases.

Short GRBs with extended emission (EE) that are composed initially of a short-hard spike and followed by a long-lasting EE, are thought to be classified as a subsection of short GRBs. The narrow energy band available during the {\em Swift} era combined with a lack of spectral information prevented discovery of the intrinsic properties of those events. In this paper, we performed a systematic search of short GRBs with EE by using all available {\em Fermi}/GBM data. The search identified 26 GBM-detected short GRBs with EE that are similar to GRB 060614 observed by {\em Swift}/BAT. We focus on investigating the spectral and temporal properties for both the hard spike and the EE components of all 26 GRBs, and explore differences and possible correlations between them. We find that while the peak energy ($E_{\rm p}$) of the hard spikes is a little bit harder than that of the EE, but their fluences are comparable. The harder $E_{\rm p}$ seems to correspond to a larger fluence and peak flux with a large scatter for both the hard spike and EE components. Moreover, the $E_{\rm p}$ of both the hard spikes and EE are compared to other short GRBs. Finally, we also compare the properties of GRB 170817A with those short GRBs with EE and find no significant statistical differences between them. We find that GRB 170817A has the lowest $E_{\rm p}$, likely because it was off-axis.

Based on the early-year observations from Neil Gehrels \emph{Swift} Observatory, \cite{liang07} performed a systematic analysis for the shallow decay component of gamma-ray bursts (GRBs) X-ray afterglow, in order to explore its physical origin. Here we revisit the analysis with an updated sample (with Swift/XRT GRBs between February 2004 and July 2017). We find that with a larger sample, 1) the distributions of the characteristic properties of the shallow decay phase (e.g. $t_{b}$ , $S_{X}$, $\Gamma_{X,1}$, and $\alpha_{X,1}$) still accords with normal or lognormal distribution; 2) $\Gamma_{X,1}$ and $\Gamma_{\gamma}$ still show no correlation, but the tentative correlations of durations, energy fluences, and isotropic energies between the gamma-ray and X-ray phases still exist; 3) for most GRBs, there is no significant spectral evolution between the shallow decay segment and its follow-up segment, and the latter is usually consistent with the external-shock models; 4) assuming that the central engine has a power-law luminosity release history as $L(t)=L_{0}(\frac{t}{t_{0}})^{-q}$, we find that the value $q$ is mainly distributed between -0.5 and 0.5, with an average value of 0.16$\pm$ 0.12; 5) the tentative correlation between $E_{\rm{iso},X}$ and $t'_{b}$ disappears, so that the global 3-parameter correlation ($E_{\rm{iso},X}-E'_{p}-t'_{b}$) becomes less significant; 6) the anti-correlation between $L_{X}$ and $t'_{b}$ and the three-parameter correlation ($E_{\rm{iso},\gamma}-L_{X}-t_{b}$) indeed exist with a high confidence level. Overall, our results are generally consistent with \cite{liang07}, confirming their suggestion that the shallow decay segment in most bursts is consistent with an external forward shock origin, probably due to a continuous energy injection from a long-lived central engine.

The event GW170817/GRB 170817A, discovered via the successful joint observation of its gravitational wave radiation and its multi-wavelength electromagnetic counterparts, was the first definite "smoking-gun" from the merger of two neutron stars (NSs). However, the remnant of the merger remains unknown. Piro et al. recently claimed that a low-significance X-ray variability in GRB 170817A. By systematically comparing the properties of variability in the afterglow of GRB 170817A and X-ray flares in GRB afterglows, we find that this X-ray variability seems to share similar statistical correlations with X-ray flares in GRB afterglows. We further investigate several possible merger product scenarios to see whether they can produce the observed X-ray variability in GRB 170817A. The first scenario invokes a stable magnetar as the central engine producing the later X-ray variability via differential rotation or fall-back accretion onto the NS. The second scenario invokes a black hole as the central engine with a fall-back accretion process. The final scenario is a central engine with a long-lived supra-massive NS. We find that the first two scenarios have difficulty producing the later X-ray variability, which requires either an impractical NS magnetic field or an extraordinarily large stellar envelope and an extremely long accretion timescale. However, the third scenario seems to be consistent with observations, and the later X-ray variability can be produced by the magnetosphere which is expelled following the collapse of the NS with a $B_p\in(3.6, 13.5)\times10^{13}$ G.

The merging process of binary neutron stars is the source of the quadrupole gravitational-wave (GW) radiation. Its remnant can be either a supra-massive or a stable NS. Such evidence of magnetar signature has been supported indirectly by observed X-ray plateau of some gamma-ray bursts (GRBs) afterglow. Recently, Xue et al. (2019) discovered an X-ray transient CDF-S XT2 that is claimed to be powered by a stable magnetar from the merger of double NS. In this paper, we revisit the X-ray emission of CDF-S XT2 and find that it is more consistent with a supra-massive magnetar central engine, surviving thousands of seconds to collapse black hole. We present the comparisons of the X-ray plateau luminosity, break time, and the parameters of magnetar between CDF-S XT2 and other short GRBs with internal plateau samples. By adapting the collapse time to constrain the equation of state (EOS), we find that three EOSs (GM1, DD2, and DDME2) are consistent with the observational data. On the other hand, if the most released rotation energy of magnetar is dominated by GW radiation, we also constrain the upper limit of ellipticity of NS for given EOS, and it is in a range in $[0.89-1.8]\times 10^{-3}$. Its GW signal cannot be detected by aLIGO or even for more sensitive Einstein Telescope in the future.

The long-lasting X-ray plateau emission in long gamma-ray bursts (GRBs) shows observational evidence for ongoing energy injection, which may be from magnetar spindown due to energy released via either magnetic dipole (MD) or gravitational wave (GW) radiation. In this paper, by systematically analyzing the {\em Swift}/XRT light curves detected before 2018 July, we find 45 light curves with a measured redshift that monotonically decay as a smooth broken power law. By assuming that the central engines of these GRBs are newly born magnetars, we measure the braking index $n$ of putative millisecond magnetars, due to MD and GW radiations. The inferred braking indices are not close to 3 or 5, but range between them with a normal distribution ($n_{\rm c}=4.02\pm 0.11$). We define a dimensionless parameter $\Re$, which is the ratio between the MD and GW components, and find that the energy released via magnetar spindown in most GRBs of our sample is dominated by GW radiation for $P_0=3$ ms and $\epsilon=0.005$ and 0.01. On the other hand, we find that $\Re$ and the braking index $n$ seem to be anticorrelated within a large systematic error at $t=0$, but depend on the values of the parameters $P_0$ and $\epsilon$. These results suggest that the contribution of GW radiation cannot be ignored, and that a larger braking index leads to GWs dominating the energy released during magnetar spindown if indeed magnetars are operating in some long GRBs.

The radiation mechanism of the prompt emission of gamma-ray bursts (GRBs) remains an open question. Although their spectra are usually well fitted with the empirical Band function, which is widely believed to be fully non-thermal and interpreted as an optically thin synchrotron emission, accumulating evidence shows that a thermal component actually exists. In this paper, a multi-color blackbody (mBB) model is proposed for the time-integrated spectrum of GRB 081221 by assuming a power-law distribution of the thermal luminosities with temperature, which manifests photospheric emissions from a different radius and/or angle. The effects of the minimum temperature $kT_{\min }$, the maximum temperature $kT_{\max }$ and the power law index $m$ of the luminosity distribution of an mBB are discussed. The fitting to the time-integrated spectrum during the bright phase (from 20s to 30s since the trigger) of GRB 081221 by the mBB model yields $kT_{\min }$ = $4.4\pm 0.3$ keV, $kT_{\max }$ = 57.0 $_{-1.4}^{+1.6}$ keV, and $m=-0.46_{-0.06}^{+0.05}$. When the time bin is small enough, the time-resolved spectra of GRB 081221 are well fitted with a series of single-temperature blackbodies. Our results imply the prompt emission of GRB 081221 is dominated by the photosphere emission and its time-integrated spectrum is a superposition of pure blackbody components at different times, indicating that some empirical Band spectra may be interpreted as mBB if the temperature is widely distributed.

A newly-born magnetar is thought to be central engine of some long gamma-ray bursts (GRBs). We investigate the evolution of the electromagnetic (EM) emission from the magnetic dipole (MD) radiation wind injected by spin-down of a newly-born magnetar via both quadrupole gravitational-wave (GW) and MD radiations. We show that the EM luminosity evolves as $L_{\rm em}\propto (1+t/\tau_c)^{\alpha}$, and $\alpha$ is $-1$ and $-2$ in the GW and MD radiation dominated scenarios, respectively. Transition from the GW to MD radiation dominated epoch may show up as a smooth break with slope changing from $-1$ to $-2$. If the magnetar collapses to a black hole before $\tau_c$, the MD radiation should be shut down, then the EM light curve should be a plateau followed by a sharp drop. The expected generic light curve in this paradigm is consistent with the canonical X-ray light curve of {\em Swift} long GRBs. The X-ray emission of several long GRBs are identified and interpreted as magnetar spin-down via GW or MD, as well as constrain the physical parameters of magnetar. The combination of MD emission and GRB afterglows may make the diversity of the observed X-ray light curves. This may interpret the observed chromatic behaviors of the X-ray and optical afterglow light curves and the extremely low detection rate of a jet-like break in the X-ray afterglow light curves of long GRBs.

Two-episode emission components separated by quiescent gaps in the prompt emission of gamma-ray bursts (GRBs) have been observed in the {\em Swift} era, but there is a lack of spectral information due to the narrow energy band of the {\em Swift}/Burst Alert Telescope. In this paper, a systematic analysis of the spectral and temporal properties of the prompt emission of 101 {\em Fermi}/Gamma-ray Burst Monitor detected long GRBs show the existence of two-episode emission components in the light curves, with quiescent times of up to hundreds of seconds. We focus on investigating the differences of those two emission episodes. We find that the light curves of the two emission components exhibit different behavior, e.g., a soft emission component that either precedes or follows the main prompt emission or that the intensity of the two emission episodes are comparable with each other. No statistically significant correlation in the duration of the two emission episodes can be claimed. We define a new parameter $\varepsilon$ as the ratio of the peak flux of the first and second emission episodes and find that a higher $\varepsilon$ corresponds to a larger fluence. The preferred spectral model in our analysis is a cutoff power-law model for most GRBs. The distribution of $E_p$ for episodes I and II range from tens of keV to 1000 keV with a lognormal fit and there are no significant differences between them. Moreover, we do not find significant relationships between $\varepsilon$ and $E_p$ for the two emission episodes. Those results suggest that these two-episode emission components likely share the same physical origin.

The favored progenitor model for Gamma-ray Bursts (GRBs) with Supernova (SN) association is the core collapse of massive stars. One possible outcome of such a collapse is a rapidly spinning, strongly magnetized neutron star ("magnetar"). We systematically analyze the multi-wavelength data of GRB/SN associations detected by several instruments before 2017 June. Twenty GRB/SN systems have been confirmed via direct spectroscopic evidence or a clear light curve bump, as well as some spectroscopic evidence resembling a GRB-SN. We derive/collect the basic physical parameters of the GRBs and the SNe, and look for correlations among these parameters. We find that the peak brightness, $\rm ^{56}Ni$ mass, and explosion energy of SNe associated with GRBs are statistically higher than other Type Ib/c SNe. A statistically significant relation between the peak energy of GRBs and the peak brightness of their associated SNe is confirmed. No significant correlations are found between the GRB energies (either isotropic or beaming-corrected) and the supernova energy. We investigate the energy partition within these systems and find that the beaming-corrected GRB energy of most systems is smaller than the SN energy, with less than 30\% of the total energy distributed in the relativistic jet. The total energy of the systems is typically smaller than the maximum available energy of a millisecond magnetar ($2\times 10^{52}$ erg), especially if aspherical SN explosions are considered. The data are consistent with-though not a proof of-the hypothesis that most, but not all, GRB/SN systems are powered by millisecond magnetars.

Double neutron star (DNS) merger events are promosing candidates of short Gamma-ray Burst (sGRB) progenitors as well as high-frequecy gravitational wave (GW) emitters. On August 17, 2017, such a coinciding event was detected by both the LIGO-Virgo gravitational wave detector network as GW170817 and Gamma-Ray Monitor on board NASA's {\it Fermi} Space Telescope as GRB 170817A. Here we show that the fluence and spectral peak energy of this sGRB fall into the lower portion of the distributions of known sGRBs. Its peak isotropic luminosity is abnormally low. The estimated event rate density above this luminosity is at least $190^{+440}_{-160} {\rm Gpc^{-3} \ yr^{-1}}$, which is close to but still below the DNS merger event rate density. This event likely originates from a structured jet viewed from a large viewing angle. There are similar faint soft GRBs in the {\it Fermi} archival data, a small fraction of which might belong to this new population of nearby, low-luminosity sGRBs.

Motivating by the discovery of association between GW 170817 and sGRB 170817A, we present a comprehensive analysis for sGRBs observed with Fermi/GBM in 9 operation years and study the properties of sGRB 170817A -like events. We derive a catalog of 275 typical sGRBs and 48 sGRB 170817A-like weak events from the GBM data of 2217 GRBs. We visibly identify two patterns of their light curve, single episode (Pattern I, 61\% of the SGRBs) and multiple episodes (Pattern II, 39\% of the SGRBs). Their duration distribution shows a tentative bimodal feature. Their spectra can be fitted with a cutoff power-law model, except for 4 sGRBs, and the spectral indices normally distribute at $\Gamma=0.69\pm 0.40$. Their $E_p$ values show a tentative bimodal distribution with peaks at 145 keV and 500 keV. No correlation among $T_{90}$, $E_p$, and $\Gamma$ is found. GRB 170817A is a soft, weak sGRB with $ E_{p}=124\pm 106$ keV, $L_{\rm iso}=(5.67\pm4.65)\times10^{46}\rm ~erg~s^{-1}$, and $E_{\rm iso}=(3.23\pm2.65)\times10^{46}\rm ~erg$. It follows the $E_{\rm iso}-E_{\rm p}$ relation of typical short GRBs. Its lightcurve is of Pattern II. Two lightcurve patterns, together with the potential two components in the $E_{\rm p}$ and $T_{90}$ distributions, we suspect that the current sample may include two distinct types of sGRBs from different progenitors. sGRB 170817A-like events may be from NS-NS mergers and those sGRBs with a Pattern I lightcurve may be from another distinct type of compact binary.

GRB 160625B is an extremely bright GRB with three distinct emission episodes. By analyzing its data observed with the GBM and LAT on board the {\em Fermi} mission, we find that a multi-color black body (mBB) model can be used to fit the spectra of initial short episode (Episode I) very well within the hypothesis of photosphere emission of a fireball model. The time-resolved spectra of its main episode (Episode II), which was detected with both GBM and LAT after a long quiet stage ($\sim 180$ seconds) of the initial episode, can be fitted with a model composing of an mBB component plus a cutoff power-law (CPL) component. This GRB was detected again in the GBM and LAT bands with a long extended emission (Episode III) after a quiet period of $\sim 300$ seconds. The spectrum of Episode III is adequately fitted with a CPL plus a single power-law models, and no mBB component is required. These features may imply that the emission of three episodes are dominated by distinct physics process, i.e., Episode I is possible from cocoon emission surrounding the relativistic jet, Episode II may be from photosphere emission and internal shock of relativistic jet, and Episode III is contributed by internal and external shocks of relativistic jet. On the other hand, both X-ray and optical afterglows are consistent with standard external shocks model.

The X-ray plateau followed by a steep decay ("internal plateau") has been observed in both long and short gamma-ray burst (GRBs), implying a millisecond magnetar operating in some GRBs. The sharp decay at the end of plateau, marking the abrupt cessation of the magnetar central engine, has been considered as the collapse of a supra-massive magnetar to a black hole (BH) when it spins down. If "internal plateau" is indeed the evidence of a magnetar central engine, a natural expectation is a signature from the new-born BH in some candidates. In this work, we find that GRB 070110 is a particular case, which shows a small X-ray bump following its "internal plateau". We interpret the plateau with a spin-down supra-massive magnetar and the X-ray bump with a fall-back BH accretion. This indicates that the new-born BH is likely active in some GRBs. Therefore, GRB 070110-like events may provide a further support to the magnetar central engine model and enable us to investigate the properties of the magnetar as well as the new-born BH.

The jet compositions, central engines, and progenitors of gamma-ray bursts (GRBs) remain open questions in GRB physics. Applying broadband observations, including GRB prompt emission and afterglow properties derived from {\em Fermi} and {\em Swift} data, as well as from Keck host-galaxy observations, we address these questions for the peculiar, bright GRB 110731A. By using the pair-opacity method, we derive $\Gamma_{0}>190$ during the prompt emission phase. Alternatively, we derive $\Gamma_{0} \approx 580$ and $\Gamma_{0} \approx 154$ by invoking the early-afterglow phase within the homogeneous density and wind cases, respectively. On the other hand, nondetection of a thermal component in the spectra suggests that the prompt emission is likely powered by dissipation of a Poynting-flux-dominated jet leading to synchrotron radiation in an optically thin region. The nondetection of a jet break in the X-ray and optical bands allows us to place a lower limit on the jet opening angle $\theta_{\rm j}>5.5^{\circ}$. Within a millisecond magnetar central engine scenario, we derive the period $P_{0}$ and polar magnetic field strength $B_{\rm p}$, which have extreme (but still allowed) values. The moderately short observed duration (7.3\,s) and relatively large redshift ($z=2.83$) places the burst as a "rest-frame short" GRB, so the progenitor of the burst is subject to debate. Its relatively large $f_{{\rm eff}, z}$ parameter (ratio of the 1\,s peak flux of a pseudo-GRB and the background flux) and a large physical offset from a potential host galaxy suggest that the progenitor of GRB 110731A may be a compact-star merger.

The merger of a double neutron star (NS-NS) binary may result in a rapidly rotating massive NS with an extremely strong magnetic field (i.e., a millisecond magnetar). In this case, the magnetic spin-down of the NS remnant provides an additional source of sustained energy injection, which would continuously power the merger ejecta. The thermal emission from the merger ejecta would give rise to a bright optical "magnetar-powered merger-nova". In this work, we carry out a complete search for magnetar-powered merger-nova from Swift short gamma-ray burst (SGRB) sample. We focus on short GRBs with extended emission or internal plateau, which may signify the presence of magnetars as the central engine. We eventually find three candidates of magnetar-powered merger-nova from the late observations of GRB 050724, GRB 070714B and GRB 061006. With standard parameter values, the magnetar remnant scenario could well interpret the multi-band data of all three bursts, including the extended emission and their late chromatic features in the optical and X-ray data. The peak luminosities of these merger-novae reach several times $10^{42}~{\rm erg~s^{-1}}$, more than one order of magnitude brighter than the traditional "kilo-novae" with peak luminosity of $\sim 10^{41}~{\rm erg~s^{-1}}$. Intense, multi-color late time observations of short GRBs are encouraged to identify more merger-novae in the future.

[Zhang 2014] propose to redefine the true GRB central engine activity duration, $t_{\rm burst}$, by considering the contributions from the prompt $\gamma$-ray emission, X-ray flare and internal plateau features. With a comprehensive study on a large sample of Swift GRBs, it is shown that the $t_{\rm burst}$ distribution in the observer frame consists of a bimodal feature, suggesting the existence of a new population of ultra-long GRBs. In this work, we make a series of further studies on $t_{\rm burst}$: we update the Swift GRB sample up to June 2016; we investigate the properties of $t_{\rm burst}$ distribution in the rest frame; we redefine $t_{\rm burst}$ by involving external plateau contributions; we make a multiple parameter analysis to investigate whether the bursts within the ultra-long population being statistically different in sense of other features besides the duration distribution. We find that for all situations, the distribution of $t_{\rm burst}$ requires two normal distributions in logarithmic space to provide a good fit, both in observer frame and rest frame. Considering the observational gap effect would not completely erase the bimodal distribution feature. However the bursts within the ultra-long population may have no statistically different in sense of other features besides the duration term. We thus suggest that if the ultra-long population of GRBs indeed exists, their central engine mechanism and radiation mechanism should be similar to the normal population, but they have longer central engine activity timescale.

By fitting the bolometric light curves of 31 super-luminous supernovae (SLSNe) with the magnetar engine model, we derive the ejecta masses and magnetar parameters for these SLSNe. The lower boundary of magnetic field strengths of SLSN magnetars can be set just around the critical field strength $B_{\rm c}$ of electron Landau quantization. In more details, SLSN magnetars can further be divided into two subclasses of magnetic fields of $\sim(1-5)B_{\rm c}$ and $\sim(5-10)B_{\rm c}$, respectively. It is revealed that these two subclasses of magnetars are just associated with the slow-evolving and fast-evolving bolometric light curves of SLSNe. In comparison, the magnetars harbored in gamma-ray bursts (GRBs) and associated hypernovae are usually inferred to have much higher magnetic fields with a lower boundary about $\sim10B_{\rm c}$. This robustly suggests that it is the magnetic fields that play the crucial role in distinguishing SLSNe from GRBs/hypernovae. The rotational energy of SLSN magnetars are found to be correlated with the masses of supernova ejecta, which provides a clue to explore the nature of their progenitors. Moreover, the distribution of ejecta masses of SLSNe is basically intermediate between those of normal core-collapse supernovae and hypernovae. This could indicate an intrinsic connection among these different stellar explosions.

GRB 160821B is a short gamma-ray burst (GRB) at redshift $z=0.16$, with a duration less than 1 second and without detection of any "extended emission" up to more than 100 seconds in both {\em Swift}/BAT and {\em Fermi}/GBM bands. An X-ray plateau with a sharp drop 180 seconds after the BAT trigger was observed with {\em Swift}/XRT. No supernova or kilo-nova signature was detected. Assuming the central engine of this SGRB is a recently born supra-massive magnetar, we can explain the SGRB as jet radiation and its X-ray plateau as the internal energy dissipation of the pulsar wind as it spins down. We constrain its surface magnetic field as $B_{\rm p}<3.12\times 10^{16}$ G and initial spin period as $P_0< 8.5\times 10^{-3}$ seconds. Its equation of state is consistent with the GM1 model with $M_{\rm TOV} \sim 2.37 M_\odot$ and ellipticity $\epsilon<0.07$. Its gravitational wave (GW) radiation may be detectable with the future Einstein Telescope, but is much weaker than the current detectability limit of advanced-LIGO. The GW radiation of such an event would be detectable by advanced-LIGO if it occurred at a distance of 100 Mpc ($z=0.023$).

We summarize basic observational properties of gamma-ray bursts (GRBs), including prompt emission properties, afterglow properties, and classification schemes. We also briefly comment on the current physical understanding of these properties.

Gamma ray bursts (GRBs) are classified into long and short categories based on their durations. Broad band studies suggest that these two categories of objects roughly correspond to two different classes of progenitor systems, i.e. compact star mergers (Type I) vs. massive star core collapse (Type II). However, the duration criterion sometimes leads to mis-identification of the progenitor systems. We perform a comprehensive multi-wavelength comparative study between duration-defined long GRBs and short GRBs as well as the so-called "consensus" long GRBs and short GRBs (which are believed to be more closely related to the two types of progenitor systems). The parameters we study include two parts: the prompt emission properties including duration ($T_{90}$), spectral peak energy ($E_{\rm p}$), low energy photon index ($\alpha$), isotropic $\gamma$-ray energy ($E_{\rm \gamma, iso}$), isotropic peak luminosity ($L_{\rm p,iso}$), and the amplitude parameters ($f$ and $f_{\rm eff}$); and the host galaxy properties including stellar mass ($M_*$), star formation rate (SFR), metallicity ([X/H]), half light radius ($R_{50}$), angular and physical ($R_{\rm off}$) offset of the afterglow from the center of the host galaxy, the normalized offset ($r_{\rm off}=R_{\rm off}/R_{50}$), and the brightness fraction $F_{\rm light}$. For most parameters, we find interesting overlapping properties between the two populations in both 1D and 2D distribution plots. The three best parameters for the classification purpose are $T_{90}$, $f_{\rm eff}$, and $F_{\rm light}$. However, no single parameter alone is good enough to place a particular burst into the right physical category, suggesting a need of multiple criteria for physical classification.

A rapidly spinning, strongly magnetized neutron star is invoked as the central engine for some Gamma-ray bursts (GRBs), especially, the $"$internal plateau$"$ feature of X-ray afterglow. However, for these $"$internal plateau$"$ GRBs, how to produce their prompt emission remains an open question. Two different physical process have been proposed in the literature, (1) a new-born neutron star is surrounded by a hyper-accreting and neutrino cooling disk, the GRB jet can be powered by neutrino annihilation aligning the spin axis; (2) a differentially rotating millisecond pulsar was formed due to different angular velocity between the interior core and outer shell parts of the neutron star, which can power an episodic GRB jet. In this paper, by analyzing the data of one peculiar GRB 070110 (with internal plateau), we try to test which model being favored. By deriving the physical parameters of magnetar with observational data, the parameter regime for initial period ($P_{0\rm }$) and surface polar cap magnetic field ($B_{\rm p}$) of the central NS are $(0.96\sim 1.2 )~\rm ms$ and $(2.4\sim 3.7)\times 10^{14}~\rm G$, respectively. The radiative efficiency of prompt emission is about $\eta_{\gamma} \sim 6\%$. However, the radiative efficiency of internal plateau ($\eta_{\rm X}$) is larger than $31\%$ assuming the $M_{\rm NS}\sim1.4 M_{\odot}$ and $P_{0\rm }\sim1.2 ~\rm ms$. The clear difference between the radiation efficiencies of prompt emission and internal plateau implies that they maybe originated from different components (e.g. prompt emission from the relativistic jet powered by neutrino annihilation, while the internal plateau from the magnetic outflow wind).

Binary neutron star mergers are strong gravitational wave (GW) sources and the leading candidates to interpret short duration gamma-ray bursts (SGRBs). Under the assumptions that SGRBs are produced by double neutron star mergers and that the X-ray plateau followed by a steep decay as observed in SGRB X-ray light curves marks the collapse of a supra-massive neutron star to a black hole (BH), we use the statistical observational properties of {\em Swift} SGRBs and the mass distribution of Galactic double neutron star systems to place constraints on the neutron star equation of state (EoS) and the properties of the post-merger product. We show that current observations already put following interesting constraints: 1) A neutron star EoS with a maximum mass close to a parameterization of $M_{\rm max} = 2.37\,M_\odot (1+1.58\times10^{-10} P^{-2.84})$ is favored; 2) The fractions for the several outcomes of NS-NS mergers are as follows: $\sim40\%$ prompt BHs, $\sim30\%$ supra-massive NSs that collapse to BHs in a range of delay time scales, and $\sim30\%$ stable NSs that never collapse; 3) The initial spin of the newly born supra-massive NSs should be near the breakup limit ($P_i\sim1 {\rm ms}$), which is consistent with the merger scenario; 4) The surface magnetic field of the merger products is typically $\sim 10^{15}$ G; 5) The ellipticity of the supra-massive NSs is $\epsilon \sim (0.004 - 0.007)$, so that strong GW radiation is released post the merger; 6) Even though the initial spin energy of the merger product is similar, the final energy output of the merger product that goes into the electromagnetic channel varies in a wide range from several $10^{49}$ erg to several $10^{52}$ erg, since a good fraction of spin energy is either released in the form of GW or falls into the black hole as the supra-massive NS collapses.

Observations of very early multi-wavelength afterglows are critical to reveal the properties of the radiating fireball and its environment as well as the central engine of gamma-ray bursts (GRBs). We report our optical observations of GRB 111228A from 95 sec to about 50 hours after the burst trigger and investigate its properties of the prompt gamma-rays and the ambient medium using our data and the data observed with {\em Swift} and {\em Fermi} missions. Our joint optical and X-ray spectral fits to the afterglow data show that the ambient medium features as low dust-to-gas ratio. Incorporating the energy injection effect, our best fit to the afterglow lightcurves with the standard afterglow model via the Markov Chain Monte Carlo (MCMC) technique shows that $\epsilon_e=(6.9\pm 0.3)\times 10^{-2}$, $\epsilon_B=(7.73\pm 0.62)\times 10^{-6}$, $E_{\rm K}=(6.32\pm 0.86)\times 10^{53}\rm erg$, $n=0.100\pm 0.014$ cm$^{-3}$. The low medium density likely implies that the afterglow jet may be in a halo or in a hot ISM. Achromatic shallow decay segment observed in the optical and X-ray bands is well explained with the long-lasting energy injection from the central engine, which would be a magnetar with a period of about 1.92 ms inferred from the data. The $E_p$ of its time-integrated prompt gamma-ray spectrum is $\sim 26$ KeV. Using the initial Lorentz factor ($\Gamma_0=476^{+225}_{-237}$) derived from our afterglow model fit, it is found that GRB 111228A satisfies the $L_{\rm iso}-E_{\rm p,z}-\Gamma_0$ relation and bridges the typical GRBs and low luminosity GRBs in this relation.

One favored progenitor model for short duration gamma-ray bursts (SGRBs) is the coalescence of two neutron stars (NS-NS). One possible outcome of such a merger would be a rapidly spinning, strongly magnetized neutron star (known as a millisecond magnetar). These magnetars may be "supra-massive", implying they would collapse to black holes after losing centrifugal support due to magnetic dipole spindown. By systematically analyzing the BAT-XRT light curves of all short GRBs detected by {\em swift}, we test how well the data are consistent with this central engine model of short GRBs. We find that the so-called "extended emission" observed with BAT in some short GRBs are fundamentally the same component as the "internal X-ray plateau" as observed in many short GRBs, which is defined as a plateau in the lightcurve followed by a very rapid drop. Based on how likely a short GRB hosts a magnetar, we characterize the entire {\em Swift} short GRB sample into three categories: the "internal plateau" sample, the "external plateau" sample, and the "no plateau" sample. The derived magnetar surface magnetic field $B_{\rm p}$ and the initial spin period $P_0$ fall into the reasonable range. No GRBs in the internal plateau sample have the total energy exceeding the maximum energy budget of a millisecond magnetar. Assuming that the rapid fall time at the end of the internal plateau is the collapse time of a supra-massive magnetar to a black hole, and applying the measured mass distribution of NS-NS systems in our Galaxy, we constrain the neutron star equation of state (EOS). The data suggest that the NS equation of state is close to the GM1 model, which has a maximum non-rotating NS mass $M_{\rm TOV} \sim 2.37 M_\odot$.

Traditionally gamma-ray bursts (GRBs) are classified in the $T_{90}$-hardness ratio two-dimensional plane into long/soft and short/hard GRBs. In this paper, we suggest to add the "amplitude" of GRB prompt emission as the third dimension as a complementary criterion to classify GRBs, especially those of short durations. We define three new parameters $f$, $f_{\rm eff}$ and $f_{\rm eff,z}$ as ratios between the measured/simulated peak flux of a GRB/pseudo-GRB and the flux background, and discuss the applications of these parameters to GRB classification. We systematically derive these parameters to find that most short GRBs are likely not "tip-of-iceberg" of long GRBs. However, one needs to be cautious if a short GRB has a relatively small $f$ (e.g. $f<1.5$), since the chance for an intrinsically long GRB to appear as a "disguised" short GRB is higher. Based on avaialble data, we quantify the probability of a disguised short GRB below a certain $f$ value is as $P (<f)\sim 0.78^{+0.71}_{-0.4} f^{-4.33\pm 1.84}$. By progressively "moving" a long GRB to higher redshifts through simulations, we also find that most long GRBs would show up as rest-frame short GRBs above a certain redshift.

A rapidly spinning, strongly magnetized neutron star (magnetar) has been proposed as one possible candidate of the central engine of gamma-ray bursts (GRBs). We systematically analyze the Swift/XRT light curves of long GRBs detected before 2013 August, and characterize them into four categories based on how likely they may harbor a magnetar central engine: Gold, Silver, Aluminum, and Non-magnetar. We also independently analyze the data of short GRBs with a putative magnetar central engine. We then perform a statistical study of various properties of the magnetar samples and the non-magnetar sample, and investigate whether the data are consistent with the hypothesis that there exist two types of central engines. By deriving the physical parameters of the putative magnetars, we find that the observations of the Gold and Silver samples are generally consistent with the predictions of the magnetar model. For a reasonable beaming factor for long GRBs, the derived magnetar surface magnetic field $B_p$ and initial spin period $P_0$ fall into the reasonable range. Magnetar winds in short GRBs, on the other hand, are consistent with being isotropic. No GRB in the magnetar sample has a beam-corrected total energy exceeding the maximum energy budget defined by the initial spin energy of the magnetar, while some non-magnetar GRBs do violate such a limit. With beaming correction, on average the non-magnetar sample is more energetic and luminous than the magnetar samples. Our analysis hints that millisecond magnetars are likely operating in a good fraction, but probably not all, GRBs.

We present a time-resolved spectral analysis of 51 long and 11 short bright GRBs observed with the {\em Femri}/GBM, paying special attention to $E_{\rm p}$ evolution within a same burst. Among 8 single-pulse long GRBs, 5 show hard-to-soft evolution, while 3 show intensity-tracking. The multi-pulse long GRBs have more complicated patterns. Among the GRBs whose time-resolved spectrum is available for the first pulse, almost half (15/32 GRBs) show clear hard-to-soft evolution, and the other half (17/32 GRBs) show clear intensity-tracking. Later pulses typically show the tracking behavior, although a hard-to-soft evolution pattern was identified in the 2nd pulse of 2 GRBs whose pulses are well separated. Statistically, the hard-to-soft evolution pulses tend to be more asymmetric than the intensity-tracking ones, with a steeper rising wing than the falling wing. Short GRBs have $E_{\rm p}$ tracking intensity exclusively with the 16ms time resolution analysis. We performed a simulation analysis, and suggest that at least for some bursts, the late intensity-tracking pulses could be a consequence of overlapping hard-to-soft pulses. However, the fact that the intensity-tracking pattern exists in the first pulse of multi-pulse long GRBs and some single-pulse GRBs suggest that intensity tracking is an independent component, which may operate in some late pulses as well. For the GRBs with measured redshifts, we present a time-resolved $E_{\rm p}-L_{\gamma, \rm iso}$ correlation analysis and show that the scatter of the correlation is comparable to that of the global Amati/Yonetoku relation. We discuss the predictions of various radiation models regarding $E_{\rm p}$ evolution, as well as the possibility of a precession jet in GRBs. It seems that the data pose great challenge to all these models, and hold the key to unveil the physics of GRB prompt emission.