Search

Rapidly rotating newborn magnetars, which originate from binary neutron star (NS) mergers and serve as the central engines of short gamma-ray bursts (GRBs), may leave some imprints on their prompt gamma-ray light curves even though they are far from their radiating fireballs. A high-frequency quasi-periodic oscillation (QPO) would be a unique feature for the magnetar central engine, especially a hypermassive magnetar. By conducting a systematic analysis of the prompt gamma-ray light curves from 605 short GRBs observed by {\em Fermi}/Gamma-ray Burst Monitor, we have identified such QPO signals in three GRBs (e.g. GRB 120323A, GRB 181222B, and GRB 190606A). The QPOs that peaked at $1258^{+6}_{-6}$ Hz for GRB 120323A, $623^{+4}_{-4}$ Hz for GRB 181222B, and $1410^{+4}_{-5}$ Hz for GRB 190606A are all with a confidence level above 5.2 $\sigma$. The high-frequency QPO signals of those three short GRBs may be caused by a hypermassive magnetar acting as the central engine in a binary NS merger of a binary NS., Comment: 11 pages, 1 Table, 7 Figures, accepted for publication in MNRAS, and matched with the published verison

The determination of the equation of state (EOS) of a neutron star (NS) and its maximum mass is very important for understanding the formation and properties of NSs under extreme conditions, but they remain open questions. Short-duration gamma-ray bursts (GRBs) are believed to originate from the merger of binary NSs or giant flares (GFs) of soft gamma repeaters (SGRs). Recently, the high-frequency quasi-periodic oscillations (QPOs) have been claimed to be identified from two short GRBs (GRB 931101B and GRB 910711). In this paper, we propose that the observed high-frequency QPOs in these two short GRBs result from torsional oscillations in the GFs of SGRs associated with cold NSs, or from radial oscillations of hypermassive NSs as the hot remnants of binary NS mergers, and then to constrain the EOS of NSs. For torsional oscillations, the six selected EOSs (TM1, NL3, APR, SLy4, DDME2, and GM1) of NSs suitable for the zero-temperature condition exhibit significant overlap in mass ranges, suggesting that we cannot constrain the EOS of NSs. For radial oscillations, the six selected EOSs (IUF, TM1, TMA, FSG, BHBLp, and NL3) of NSs suitable for the high-temperature condition cannot be ruled out when redshift is considered. However, it is found that the EOS can only be constrained if the redshift and temperature of the remnant can be measured., Comment: 16 pages, 2 Tables, 10 Figures, accepted for publication in ApJ, and matched with the published verison

Gravitational waves (GWs) from core-collapse supernovae (CCSNe) have been proposed as a means to probe the internal physical properties of supernovae. However, due to their complex time-frequency structure, effectively searching for and extracting GW signals from CCSNe remains an unsolved challenge. In this paper, we apply the improved multisynchrosqueezing transform (IMSST) method to reconstruct simulated GW data based on the advanced LIGO (aLIGO) and Einstein Telescope (ET) detectors. These data are generated by the magnetorotational and neutrino-driven mechanisms, and we use the match score as the criterion for evaluating the quality of the reconstruction. To assess whether the reconstructed waveforms correspond to true GW signals, we calculate the false alarm probability of reconstruction (FAPR). For GW sources located at 10 kpc and datasets where the waveform amplitudes are normalized to $5 \times 10^{-21}$ observed by aLIGO, FAPR are $2.1 \times 10^{-2}$ and $6.2 \times 10^{-3}$, respectively. For GW sources at 100 kpc and with waveform amplitudes normalized to $5 \times 10^{-21}$ observed by ET, FAPR are $1.3 \times 10^{-1}$ and $1.5 \times 10^{-2}$, respectively. When the gravitational wave strain reaches $7 \times 10^{-21}$ and the match score threshold is set to 0.75, the IMSST method achieves maximum reconstruction distances of approximately 37 kpc and 317 kpc for aLIGO and ET, respectively. Finally, we compared the performance of IMSST and STFT in waveform reconstruction based on the ET. The results show that the maximum reconstructable distance using STFT is 186 kpc.

Long-duration gamma-ray bursts (GRBs) are thought to be from core collapse of massive stars, and a rapidly spinning magnetar or black hole may be formed as the central engine. The extended emission in the prompt emission, flares and plateaus in X-ray afterglow, are proposed to be as the signature of central engine re-activity. However, the directly evidence from observations of identifying the central engines remain an open question. In this paper, we systemically search for long-duration GRBs that consist of bumps in X-ray afterglow detected by Swift/XRT, and find that the peak time of the X-ray bumps exhibit bimodal distribution (defined as early and late bumps) with division line at $t=7190$ s. Although we cannot rule out that such a bimodality arises from selection effects. We proposed that the long-duration GRBs with an early (or late) bumps may be originated from the fall-back accretion onto a new-born magnetar (or black hole). By adopting MCMC method to fit the early (or late) bumps of X-ray afterglow with the fall-back accretion of magnetar (or black hole), it is found that the initial surface magnetic filed and period of magnetars for most early bumps are clustered around $5.88\times10^{13}$ G and $1.04$ ms, respectively. Meanwhile, the derived accretion mass of black hole for late bumps is range of $[4\times10^{-4}, 1.8\times10^{-2}]~M_{\odot}$, and the typical fall-back radius is distributed range of $[1.04, 4.23]\times 10^{11}$ cm which is consistent with the typical radius of a Wolf-Rayet star. However, we also find that the fall-back accretion magnetar model is disfavored by the late bumps, but the fall-back accretion of black hole model can not be ruled out to interpret the early bumps of X-ray afterglow., Comment: 14 pages, 5 figures, accepted for publication in MNRAS

The prompt emission mechanism of gamma-ray bursts (GRBs) is a long-standing open question, and GRBs have been considered as potential sources of high-energy neutrinos. Despite many years of search for the neutrino events associated with GRBs from IceCube, there were no results. However, the absence of search results for neutrinos provides a unique opportunity to constrain the parameter space of GRB jet models. In this paper, we chose four peculiar GRBs with two different types of jet composition to investigate neutrino emission. It is found that only GRB 211211A could be well constrained within the dissipative photosphere model. By adopting the specific parameters of the photosphere, one can obtain $\varepsilon _{\rm p}/\varepsilon_{\rm e}<8$ for $f_{\rm p}>0.2$ from GRB 211211A. For the Internal-collision-induced Magnetic Reconnection and Turbulence (ICMART) model, we can effectively constrain neither GRB 230307A nor GRB 080916C. Moreover, we also investigate the detection prospects of high-energy neutrinos from GRBs and find that it is difficult to detect at least one high-energy neutrino associated with GRBs from the ICMART model even during the IceCube-Gen2 operation. For the GRB 211211A-like events, it is possible to detect at least one neutrino coincident with the gravitational wave during the IceCube-Gen2 operation, if such an event originated from mergers of compact stars within the photosphere dissipation., Comment: 15 pages, 10 figures, matched with the published verison

Short gamma-ray bursts (SGRBs) with extended emission (EE) are composed of initial main emission (ME) with a short-hard spike, followed by a long-lasting EE. Whether the ME and EE originated from the same origin or not, as well as the jet composition, remains an open question. In this paper, we present a systematic analysis of 36 gamma ray bursts (GRBs) in our sample, which are identified as the category of SGRBs with EE as observed by Fermi/Gamma-ray Burst Monitor. By extracting time-integrated spectra of ME and EE with cutoff power-law or Band models for our sample, we find that 20 out of 36 SGRBs have $\alpha$ values that exceed the death line (e.g., -2/3) of synchrotron emission within either ME or EE phases, and we suggest that the quasi-thermal component should exist in the prompt emission. Then, we extract the time-resolved spectra of our samples, but only four GRBs are bright enough to extract the time-resolved spectra. We find that both thermal and nonthermal emissions do exist in the prompt emission of those four bright GRBs, which suggests that a hybrid jet (e.g., matter and Poynting-flux outflow) in GRBs should exist. Moreover, strong positive correlations (e.g., $F_{\rm tot}-\Gamma$ and $F_{\rm tot}-kT$) in the time-resolved spectra of ME and EE for those four GRBs have been discovered. This indicates that the spectral evolution of both ME and EE seem to share similar behavior, possibly from the same physical origin., Comment: 30 pages, 9 figures and 11 tables, matched with the published verison

We present optical photometry for the afterglow of GRB 221009A, in some respects the most extraordinary gamma-ray burst (GRB) ever observed. Good quality in the R-band light curve is obtained, covering 0.32-19.57 days since the Fermi-GBM trigger. We find that a weak bump emerges fromthe declining afterglow at $t \approx 11$ days; a supernova (SN) may be responsible. We use a smooth broken power-law and $^{56}\mathrm{Ni}$ model to fit the light curve. The best-fitting results reveal that the SN ejected a total mass of $M_\mathrm{ej} = 3.70 M_\odot$, a $^{56}\mathrm{Ni}$ mass of $M_\mathrm{Ni} = 0.23 M_\odot$, and a kinetic energy of $E_\mathrm{SN,K} = 2.35 \times 10^{52} \mathrm{erg}$. We also compare GRB 221009A with other GRB-SN events based on a GRB-associated SN sample, and find that only SN 2003lw and SN 2011kl can be obviously revealed in the afterglow of GRB 221009A by setting these objects at its distance. This suggests that a supernova (SN 2022xiw) is possibly obscured by the brighter afterglow emission from GRB 221009A.

Magnetars are potential energy sources or central engines for numerous transient phenomena in the Universe. How newborn magnetars evolve in different environments remains an open question. Based on both observed and candidate magnetars, it is found that the periods of all magnetars or candidates appear as a bimodal distribution, and are defined as the ``long-P'' and ``short-P'' magnetar subclasses, respectively. We find that for the ``short-P'' subclass of magnetars, the $\dot{P}$ values also appear as a bimodal distribution, and therefore can be classified as ``high-$\dot{P}$ short-P'' and ``low-$\dot{P}$ short-P'' magnetar subclasses. In this paper, we use Monte Carlo simulations to generate synthetic magnetar populations and investigate the evolution of the ``high-$\dot{P}$ short-P'' and ``low-$\dot{P}$ short-P'' magnetar subclasses by considering both the magnetar spin and inclination, as well as the decay of their magnetic field within their evolution in both vacuum and plasma-filled magnetospheres. We find that the magnetar evolution is dependent on both spin and magnetic field, but seems to be insensitive to inclination evolution and magnetospheric environment for the ``high-$\dot{P}$ short-P'' sub-class. In comparison for the case of ``high-$\dot{P}$ short-P'', the magnetar evolution is dependent on spin, magnetic field, and inclination evolution, as well as the magnetospheric environment. The best evolution model should be the case of inclination evolution in vacuum with a small value of $\overline{\mathrm{FOM}}$. The differences in the best-fit parameters also suggest that the ``high-$\dot{P}$ short-P'' and ``low-$\dot{P}$ short-P'' magnetar subclasses may be tracking with different evolution channels., Comment: 18 pages, 2 Tables, and 9 Figures. PRD in press, and matched with the published verison

The nearby long-duration GRB 191019A recently detected by Swift lacks an associated supernova and belongs to a host galaxy with little star formation activity, suggesting that the origin of this burst is the result of a merger of two compact objects with dynamical interactions in a high-density medium of an active galactic nucleus. Given the potential motivation of this event, and given that it occurs in such a high-density environment, the ejecta-circumstellar medium (CSM) interaction cannot be ignored as possibly contributing to the kilonova emission. Here, we theoretically calculate the kilonova emission by considering the contribution of the ejecta-CSM interaction in a high-density environment. We find that the contribution to the kilonova emission from the ejecta-CSM interaction will dominate at a later time, and a smaller ejecta mass will have a stronger kilonova emission from the ejecta-CSM interaction. Moreover, we try to apply it to GRB 191019A, but we find that it is difficult to identify the possible kilonova emission from the observations, due to the contribution of the bright host galaxy. On the other hand, less injected mass (less than $M_{\rm ej}=2\times10^{-5}M_{\odot}$) will be required if one can detect the kilonova emission associated with a GRB 191019A-like event in the future. The {\em r}-process-powered and spin energy contributions from the magnetar are also discussed., Comment: 10 pages, 6 figures, updated the rederence list. ApJ in press, and matched with the published verison

The prompt emission, X-ray plateau, and X-ray flares of Gamma-ray bursts (GRB) are thought to be from internal dissipation, and the magnetar as the central engine with propeller fallback accretion is proposed to interpret the observed phenomena of GRBs. In this paper, by systematically searching for X-ray emission observed by Swift/Xry Telescope, we find that seven robust GRBs include both X-ray flares and plateau emissions with measured redshift. More interestingly, the X-ray flares/bumps for those seven GRBs are simultaneously observed in the gamma-ray band. By adopting the propeller fallback accretion model to fit the observed data, it is found that the free parameters of two GRBs (140512A and 180329B) can be constrained very well, while in the other five cases, more or less, they are not all sufficiently constrained. On the other hand, this requires that the conversion efficiency of the propeller is to be two or three times higher than that of the spindown dipole radiation of the magnetar. If this is the case, it is contradictory to the expectation from the propeller model: namely, a dirtier ejecta should be less efficient in producing gamma-ray emissions. Our results hint that at least the magnetar central engine with propeller fallback accretion model cannot interpret very well both the GRB X-ray flares simultaneously observed in the gamma-ray band and the X-ray flares of GRBs with a high Lorentz factor., Comment: 13 pages, 5 figures, 1 table, updated references, match with published version

Black hole-neutron star (BH-NS) mergers are expected to emit gravitational-wave (GW) and electromagnetic (EM) counterparts when the NS is tidally disrupted or plunges into the BH. Recently, GW 200105 and GW200115 were claimed as originating in BH-NS mergers, even GW 200105 remains in debate. Several optical source candidates are reported to possible associate with the two GW events, but not confirmed yet. In this work, we assume that the BH is charged (the NS is naturally charged) and try to constrain the charge of the BH by using the possible associated EM emission from the charged BH and NS system working in the inspiral regime. We adopt electric and magnetic dipole radiations for the binaries which power a Poynting-flux-dominated outflow to accelerate electrons. Then, it produces the observed EM radiation via synchrotron radiation. We find that the conversion efficiency in the X-ray band is much higher than that of the ultraviolet (UV), near-infrared, and radio bands. The estimated maximum charge-to-mass ratio (the charge for unit mass) of the BH is $1.12\times 10^{-6}$ and $1.53\times 10^{-6}$ esu for the binary systems of GW200105 and GW200115, respectively, if magnetic field strength $B_{p}\lesssim ~10^{16}$ G and period $P>~1$ ms for the NS spin., Comment: 13 pages, 4 Tables, and 2 Figures. Accepted for publication in PRD

The gravitational waves (GW) from core-collapse supernovae (CCSN) have been proposed as a probe to investigate physical properties inside of the supernova. However, how to search and extract the GW signals from core-collapse supernovae remains an open question due to its complicated time-frequency structure. In this paper, we apply the Ensemble Empirical Mode Decomposition (EEMD) method to decompose and reconstruct simulated GW data generated by magnetorotational mechanism and neutrino-driven mechanism within the advanced LIGO, using the match score as the criterion for assessing the quality of the reconstruction. The results indicate that by decomposing the data, the sum of the first six intrinsic mode functions (IMFs) can be used as the reconstructed waveform. To determine the probability that our reconstructed waveform corresponds to a real GW waveform, we calculate the false alarm probability of reconstruction (FAPR). By setting the threshold of the match score to be 0.75, we obtain FAPR of GW sources at a distance of 5 kpc and 10 kpc to be $6\times10^{-3}$ and $1\times10^{-2}$ respectively. If we normalize the maximum amplitude of the GW signal to $5\times10^{-21}$, the FAPR at this threshold is $4\times10^{-3}$. Furthermore, in our study, the reconstruction distance is not equivalent to the detection distance. When the strain of GW reaches $7 \times 10^{-21}$, and the match score threshold is set at 0.75, we can reconstruct GW waveform up to approximately 36 kpc., Comment: 9 pages, 6 figures. Accepted by MNRAS

The nature of the merger remnant of binary neutron star (BNS) remains an open question. From the theoretical point of view, one possible outcome is a supra-massive neutron star (SMNS), which is supported by rigid rotation and through its survival of hundreds of seconds before collapsing into a black hole (BH). If this is the case, the SMNS can emit continuous gravitational waves (GW) and electromagnetic (EM) radiation, particularly in the X-ray band. In this work, the ellipticity and initial frequency of SMNS are constrained with a Bayesian framework using simulated X-ray and GW signals, which could be detected by The Transient High Energy Sky and Early Universe Surveyor (THESEUS) and Einstein Telescope (ET), respectively. We found that only considering the X-ray emission can not completely constrain the initial frequency and ellipticity of the SMNS, but it can reduce the ranges of the parameters. Afterwards, we can use the posterior distribution of the X-ray parameter estimates as a prior for the GW parameter estimates. It was found that the 95$\%$ credible region of the joint X-ray-GW analysis was about $10^5$ times smaller than that of the X-ray analysis alone., Comment: Accepted by MNRAS

Long-duration GRB 211211A that lacks a supernova emission even down to very stringent limits at such a low redshift $z=0.076$ and is associated with kilonova emission, suggests that its physical origin is from a binary compact star merger. By reanalyzing its data observed with the Gamma-Ray Burst Monitor on board the Fermi mission, we find that both time-integrated and time-resolved spectra can be fitted well by using a 2SBPL plus blackbody (2SBPL+BB) model in the prompt emission. The bulk Lorentz factors ($\Gamma_{\rm ph}$) of the outflow can be inferred by invoking the observed thermal emission at the photosphere radius within a pure fireball model, and we find out that the temporal evolution of $\Gamma_{\rm ph}$ seems to be tracking with the light curve. The derived values of $\Gamma_{\rm ph}$ are also consistent with the $\Gamma_{\rm ph}$-$L_{\gamma, \rm iso}$/$E_{\gamma, \rm iso}$ correlations that had been found in other bursts. Moreover, we also calculate the magnetization factor $\sigma_{0}$ in the central engine and $\sigma_{\rm ph}$ at the photosphere radius within the framework of a hybrid jet model, and find that the values of both $1+\sigma_{\rm 0}$ and $1+\sigma_{\rm ph}$ are larger than 1 for different time slices. It suggests that at least the Poynting-flux component is indeed existent in the outflow. If this is the case, one possible physical interpretation of thermal and nonthermal emissions in GRB 211211A is from the contributions of both $\nu\bar{\nu}$ annihilation and the Blandford-Znajek mechanisms in the relativistic jet when a stellar mass black hole resides in the central engine., Comment: 12 pages, 6 figures, and 4 tables. ApJ in press, and matched with the published verison

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 $\sim 0.1$ s without extended emission (EE) 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 $E_{\rm p}=113^{+9}_{-7}$ keV, and isotropic energy $E_{\gamma,iso} =1.36^{+0.17}_{-0.14}\times 10^{51}~\rm erg$. In order to reveal the possible physical origin of GRB 201221D, we adopted multi-wavelength criteria (e.g., Amati relation, $\varepsilon$-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 $\hat{\alpha}$ is larger than $2+\hat{\beta}$ 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., Comment: 15 pages, 3 tables, 11 Figures. Accepted for published in RAA

Weak and continuous gravitational-wave (GW) radiation can be produced by newborn magnetars with deformed structure and is expected to be detected by the Einstein telescope in the near future. In this work we assume that the deformed structure of a nascent magnetar is not caused by a single mechanism but by multiple time-varying quadrupole moments such as those present in magnetically induced deformation, starquake-induced ellipticity, and accretion column-induced deformation. The magnetar loses its angular momentum through accretion, magnetic dipole radiation, and GW radiation. Within this scenario, we calculate the evolution of GWs from a newborn magnetar by considering the above three deformations. We find that the GW evolution depends on the physical parameters of the magnetar (e.g., period and surface magnetic field), the adiabatic index, and the fraction of poloidal magnetic energy to the total magnetic energy. In general the GW radiation from a magnetically induced deformation is dominant if the surface magnetic field of the magnetar is large, but the GW radiation from magnetar starquakes is more efficient when there is a larger adiabatic index if all other magnetar parameters remain the same. We also find that the GW radiation is not very sensitive to different magnetar equations of state., Comment: 12 pages, 6 figures, Physical Review D in press, and matched with the published version

It is generally believed that long-duration gamma-ray bursts (GRBs) are associated with massive star core-collapse, whereas short-duration GRBs are associated with mergers of compact star binaries. However, growing observations have suggested that oddball GRBs do exist, and multiple criteria (prompt emission properties, supernova/kilonova associations, and host galaxy properties) rather than burst duration only are needed to classify GRBs physically. A previously reported long-duration burst, GRB 060614, could be viewed as a short GRB with extended emission if it were observed at a larger distance and was associated with a kilonova-like feature. As a result, it belongs to the Type-I (compact star merger) GRB category and is likely of the binary neutron star merger origin. Here we report a peculiar long-duration gamma-ray burst, GRB 211211A, whose prompt emission properties in many aspects differ from all known Type-I GRBs, yet its multi-band observations suggest a non-massive-star origin. In particular, significant excess emission in both optical and near-infrared wavelengths has been discovered, which resembles kilonova emission as observed in some Type-I GRBs. These observations point towards a new progenitor type of GRBs. A scenario invoking a white dwarf-neutron star merger with a post-merger magnetar engine provides a self-consistent interpretation for all the observations, including prompt gamma-rays, early X-ray afterglow, as well as the engine-fed kilonova emission., Comment: Updated; Authors' version. 44 pages, 8 figures, 6 tables

Long-duration gamma-ray bursts (GRBs) associated with supernovae (SNe) are believed to originate from massive star core-collapse events, whereas short-duration GRBs that are related to compact star mergers are expected to be accompanied by kilonovae. GRB 211227A, which lasted about 84 s, had an initial short/hard spike followed by a series of soft gamma-ray extended emission at redshift $z=$0.228. We performed follow-up observations of the optical emission using BOOTES, LCOGT, and the Lijiang 2.4m telescope, but we detected no associated supernova signature, even down to very stringent limits at such a low redshift. We observed the host galaxy within a large error-circle and roughly estimate the physical offset of GRB 211227A as $20.47\pm14.47$ kpc from the galaxy center. These properties are similar to those of GRB 060614, and suggest that the progenitor of GRB 211227A is not favored to be associated with the death of massive stars. Hence, we propose that GRB 211227A originates from a compact star merger. Calculating pseudo-kilonova emission for this case by adopting the typical parameters, we find that any associated pseudo-kilonova is too faint to be detected. If this is the case, it explains naturally the characteristics of the prompt emission, the lack of SN and kilonova emission, and the large physical offset from the galaxy center., Comment: 12 pages, 7 figures, and 1 Table. Matched with the published verison

Assuming that the shallow-decaying phase in the early X-ray lightcurves of gamma-ray bursts (GRBs) is attributed to the dipole radiations (DRs) of a newborn magnetar, we present a comparative analysis for the magnetars born in death of massive stars and merger of compact binaries with long and short GRB (lGRB and sGRB) data observed with the {\em Swift} mission. We show that the typical braking index ($n$) of the magnetars is $\sim 3$ in the sGRB sample, and it is $\sim 4$ for the magnetars in the lGRB sample. Selecting a sub-sample of the magnetars whose spin-down is dominated by DRs ($n\lesssim 3$) and adopting a universal radiation efficiency of $0.3$, we find that the typical magnetic field strength ($B_p$) is $10^{16}$ G {\em vs.} $10^{15}$ G and the typical initial period ($P_0$) is $\sim 20$ ms {\em vs.} $2$ ms for the magnetars in the sGRBs {\em vs.} lGRBs. They follow the same relation between $P_0$ and the isotropic GRB energy as $ P_0\propto E_{\rm jet}^{-0.4}$. We also extend our comparison analysis to superluminous supernovae (SLSNe) and stable pulsars. Our results show that a magnetar born in merger of compact stars tends to have a stronger $B_p$ and a longer $P_0$ by about one order of magnitude than that born in collapse of massive stars. Its spin-down is dominated by the magnetic DRs as old pulsars, being due to its strong magnetic field strength, whereas the early spin-down of magnetars born in massive star collapse is governed by both the DRs and gravitational wave (GW) emission. A magnetar with a faster rotation speed should power a more energetic jet, being independent of its formation approach., Comment: 10 pages, 2 tables, 7 figures. Published in MNRAS

Observationally, there is a small fraction of Gamma-ray bursts (GRBs) with prompt emission observed by Fermi/GBM that are composed of two pulses. Occasionally, the distance to a GRB may be lensed when a high mass astrophysical object resides in the path between the GRB source and observer. In this paper, we describe GRB 200716C, which has a two-pulse emission and duration of a few seconds. We present a Bayesian analysis identifying gravitational lensing in both temporal and spectral properties, and calculate the time delay ($\Delta t\sim 1.92$ s) and magnification ($\gamma\sim 1.5$) between those two pulses based on the temporal fits. One can roughly estimate the lens mass to be about $2.4\times 10^{5}~M_{\odot}$ in the rest frame. We also calculate the false alarm probability for this detection to be about 0.07\% with trial factors, and a present-day number density of about $808 \rm~Mpc^{-3}$ with an energy density $\Omega\sim 1.4\times 10^{-3}$. If the first pulse of this GRB near the trigger time is indeed gravitationally echoed by a second pulse, GRB 200716C may be a short GRB candidate with extended emission., Comment: 22 pages, 2 tables, 7 figures, ApJL in press, added two references

One possible progenitor of short gamma-ray bursts (GRBs) is thought to be from a double neutron star (NS) merger, and the remnant of such a merger may be a supramassive NS, which is supported by rigid rotation and through its survival of hundreds of seconds before collapsing into a black hole (BH). If this is the case, an optical/infrared transient (namely merger-nova) is generated from the ejected materials and it is powered by radioactive decay from $r$-process, spin-down energy from a supramassive NS, as well as the magnetic wind from a newborn BH. In this paper, we systematically search for the signature of a supramassive NS central engine by analyzing the X-ray emission of short GRBs with internal plateau observed by {\em Swift}, and we find that five candidates of short GRBs have such feature with redshift measurement. Then, we calculate the possible merger-nova emission from those candidates for given the typical model parameters by considering the above three energy sources, and compare its brightness with the sensitivity of some optical telescopes. We find that the merger-nova emission of GRB 060801 in K-, r-, and U-bands with variations of $M_{\rm ej}$ ($10^{-4}-10^{-2} M_{\odot}$), $\kappa$ ($0.1-10 ~\rm cm^{2}~g^{-1}$), and $\beta$ ($0.1-0.3$) is very difficult to detect using the Vera C. Rubin, Panoramic Survey Telescope and Rapid Response System (Pan-STARRS), the Zwicky Transient Facility, and the Roman Space Telescope (Roman), except for the case of large ejecta mass $M_{\rm ej}=10^{-2} M_{\odot}$. However, we are very hopeful to detect the merger-nova emission of GRBs 090515, 100625A, and 101219A using more sensitive instruments, such as Vera C. Rubin, Pan-STARRS, and Roman. Moreover, the merger-nova emission of GRB 160821B is bright enough to detect in our calculations, and it is also consistent with current real observations of merger-nova emission., Comment: 21 pages, 1 table, 6 figures. ApJ in press, two more references added, matched the published version

The merger of binary neutron stars (NS-NS) as the progenitor of short Gamma-ray bursts (GRBs) has been confirmed by the discovery of the association of the gravitational wave (GW) event GW170817 with GRB 170817A. However, the merger product of binary NS remains an open question. An X-ray plateau followed by a steep decay ("internal plateau") has been found in some short GRBs, implying that a supra-massive magnetar operates as the merger remnant and then collapses into a newborn black hole (BH) at the end of the plateau. X-ray bump or second-plateau following the "internal plateau" is considered as the expected signature from the fallback accretion onto this newborn BH through Blandford-Znajek mechanism (BZ). At the same time, a nearly isotropic wind driven by Blandford-Paynemechanism (BP) from the newborn BH's disk can produce a bright kilonova. Therefore, the bright kilonova observation for a short GRB with "internal plateau" (and followed by X-ray bump or second-plateau) provides further evidence for this scenario. In this paper, we find that GRB 160821B is a candidate of such a case, and the kilonova emission of GRB 160821B is possibly powered by the BP wind from a newborn BH. Future GW detection of GRB 160821B-like events may provide further support to this scenario, enable us to investigate the properties of the magnetar and the newborn BH, and constrain the equation of state of neutron stars., Comment: 6 pages, 3 figures. Accepted for publication in ApJ

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 model can adequately 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 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 a millisecond magnetar, which loses its rotation energy via GW quadrupole radiation. Then, the abrupt drop decay is caused by the magnetar collapsing into a black hole before switching to EM-dominated emission. This is the first short GRB for which the X-ray emission has such an intriguing feature powered by a magnetar via GW-dominated radiation. If this is the case, one can estimate the physical parameters of a magnetar, the GW signal powered by a magnetar and the merger-nova emission are also discussed., Comment: 21 pages, 2 Tables, 7 figures, ApJL in press, two more references added, matched the published version

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., Comment: 7 pages, 1 table, and 4 figures. ApJ in press, matched the published version

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., Comment: 12 pages, 6 figures, and 1 table. MNRAS in press, matched the published version

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., Comment: 23 pages, 12 figures. Accepted by ApJ

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., Comment: 7 pages, 2 figures. Accepted for publication in MNRAS

Binary neutron star (NS) mergers may result in remnants of supra-massive or even stable NS, which have been supported indirectly by observed X-ray plateau of some gamma-ray bursts (GRBs) afterglow. Recently, Xue et al. (2019) discovered a X-ray transient CDF-S XT2 that is powered by a magnetar from merger of double NS via X-ray plateau and following stepper phase. However, the decay slope after the plateau emission is a little bit larger than the theoretical value of spin-down in electromagnetic (EM) dominated by losing its rotation energy. In this paper, we assume that the feature of X-ray emission is caused by a supra-massive magnetar central engine for surviving thousands of seconds to collapse black hole. Within this scenario, 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 adopting 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 range in $[0.32-1.3]\times 10^{-3}$. Its GW signal can not be detected by aLIGO or even for more sensitive Einstein Telescope in the future., Comment: 13 pages, 5 figures,1 table. Accepted for publication by Research in Astronomy and Astrophysics

Gamma-ray burst (GRB) central engines and jet production mechanisms are still open questions. Assuming that the shallow decay segments of canonical X-ray afterglow lightcurves of {\em Swift} GRBs are attributed to the magnetic dipole (MD) radiations of newly-born magnetars, we derive the parameters of the magnetars and explore their possible relations to jet and MD wind emission. We show that the magnetar initial spin period ($P_0$) are tightly correlated with the jet energy ($E_{\rm jet}$), which is almost proportional to the wind energy ($E_{\rm wind}$). Our least square fits yield $P_0\propto E^{-0.36\pm 0.03}_{\rm jet}$ and $E_{\rm wind}\propto E^{0.91\pm 0.07}_{\rm jet}$. These relations may imply that a magnetar with faster rotating speed can power a more energetic GRB, and energy partition between the jet and wind may be quasi-universal. Although the $P_0-E_{\rm jet}$ relation is driven by a few sub-energetic GRBs in our sample, our Monte Carlo simulation analysis shows that sample selection biases from instrumental flux limits and contaminations of the bright jet afterglows cannot make this correlation. Within this jet-wind paradigm, we propose that GRB 101225A-like ultra-long GRBs, whose prompt gamma-ray/X-ray lightcurves are featured as a long-lasting plateau with a sharp drop, may be the orphan MD wind emission being due to misalignment of their jet axis to the light of sight. Brief discussion on the orphan MD wind emission and its association with the gravitational wave radiation of newly-born magnetars is presented., Comment: 19 pages, 2 tables, 12 figures, accepted for publication in ApJ

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., Comment: 8 pages, 3 figures. ApJ in press, matched the published version

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., Comment: Accepted for publication in the ApJ, 23 pages, 6 figures

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., Comment: 7 pages, 2 figures, Accepted for publication in MNRAS

The spectral evolution and spectral lag behavior of 92 bright pulses from 84 gamma-ray bursts (GRBs) observed by the Fermi GBM telescope are studied. These pulses can be classified into hard-to-soft pulses (H2S, 64/92), H2S-dominated-tracking pulses (21/92), and other tracking pulses (7/92). We focus on the relationship between spectral evolution and spectral lags of H2S and H2S-dominated-tracking pulses. %in hard-to-soft pulses (H2S, 64/92) and H2S-dominating-tracking (21/92) pulses. The main trend of spectral evolution (lag behavior) is estimated with $\log E_p\propto k_E\log(t+t_0)$ (${\hat{\tau}} \propto k_{\hat{\tau}}\log E$), where $E_p$ is the peak photon energy in the radiation spectrum, $t+t_0$ is the observer time relative to the beginning of pulse $-t_0$, and ${\hat{\tau}}$ is the spectral lag of photons with energy $E$ with respect to the energy band $8$-$25$ keV. For H2S and H2S-dominated-tracking pulses, a weak correlation between $k_{{\hat{\tau}}}/W$ and $k_E$ is found, where $W$ is the pulse width. We also study the spectral lag behavior with peak time $t_{\rm p_E}$ of pulses for 30 well-shaped pulses and estimate the main trend of the spectral lag behavior with $\log t_{\rm p_E}\propto k_{t_p}\log E$. It is found that $k_{t_p}$ is correlated with $k_E$. We perform simulations under a phenomenological model of spectral evolution, and find that these correlations are reproduced. We then conclude that spectral lags are closely related to spectral evolution within the pulse. The most natural explanation of these observations is that the emission is from the electrons in the same fluid unit at an emission site moving away from the central engine, as expected in the models invoking magnetic dissipation in a moderately-high-$\sigma$ outflow., Comment: 58 pages, 11 figures, 3 tables. ApJ in press

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., Comment: 22 pages, 1 Table, 5 Figures. updated the reference list and matched the published version

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., Comment: 33 pages, 3 Tables, 11 Figures, updated the reference list and matched the published version

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., Comment: 37 pages, 3 tables, 12 figures

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., Comment: 6 pages, 2 figures, accepted for publication in ApJ

[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., Comment: Accepted for publication in ApJ

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., Comment: 18 pages, 10 figures, 3 tables, accepted for publication in ApJ

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., Comment: 28 pages, 14 figures, and 2 Tables. Accepted for publication in ApJ

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$)., Comment: 10 pages, 7 figures, Accepted for publication in the ApJ

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., Comment: 36 pages, 6 figures. To be published in Space Science Review and Chapter 1 of ISSI-Beijing book "Gamma-ray bursts: a tool to explore the young universe"

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., Comment: Accepted by ApJS

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., Comment: Match with published version: 2017, ApJ, 873, 50

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)., Comment: 5 pages, 2 figures. 1 Table, Accepted for publication in MNRAS

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., Comment: 33 pages, 18 figures, 2 tables. Accepted by ApJ

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., Comment: Accepted for publication in Physics Review D

Relations between radio surface brightness ($\Sigma$) and diameter ($D$) of supernova remnants (SNRs) are important in astronomy. In this paper, following the work Duric \& Seaquist (1986) at adiabatic phase, we carefully investigate shell-type supernova remnants at radiative phase, and obtain theoretical $\Sigma$-$D$ relation at radiative phase of shell-type supernova remnants at 1 GHz. By using these theoretical $\Sigma$-$D$ relations at adiabatic phase and radiative phase, we also roughly determine phases of some supernova remnant from observation data., Comment: 12 pages, 2 tables, major revision, version accepted by Journal of Astrophysics and Astronomy

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$., Comment: 21 pages, 4 tables, 11 figures. Version published in ApJ, 805, 89

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., Comment: 19 pages, 41 figures, 3 Tables. ApJ, in press