Using archival {\it Chandra} observations with a total effective exposure of 323 ks, we derive an updated catalog of point sources in the bulge globular cluster M28. The catalog contains 502 X-ray sources within an area of $\sim475\, \rm arcmin^{2}$, and more than $90\%$ of these sources are first detected in this cluster. We find significant dips in the radial distribution profiles of X-ray sources in M28, with the projected distance and width of the distribution dip for bright ($L_{X} \gtrsim 4.5\times 10^{30} {\rm\ erg\ \,s^{-1}}$) X-ray sources are larger than the faint ($L_{X} \lesssim 4.5\times 10^{30} {\rm\ erg\ \,s^{-1}}$) sources. The "generalized King model" fitting give a slightly larger average mass for the bright sources ($1.30\pm0.15\,M_{\odot}$) than the faint sources ($1.09\pm0.14\,M_{\odot}$), which support a universal mass segregation delay between heavy objects in GCs. Compared with 47 Tuc and Terzan 5, we show that the dynamical age of M28 is comparable to Terzan 5 and much smaller than 47 Tuc, but it is evolving more fast (i.e., with smaller two-body relaxation timescale) than 47 Tuc. These features may suggest an acceleration effect of cluster dynamical evolution by tidal shock in M28. Besides, we find an abnormal deficiency of X-ray sources in the central region ($R \lesssim 1.5 \rm~arcmin$) of M28 than its outskirts, which indicate that M28 may have suffered an early phase of primordial binary disruption within its central region, and mass segregation effect will erase such a phenomenon as cluster evolve to older dynamical age.

Most dynamically confirmed stellar-mass black holes and the candidates were originally selected from X-ray outbursts. In the present work, we search for black hole candidates in the LAMOST survey by using the spectra along with photometry from the ASAS-SN survey, where the orbital period of the binary may be revealed by the periodic light curve, such as the ellipsoidal modulation type. Our sample consists of 9 binaries, where each source contains a giant star with large radial velocity variation ($\Delta V_{\rm R} > 70~{\rm km~s^{-1}}$) and periods known from light curves. We focus on the 9 sources with long periods ($T_{\rm ph} > 5$ days) and evaluate the mass $M_2$ of the optically invisible companion. Since the observed $\Delta V_{\rm R}$ from only a few repeating spectroscopic observations is a lower limit of the real amplitude, the real mass $M_2$ can be significantly higher than the current evaluation. It is likely an efficient method to place constraints on $M_2$ by combining $\Delta V_{\rm R}$ from LAMOST and $T_{\rm ph}$ from ASAS-SN, particularly by the ongoing LAMOST Medium Resolution Survey.

We propose a method to search for stellar-mass black hole (BH) candidates with giant companions from spectroscopic observations. Based on the stellar spectra of LAMOST Data Release 6, we obtain a sample of seven giants in binaries with large radial velocity variation $\Delta V_R > 80~{\rm km~s^{-1}}$. With the effective temperature, surface gravity, and metallicity provided by LAMOST, and the parallax given by {\it Gaia}, we can estimate the mass and radius of the giant, and therefore evaluate the possible mass of the optically invisible star in the binary. We show that the sources in our sample are potential BH candidates, and are worthy of dynamical measurement by further spectroscopic observations. Our method may be particularly valid for the selection of BH candidates in binaries with unknown orbital periods.

The occurrence rates of bright X-ray flares in z<1 gamma-ray bursts (GRBs) with or without observed supernovae (SNe) association were compared. Our Sample I: the z<1 long GRBs (LGRBs) with SNe association (SN-GRBs) and with early Swift/X-Ray Telescope (XRT) observations, consists of 18 GRBs, among which only two GRBs have bright X-ray flares. Our Sample II: for comparison, all the z<1 LGRBs without observed SNe association and with early Swift/XRT observations, consists of 45 GRBs, among which 16 GRBs present bright X-ray flares. Thus, the study indicates a lower occurrence rate of bright X-ray flares in Sample I (11.1%) than in Sample II (35.6%). In addition, if dim X-ray fluctuations are included as flares, then 16.7% of Sample I and 55.6% of Sample II are found to have flares, again showing the discrepancy between these two samples. We examined the physical origin of these bright X-ray flares and found that most of them are probably related to the central engine reactivity. To understand the discrepancy, we propose that such a lower occurrence rate of flares in the SN-GRB sample may hint at an energy partition among the GRB, SNe, and X-ray flares under a saturated energy budget of massive star explosion.

Short gamma-ray bursts may originate from the merger of double neutron stars (NS) or that of a black hole (BH) and an NS. We propose that the bright X-ray flare related to the central engine reactivity may hint a BH-NS merger, since such a merger can provide more fall-back materials and therefore a more massive accretion disk than the NS-NS merger. Based on the observed 49 short bursts with Swift/X-ray Telescope follow-up observations, we find that three bursts have bright X-ray flares, among which three flares from two bursts are probably related to the central engine reactivity. We argue that these two bursts may originate from the BH-NS merger rather than the NS-NS merger. Our suggested link between the central engine-powered bright X-ray flare and the BH-NS merger event can be checked by the future gravitational wave detections from advanced LIGO and Virgo.

Outflows launched from a remnant disc of compact binary merger may have essential contribution to the kilonova emission. Numerical calculations are conducted in this work to study the structure of accretion flows and outflows. By the incorporation of limited-energy advection in the hyper-accretion discs, outflows occur naturally from accretion flows due to imbalance between the viscous heating and the sum of the advective and radiative cooling. Following this spirit, we revisit the properties of the merger outflow ejecta. Our results show that around $10^{-3} \sim 10^{-1} M_\odot$ of the disc mass can be launched as powerful outflows. The amount of unbound mass varies with the disc mass and the viscosity. The outflow-contributed peak luminosity is around $10^{40} \sim 10^{41}$ erg/s. Such a scenario can account for the observed kilonovae associated with short gamma-ray bursts, including the recent event AT2017gfo (GW170817).

Two plateaus and one following bump in the X-ray light curve of GRB 170714A have been detected by the \textit{Swift}/X-Ray Telescope, which could be very meaningful for the central engine of gamma-ray bursts (GRBs), implying that the origin of this burst might be different from that of other ultra-long GRBs. We propose that merging two neutron stars into a hyper-massive quark star (QS) and then collapsing into a black hole (BH), with a delay time around $10^4$~s, could be responsible for those X-ray components. The hyper-massive QS is initially in a fluid state, being turbulent and differentially rotating, but would be solidified and release its latent heat injected into the GRB fireball (lasting about $10^3$~s during the liquid-solid phase transition). Magnetic field as high as $\sim 10^{15}$~G could be created by dynamo action of the newborn liquid QS, and a magnetar-like central engine (after solidification) supplies significant energy for the second plateau. More energy could be released during a fall-back accretion after the post-merger QS collapses to a BH, and the X-ray bump forms. This post-merger QS model might be tested by future observations, with either advanced gravitational wave detectors (e.g., advanced LIGO and VIRGO) or X-ray/optical telescopes.

Central engine of gamma-ray bursts (GRBs) may be intermittent and launch several episodes of ejecta separated by a long quiescent interval. In this scenario, an external shock is formed due to the propagation of the first launched ejecta into the circum-burst medium and the later launched ejecta may interact with the external shock at later period. Owing to the internal dissipation, the later launched ejecta may be observed at a later time ($t_{\rm{jet}}$). In this paper, we study the relation of $t_{\rm{b}}$ and $t_{\rm{jet}}$, where $t_{\rm{b}}$ is the collision time of the later launched ejecta with the formed external shock. It is found that the relation of $t_{\rm{b}}$ and $t_{\rm{jet}}$ depends on the bulk Lorentz factor ($\Gamma_{\rm{jet}}$) of the later launched ejecta and the density ($\rho$) of the circum-burst medium. If the value of $\Gamma_{\rm{jet}}$ or $\rho$ is low, the $t_{\rm{b}}$ would be significantly larger than $t_{\rm{jet}}$. However, the $t_{\rm{b}}\sim t_{\rm{jet}}$ can be found if the value of $\Gamma_{\rm{jet}}$ or $\rho$ is significantly large. Our results can explain the large lag of the optical emission relative to the $\gamma$-ray/X-ray emission in GRBs, e.g., GRB~111209A. For GRBs with a precursor, our results suggest that the energy injection into the external shock and thus more than one external-reverse shock may appear in the main prompt emission phase. According to our model, we estimate the Lorentz factor of the second launched ejecta in GRB~160625B.

We derive a simple analytical formula to describe the evolution of spectral index $\beta$ in the steep decay phase shaped by the curvature effect with assumption that the spectral parameters and Lorentz factor of jet shell is the same for different latitude. Here, the value of $\beta$ is estimated in 0.3$-$10keV energy band. For a spherical thin shell with a cutoff power law (CPL) intrinsic radiation spectrum, the spectral evolution can be read as a linear function of observer time. For the situation with Band function intrinsic radiation spectrum, the spectral evolution may be complex. If the observed break energy of radiation spectrum is larger than 10keV, the spectral evolution is the same as that shaped by jet shells with a CPL spectrum. If the observed break energy is less than 0.3keV, the value of $\beta$ would be a constant. Others, the spectral evolution can be approximated as a logarithmal function of the observer time in generally.

The curvature effect may be responsible for the steep decay phase observed in gamma-ray bursts. For testing the curvature effect with observations, the zero time point $t_0$ adopted to plot observer time and flux on a logarithmic scale should be appropriately selected. In practice, however, the true $t_0$ cannot be directly constrained from the data. Then, we move $t_0$ to a certain time in the steep decay phase, which can be easily identified. In this situation, we derive an analytical formula to describe the flux evolution of the steep decay phase. The analytical formula is read as $F_\nu\propto (1+\tilde t_{\rm obs}/{\tilde t_c})^{-\alpha}$ with $\alpha(\tilde{t}_{\rm obs})=2+{\int_{0}^{\log (1+\tilde{t}_{\rm obs}/{\tilde{t}_c})} {\beta(\tau)d[\log(1+\tau/{\tilde{t}_c})]}}/{\log (1 + {\tilde t}_{\rm obs}/{{\tilde t}_c})}$, where $F_\nu$ is the flux observed at frequency $\nu$, $\tilde t_{\rm obs}$ is the observer time by setting zero time point $t_0$ at a certain time in the steep decay phase, $\beta$ is the spectral index estimated around $\nu$, and ${\tilde t}_c$ is the decay timescale of the phase with $\tilde{t}_{\rm obs}{\geqslant}0$. We test the analytical formula with the data from numerical calculations. It is found that the analytical formula presents a well estimation about the evolution of flux shaped by the curvature effect. Our analytical formula can be used to confront the curvature effect with observations and estimate the decay timescale of the steep decay phase.

Recent simulations on super-Eddington accretion flows have shown that, apart from the diffusion process, the vertical advection based on magnetic buoyancy can be a more efficient process to release the trapped photons in the optically thick disk. As a consequence, the radiative luminosity from the accretion disk can be far beyond the Eddington value. Following this spirit, we revisit the structure and radiation of hyper-accretion disks with mass accretion rates in the range $10^{-3}\sim 10~M_{\sun}~{\rm s}^{-1}$. Our results show that, due to the strong cooling through the vertical advection, the disk temperature becomes lower than that in the classic model without the vertical advection process, and therefore the neutrino luminosity from the disk is lower. On the other hand, the gamma-ray photons released through the vertical advection can be extremely super-Eddington. We argue that the large amount of escaped gamma-ray photons may have more significant contribution to the primordial fireball than the neutrino annihilation, and may hint a link between gamma-ray bursts and kilonovae in the black hole hyper-accretion scenario.

We present time-resolved spectral analysis of the steep decay segments of 29 bright X-ray flares of gamma-ray bursts (GRBs) observed with the Swift/X-ray telescope, and model their lightcurves and spectral index evolution behaviors with the curvature effect model. Our results show that the observed rapid flux decay and strong spectral index evolution with time can be well fit with this model, and the derived characteristic timescales ($t_c$) are in the range of $33\sim 264$ seconds. Using an empirical relation between the peak luminosity and the Lorentz factor derived from the prompt gamma-rays, we estimate the Lorentz factors of the flares ($\Gamma_{\rm X}$). We obtain $\Gamma_{\rm X}=17\sim 87$ with a median value of $52$, which is smaller than the initial Lorentz factors of prompt gamma-ray fireballs. With the derived $t_c$ and $\Gamma_{\rm X}$, we constrain the radiating regions of 13 X-ray flares, yielding $R_{\rm X}=(0.2\sim 1.1)\times 10^{16}$ cm, which are smaller than the radii of the afterglow fireballs at the peak times of the flares. A long evolution feature from prompt gamma-ray phase to the X-ray epoch is found by incorporating our results with a sample of GRBs whose initial Lorentz factors are available in literatures, i.e., $\Gamma\propto [t_{p}/(1+z)]^{-0.69\pm 0.06}$. These results may shed lights on the long term evolution of GRB central engines.

This work focuses on a sample of seven extremely late-time X-ray flares with peak time $t_{\rm p} > 10^4 {\rm s}$, among which two flares can be confirmed as the late-time activity of central engine. The main purpose is to investigate the mechanism of such late-time flares based on the internal origin assumption. In the hyper-accreting black hole (BH) scenario, we study the possibility of two well-known mechanisms as the central engine to power such X-ray flares, i.e., the neutrino-antineutrino annihilation and the Blandford-Znajek (BZ) process. Our results show that the annihilation luminosity is far below the observational data. Thus, the annihilation mechanism cannot account for such late-time flares. For the BZ process, if the role of outflows is taken into consideration, the inflow mass rate near the horizon will be quite low such that the magnetic field will probably be too weak to power the observed X-ray flares. We therefore argue that, for the late-time flares with internal origin, the central engine is unlikely to be associated with BHs. On the contrary, a fast rotating neutron star with strong bipolar magnetic fields may be responsible for such flares.

The emission from black hole binaries (BHBs) and active galactic nuclei (AGNs) displays significant aperiodic variabilities. The most promising explanation for these variabilities is the propagating fluctuations in the accretion flow. It is natural to expect that the mechanism driving variabilities in BHBs and AGNs may operate in a black hole hyper-accretion disk, which is believed to power gamma-ray bursts (GRBs). We study the variabilities of jet power in GRBs based on the model of propagating fluctuations. It is found that the variabilities of jet power and the temporal profile of erratic spikes in this scenario are similar to those in observed light curves of prompt gamma-ray emission of GRBs. Our results show that the mechanism driving X-ray variabilities in BHBs and AGNs may operate in the central engine to drive the variabilities of GRBs.

X-ray flares are generally supposed to be produced by the later central engine activities, and may share the similar physical origin with prompt emission of gamma-ray bursts (GRBs). In this paper, we have analyzed all significant X-ray flares from the GRBs observed by {\em Swift} from April 2005 to March 2015. The catalog contains 468 bright X-ray flares, including 200 flares with redshifts. We obtain the fitting results of X-ray flares, such as start time, peak time, duration, peak flux, fluence, peak luminosity, and mean luminosity. The peak luminosity decreases with peak time, following a power-law behavior $L_p \propto T_{peak,z}^{-1.27}$. The flare duration increases with peak time. The 0.3-10 keV isotropic energy of X-ray flares distribution is a lognormal peaked at $10^{51.2}$ erg. We also study the frequency distributions of flare parameters, including energies, durations, peak fluxes, rise times, decay times and waiting times. Power-law distributions of energies, durations, peak fluxes, and waiting times are found in GRB X-ray flares and solar flares. These distributions could be well explained by a fractal-diffusive, self-organized criticality model. Some theoretical models basing on magnetic reconnection have been proposed to explain X-ray flares. Our result shows that the relativistic jets of GRBs may be Poynting-flux dominated.

The variability of accretion rate, which is believed to induce the aperiodic variability of X-ray emission from disk, may affect the energy injection into the jet. In this spirit, a correlation between disk emission and jet emission can be formed even if the mean luminosity of disk emission keeps constant. In this work, these correlations are found in the situation that the luminosity of disk emission is variable and kept with a constant mean value. The obtained correlations may be shallower than that of the fundamental plane of black hole activity. In addition, the slope of correlation may increase with increasing observed frequency of jet emission. For the luminosities spacing with three days, the slope of correlation decreases with increasing black hole mass. The deviation of our found correlations from that of the fundamental plane is related to the suppression of variability in the jet emission in comparison with that in the disk emission. This mechanism may work in some sources in which shallower correlations have been reported. Moreover, it implies that luminosities used to estimate the relation of fundamental plane should cover an appropriate timescale, in which the variability of jet emission is not significantly suppressed.