Evidence for the central engine of gamma-ray bursts (GRBs) has been collected in the Neil Gehrels \emph{Swift} data. For instance, some GRBs show an internal X-ray plateau followed by very steep decay, which is difficult to be interpreted within the framework of a black hole (BH) central engine, but are consistent within a rapidly spinning magnetar engine picture. The very steep decay at the end of the plateau suggests a sudden cessation of the central engine, which is explained as the collapse of a supra-massive magnetar into a black hole when it spins down. Here we propose that some additional evidence, such as a second X-ray plateau feature would show up, if the fall-back accretion could activate the newborn BH and sufficient energy could be transferred from the newborn BH to the GRB blast wave. With a systematic data analysis for all long GRBs, we find three candidates in \emph{Swift} sample, i.e., GRBs 070802, 090111, and 120213A, whose X-ray afterglow light curves contain two plateaus, with the first one being an internal plateau. We find that in a fairly loose and reasonable parameter space, the second X-ray plateau data for all 3 GRBs could be well interpreted with our proposed model. Future observations are likely to discover more similar events, which could offer more information of the properties of the magnetar as well as the newborn BH.

A central compact object (CCO, e.g. a black hole) with an accretion disk has been suggested as the common central engine of various astrophysical phenomena, such as gamma-ray bursts (GRBs), tidal disruption events (TDEs) and active galactic nuclei (AGNs). A jet powered by such a system might precess due to the misalignment of the angular momenta of the CCO and accretion disk. Some quasi-periodic behaviors observed in the light curves of these phenomena can be well interpreted within the framework of a precessing jet model. In this paper, we study the emission polarization of precessing jets in the three kinds of phenomena. The polarization angle also shows a gradual change for the synchrotron emission in both the random and toroidal magnetic field configurations with the precessing jet, while it can only change abruptly by $90^\circ$ for the non-precessing top-hat jet. Polarization properties are periodic due to the assumptions made in our model. The polarization observations are crucial to confirm the precession nature of jets in GRBs, TDEs and AGNs.

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 central engine of gamma-ray bursts (GRBs) is poorly constrained. There exist two main candidates: a fast-rotating black hole and a rapidly spinning magnetar. Furthermore, the X-ray plateaus are widely accepted by the energy injection into the external shock. In this paper, we systematically analyze the \emph{Swift}/XRT light curves of 101 GRBs having plateau phases and known redshifts (before May 2017). Since a maximum energy budget ($\sim2\times10^{52}$ erg) exists for magnetars but not for a black hole, this provides a good clue to identify the type of GRB central engine. We calculate the isotropic kinetic energy $E_{\rm K,iso}$ and the isotropic X-ray energy release $E_{\rm X,iso}$ for individual GRB. We identify three categories based on how likely a black hole harbor at central engine: 'Gold' (9 out of 101, both $E_{\rm X,iso}$ and $E_{\rm K,iso}$ exceed the energy budget), 'Silver' (69 out of 101, $E_{\rm X,iso}$ less than the limit but $E_{\rm K,iso}$ is greater than the limit), and 'Bronze' (23 out of 101, the energies are not above the limit). We then derive and test the black hole parameters with the Blandford Znajek mechanism, and find that the observations of the black hole candidate ('Gold'+'Silver') samples are consistent with the expectations of the black hole model. Furthermore, we also test the magnetar candidate ('Bronze') sample with the magnetar model, and find that the magnetar surface magnetic led ($B_{p}$) and initial spin period ($P_{0}$) are consistent with the expectations of the magnetar model. Our analysis indicates that, if the magnetar wind is isotropic, a magnetar central engine is possible for 20\% of the analyzed GRBs. For most GRBs a black hole is most likely operating.

We offer an explanation for the anti-correlation between the minimum variability timescale ($MTS$) in the prompt emission light curve of gamma-ray bursts (GRBs) and the estimated bulk Lorentz factor of these GRBs, in the context of a magnetically arrested disk (MAD) model. In particular, we show that previously derived limits on the maximum available energy per baryon in a Blandford-Znajek jet leads to a relationship between the characteristic MAD timescale, $t_{MAD}$, in GRBs and the maximum bulk Lorentz factor: $t_{MAD} \propto \Gamma^{-6}$, somewhat steeper than (although within the error bars of) the fitted relationship found in the GRB data. Similarly, the MAD model also naturally accounts for the observed anti-correlation between $MTS$ and gamma-ray luminosity $L$ in the GRB data, and we estimate the accretion rates of the GRB disks (given these luminosities) in the context of this model. Both of these correlations ($MTS-\Gamma$ and $MTS-L$) are also observed in the AGN data, and we discuss the implications of our results in the context of both GRB and blazar systems.

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.

GRB 170817A is a weak short gamma-ray burst (GRB) accompanied by the gravitational wave (GW) event GW170817. It is believed, that an off beaming relativistic jet, produces this weak GRB. Here we use the $E_{\rm p,i}- E_{\rm iso}$ and $\Gamma - E_{\rm iso}$ relations to determine the Lorentz factor $\Gamma$ and the viewing angle from the edge of the jet $\theta'_{\rm obs}$ of GRB 170817A. Our result shows, $\Gamma = 13.4 ^{+9.8}_{-5.5}$ and $\theta'_{\rm obs} = 4.3^{+1.8}_{-1.5} \,^{\circ}$. This corresponds to the on-axis $E_{\rm p,i} = 415^{+361}_{-167}$ keV and $E_{\rm iso} = (2.4^{+1.6}_{-1.9}) \times 10^{47}$ ergs of the GRB, which is an intrinsically weak short GRB. Interestingly, the Doppler factor and the luminosity follow a universal relation from GRBs and blazars, which indicates they may share similar radiation mechanism.

Mergers of neutron star-neutron star (NS-NS) or neutron star-black hole (NS-BH) binaries are candidate sources of gravitational waves (GWs). At least a fraction of the merger remnant should be a stellar mass BH with a sub-relativistic ejecta. A collimated jet is launched via Blandford-Znajek mechanism from the central BH to trigger a short gamma-ray burst (SGRB). At the same time, a near-isotropic wind can be driven by the Blandford-Payne mechanism (BP). In previous work, additional energy injection to the ejecta from the BP mechanism was ignored, and radioactive decay has long been thought as the main source of energy for "kilonova" in the BH scenario. In this letter, we propose that the wind driven by BP mechanism from the new-born BH-disk can heat up and push the ejecta during the prompt emission phase or even at late time when there is mass fall-back. Such a BP-powered merger-nova could be bright at optical band even for a low-luminosity SGRB. The detection of a GW event with a merger product of BH, and accompanied by a bright merger-nova, would be a robust test of our model.

`Internal plateau' followed by a sharp decay is commonly seen in short gamma-ray burst (GRB) light curves. The plateau component is usually interpreted as the dipole emission from a supra-massive magnetar, and the sharp decay may imply the collapse of the magnetar to a black hole (BH). Fall-back accretion onto the new-born BH could produce long-lasting activities via the Blandford-Znajek (BZ) process. The magnetic flux accumulated near the BH would be confined by the accretion disks for a period of time. As the accretion rate decreases, the magnetic flux is strong enough to obstruct gas infall, leading to a magnetically-arrested disk (MAD). Within this scenario, we show that the BZ process could produce two types of typical X-ray light curves: type I exhibits a long-lasting plateau, followed by a power-law decay with slopes ranging from 5/3 to 40/9; type II shows roughly a single power-law decay with slope of 5/3. The former requires low magnetic filed strength, while the latter corresponds to relatively high values. We search for such signatures of the new-born BH from a sample of short GRBs with an internal plateau, and find two candidates: GRB 101219A and GRB 160821B, corresponding to type II and type I light curve, respectively. It is shown that our model can explain the data very well.

A hyperaccreting stellar-mass black hole has been proposed as the candidate central engine of gamma-ray bursts (GRBs). The rich observations of GRBs by \textit{Fermi} and \textit{Swift} make it possible to constrain the central engine model by comparing the model predications against data. This paper is dedicated to studying the temporal evolution of central engine parameters for both prompt emission and afterglow phases. We consider two jet launching mechanisms, i.e., $\nu\bar{\nu}$ annihilations and the Blandford-Znajek (BZ) processe, and obtain analytical solutions to these two models. We then investigate the black hole central engine parameters, such as the jet power, the dimensionless entropy $\eta$, and the central engine parameter $\mu_0=\eta (1+\sigma_0)$ (where $\sigma_0$ is the initial magnetization of the engine) at the base of the jet. The black hole may be spun up by accretion, or spun down by the BZ process, leaving imprints in GRB lightcurves. Usually, a BZ jet is more powerful and is likely responsible for the late time central engine activities. However, an initially non-spinning black hole central engine may first launch a thermal "fireball" via neutrino annihilations, and then launch a Poynting-flux-dominated jet via the BZ process. Multiple flares, giant bumps, and plateaus in GRB afterglows can be well produced as the result of late time accretion onto the black hole.

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.

[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.

The Einstein's equivalence principle can be tested by using parameterized post-Newtonian parameters, of which the parameter $\gamma$ has been constrained by comparing the arrival times of photons with different energies. It has been constrained by a variety of astronomical transient events, such as gamma-ray bursts (GRBs), fast radio bursts as well as pulses of pulsars, with the most stringent constraint of $\Delta \gamma \lesssim 10^{-15}$. In this letter, we consider the arrival times of lights with different circular polarization. For a linearly polarized light, it is the combination of two circularly polarized lights. If the arrival time difference between the two circularly polarized lights is too large, their combination may lose the linear polarization. We constrain the value of $\Delta \gamma_{\rm p} < 1.6 \times 10^{-27}$ by the measurement of the polarization of GRB 110721A, which is the most stringent constraint ever achieved.

Recently, two empirical correlations related to the minimum variability timescale ($\rm MTS$) of the lightcures are discovered in gamma-ray bursts (GRBs). One is the anti-correlation between $\rm MTS$ and Lorentz factor $\Gamma$, the other is the anti-correlation between the $\rm MTS$ and gamma-ray luminosity $L_\gamma$. Both the two correlations might be used to explore the activity of the central engine of GRBs. In this paper we try to understand these empirical correlations by combining two popular black hole (BH) central engine models (namely, Blandford \& Znajek mechanism and neutrino-dominated accretion flow). By taking the $\rm MTS$ as the timescale of viscous instability of the neutrino-dominated accretion flow (NDAF), we find that these correlations favor the scenario in which the jet is driven by Blandford-Znajek (BZ) mechanism.

We work on a GRB sample whose initial Lorentz factors ($\Gamma_0$) are constrained with the afterglow onset method and the jet opening angles ($\theta_{\rm j}$) are determined by the jet break time. We confirm the $\Gamma_0$ - $E_{\gamma,\rm iso}$ correlation by Liang et al. (2010), and the $\Gamma_0$ - $L_{\gamma,\rm iso}$ correlation by L{\"u} et al. (2012). Furthermore, we find correlations between $\Gamma_0$ and the beaming corrected $\gamma$-ray energy ($E_{\gamma}$) and mean $\gamma$-ray luminosity ($L_{\gamma}$). By also including the kinetic energy of the afterglow, we find rough correlations (with larger scatter) between $\Gamma_0$ and the total ($\gamma$-ray plus kinetic) energy and the total mean luminosity, both for isotropic values and beaming corrected values: these correlations allow us to test the data with GRB central engine models. Limiting our sample to the GRBs that likely have a black hole central engine, we compare the data with theoretical predictions of two types of jet launching mechanisms from BHs, i.e. the non-magnetized $\nu \bar\nu$ - annihilation mechanism, and the strongly magnetized Blandford-Znajek (BZ) mechanism. We find that the data are more consistent with the latter mechanism, and discuss the implications of our findings for GRB jet composition.

A stellar mass black hole (BH) surrounded by a neutrino-dominated accretion flow (NDAF) has been discussed in a number of works as the central engine of gamma-ray bursts (GRBs). It is widely believed that NDAF cannot liberate enough energy for bright GRBs. However, these works have been based on the assumption of "no torque" boundary condition, which is invalid when the disk is magnetized. In this paper, we present both numerical and analytical solutions for NDAFs with non-zero boundary stresses, and reexamine their properties. We find that NDAF with such boundary torque can be powerful enough to account for those bright short GRBs, energetic long GRBs and ultra-long GRBs. The disk becomes viscously unstable, which makes it possible to interpret the variability of GRB prompt emission and the steep decay phase in the early X-ray afterglow. Finally, we study the gravitational waves radiated from a processing BH-NDAF. We find that the effects of the boundary torque on the strength of the gravitational waves can be ignored.

Relativistic jets can form from at least some tidal disruption events (TDEs) of (sub-)stellar objects around supermassive black holes. We detect the millimeter (MM) emission of IGR J12580+0134 --- the nearest TDE known in the galaxy NGC 4845 at the distance of only 17 Mpc, based on Planck all-sky survey data. The data show significant flux jumps after the event, followed by substantial declines, in all six high frequency Planck bands from 100 GHz to 857 GHz. We further show that the evolution of the MM flux densities are well consistent with our model prediction from an off-axis jet, as was initially suggested from radio and X-ray observations. This detection represents the second TDE with MM detections; the other is Sw J1644+57, an on-axis jetted TDE at redshift of 0.35. Using the on- and off-axis jet models developed for these two TDEs as templates, we estimate the detection potential of similar events with the Large Millimeter Telescope (LMT) and the Atacama Large Millimeter/submillimeter Array (ALMA). Assuming an exposure of one hour, we find that the LMT (ALMA) can detect jetted TDEs up to redshifts $z\sim1$ (2), for a typical disrupted star mass of $\sim1$ M$_\odot$. The detection rates of on- and off-axis TDEs can be as high as $\sim0.6$ (13) and 10 (220) per year, respectively, for the LMT (ALMA). We briefly discuss how such observations, together with follow-up radio monitoring, may lead to major advances in understanding the jetted TDEs themselves and the ambient environment of the CNM.

Recently, the first association between an ultra-long gamma-ray burst (GRB) and a supernova is reported, i.e., GRB 111209A/SN 2011kl, which enables us to investigate the physics of central engines or even progenitors for ultra-long GRBs. In this paper, we inspect the broad-band data of GRB 111209A/SN 2011kl. The late-time X-ray lightcurve exhibits a GRB 121027A-like fall-back bump, suggesting a black hole central engine. We thus propose a collapsar model with fall-back accretion for GRB 111209A/SN 2011kl. The required model parameters, such as the total mass and radius of the progenitor star, suggest that the progenitor of GRB 111209A is more likely a Wolf-Rayet star instead of blue supergiant, and the central engine of this ultra-long burst is a black hole. The implications of our results is discussed.

Fast radio bursts (FRBs) are radio bursts characterized by millisecond durations, high Galactic latitude positions, and high dispersion measures. Very recently, the cosmological origin of FRB 150418 has been confirmed by \cite{kea16}, and FRBs are now strong competitors as cosmological probes. The simple sharp feature of the FRB signal is ideal for them to probe some of the fundamental laws of physics. Here we show that by analyzing the delay time between different frequencies, the FRB data can place stringent upper limits on the rest mass of the photon. For FRB 150418 at $z=0.492$, one can potentially reach $m_{\gamma} \leq 5.2 \times 10^{-47}$ g, which is $10^{20}$ times smaller than the rest mass of electron, and is about $10^3$ times smaller than that obtained using other astrophysical sources with the same method.

Supermassive black holes (SMBHs) can capture and tidally disrupt stars or sub-stellar objects orbiting nearby. The detections of Sw J1644+57-like events suggest that at least some TDEs can launch a relativistic jet beaming towards Earth. A natural expectation would be the existence of TDEs with a relativistic jet beaming away from Earth. The nearby TDE candidate IGR J12580+0134 provides new insights into the jet phenomenon. Combining several constraints, we find that the event invokes a $8-40$ Jupiter mass object tidally disrupted by a $3 \times 10^5 - 1.8 \times 10^7 M_\sun$ SMBH. Recently, a bright radio transient was discovered by Irwin et al. in association with IGR J12580+0134. We perform detailed modeling of the event based on a numerical jet model previously developed for the radio emission of Sw J1644+57. We find that the radio data of IGR J12580+0134 can be interpreted within an external forward shock model in the Newtonian regime. Using Sw J1644+57 as a template and properly correcting for its luminosity, we argue that the observed X-ray flux in early times is too faint to allow an on-beam relativistic jet unless the Lorentz factor is very small. Rather, the X-ray emission is likely from the disk or corona near the black hole. From various constraints, we find that the data are consistent with an off-beam relativistic jet with a viewing angle $\theta_{\rm obs} \gtrsim 30^{\rm o}$, and an initial Lorentz factor $\Gamma_j \gtrsim $ a few.This scenario can readily be tested in the upcoming VLBI observations.

The continued observations of Sw J1644+57 in X-ray and radio bands accumulated a rich data set to study the relativistic jet launched in this tidal disruption event. We find that the re-brightening feature in the radio light curve can be naturally explained by the two-component jet model. The possible origin of this structured jet are the Blandford-Znajek and Blandford-Payne mechanisms. We also show that this two-component jet model can interpret the two kinds of quasi-periodic variations in the X-ray light curve: a 200 second quasi-periodic oscillation (QPO) and a 2.7-day quasi-periodic variation. The latter is interpreted by a precessing outer jet launched near the Bardeen-Petterson radius of a warped disk. The $\sim$ 200s QPO could be associated with a second, narrower jet sweeping the observer line-of-sight periodically, which is launched from a spinning black hole in the misaligned direction with respect to the black hole's angular momentum.

Both gamma-ray bursts (GRBs) and blazars have relativistic jets pointing at a small angle from our line of sight. Several recent studies suggested that these two kinds of sources may share similar jet physics. In this work, we explore the variability properties for GRBs and blazars as a whole. We find that the correlation between minimum variability timescale (MTS) and Lorentz factor, $\Gamma$, as found only in GRBs by Sonbas et al. can be extended to blazars with a joint correlation of $\rm MTS\propto\Gamma^{-4.7\pm0.3}$. The same applies to the $\rm MTS\propto \it L_{\gamma}^{\rm -1.0\pm0.1}$ correlation as found in GRBs, which can be well extended into blazars as well. These results provide further evidence that the jets in these two kinds of sources are similar despite of the very different mass scale of their central engines. Further investigations of the physical origin of these correlations are needed, which can shed light on the nature of the jet physics.

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$.

The well-known correlation between the radio luminosity ($L_R$) and the X-ray luminosity ($L_X$) $L_R / L_X \simeq 10^{-5}$ holds for a variety of objects like active galactic nuclei, galactic black holes, solar flares and cool stars. Here we extend the relation to gamma-ray bursts (GRBs), and find the GRBs also lay on the same $L_R-L_X$ relation, with a slightly different slope as $L_R \propto L_X^{1.1}$. This relation implies the explosions in different scales may have a common underlying origin.

The continued observations of Sw J1644+57 in X-ray and radio bands accumulated a rich data set to study the relativistic jet launched in this tidal disruption event. The X-ray light curve of Sw J1644+57 from 5-30 days presents two kinds of quasi-periodic variations: a 200 second quasi-periodic oscillation (QPO) and a 2.7-day quasi-periodic variation. The latter has been interpreted by a precessing jet launched near the Bardeen-Petterson radius of a warped disk. Here we suggest that the $\sim$ 200s QPO could be associated with a second, narrower jet sweeping the observer line-of-sight periodically, which is launched from a spinning black hole in the misaligned direction with respect to the black hole's angular momentum. In addition, we show that this two-component jet model can interpret the radio light curve of the event, especially the re-brightening feature starting $\sim 100$ days after the trigger. From the data we infer that inner jet may have a Lorentz factor of $\Gamma_{\rm j} \sim 5.5$ and a kinetic energy of $E_{\rm k,iso} \sim 3.0 \times 10^{52} {\rm erg}$, while the outer jet may have a Lorentz factor of $\Gamma_{\rm j} \sim 2.5$ and a kinetic energy of $E_{\rm k,iso} \sim 3.0 \times 10^{53} {\rm erg}$.

The particular giant X-ray bump of GRB 121027A triggered by \emph{Swift} is quite different from the typical X-ray flares in gamma-ray bursts. There exhibit four parts of the observed structural variabilities in the rise and decay phase of the bump. Considering the quality of four parts of the data, we can only analyze the data from about 5300 s to about 6100 s in the bump using the stepwise filter correlation method (Gao et al. 2012), and find that the $86^{+5.9}_{-9.4}~\rm s$ periodic oscillation may exist, which is confirmed by the Lomb-Scargle method (Scargle 1982). Furthermore, a jet precession model (Liu et al. 2010) is proposed to account for such a variability.

Gamma-ray bursts are most luminous explosions in the universe. Their ejecta are believed to move towards Earth with a relativistic speed. The interaction between this "relativistic jet" and a circum burst medium drives a pair of (forward and reverse) shocks. The electrons accelerated in these shocks radiate synchrotron emission to power the broad-band afterglow of GRBs. The external shock theory is an elegant theory, since it invokes a limit number of model parameters, and has well predicted spectral and temporal properties. On the other hand, depending on many factors (e.g. the energy content, ambient density profile, collimation of the ejecta, forward vs. reverse shock dynamics, and synchrotron spectral regimes), there is a wide variety of the models. These models have distinct predictions on the afterglow decaying indices, the spectral indices, and the relations between them (the so-called "closure relations"), which have been widely used to interpret the rich multi-wavelength afterglow observations. This review article provides a complete reference of all the analytical synchrotron external shock afterglow models by deriving the temporal and spectral indices of all the models in all spectral regimes, including some regimes that have not been published before. The review article is designated to serve as a useful tool for afterglow observers to quickly identify relevant models to interpret their data. The limitations of the analytical models are reviewed, with a list of situations summarized when numerical treatments are needed.

Synchrotron self-Compton (SSC) scattering is an important emission mechanism in many astronomical sources, such as gamma-ray bursts (GRBs) and active galactic nuclei (AGNs). We give a complete presentation of the analytical approximations for the Compton scattering of synchrotron emission with both weak and strong synchrotron self-absorption. All possible orders of the characteristic synchrotron spectral breaks ($\nu_{\rm a}$, $\nu_{\rm m}$, and $\nu_{\rm c}$) are studied. In the weak self-absorption regime, i.e., $\nu_{\rm a} < \nu_c$, the electron energy distribution is not modified by the self-absorption process. The shape of the SSC component broadly resembles that of synchrotron, but with the following features: The SSC flux increases linearly with frequency up to the SSC break frequency corresponding to the self-absorption frequency $\nu_{\rm a}$; and the presence of a logarithmic term in the high-frequency range of the SSC spectra makes it harder than the power-law approximation. In the strong absorption regime, i.e. $\nu_{\rm a} > \nu_{\rm c}$, heating of low energy electrons due to synchrotron absorption leads to pile-up of electrons, and form a thermal component besides the broken power-law component. This leads to two-component (thermal + non-thermal) spectra for both the synchrotron and SSC spectral components. For $\nu_{\rm c} < \nu_{\rm a} < \nu_{\rm m}$, the spectrum is thermal (non-thermal) -dominated if $\nu_a > \sqrt{\nu_m \nu_c}$ ($\nu_a < \sqrt{\nu_m \nu_c}$). Similar to the weak-absorption regime, the SSC spectral component is broader than the simple broken power law approximation. We derive the critical condition for strong absorption (electron pile-up), and discuss a case of GRB reverse shock emission in a wind medium, which invokes $\nu_{\rm a} > {\rm max} (\nu_{\rm m}, \nu_{\rm c})$.

A particularly interesting discovery in observations of GRB 121027A is that of a giant X-ray bump detected by the Swift/X-Ray Telescope. The X-ray afterglow re-brightens sharply at about 1000 s after the trigger by more than two orders of magnitude in less than 200 s. This X-ray bump lasts for more than 10 ks. It is quite different from typical X-ray flares. In this Letter we propose a fall-back accretion model to interpret this X-ray bump within the context of the collapse of a massive star for a long-duration gamma-ray burst. The required fall-back radius of about 3.5e10 cm and mass of about 0.9-2.6 solar masses imply that a significant part of the helium envelope should survive through the mass loss during the last stage of the massive progenitor of GRB 121027A.

A hyper-accreting stellar-mass black hole has been long speculated as the best candidate of central engine of gamma-ray bursts (GRBs). Recent rich observations of GRBs by space missions such as Swift and Fermi pose new constraints on GRB central engine models. In this paper, we study the baryon loading processes of a GRB jet launched from a black hole central engine. We consider a relativistic jet powered by $\nu \bar\nu$-annihilation or by the Blandford-Znajek (BZ) mechanism. We consider baryon loading from a neutrino-driven wind from a neutrino-cooling-dominated accretion flow. For a magnetically dominated BZ jet, we consider neutron-drifting from the magnetic wall surrounding the jet and subsequent positron capture and proton-neutron inelastic collisions. The minumim baryon loads in both types of jet are calculated. We find that in both cases, a more luminous jet tends to be more baryon poor. A neutrino-driven "fireball" is typically "dirtier" than a magnetically dominated jet, while a magnetically dominated jet can be much cleaner. Both models have the right scaling to interpret the empirical $\Gamma-L_{\rm iso}$ relation discovered recently. Since some neutrino-driven jets have too much baryon loading as compared with the data, we suggest that at least a good fraction of GRBs should have a magnetically dominated central engine.

The X-ray transient source Sw J1644+57 recently discovered by Swift is believed to be triggered by tidal disruption of a star by a rapidly spinning supermassive black hole (SMBH). For such events, the outer disk is very likely misaligned with respect to the equatorial plane of the spinning SMBH, since the incoming star before disruption most likely has an inclined orbital plane. The tilted disk is subject to the Lense-Thirring torque, which tends to twist and warp due to the Bardeen-Petterson effect. The inner disk tends to align with the SMBH spin, while the outer region tends to remain in the stellar orbital plane, with a transition zone around the Bardeen-Petterson radius. The relativistic jet launched from the spinning SMBH would undergo precession. The 5-30 day X-ray lightcurve of Sw J1644+57 shows a quasi-periodic (2.7-day) variation with noticeable narrow dips. We numerically solve a warped disk and propose a jet-precessing model by invoking a Blandford-Znajek jet collimated by a wind launched near the Bardeen-Petterson radius. Through simulations, we show that the narrow dips in the X-ray lightcurve can be reproduced for a range of geometric configurations. From data we infer that the inclination angle of the initial stellar orbit is in the range of $10^{\circ}-20^{\circ}$ from the SMBH equatorial plane, that the jet should have a moderately high Lorentz factor, and that the inclination angle, jet opening angle, and observer's viewing angle are such that the duty cycle of the line-of-sight sweeping the jet cone is somewhat less than 0.5.

Soft extended emission (EE) following initial hard spikes up to 100 seconds was observed with {\em Swift}/BAT for about half of short-type gamma-ray bursts (SGRBs). This challenges the conversional central engine models of SGRBs, i.e., compact star merger models. In the framework of the black hole-neutron star merger models, we study the roles of the radial angular momentum transfer in the disk and the magnetic barrier around the black hole for the activity of SGRB central engines. We show that the radial angular momentum transfer may significantly prolong the lifetime of the accretion process and multiple episodes may be switched by the magnetic barrier. Our numerical calculations based on the models of the neutrino-dominated accretion flows suggest that the disk mass is critical for producing the observed EE. In case of the mass being $\sim 0.8M_{\odot}$, our model can reproduce the observed timescale and luminosity of both the main and EE episodes in a reasonable parameter set. The predicted luminosity of the EE component is lower than the observed EE with about one order of magnitude and the timescale is shorter than 20 seconds if the disk mass being $\sim 0.2M_{\odot}$. {\em Swift}/BAT-like instruments may be not sensitive enough to detect the EE component in this case. We argue that the EE component would be a probe for merger process and disk formation for compact star mergers.

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.

The bulk Lorentz factor of the gamma-ray burst (GRB) ejecta (Gamma_0) is a key parameter to understand the GRB physics. Liang et al. have discovered a correlation between Gamma_0 and isotropic gamma-ray energy: Gamma_0\proptoE_{gamma,iso,52}^{0.25}. By including more GRBs with updated data and more methods to derive Gamma_0, we confirm this correlation and obtain Gamma_0~91E_{gamma,iso,52}^{0.29}. Evaluating the mean isotropic gamma-ray luminosities L_{gamma,iso} of the GRBs in the same sample, we discover an even tighter correlation Gamma_0~249 L_{gamma,iso,52}^{0.30}. We propose an interpretation to this later correlation. Invoking a neutrino-cooled hyperaccretion disk around a stellar mass black hole as the central engine of GRBs, we derive jet luminosity powered by neutrino annihilation and baryon loading from a neutrino-driven wind. Applying beaming correction, we finally derive Gamma_0\proptoL_{gamma,iso}^{0.22}, which is well consistent with the data. This suggests that the central engine of long GRBs is likely a stellar mass black hole surrounded by a hyper-accreting disk.

We propose a two-component jet model consistent with the observations of several gamma ray bursts (GRBs) and active galactic nuclei (AGNs). The jet consists of inner and outer components, and they are supposed to be driven by the Blandford-Znajek (BZ) and Blandford-Payne (BP) processes, respectively. The baryons in the BP jet is accelerated centrifugally via the magnetic field anchored in the accretion disk. The BZ jet is assumed to be entrained a fraction of accreting matter leaving the inner edge of the accretion disk, and the baryons are accelerated in the conversion from electromagnetic energy to the kinetic energy. By fitting the Lorentz factors of some GRBs (GRB 030329, GRB 051221A, GRB 080413B) and AGNs (Cen A, Mkn 501 and Mkn 421) with this model, we constrain the physical parameters related to the accretion and outflow of these two kind of objects. We conclude that the spine/sheath structure of the jet from these sources can be interpreted naturally by the BZ and BP processes.

Recently a hard X-ray transient event, Sw J1644+57, was discovered by the Swift satellite. It likely marks the onset of a relativistic jet from a supermassive black hole, possibly triggered by a tidal disruption event. Another candidate in the same category, Sw J2058+05, was also reported. The low event rate suggests that only a small fraction of TDEs launch relativistic jets. A common speculation is that these rare events are related to rapidly spinning black holes. We attribute jet launching to the Blandford-Znajek mechanism, and use the available data to constrain the black hole spin parameter for the two events. It is found that the two black holes indeed carry a moderate to high spin, suggesting that black hole spin is likely the crucial factor behind the Sw J1644+57 - like events.

We discuss transportation and redistribution of energy and angular momentum in the magnetic connection (MC) process and Blandford-Payne (BP) process. The MC results in readjusting the interior viscous torque, and its effects are operative not only in but also beyond the MC region. The BP process is invoked to transfer the "excessive" angular momentum from an accretion disc. In addition, we derive a criterion for the interior viscous torque to resolve the puzzle of the overall equilibrium of angular momentum in disc accretion. It turns out that the BP efficiency of extracting angular momentum and the intensity of the outflow are required to be greater than some critical values.

Recent observations and theoretical work on gamma-ray bursts (GRBs) favor the central engine model of a Kerr black hole (BH) surrounded by a magnetized neutrino-dominated accretion flow (NDAF). The magnetic coupling between the BH and disk through a large-scale closed magnetic field exerts a torque on the disk, and transports the rotational energy from the BH to the disk. We investigate the properties of the NDAF with this magnetic torque. For a rapid spinning BH, the magnetic torque transfers enormous rotational energy from BH into the inner disk. There are two consequences: (i) the luminosity of neutrino annihilation is greatly augmented; (ii) the disk becomes thermally and viscously unstable in the inner region, and behaves S-Shape of the surface density versus accretion rate. It turns out that magnetically torqued NDAF can be invoked to interpret the variability of gamma-ray luminosity. In addition, we discuss the possibility of restarting the central engine to produce the X-ray flares with required energy.

We propose a model of magnetic connection (MC) of a black hole with its surrounding accretion disc based on large-scale magnetic field. The MC gives rise to transport of energy and angular momentum between the black hole and the disc, and the closed field lines pipe the hot matter evaporated from the disc, and shape it in the corona above the disc to form a magnetically induced disc-corona system, in which the corona has the same configuration as the large-scale magnetic field. We numerically solve the dynamic equations in the context of the Kerr metric, in which the large-scale magnetic field is determined by dynamo process and equipartition between magnetic pressure and gas pressure. Thus we can obtain a global solution rather than assuming the distribution of large-scale magnetic field beforehand. The main MC effects lie in three aspects. (1) The rotational energy of a fast-spinning black hole can be extracted, enhancing the dissipation in the accretion disc, (2) the closed field lines provide a natural channel for corona matter escaping from disc and finally falling into black hole, and (3) the scope of the corona can be bounded by the conservation of magnetic flux. We simulate the high-energy spectra of this system by using Monte-Carlo method, and find that the relative hardness of the spectra decreases as accretion rate or black hole spin $a_*$ increases. We fit the typical X-ray spectra of three black-hole binaries (GRO J1655-40, XTE 1118+480 and GX 339-4) in the low/hard or very high state.

We propose a three-stage model with Blandford-Znajek (BZ) and hyperaccretion process to interpret the recent observations of early afterglows of Gamma-Ray Bursts (GRBs). In the first stage, the prompt GRB is powered by a rotating black hole (BH) invoking the BZ process. The second stage is a quiet stage, in which the BZ process is shut off, and the accretion onto the BH is depressed by the torque exerted by the magnetic coupling (MC) process. Part of the rotational energy transported by the MC process from the BH is stored in the disk as magnetic energy. In the third stage, the MC process is shut off when the magnetic energy in the disk accumulates and triggers the magnetic instability. At this moment, the hyperaccretion process may onset, and the jet launched in this restarted central engine generates the observed X-ray flares. This model can account for energies and timescales of GRBs with X-ray flares observed in early afterglows.

Many models of gamma-ray bursts suggest a common central engine: a black hole of several solar masses accreting matter from a disk at an accretion rate from 0.01 to 10 M_\sun / s. The inner region of the disk is cooled by neutrino emission and large amounts of its binding energy were liberated, which could trigger the fireball. We improve the neutrino-dominated accreting flows by considering the effects of the magnetic fields, and find that more than half of the liberating energy can be extracted directly by the large-scale magnetic fields on the disk. And it turns out that the temperature of the disk is a bit lower than the neutrino-dominated accreting flows without magnetic field. Therefore, The outflows are magneticallydominated rather than neutrino dominated. In our model, neutrino mechanism can fuel some GRBs (not the brightest ones), but cannot fuel X-ray flares. However, the magnetic processes (both BZ and electromagnetic luminosity from a disk) are viable mechanisms for most of GRBs and the following X-ray flares.

An extreme Kerr black hole (BH) surrounded by a precessing disk is invoked to explain the light curves of gamma-ray bursts (GRBs) based on the coexistence of the Blandford-Znajek (BZ) and the magnetic coupling (MC) processes. The overall shape of the light curves and the duration of GRBs are interpreted by the evolution of the half-opening angle of the magnetic flux on the BH horizon, and the complex temporal structures are modulated by the precession and nutation of the jet powered by the BZ process. The time profile of the emission exhibits a fast rise and a slow decay due to the effect of the evolution of the half-opening angle. The light curves of several GRBs are well fitted by this model with only six free parameters.

A toy model for gamma-ray burst supernovae (GRB-SNe) is discussed by considering the effects of screw instability of magnetic field in black hole (BH) magnetosphere. The screw instability in the Blandford-Znajek (BZ) process (henceforth SIBZ) can coexist with the screw instability in the magnetic coupling (MC) process (henceforth SIMC). It turns out that both SIBZ and SIMC occur inevitably, provided that the following parameters are greater than some critical values, i.e., (i) the BH spin, (ii) the power-law index describing the magnetic field at the disk, and (iii) the vertical height of the astrophysical load above the equatorial plane of the rotating BH. The features of several GRBs are well fitted. In our model the durations of the long GRBs depend on the evolve time of the half-opening angle. A small fraction of energy is extracted from the BH via the BZ process to power a GRB, while a large fraction of energy is extracted from the BH via the MC process to power an associated supernova. In addition, the variability time scales of tens of msec in the light curves of the GRBs are fitted by two successive flares due to SIBZ.

A toy model for gamma-ray burst-supernovae (GRB-SN) is discussed by considering the coexistence of baryon poor outflows from black holes (BHs) and a powerful spin-connection to surrounding disk, giving rise to consistent calorimetry as described by van Putten (2001a) in a variant of the Blandford-Znajek process (BZ, 1977). In this model the half-opening angle of the magnetic flux tube on the horizon is determined by the mapping relation between the angular coordinate on the BH horizon and the radial coordinate on the surrounding accretion disk. The GRB is powered by the baryon poor outflows in the BZ process, and the associated SN is powered by very small fraction of the spin energy transferred from the BH to the disk in the magnetic coupling (MC) process. The timescale of the GRB is fitted by the duration of the open magnetic flux on the horizon. It turns out that the data of several GRB-SNe are well fitted with our model.

Screw instability of the magnetic field connecting a rotating black hole (BH) with its surrounding disk is discussed based on the model of the coexistence of the Blandford-Znajek (BZ) process and the magnetic coupling (MC) process (CEBZMC). A criterion for the screw instability with the state of CEBZMC is derived based on the calculations of the poloidal and toroidal components of the magnetic field on the disk. It is shown by the criterion that the screw instability will occur, if the BH spin and the power-law index for the variation of the magnetic field on the disk are greater than some critical values. It turns out that the instability occurs outside some critical radii on the disk. It is argued that the state of CEBZMC always accompanies the screw instability. In addtition, we show that the screw instability contributes only a small fraction of magnetic extraction of energy from a rotating BH.

Effects of magnetic coupling (MC) of a rotating black hole (BH) with its surrounding accretion disk are discussed in detail in the following aspects: (i) The mapping relation between the angular coordinate on the BH horizon and the radial coordinate on the disk is modified based on a more reasonable configuration of magnetic field, and a condition for coexistence of the Blandford-Znajek (BZ) and the MC process is derived. (ii) The transfer direction of energy and angular momentum in MC process is described equivalently by the co-rotation radius and by the flow of electromagnetic angular momentum and redshifted energy, where the latter is based on an assumption that the theory of BH magnetosphere is applicable to both the BZ and MC processes. (iii) The profile of the current on the BH horizon and that of the current density flowing from the magnetosphere onto the horizon are given in terms of the angular coordinate of the horizon. It is shown that the current on the BH horizon varies with the latitude of the horizon and is not continuous at the angular boundary between the open and closed magnetic field lines. (iv) The MC effects on disk radiation are discussed, and a very steep emissivity is produced by MC process, which is consistent with the recent XMM-Newton observation of the nearby bright Seyfert 1 galaxy MCG-6-30-15 by a variety of parameters of the BH-disk system.

An analytic model of rotating hot spot is proposed to explain kilohertz quasi-periodic oscillations (kHz QPO) in X-ray binaries, which is based on the magnetic coupling (MC) of a rotating black hole (BH) with its surrounding accretion disc. The hot spot in the inner region of the disc is produced by energy transferred from a spinning BH with non-axisymmetric magnetic field. The frequency and energy band of the hot spot turns out to be related to six parameters, by which the strength and the azimuthal angular region of the bulging magnetic field on the BH horizon, the mass and spin of the BH, and the power-law index of the magnetic field varying with the radial coordinate on the disc are described. In addition, the correlation of the fluctuation of the bulging magnetic field with the widths of QPO frequency is discussed.

The transfer of energy and angular momentum in the magnetic coupling (MC) of a rotating black hole (BH) with its surrounding accretion disc is discussed based on a mapping relation derived by considering the conservation of magnetic flux with two basic assumptions: (i)the magnetic field on the horizon is constant, (ii) the magnetic field on the disc surface varies as a power law with the radial coordinate of the disc. The following results are obtained: (i) the transfer direction of energy and angular momentum between the BH and the disc depends on the position of a co-rotation radius relative to the MC region on the disc, which is eventually determined by the BH spin; (ii) the evolution characteristics of a rotating BH in the MC process without disc accretion are depicted in a parameter space, and a series of values of the BH spin are given to indicate the evolution characteristics; (iii) the efficiency of converting accreted mass into radiation energy of a BH-disc system is discussed by considering coexistence of disc accretion and the MC process; (iv) the MC effects on disc radiation and emissivity index are discussed and it is concluded that they are consistent with the recent XMM-Newton observation of the nearby bright Seyfert 1 galaxy MCG-6-30-15 with reference to a variety of parameters of the BH-disc system.

An analytic model of the evolution of a rotating black hole (BH) is proposed by considering the coexistence of disk accretion with the Blandford-Znajek process. The evolutionary characteristics of the BH are described in terms of three parameters: the BH spin $a_*$, the ratio $k$ of the angular velocity of the magnetic field lines to that of the BH horizon and the parameter $\lambda$ indicating the position of the inner edge of the disk. It is shown that the ratio $k$ being a little greater than 0.5 affects the evolutionary characteristics of the BH significantly, and the BH spin increases rather than decreases in its evolutionary process provided that the initial value of the BH spin is located in an appropriate value range determined by ratio $k$. Our calculations show that the system of a BH accretion disk with $k=0.6$ might provide a much higher output energy in a shorter timescale for gamma-ray bursts than the same system with $k=0.5$.

Evolution characteristics of a Kerr black hole (BH) are investigated by considering coexistence of disc accretion with the Blandford-Znajek process (the BZ process). (i) The rate of extracting energy from the rotating BH in the BZ process and that in MC process are expressed by a unified formula, which is derived by using an improved equivalent cicuit. (ii) The mapping relation between the angular coordinate on the BH horizon and the radial coordinate on the disc is given in the context of general relativity and conservation of magnetic flux. (iii) The power and torque in the BZ process are compared with those in MC process in detail. (iv)Evolution characteristics of the BH and energy extracting efficiency are discussed by using the characteristics functions of BH evolution in the corresponding parameter space. (v) Power dissipation on the BH horizon and BH entropy increase are discussed by considering the coexistence of the above energy mechanisms.