Observationally, fast radio bursts (FRBs) can be divided into repeating and apparently non-repeating (one-off) ones. It is unclear whether all FRBs repeat and whether there are genuine non-repeating FRBs. We attempt to address these questions using Monte Carlo simulations. We define a parameter $T_c$ at which the accumulated number of non-repeating sources becomes comparable to the total number of the repeating sources, which is a good proxy to denote the intrinsic repeater fraction among FRBs. Assuming that both types of sources exist, we investigate how the {\em observed} repeater fraction evolves with time for different parameters, including $T_c$, the repeating rate distribution power-law index $q$, the minimum repetition rate $r_{\rm 0,min}$, as well as the Weibull distribution index $k$ for repeating bursts. We find that unless $T_c \rightarrow \infty$ (i.e. there is no genuine non-repeating FRB source), the observed repeater fraction should increase with time first, reaching a peak, and then declines. The peak time $T_p$ and the peak fraction $F_{\rm r,obs,p}$ depend on $T_c$ and other repeating rate parameters. With the current data, one may pose a lower limit $T_c > (0.5-25)$ d for reasonable parameter values. Future continuous monitoring of FRBs using wide-field radio telescopes such as CHIME would measure or set a more stringent lower limit on $T_c$. The detection of a peak in the observed repeater fraction would disfavor the ansatz that ``all FRB sources repeat''.

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

The newly discovered 16.35 days period for repeating FRB 180916.J0158+65 provides an essential clue for understanding the sources and emission mechanism of repeating FRBs. Many models propose that the periodically repeating FRBs might be related to binary star systems that contain at least one neutron star (NSC-FRB system). It has been suggested that the NS "combed" by the strong wind from a companion star might provide a solution. Following the "Binary Comb" model, we use the population synthesis method to study in detail the properties of the companion stars and the nature of NSC-FRB systems. Our main findings are: 1) the companion star is most likely to be a B-type star; 2) the period of 16 days of FRB 180916 happens to fall in the most probable period range, which may explain why FRB 180916 was the first detected periodically repeating FRB, and we expect to observe more periodically repeating FRBs with periods around 10-30 days; 3) the birth rate for NSC-FRB system is large enough to fulfill the event rate requirement set by the observation of FRB 180916, which supports the proposal that NSC-FRB can provide one significant channel for producing periodically repeating FRBs.

Fast radio bursts (FRBs) are a mysterious astrophysical phenomenon of bright pulses emitted at radio frequencies, and it is expected for them to be frequently detected in the future. The dispersion measures of FRBs are related to cosmological parameters, and thus FRBs have a potential to be developed into a new cosmological probe if their data can be largely accumulated in the future. In this work, we study the capability of future FRB data for improving cosmological parameter estimation in two dynamical dark energy models. We find that the simulated FRB data can break the parameter degeneracies inherent in the current cosmic microwave background (CMB) data. Therefore, the combination of the CMB and FRB data can significantly improve the constraints on the Hubble constant and dark energy parameters, compared to those using CMB or FRB alone. If ten thousand events of FRBs with known redshifts are detected in the future, they would behave better than the baryon acoustic oscillation (BAO) data in breaking the parameter degeneracies inherent in the CMB data. We also find that the combination of FRB and gravitational wave (GW) standard siren data provides an independent low-redshift probe to verify the results from CMB and BAO data. For the data combination of CMB, GW, and FRB, it is found that the main contribution to the constraints comes from the CMB and GW data, but the inclusion of FRB data still can evidently improve the constraint on the baryon density.

Fast radio bursts (FRBs) are mysterious millisecond-duration radio transients observed from extragalactic distances. Magnetars have been long speculated as the possible engine to power repeating bursts from FRB sources, but no convincing evidence has been collected so far. Recently, a Galactic magnetar dubbed Soft Gamma-ray Repeater (SGR) J1935+2154 entered an active phase by emitting intense soft gamma-ray bursts. One fast radio burst with two peaks (hereafter FRB 200428) and a luminosity comparable to those of extragalactic FRBs was detected from the source in association with a soft $\gamma$-ray / hard X-ray flare. Here we report an eight-hour targeted observational campaign with the Five-hundred-meter Aperture Spherical radio Telescope (FAST) in four observing sessions, assisted by multi-wavelength (optical and hard X-rays) observations. During the third session, twenty-nine SGR bursts were detected by the {\it Fermi} Gamma-ray Burst Monitor (GBM). Throughout the observing period, especially during the scrutinized epochs when SGR bursts reached the Earth, no single dispersed pulsed emission was detected by FAST. This places a fluence upper limit eight orders of magnitude deeper than the fluence of FRB 200428. Our results suggest that FRB -- SGR burst associations are rather rare. FRBs may be highly relativistic and geometrically beamed or may carry narrow spectra with the characteristic frequencies of the most bursts far outside the FAST band. More intriguingly, it is also possible that the physical condition to achieve coherent radiation in SGR bursts is contrived. Only under extreme conditions could an FRB be made in association with an SGR burst.

Fast Radio Bursts (FRBs) are bright radio transients with millisecond duration at cosmological distances. Since compact dark matter/objects (COs) could act as lenses and cause split of this kind of very short duration signals, Mu$\rm{\tilde{n}}$oz et al. (2016) has proposed a novel method to probe COs with lensing of FRBs. In this Letter, we for the first time apply this method to real data and give constraints of the nature of COs with currently available FRB observations. We emphasize the information from dynamic spectra of FRBs is quite necessary for identifying any lensed signals and find no echoes in the existing data. The null search gives a constraint comparable to that from galactic wide binaries, though the methods of redshift inference from dispersion measure would impact a little. Furthermore, we make an improved forecast basing on the distributions of real data for the ongoing and upcoming telescopes. Finally, we discuss the situation where one or more lensed signals will be detected. In such a case, the parameter space of COs can be pinned down very well since the lens mass can be directly determined through the observed flux ratio and time delay between split images.

Fast radio bursts (FRBs) are short pulses observed in radio band from cosmological distances, some of which emit repeating bursts. The physical origins of these mysterious events have been subject to wide speculations and heated debates. One class of models invoke soft gamma-ray repeaters (SGRs), or magnetars, as the sources of FRBs. Magnetars are rotating neutron stars with extremely strong magnetic field and can sporadically emit bursts from X-ray ($\sim$keV) to soft gamma-ray ($\sim$sub-MeV) with duration from $10^{-2}$~s to $10^2$~s. However, even though some bright radio bursts have been observed from some magnetars, no FRB-like events had been detected to be associated with any magnetar burst, including one giant flare, and no radio burst has been associated with any X-ray event from any magnetar. Therefore, there is still no observational evidence for magnetar-FRB association up to today. Recently, a pair of FRB-like bursts (FRB~200428 hereafter) separated by 30 milliseconds (ms) were detected from the general direction of the Galactic magnetar SGR~J1935+2154. Here we report the detection of a non-thermal X-ray burst in the 1-250 keV energy band with the Insight-HXMT satellite, which we identify as emitted from SGR~J1935+2154. The burst showed two hard peaks with a separation of $\sim30$ ms, consistent with the separation between the two bursts in FRB~200428. The delay time between the double radio and X-ray peaks is $\sim8.57$ s, fully consistent with the dispersion delay of FRB~200428. We thus identify the non-thermal X-ray burst is associated with FRB~200428 whose high energy counterpart is the two hard peaks in X-ray. Our results suggest that the non-thermal X-ray burst and FRB~200428 share the same physical origin in an explosive event from SGR~J1935+2154.

In this paper, we present a numerical method to study the predicted lightcurves as a function of viewing angle. We extrapolate the fitting formulae for the mass and velocity of tidal dynamical ejecta across a wide range of mass ratio validated with 66 simulations and use them in the calculations of the kilonova lightcurves. The calculated peak luminosity of a BH-NS merger kilonova is typically about a few times $10^{41}\ {\rm erg\ s^{-1}}$, which is always $\lesssim4.5\times10^{41}\ {\rm erg\ s^{-1}}$. This corresponds to the AB absolute magnitudes fainter than $\sim -15\ {\rm mag}$ in optical and $\sim -16\ {\rm mag}$ in infrared. Since the projected photosphere area of the dynamical ejecta is much larger than that of the disk wind outflows, the dynamical ejecta usually contribute to the majority of the kilonova emission from BH-NS mergers. The fitted blackbody temperature and the shape of the observed multi-band lightcurves are insensitive to the line of sight. The peak time of the observed multi-band lightcurves, affected by the light propagation effect, is related to the relative motion direction between the dynamical ejecta and the observer. The observed luminosity varies with the projected photosphere area determined by the viewing angles. However, the predicted peak luminosity only varies by a factor of $\sim (2 - 3)$ (or by $\sim1\ {\rm mag}$) for different viewing angles. When the short-duration gamma-ray burst afterglow is taken into account, for an on-axis geometry, the kilonova emission is usually outshone by the afterglow emission and can be only observed in the redder bands, especially in the $K$-band at late times. Compared with GW170817/AT2017gfo, the BH-NS merger kilonovae are optically dim but possibly infrared bright. At the same epoch after the merger, the blackbody fitting temperature of the BH-NS merger kilonovae is lower than that of GW170817/AT2017gfo.

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.

We use the In data collected during the 2019 outburst from X-ray pulsar 4U 1901+03 to complement the orbital parameters reported by Fermi/GBM. Using the Insight-HXMT, we examine the correlation between the derivative of the intrinsic spin frequency and bolometric flux based on accretion torque models. It was found that the pulse profiles significantly evolve during the outburst. The existence of two types of the profile's pattern discovered in the Insight-HXMT data indicates that this source experienced transition between a super-critical and a sub-critical accretion regime during its 2019 outburst. Based on the evolution of the pulse profiles and the torque model, we derive the distance to 4U 1901+03 as 12.4+-0.2 kpc.

We study the evolution of the temporal properties of MAXI 1820+070 during the 2018 outburst in its hard state from MJD 58190 to 58289 with Insight-HXMT in a broad energy band 1-150 keV. We find different behaviors of the hardness ratio, the fractional rms and time lag before and after MJD 58257, suggesting a transition occurred at around this point. The observed time lags between the soft photons in the 1-5 keV band and the hard photons in higher energy bands, up to 150 keV, are frequency-dependent: the time lags in the low-frequency range, 2-10 mHz, are both soft and hard lags with a timescale of dozens of seconds but without a clear trend along the outburst; the time lags in the high-frequency range, 1-10 Hz, are only hard lags with a timescale of tens of milliseconds; first increase until around MJD 58257 and decrease after this date. The high-frequency time lags are significantly correlated to the photon index derived from the fit to the quasi-simultaneous NICER spectrum in the 1-10 keV band. This result is qualitatively consistent with a model in which the high-frequency time lags are produced by Comptonization in a jet.

Five fast radio bursts (FRBs), including three apparently non-repeating ones FRB 180924, FRB 181112, and FRB 190523, and two repeaters, FRB 121102 and FRB 180916.J0158+65, have already been localized so far. We apply a method developed recently by us (Li et al. 2019) to these five localized FRBs to give a cosmology-insensitive estimate of the fraction of baryon mass in the IGM, $f_{\rm IGM}$. Using the measured dispersion measure (DM) and luminosity distance $d_{\rm L}$ data (inferred from the FRB redshifts and $d_{\rm L}$ of type Ia supernovae at the same redshifts) of the five FRBs, we constrain the local $f_{\rm IGM} = 0.84^{+0.16}_{-0.22}$ with no evidence of redshift dependence. This cosmology-insensitive estimate of $f_{\rm IGM}$ from FRB observations is in excellent agreement with previous constraints using other probes. Moreover, using the three apparently non-repeating FRBs only we get a little looser but consistent result $f_{\rm IGM} = 0.74^{+0.24}_{-0.18}$. In these two cases, reasonable estimations for the host galaxy DM contribution (${\rm DM_{host}}$) can be achieved by modelling it as a function of star formation rate. The constraints on both $f_{\rm IGM}$ and ${\rm DM_{host}}$ are expected to be significantly improved with the rapid progress in localizing FRBs.

Glitches correspond to sudden jumps of rotation frequency ($\nu$) and its derivative ($\dot{\nu}$) of pulsars, the origin of which remains not well understood yet, partly because the jump processes of most glitches are not well time-resolved. There are three large glitches of the Crab pulsar, detected in 1989, 1996 and 2017, which were found to have delayed spin-up processes before the normal recovery processes. Here we report two additional glitches of the Crab pulsar occurred in 2004 and 2011 for which we discovered delayed spin up processes, and present refined parameters of the largest glitch occurred in 2017. The initial rising time of the glitch is determined as $<0.48$ hour. We also carried out a statistical study of these five glitches with observed spin-up processes. The two glitches occurred in 2004 and 2011 have delayed spin-up time scales ($\tau_{1}$) of $1.7\pm0.8$\,days and $1.6\pm0.4$\,days, respectively. We find that the $\Delta{\nu}$ vs. $|\Delta{\dot\nu}|$ relation of these five glitches is similar to those with no detected delayed spin-up process, indicating that they are similar to the others in nature except that they have larger amplitudes. For these five glitches, the amplitudes of the delayed spin-up process ($|\Delta{\nu}_{\rm d1}|$) and recovery process ($\Delta{\nu}_{\rm d2}$), their time scales ($\tau_{1}$, $\tau_{2}$), and permanent changes in spin frequency ($\Delta{\nu}_{\rm p}$) and total frequency step ($\Delta{\nu}_{\rm g}$) have positive correlations. From these correlations, we suggest that the delayed spin-up processes are common for all glitches, but are too short and thus difficult to be detected for most glitches.

No robust detection of prompt electromagnetic counterparts to fast radio bursts (FRBs) has yet been obtained, in spite of several multi-wavelength searches carried out so far. Specifically, X/gamma-ray counterparts are predicted by some models. We planned on searching for prompt gamma-ray counterparts in the Insight-Hard X-ray Modulation Telescope (Insight-HXMT) data, taking advantage of the unique combination of large effective area in the keV-MeV energy range and of sub-ms time resolution. We selected 39 FRBs that were promptly visible from the High-Energy (HE) instrument aboard Insight-HXMT. After calculating the expected arrival times at the location of the spacecraft, we searched for a significant excess in both individual and cumulative time profiles over a wide range of time resolutions, from several seconds down to sub-ms scales. Using the dispersion measures in excess of the Galactic terms, we estimated the upper limits on the redshifts. No convincing signal was found and for each FRB we constrained the gamma-ray isotropic-equivalent luminosity and the released energy as a function of emission timescale. For the nearest FRB source, the periodic repeater FRB180916.J0158+65, we find $L_{\gamma,iso}<5.5\times 10^{47}$ erg/s over 1 s, whereas $L_{\gamma,iso}<10^{49}-10^{51}$ erg/s for the bulk of FRBs. The same values scale up by a factor of ~100 for a ms-long emission. Even on a timescale comparable with that of the radio pulse itself no keV-MeV emission is observed. A systematic association with either long or short GRBs is ruled out with high confidence, except for subluminous events, as is the case for core-collapse of massive stars (long) or binary neutron star mergers (short) viewed off axis. Only giant flares from extra-galactic magnetars at least ten times more energetic than Galactic siblings are ruled out for the nearest FRB.

The post-merger product of the first binary neutron star merger event detected in gravitational waves, GW170817, depends on neutron star equation of state (EoS) and is not well determined. We generally discuss the constraints one may pose on the maximum mass of a non-spinning neutron star, $M_{\rm TOV}$, based on the observations and some EoS-independent universal relations of rapidly-spinning neutron stars. If the merger product is a black hole after a brief hypermassive neutron star (HMNS) phase, we derive $M_{\rm TOV} < 2.09^{+0.11}_{-0.09}(^{+0.06}_{-0.04}) M_{\odot}$ at the 2$\sigma$ (1$\sigma$) level. The cases for a massive neutron star (MNS), either a supra-massive neutron star (SMNS) or even a stable neutron star (SNS), are also allowed by the data. We derive $2.09^{+0.11}_{-0.09}(^{+0.06}_{-0.04} M_{\odot}) \leq M_{\rm TOV}< 2.43^{+0.10}_{-0.08}(^{+0.06}_{-0.04}) M_{\odot}$ for the SMNS case and $M_{\rm TOV} \geq 2.43^{+0.10}_{-0.08}(^{+0.06}_{-0.04})M_{\odot}$ for the SNS case, at the $2\sigma$ ($1\sigma$) confidence level. In the MNS cases, we also discuss the constraints on the neutron star parameters (the dipolar magnetic field strength at the surface $B_p$ and the ellipticity $\epsilon$) that affect the spindown history, by considering different MNS survival times, e.g. 300 s, 1 d, and 155 d after the merger, as suggested by various observational arguments. We find that once an SMNS is formed, without violating the EM observational constraints, there always exist a set of ($B_p, \epsilon$) parameters that allow the SMNS to survive for 300s, 1 d, 155 d, or even longer.

Bars in disc-dominated galaxies are able to drive gas inflow inside the corotation radius, thus enhancing the central star formation rate (SFR). Previous work, however, has found that disc-dominated galaxies with centrally suppressed SFR frequently host a bar. Here we investigate possible causes for the suppression of central SFR in such cases. We compare physical properties of a sample of disc-dominated barred galaxies with high central SFR (HC galaxies) with those of a sample of disc-dominated barred galaxies with low central SFR (LC galaxies). We find that the two samples have on average similar HI content and bars of similar strength. But we also find that the HCs have bluer colors than LCs, and that outside the bar region they host stronger spiral arms than the LCs where closed rings are more often seen. We discuss and evaluate the possible causes for the suppression of the central SFR in the LC galaxies as opposed to its enhancement in the HC galaxies.

We study the statistical properties of 320 bulges of disk galaxies in the Carnegie-Irvine Galaxy Survey, using robust structural parameters of galaxies derived from image fitting. We apply the Kormendy relation to classify classical and pseudo bulges and characterize bulge dichotomy with respect to bulge structural properties and physical properties of host galaxies. We confirm previous findings that pseudo bulges on average have smaller S\'{e}rsic indices, smaller bulge-to-total ratios, and fainter surface brightnesses when compared with classical bulges. Our sizable sample statistically shows that pseudo bulges are more intrinsically flattened than classical bulges. Pseudo bulges are most frequent (incidence $\gtrsim 80\%$) in late-type spirals (later than Sc). Our measurements support the picture in which pseudo bulges arose from star formation induced by inflowing gas, while classical bulges were born out of violent processes such as mergers and coalescence of clumps. We reveal differences with the literature that warrant attention: (1) the bimodal distribution of S\'{e}rsic indices presented by previous studies is not reproduced in our study; (2) classical and pseudo bulges have similar relative bulge sizes; and (3) the pseudo bulge fraction is considerably smaller in early-type disks compared with previous studies based on one-dimensional surface brightness profile fitting. We attribute the above differences to our improved image quality, more robust bulge-to-disk decomposition technique, and different classification criteria applied. Moreover, we find that barred galaxies do not host more pseudo bulges or more prominent pseudo bulges than unbarred galaxies. Various implications of these findings are discussed.

In the literature, compact binary coalescences (CBCs) have been proposed as one of the main scenarios to explain the origin of some non-repeating fast radio bursts (FRBs). The large discrepancy between the FRB and CBC event rate densities suggest their associations, if any, should only apply at most for a small fraction of FRBs. Through a Bayesian estimation method, we show how a statistical analysis of the coincident associations of FRBs with CBC gravitational wave (GW) events may test the hypothesis of these associations. We show that during the operation period of advanced LIGO, the detection of $\sim100$ ($\sim1000$) GW-less FRBs with dispersion measure (DM) values smaller than 500 ${\rm pc~cm^{-3}}$ could reach the constraint that less than $10\%$ (or $1\%$) FRBs are related to binary black hole (BBH) mergers. The same number of FRBs with DM values smaller than 100 ${\rm pc~cm^{-3}}$ is required to reach the same constraint for binary neutron star (BNS) mergers. With the upgrade of GW detectors, the same constraints for BBH and BNS mergers can be reached with less FRBs or looser requirements for the DM values. It is also possible to pose constraints on the fraction of each type of CBCs that are able to produce observable FRBs based on the event density of FRBs and CBCs. This would further constrain the dimensionless charge of black holes in binary BH systems.

We report on our analysis of the 2019 outburst of the X-ray accreting pulsar 4U 1901+03 observed with Insight-HXMT and NICER. Both spectra and pulse profiles evolve significantly in the decaying phase of the outburst. Dozens of flares are observed throughout the outburst. They are more frequent and brighter at the outburst peak. We find that the flares, which have a duration from tens to hundreds of seconds, are generally brighter than the persistent emission by a factor of $\sim$ 1.5. The pulse profile shape during the flares can be significantly different than that of the persistent emission. In particular, a phase shift is clearly observed in many cases. We interpret these findings as direct evidence of changes of the pulsed beam pattern, due to transitions between the sub- and super-critical accretion regimes on a short time scale. We also observe that at comparable luminosities the flares' pulse profiles are rather similar to those of the persistent emission. This indicates that the accretion on the polar cap of the neutron star is mainly determined by the luminosity, i.e., the mass accretion rate.

A new black hole X-ray binary (BHXRB) MAXI J1535-571 was discovered by MAXI during its outburst in 2017. Using observations taken by the first Chinese X-ray satellite, the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT), we perform a joint spectral analysis (2-150 keV) in both energy and time domains. The energy spectra provide the essential input for probing the intrinsic Quasi-Periodic Oscillation (QPO) fractional rms spectra (FRS). Our results show that during the intermediate state, the energy spectra are in general consistent with those reported by Swift/XRT and NuSTAR. However, the QPO FRS become harder and the FRS residuals may suggest the presence of either an additional power-law component in the energy spectrum or a turn-over in the intrinsic QPO FRS at high energies.

The joint detection of GW 170817 and GRB 170817A indicated that at least a fraction of short gamma ray bursts (SGRBs) originate from binary neutron star (BNS) mergers. One possible remnant of a BNS merger is a rapidly rotating, strongly magnetized neutron star, which has been discussed as one possible central engine for GRBs. For a rapidly rotating magnetar central engine, the deposition of the rotation energy into the ejecta launched from the merger could lead to bright radio emission. The brightness of radio emission years after an SGRB would provide an estimation of the kinetic energy of ejecta and hence, a possible constraint on the BNS merger product. We perform a more detailed calculation on the brightness of radio emission from the interaction between the merger ejecta and circumburst medium in the magnetar scenario, invoking several important physical processes such as generic hydrodynamics, relativistic effects, and the deep Newtonian phase. We use the model to constrain the allowed parameter space for 15 SGRBs that have late radio observations. Our results show that an injected energy $E_{\rm inj} \sim 10^{52}$ erg is allowed for all the cases, which suggests that the possibility of a supra-massive or hyper-massive neutron star remnant is not disfavored by the available radio data.

We study the 2018 outburst of Aql X-1 via the monitor of all sky X-ray image (MAXI) data. We show that the outburst starting in February 2018 is a member of short-low class in the frame of outburst duration and the peak count rate although the outburst morphology is slightly different from the other fast-rise-exponential-decay (FRED) type outbursts with a milder rising stage. We study the partial accretion in the weak propeller stage of Aql X-1 via the MAXI data of the 2018 outburst. We report on the spectral analysis of 3 observations of Aquila X-1 obtained by Insight - hard X-ray modulation telescope (Insight-HXMT) during the late decay stage of the 2018 outburst. We discuss that the data taken by Insight-HXMT is just after the transition to the weak propeller stage. Our analysis shows the necessity of a comptonization component to take into account the existence of an electron cloud resulting photons partly up-scattered.

As China's first X-ray astronomical satellite, the Hard X-ray Modulation Telescope (HXMT), which was dubbed as Insight-HXMT after the launch on June 15, 2017, is a wide-band (1-250 keV) slat-collimator-based X-ray astronomy satellite with the capability of all-sky monitoring in 0.2-3 MeV. It was designed to perform pointing, scanning and gamma-ray burst (GRB) observations and, based on the Direct Demodulation Method (DDM), the image of the scanned sky region can be reconstructed. Here we give an overview of the mission and its progresses, including payload, core sciences, ground calibration/facility, ground segment, data archive, software, in-orbit performance, calibration, background model, observations and some preliminary results.

We present a detailed timing study of the brightest persistent X-ray source Sco X-1 using the data collected by the Hard X-ray Modulation Telescope ($Insight$-HXMT) from July 2017 to August 2018. A complete $Z$-track hardness-intensity diagram (HID) is obtained. The normal branch oscillations (NBOs) at $\sim$ 6 Hz in the lower part of the normal branch (NB) and the flare branch oscillations (FBOs) at $\sim$ 16 Hz in the beginning part of the flaring branch (FB) are found in observations with the Low Energy X-ray Telescope (LE) and the Medium Energy X-ray Telescope (ME) of $Insight$-HXMT, while the horizontal branch oscillations (HBOs) at $\sim$ 40 Hz and the kilohertz quasi-periodic oscillations (kHz QPOs) at $\sim$ 800 Hz are found simultaneously up to 60 keV for the first time on the horizontal branch (HB) by the High Energy X-ray Telescope (HE) and ME. We find that for all types of the observed QPOs, the centroid frequencies are independent of energy, while the root mean square (rms) increases with energy; the centroid frequencies of both the HBOs and kHz QPOs increase along the $Z$-track from the top to the bottom of the HB; and the NBOs show soft phase lags increasing with energy. A continuous QPO transition from the FB to NB in $\sim$ 200 s are also detected. Our results indicate that the non-thermal emission is the origin of all types of QPOs, the innermost region of the accretion disk is non-thermal in nature, and the corona is nonhomogeneous geometrically.

The evidences for the influence of thermonuclear (type-I) X-ray bursts upon the surrounding environments in neutron star low-mass X-ray binaries (LMXB) were detected previously via spectral and timing analyses. Benefitting from a broad energy coverage of Insight-HXMT, we analyze one photospheric radius expansion (PRE) burst, and find an emission excess at soft X-rays. Our spectral analysis shows that, such an excess is not likely relevant to the disk reflection induced by the burst emission and can be attributed to an enhanced pre-burst/persistent emission. We find that the burst and enhanced persistent emissions sum up to exceed Eddington luminosity by $\sim$ 40 percentages. We speculate that the enhanced emission is from a region beyond the PRE radius, or through the Comptonization of the corona.

Corona cooling was detected previously from stacking a series of short type-I bursts occurred during the low/had state of atoll outburst. Type-I bursts are hence regarded as sharp probe to our better understanding on the basic property of the corona. The launch of the first Chinese X-ray satellite Insight-HXMT has large detection area at hard X-rays which provide almost unique chance to move further in this research field. We report the first detection of the corona cooling by Insight-HXMT from single short type-I burst showing up during {\bf flare} of 4U 1636-536. This type-I X-ray burst has a duration of $\sim$13 seconds and hard X-ray shortage is detected with significance 6.2~$\sigma$ in 40-70 keV. A cross-correlation analysis between the lightcurves of soft and hard X-ray band, shows that the corona shortage lag the burst emission by 1.6 $\pm$1.2~s. These results are consistent with those derived previously from stacking a large amount of bursts detected by RXTE/PCA within a series of {\bf flares} of 4U 1636-536. Moreover, the broad bandwidth of Insight-HXMT allows as well for the first time to infer the burst influence upon the continuum spectrum via performing the spectral fitting of the burst, which ends up with the finding that hard X-ray shortage appears at around 40 keV in the continuum spectrum. These results suggest that the evolution of the corona along with the outburst{\bf /flare} of NS XRB may be traced via looking into a series of embedded type-I bursts by using Insight-HXMT.

The Hard X-ray Modulation Telescope (HXMT) named Insight is China's first X-ray astronomical satellite. The Low Energy X-ray Telescope (LE) is one of its main payloads onboard. The detectors of LE adopt swept charge device CCD236 with L-shaped transfer electrodes. Charges in detection area are read out continuously along specific paths, which leads to a time response distribution of photons readout time. We designed a long exposure readout mode to measure the time response distribution. In this mode, CCD236 firstly performs exposure without readout, then all charges generated in preceding exposure phase are read out completely. Through analysis of the photons readout time in this mode, we obtained the probability distribution of photons readout time.

We present a study of timing properties of the accreting pulsar 2S 1417-624 observed during its 2018 outburst, based on Swift/BAT, Fermi/GBM, Insight-HXMT and NICER observations. We report a dramatic change of the pulse profiles with luminosity. The morphology of the profile in the range 0.2-10.0keV switches from double to triple peaks at $\sim2.5$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$ and from triple to quadruple peaks at $\sim7$ $\rm \times 10^{37}{\it D}_{10}^2\ erg\ s^{-1}$. The profile at high energies (25-100keV) shows significant evolutions as well. We explain this phenomenon according to existing theoretical models. We argue that the first change is related to the transition from the sub to the super-critical accretion regime, while the second to the transition of the accretion disc from the gas-dominated to the radiation pressure-dominated state. Considering the spin-up as well due to the accretion torque, this interpretation allows to estimate the magnetic field self-consistently at $\sim7\times 10^{12}$G.

The long-term evolution of the centroid energy of the CRSF in Her X-1 is still a mystery. We report a new measurement from a campaign between {\sl Insight}-HXMT and {\sl NuSTAR} performed in February 2018. Generally, the two satellites show well consistent results of timing and spectral properties. The joint spectral analysis confirms that the previously observed long decay phase has ended, and that the line energy instead keeps constant around 37.5 keV after flux correction.

We present multi-component decomposition of high-quality $R$-band images of 312 disk galaxies from the Carnegie-Irvine Galaxy Survey. In addition to bulges and disks, we successfully model nuclei, bars, disk breaks, nuclear/inner lenses, and inner rings. Our modeling strategy treats nuclear rings and nuclear bars as part of the bulge component, while other features such as spiral arms, outer lenses, and outer rings are omitted from the fits because they are not crucial for accurate bulge measurements. The error budget of bulge parameters includes the uncertainties from sky level measurements and model assumptions. Comparison with multi-component decomposition from the \textit{Spitzer} Survey of Stellar Structure in Galaxies reveals broad agreement for the majority of the overlapping galaxies, but for a considerable fraction of galaxies there are significant differences in bulge parameters caused by different strategies in model construction. We confirm that on average bulge prominence decreases from early to late-type disk galaxies, although the large scatter of bulge-to-total ratios in each morphological bin limits the application of Hubble type as an accurate predictor of bulge-to-total ratio. In contrast with previous studies claiming that barred galaxies host weaker bulges, we find that barred and unbarred spiral galaxies have similar bulge prominence.

We report on analysis of observations of the bright transient X-ray pulsar \src obtained during its 2017-2018 giant outburst with Insight-HXMT, \emph{NuSTAR}, and \textit{Swift} observatories. We focus on the discovery of a sharp state transition of the timing and spectral properties of the source at super-Eddington accretion rates, which we associate with the transition of the accretion disk to a radiation pressure dominated (RPD) state, the first ever directly observed for magnetized neutron star. This transition occurs at slightly higher luminosity compared to already reported transition of the source from sub- to super-critical accretion regime associate with onset of an accretion column. We argue that this scenario can only be realized for comparatively weakly magnetized neutron star, not dissimilar to other ultra-luminous X-ray pulsars (ULPs), which accrete at similar rates. Further evidence for this conclusion is provided by the non-detection of the transition to the propeller state in quiescence which strongly implies compact magnetosphere and thus rules out magnetar-like fields.

Kilonova explosions typically release $\sim 10^{50-51}$ erg in kinetic energy, which is sufficient to constitute a kilonova remnant (KNR); however, it has not yet been confirmed. In this work, we investigate the probable association between G4.8+6.2 and the guest star of AD 1163, which is recorded by the Korea ancient astronomers. Although the evidence available is insufficient to draw a definite conclusion, it is at least theoretically self-consistent that the guest star of AD 1163 was a historical kilonova associated with G4.8+6.2, considering the possible short visible timescale of AD 1163, the relatively high Galactic latitude of G4.8+6.2, and that G4.8+6.2 is spatially coincident with the guest star of AD 1163. Further observation of G4.8+6.2 is needed to test our hypothesis. If our interpretation is correct, our results indicate that young KNRs should have a large diameter and low surface brightness, unlike other young supernova remnants.

Extragalactic dispersion measures (DMs) obtained from observations of fast radio bursts (FRBs) are an excellent tool for probing intergalactic medium (IGM) and for conducting cosmography. However, the DM contribution from the IGM (${\rm DM_{IGM}}$) depends on the fraction of baryon mass in the IGM, $f_{\rm IGM}$, which is not properly constrained. As $f_{\rm IGM}(z)$ is geometrically related to the Hubble parameter $H(z)$ and ${\rm DM_{IGM}}(z)$, here we propose that combining two independent measurements of FRBs and $H(z)$ in similar redshift ranges provides a novel and cosmology-free method to constrain the evolution of $f_{\rm IGM}(z)$. Under the assumption that $f_{\rm IGM}$ is evolving with redshift in a functional form, we forecast that the evolution of $f_{\rm IGM}(z)$ can be well inferred in a combined analysis of $\sim3000$ ${\rm DM_{IGM}}(z)$ derived from FRBs and $\sim50$ $H(z)$ derived from the Hubble parameter data. Though the efficiency of our method is not as good as that of the other model-independent method involving the joint measurements of DM and luminosity distance of FRBs, our method offers a new model-independent way to constrain $f_{\rm IGM}(z)$.

The LIGO and Virgo scientific collaboration (LVC) alerted two binary neutron star (BNS) merger candidates, S190425z and S190426c. Fermi-GBM observed 55.6\% (for S190425z) and 100\% (for S190426c) of the probability regions of both events at the respective merger times, but no gamma-ray burst (GRB) was detected in either case. The derived luminosity upper limits suggest that a short GRB similar to GRB 170817A would not be detectable for both cases due to their larger distances than GW170817. Assuming that the jet profile obtained from the GW170817/GRB 170817A is quasi-universal for all BNS-GRB associations, we derive that the viewing angles of S190425z and S190426c should be $> (0.11-0.41)$ and $> (0.09-0.39)$, respectively. Through Monte Carlo simulations, we show that with the GRB 170817A-like jet structure, all sky gamma-ray detectors, such as GBM and GECAM, are expected to detect $\sim 4.6\%$, $3.9\%$, $1.7\%$ and $6.6\%$, $5.7\%$, $2.8\%$ BNS mergers triggered by aLIGO, aLIGO A+ and ET, respectively. The joint detection fraction would be largely reduced for Swift-BAT, SVOM-ECLAIRS and Einstein Probe, whose sensitivities are better but fields of view are smaller.

In this work, we report the in-orbit demonstration of X-ray pulsar navigation with Insight-Hard X-ray Modulation Telescope (Insight-HXMT), which was launched on Jun. 15th, 2017. The new pulsar navigation method 'Significance Enhancement of Pulse-profile with Orbit-dynamics' (SEPO) is adopted to determine the orbit with observations of only one pulsar. In this test, the Crab pulsar is chosen and observed by Insight-HXMT from Aug. 31th to Sept. 5th in 2017. Using the 5-day-long observation data, the orbit of Insight-HXMT is determined successfully with the three telescopes onboard - High Energy X-ray Telescope (HE), Medium Energy X-ray Telescope (ME) and Low Energy X-ray Telescope (LE) - respectively. Combining all the data, the position and velocity of the Insight-HXMT are pinpointed to within 10 km (3 sigma) and 10 m/s (3 sigma), respectively.

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

Fast radio bursts (FRBs) are highly dispersed and probably extragalactic radio flashes with millisecond-duration. Recently, the Canadian Hydrogen Intensity Mapping Experiment (using the CHIME/FRB instrument) has reported detections of 13 FRBs during a pre-commissioning phase. It is more exciting that one of the 13 FRBs is a second source of repeaters which suggests that CHIME/FRB and other wide-field sensitive radio telescopes will find a substantial population of repeating FRBs. We have proposed strongly lensed repeating FRBs as a precision cosmological probe, e.g. constraining the Hubble constant and model-independently estimating the cosmic curvature. Here, we study complementary constraints on the equation of state of dark energy from strongly lensed FRBs to currently available popular probes. It is found that, in the framework of Chevalier-Polarski-Linder parametrization, adding time delay measurement of 30 strongly lensed FRB systems to cosmic microwave background radiation and type Ia supernovae can improve the dark energy figure of merit by a factor 2. In the precision cosmology era, this improvement is of great significance for studying the nature of dark energy.

The Gamma-Ray Integrated Detectors (GRID) is a space mission concept dedicated to monitoring the transient gamma-ray sky in the energy range from 10 keV to 2 MeV using scintillation detectors onboard CubeSats in low Earth orbits. The primary targets of GRID are the gamma-ray bursts (GRBs) in the local universe. The scientific goal of GRID is, in synergy with ground-based gravitational wave (GW) detectors such as LIGO and VIRGO, to accumulate a sample of GRBs associated with the merger of two compact stars and study jets and related physics of those objects. It also involves observing and studying other gamma-ray transients such as long GRBs, soft gamma-ray repeaters, terrestrial gamma-ray flashes, and solar flares. With multiple CubeSats in various orbits, GRID is unaffected by the Earth occultation and serves as a full-time and all-sky monitor. Assuming a horizon of 200 Mpc for ground-based GW detectors, we expect to see a few associated GW-GRB events per year. With about 10 CubeSats in operation, GRID is capable of localizing a faint GRB like 170817A with a 90% error radius of about 10 degrees, through triangulation and flux modulation. GRID is proposed and developed by students, with considerable contribution from undergraduate students, and will remain operated as a student project in the future. The current GRID collaboration involves more than 20 institutes and keeps growing. On August 29th, the first GRID detector onboard a CubeSat was launched into a Sun-synchronous orbit and is currently under test.

The recently discovered neutron star transient Swift J0243.6+6124 has been monitored by {\it the Hard X-ray Modulation Telescope} ({\it Insight-\rm HXMT). Based on the obtained data, we investigate the broadband spectrum of the source throughout the outburst. We estimate the broadband flux of the source and search for possible cyclotron line in the broadband spectrum. No evidence of line-like features is, however, found up to $\rm 150~keV$. In the absence of any cyclotron line in its energy spectrum, we estimate the magnetic field of the source based on the observed spin evolution of the neutron star by applying two accretion torque models. In both cases, we get consistent results with $B\rm \sim 10^{13}~G$, $D\rm \sim 6~kpc$ and peak luminosity of $\rm >10^{39}~erg~s^{-1}$ which makes the source the first Galactic ultraluminous X-ray source hosting a neutron star.

For non-spinning, charged (Reissner-Nordstr\"om) black holes, the particles with an opposite sign of charge with respect to that of the black hole will be pulled into the black hole by the extra electromagnetic force. Such a hole will be quickly neutralized so that there should not exist significantly charged, non-spinning black holes in the universe. The case of spinning, charged (Kerr-Newmann) black holes is more complicated. For a given initial position and initial velocity of the particle, an oppositely charged particle does not always more easily fall into the black hole than a neutral particle. The possible existence of a magnetosphere further complicate the picture. One therefore cannot straightforwardly conclude that a charged spinning black hole will be neutralized. In this paper, we make the first step to investigate the neutralization of Kerr-Newmann black holes without introducing a magnetosphere. We track the particle trajectories under the influence of the curved spacetime and the electromagnetic field carried by the spinning, charged black hole. A statistical method is used to investigate the neutralization problem. We find a universal dependence of the falling probability into the black hole on the charge of the test particle, with the oppositely charged particles having a higher probability of falling. We therefore conclude that charged, spinning black holes without a magnetosphere should be quickly neutralized, consistent with people's intuition. The neutralization problem of Kerr-Newmann black holes with a co-rotating force-free magnetosphere is subject to further studies.

Cosmic backreaction as an additional source of the expansion of the universe has been a debate topic since the discovery of cosmic acceleration. The major concern is whether the self interaction of small-scale nonlinear structures would source gravity on very large scales. Gregory Ryskin argued against the additional inclusion of gravitational interaction energy of astronomical objects, whose masses are mostly inferred from gravitational effects and hence should already contain all sources with long-range gravity forces. Ryskin proposed that the backreaction contribution to the energy momentum tensor is instead from the rest of the universe beyond the observable patch. Ryskin's model elegantly solves the fine-tuning problem and is in good agreement with the Hubble diagram of Type Ia supernovae. In this article we revisit Ryskin's model and show that it is {\it inconsistent} with at least one of the following statements: (i) the universe is matter-dominated at low redshift ($z\lesssim 2$); (ii) the universe is radiation-dominated at sufficiently high redshift; (iii) matter density fluctuations are tiny ($\lesssim 10^{-4}$) at the recombination epoch.

Fast radio bursts (FRBs) are short duration ($\sim$millisecond) radio transients with cosmological origin. The simple sharp features of the FRB signal have been utilized to probe two fundamental laws of physics, namey, testing Einstein's weak equivalence principle and constraining the rest mass of the photon. Recently, \cite{hessels18} found that after correcting for dispersive delay, some of the bursts in FRB 121102 have complex time-frequency structures that include sub-pulses with a time-frequency downward drifting property. Using the delay time between sub-pulses in FRB 121102, here we show that the parameterized post-Newtonian parameter $\gamma$ is the same for photons with different energies to the level of $\left|\gamma_1-\gamma_2\right|<2.5\times 10^{-16}$, which is 1000 times better than previous constraints from FRBs using similar methods. We also obtain a stringent constraint on the photon mass, $m_{\gamma} < 5.1\times10^{-48}$ g, which is 10 times smaller than previous best limits on the photon mass derived through the velocity dispersion method.

With a selected sample of neutron star (NS) equation-of-states (EOSs) that are consistent with the current observations and have a range of maximum masses, we investigate the relations between NS gravitational mass $M_g$ and baryonic mass $M_b$, and the relations between the maximum NS mass supported through uniform rotation ($M_{\rm max}$) and that of nonrotating NSs ($M_{\rm TOV}$). We find that if one intends to apply an EOS-independent quadratic, universal transformation formula ($M_b=M_g+A\times M_{g}^2$) to all EOSs, the best fit $A$ value is 0.080 for non-rotating NSs only and 0.073 when different spin periods are considered. The residual error of the transformation is as large as $\sim0.1M_{\odot}$. For different EOSs, we find that the parameter $A$ for non-rotating NSs is proportional to $R_{1.4}^{-1}$ (where $R_{1.4}$ is NS radius for 1.4$M_\odot$ in unit of km). For a particular EOS, if one adopts the best-fit parameters for different spin periods, the residual error of the transformation is smaller, which is of the order of 0.01$M_\odot$ for the quadratic form and less than 0.01$M_\odot$ for the cubic form ($M_b=M_g+A_1\times M_{g}^2+A_2\times M_{g}^3$). We also find a very tight and general correlation between the normalized mass gain due to spin $\Delta m\equiv(M_{\rm max}-M_{\rm TOV})/M_{\rm TOV}$ and the spin period normalized to the Keplerian period ${\cal P}$, i.e. ${\rm log_{10}}\Delta m = (-2.74\pm0.05){\rm log_{10}}{\cal P}+{\rm log_{10}}(0.20\pm 0.01)$, which is independent of EOS models. Applications of our results to GW170817 is discussed.

We formulate and solve the problem of spherically symmetric, steady state, adiabatic accretion onto a Schwarzschild-like black hole obtained recently. We derive the general analytic expressions for the critical points, the critical velocity, the critical speed of sound, and subsequently the mass accretion rate. The case for polytropic gas is discussed in detail. We find the parameter characterizing the breaking of Lorentz symmetry will slow down the mass accretion rate, while has no effect on the gas compression and the temperature profile below the critical radius and at the event horizon.

The excessive dispersion measure (DM) of fast radio bursts (FRBs) has been proposed to be a powerful tool to study intergalactic medium (IGM) and to perform cosmography. One issue is that the fraction of baryons in the IGM, $f_{\rm IGM}$, is not properly constrained. Here we propose a method of estimating $f_{\rm IGM}$ using a putative sample of FRBs with the measurements of both DM and luminosity distance $d_{\rm L}$. The latter can be obtained if the FRB is associated with a distance indicator (e.g. a gamma-ray burst or a gravitational wave event), or the redshift $z$ of the FRB is measured and $d_{\rm L}$ at the corresponding $z$ is available from other distance indicators (e.g. type Ia supernovae) at the same redshift. Since $d_{\rm L}/{\rm DM}$ essentially does not depend on cosmological parameters, our method can determine $f_{\rm IGM}$ independent of cosmological parameters. We parameterize $f_{\rm IGM}$ as a function of redshift and model the DM contribution from a host galaxy as a function of star formation rate. Assuming $f_{\rm IGM}$ has a mild evolution with redshift with a functional form and by means of Monte Carlo simulations, we show that an unbiased and cosmology-independent estimate of the present value of $f_{\rm IGM}$ with a $\sim 12\%$ uncertainty can be obtained with 50 joint measurements of $d_{\rm L}$ and DM. In addition, such a method can also lead to a measurement of the mean value of DM contributed from the local host galaxy.

The detection of the high-energy ($\sim290$ TeV) neutrino coincident with the flaring blazar TXS 0506+056, the first and only $3\sigma$ neutrino-source association to date, provides new, multimessenger tests of the weak equivalence principle (WEP) and Lorentz invariance. Assuming that the flight time difference between the TeV neutrino and gamma-ray photons from the blazar flare is mainly caused by the gravitational potential of the Laniakea supercluster of galaxies, we show that the deviation from the WEP for neutrinos and photons is conservatively constrained to have an accuracy of $10^{-6}-10^{-7}$, which is 3--4 orders of magnitude better than previous results placed by MeV neutrinos from supernova 1987A. In addition, we demonstrate that the association of the TeV neutrino with the blazar flare sets limits on the energy scales of quantum gravity for both linear and quadratic violations of Lorentz invariance (LIV) to $E_{\rm QG, 1}>3.2\times10^{15}-3.7\times10^{16}$ GeV and $E_{\rm QG, 2}>4.0\times10^{10}-1.4\times10^{11}$ GeV. These improve previous limits on both linear and quadratic LIV energy scales in neutrino propagation by 5--7 orders of magnitude.

Fast Radio bursts (FRBs), bright transients with millisecond durations at $\sim$ GHz and typical redshifts probably $>0.8$, are likely to be gravitationally lensed by intervening galaxies. Since the time delay between images of strongly lensed FRB can be measured to extremely high precision because of the large ratio $\sim10^9$ between the typical galaxy-lensing delay time $\sim\mathcal{O}$(10 days) and the width of bursts $\sim\mathcal{O}$(ms), we propose strongly lensed FRBs as precision probes of the universe. We show that, within the flat $\Lambda$CDM model, the Hubble constant $H_0$ can be constrained with a $\sim0.91\%$ uncertainty from 10 such systems probably observed with the Square Kilometer Array (SKA) in $<$ 30 years. More importantly, the cosmic curvature can be model-independently constrained to a precision of $\sim0.076$. This constraint can directly test the validity of the cosmological principle and break the intractable degeneracy between the cosmic curvature and dark energy.

We present the X-ray timing results of the new black hole candidate (BHC) MAXI J1535-571 during its 2017 outburst from Hard X-ray Modulation Telescope (\emph{Insight}-HXMT) observations taken from 2017 September 6 to 23. Following the definitions given by \citet{Belloni2010}, we find that the source exhibits state transitions from Low/Hard state (LHS) to Hard Intermediate state (HIMS) and eventually to Soft Intermediate state (SIMS). Quasi-periodic oscillations (QPOs) are found in the intermediate states, which suggest different types of QPOs. With the large effective area of \emph{Insight}-HXMT at high energies, we are able to present the energy dependence of the QPO amplitude and centroid frequency up to 100 keV which is rarely explored by previous satellites. We also find that the phase lag at the type-C QPOs centroid frequency is negative (soft lags) and strongly correlated with the centroid frequency. By assuming a geometrical origin of type-C QPOs, the source is consistent with being a high inclination system.

Using high-quality optical images from the Carnegie-Irvine Galaxy Survey, we perform multi-component decompositions of S0 galaxies to derive accurate structural parameters to constrain their physical origin. Many S0s do not host prominent bulges. S0 galaxies have a broad distribution of bulge-to-total ratios ($B/T$) and S\'ersic indices ($n$), with average values of $B/T=0.34\pm0.15$ and $n=2.62\pm1.02$, qualitatively consistent with the notion that S0s define a parallel sequence with and may have evolved from spiral galaxies. This is further reinforced by the incidence of bars and lenses in S0s, which when compared with the statistics in spirals, supports the idea that lenses are demised bars. However, despite their wide range of prominence, the bulges of S0s form a surprisingly uniform sequence on both the Kormendy and fundamental plane relations. There is no evidence for population dichotomy or other indications of differentiation into classical and pseudo bulges. Most of the S0s reside in the field and in groups; cluster environment is not a necessary condition for S0 production. The properties of S0 bulges show little correlation with environmental indicators, after the dependence of galaxy stellar mass on environment is taken into account. As the bulges of late-type spirals and S0s are intrinsically different, and environmental effects that may account for such evolution appear to be minimal, we conclude that late-type spirals are not plausible progenitors of S0s. The bulges of S0s likely formed at an early epoch, after which secular processes contributed little to their subsequent evolution.

The first gravitational-wave event from the merger of a binary neutron star system (GW170817) was detected recently. The associated short gamma-ray burst (GRB 170817A) has a low isotropic luminosity ($\sim 10^{47}$ erg s$ ^{-1}$) and a peak energy $E_{p}\sim 145$ keV during the initial main emission between -0.3 and 0.4 s. The origin of this short GRB is still under debate, but a plausible interpretation is that it is due to the off-axis emission from a structured jet. We consider two possibilities. First, since the best-fit spectral model for the main pulse of GRB 170817A is a cutoff power law with a hard low-energy photon index ($\alpha =-0.62_{-0.54}^{+0.49} $), we consider an off-axis photosphere model. We develop a theory of photosphere emission in a structured jet and find that such a model can reproduce a low-energy photon index that is softer than a blackbody through enhancing high-latitude emission. The model can naturally account for the observed spectrum. The best-fit Lorentz factor along the line of sight is $\sim 20$, which demands that there is a significant delay between the merger and jet launching. Alternatively, we consider that the emission is produced via synchrotron radiation in an optically thin region in an expanding jet with decreasing magnetic fields. This model does not require a delay of jet launching but demands a larger bulk Lorentz factor along the line of sight. We perform Markov Chain Monte Carlo fitting to the data within the framework of both models and obtain good fitting results in both cases.

Recently, some theoretical models predicted that a small fraction of fast radio bursts (FRBs) could be associated with gravitational waves (GWs). In this work, we discuss the possibility of using GW/FRB association systems, if commonly detected in the future, as a complementary cosmic probe. We propose that upgraded standard sirens can be constructed from the joint measurements of luminosity distances $D_{L}$ derived from GWs and dispersion measures ${\rm DM_{IGM}}$ derived from FRBs (i.e., the combination $D_{L}\cdot{\rm DM_{IGM}}$). Moreover, unlike the traditional standard-siren approach (i.e., the $D_{L}$ method) and the $\rm DM_{IGM}$ method that rely on the optimization of the Hubble constant $H_0$, this $D_{L}\cdot{\rm DM_{IGM}}$ method has the advantage of being independent of $H_0$. Through Monte Carlo simulations, we prove that the $D_{L}\cdot{\rm DM_{IGM}}$ method is more effective to constrain cosmological parameters than $D_{L}$ or ${\rm DM_{IGM}}$ separately, and enables us to achieve accurate multimessenger cosmology from around 100 GW/FRB systems. Additionally, even if GW/FRB associations are not exist, the methodology developed here can still be applied to those GWs and FRBs that occur at the same redshifts.

Limited by the sensitivities of the current gravitational wave (GW) detectors, the central remnant of the binary neutron star (NS) merger associated with GW170817 remains an open question. Considering the relatively large total mass, it is generally proposed that the merger of GW170817 would lead to a shortly lived hypermassive NS or directly produce a black hole (BH). There is no clear evidence to support or rule out a long-lived NS as the merger remnant. Here we utilize the GW and electromagnetic (EM) signals to comprehensively investigate the parameter space that allows a long-lived NS to survive as the merger remnant of GW170817. We find that for some stiff equations of state, the merger of GW170817 could, in principle, lead to a massive NS, which has a millisecond spin period. The post-merger GW signal could hardly constrain the ellipticity of the NS. If the ellipticity reaches $10^{-3}$, in order to be compatible with the multi-band EM observations, the dipole magnetic field of the NS ($B_p$) is constrained to the magnetar level of $\sim10^{14}$ G. If the ellipticity is smaller than $10^{-4}$, $B_p$ is constrained to the level of $\sim10^{10}-10^{12}\,$G. These conclusions weakly depend on the adoption of equations of state.

Recently, the LIGO-Virgo collaboration reported their first detection of gravitational wave (GW) signals from a low mass compact binary merger GW170817, which is most likely due to a double neutron star (NS) merger. With the GW signals only, the chirp mass of the binary is precisely constrained to $1.188^{+0.004}_{-0.002}~\rm{M_{\odot}}$, but the mass ratio is loosely constrained in the range $0.4-1$, so that a very rough estimation of the individual NS masses ($0.86~{\rm M_{\odot}}<M_1<1.36~\rm{M_{\odot}}$ and $1.36~{\rm M_{\odot}}<M_2<2.26~\rm{M_{\odot}}$) was obtained. Here we propose that if one can constrain the dynamical ejecta mass through performing kilonova modeling of the optical/IR data, by utilizing an empirical relation between the dynamical ejecta mass and the mass ratio of NS binaries, one may place a more stringent constraint on the mass ratio of the system. For instance, considering that the red "kilonova" component is powered by the dynamical ejecta, we reach a tight constraint on the mass ratio in the range of $0.46-0.59$. Alternatively, if the blue "kilonova" component is powered by the dynamical ejecta, the mass ratio would be constrained in the range of $0.53-0.67$. Overall, such a multi-messenger approach could narrow down the mass ratio of GW170817 system to the range of $0.46-0.67$, which gives a more precise estimation of the individual NS mass than pure GW signal analysis, i.e. $0.90~{\rm M_{\odot}}<M_1<1.16~{\rm M_{\odot}}$ and $1.61~{\rm M_{\odot}}<M_2<2.11~{\rm M_{\odot}}$.

The recent observations of GW170817 and its electromagnetic (EM) counterparts show that double neutron star mergers could lead to rich and bright EM emissions. Recent numerical simulations suggest that neutron star and neutron star/black hole (NS-NS/BH) mergers would leave behind a central remnant surrounded by a mildly isotropic ejecta. The central remnant could launch a collimated jet and when the jet propagating through the ejecta, a mildly relativistic cocoon would be formed and the interaction between the cocoon and the ambient medium would accelerate electrons via external shock in a wide angle. So that the merger-nova photons (i.e., thermal emission from the ejecta) would be scattered into higher frequency via inverse compton (IC) process when they propagating through the cocoon shocked region. We find that the IC scattered component peaks at X-ray band and it will reach its peak luminosity in order of days (simultaneously with the merger-nova emission). With current X-ray detectors, such a late X-ray component could be detected out to 200 Mpc, depending on the merger remnant properties. It could serve as an important electromagnetic counterpart of gravitational wave signals from NS-NS/BH mergers. Nevertheless, simultaneous detection of such a late X-ray signal and the merger-nova signal could shed light on the cocoon properties and the concrete structure of the jet.

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.

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.

The coincident detection of a gravitational-wave (GW) event GW170817 with electromagnetic (EM) signals (e.g., a short gamma-ray burst SGRB 170817A or a macronova) from a binary neutron star merger within the nearby galaxy NGC 4933 provides a new, multimessenger test of the weak equivalence principle (WEP), extending the WEP test with GWs and photons. Assuming that the arrival time delay between the GW signals from GW170817 and the photons from SGRB 170817A or the macronova is mainly attributed to the gravitational potential of the Milky Way, we demonstrate that the strict upper limits on the deviation from the WEP are $\Delta \gamma<1.4\times10^{-3}$ for GW170817/macronova and $\Delta \gamma <5.9\times10^{-8}$ for GW170817/SGRB 170817A. A much more severe constraint on the WEP accuracy can be achieved ($\sim0.9\times10^{-10}$) for GW170817/SGRB 170817A when we consider the gravitational potential of the Virgo Cluster, rather than the Milky Way's gravity. This provides the tightest limit to date on the WEP through the relative differential variations of the $\gamma$ parameter for two different species of particles. Compared with other multimessenger (photons and neutrinos) results, our limit is 7 orders of magnitude tighter than that placed by the neutrinos and photons from supernova 1987A, and is almost as good as or is an improvement of 6 orders of magnitude over the limits obtained by the low-significance neutrinos correlated with GRBs and a blazar flare.

The development of two-dimensional (2D) bulge-to-disk decomposition techniques has shown their advantages over traditional one-dimensional (1D) techniques, especially for galaxies with non-axisymmetric features. However, the full potential of 2D techniques has yet to be fully exploited. Secondary morphological features in nearby disk galaxies, such as bars, lenses, rings, disk breaks, and spiral arms, are seldom accounted for in 2D image decompositions, even though some image-fitting codes, such as GALFIT, are capable of handling them. We present detailed, 2D multi-model and multi-component decomposition of high-quality $R$-band images of a representative sample of nearby disk galaxies selected from the Carnegie-Irvine Galaxy Survey, using the latest version of GALFIT. The sample consists of five barred and five unbarred galaxies, spanning Hubble types from S0 to Sc. Traditional 1D decomposition is also presented for comparison. In detailed case studies of the 10 galaxies, we successfully model the secondary morphological features. Through a comparison of best-fit parameters obtained from different input surface brightness models, we identify morphological features that significantly impact bulge measurements. We show that nuclear and inner lenses/rings and disk breaks must be properly taken into account to obtain accurate bulge parameters, whereas outer lenses/rings and spiral arms have a negligible effect. We provide an optimal strategy to measure bulge parameters of typical disk galaxies, as well as prescriptions to estimate realistic uncertainties of them, which will benefit subsequent decomposition of a larger galaxy sample.

Finding the electromagnetic (EM) counterpart of binary compact star merger, especially the binary neutron star (BNS) merger, is critically important for gravitational wave (GW) astronomy, cosmology and fundamental physics. On Aug. 17, 2017, Advanced LIGO and \textit{Fermi}/GBM independently triggered the first BNS merger, GW170817, and its high energy EM counterpart, GRB 170817A, respectively, resulting in a global observation campaign covering gamma-ray, X-ray, UV, optical, IR, radio as well as neutrinos. The High Energy X-ray telescope (HE) onboard \textit{Insight}-HXMT (Hard X-ray Modulation Telescope) is the unique high-energy gamma-ray telescope that monitored the entire GW localization area and especially the optical counterpart (SSS17a/AT2017gfo) with very large collection area ($\sim$1000 cm$^2$) and microsecond time resolution in 0.2-5 MeV. In addition, \textit{Insight}-HXMT quickly implemented a Target of Opportunity (ToO) observation to scan the GW localization area for potential X-ray emission from the GW source. Although it did not detect any significant high energy (0.2-5 MeV) radiation from GW170817, its observation helped to confirm the unexpected weak and soft nature of GRB 170817A. Meanwhile, \textit{Insight}-HXMT/HE provides one of the most stringent constraints (~10$^{-7}$ to 10$^{-6}$ erg/cm$^2$/s) for both GRB170817A and any other possible precursor or extended emissions in 0.2-5 MeV, which help us to better understand the properties of EM radiation from this BNS merger. Therefore the observation of \textit{Insight}-HXMT constitutes an important chapter in the full context of multi-wavelength and multi-messenger observation of this historical GW event.

The ejecta composition is an open question in gamma-ray bursts (GRB) physics. Some GRBs possess a quasi-thermal spectral component in the time-resolved spectral analysis, suggesting a hot fireball origin. Others show a featureless non-thermal spectrum known as the "Band" function, consistent with a synchrotron radiation origin and suggesting that the jet is Poynting-flux-dominated at the central engine and likely in the emission region as well. There are also bursts showing a sub-dominant thermal component and a dominant synchrotron component suggesting a likely hybrid jet composition. Here we report an extraordinarily bright GRB 160625B, simultaneously observed in gamma-rays and optical wavelengths, whose prompt emission consists of three isolated episodes separated by long quiescent intervals, with the durations of each "sub-burst" being $\sim$ 0.8 s, 35 s, and 212 s, respectively. Its high brightness (with isotropic peak luminosity L$_{\rm p, iso}\sim 4\times 10^{53}$ erg/s) allows us to conduct detailed time-resolved spectral analysis in each episode, from precursor to main burst and to extended emission. The spectral properties of the first two sub-bursts are distinctly different, allowing us to observe the transition from thermal to non-thermal radiation between well-separated emission episodes within a single GRB. Such a transition is a clear indication of the change of jet composition from a fireball to a Poynting-flux-dominated jet.

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

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

Neutron stars may sustain a non-axisymmetric deformation due to magnetic distortion and are potential sources of continuous gravitational waves (GWs) for ground-based interferometric detectors. With decades of searches using available GW detectors, no evidence of a GW signal from any pulsar has been observed. Progressively stringent upper limits of ellipticity have been placed on Galactic pulsars. In this work, we use the ellipticity inferred from the putative millisecond magnetars in short gamma-ray bursts (SGRBs) to estimate their detectability by current and future GW detectors. For $\sim 1$ ms magnetars inferred from the SGRB data, the detection horizon is $\sim 30$ Mpc and $\sim 600$ Mpc for advanced LIGO (aLIGO) and Einstein Telescope (ET), respectively. Using the ellipticity of SGRB millisecond magnetars as calibration, we estimate the ellipticity and gravitational wave strain of Galactic pulsars and magnetars assuming that the ellipticity is magnetic-distortion-induced. We find that the results are consistent with the null detection results of Galactic pulsars and magnetars with the aLIGO O1. We further predict that the GW signals from these pulsars/magnetars may not be detectable by the currently designed aLIGO detector. The ET detector may be able to detect some relatively low frequency signals ($<50$ Hz) from some of these pulsars. Limited by its design sensitivity, the eLISA detector seems not suitable for detecting the signals from Galactic pulsars and magnetars.

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

Zhang (2013) proposed a type of GRB-less X-ray transient associated with double neutron star (NS-NS) mergers under the conjecture of a rapidly-spinning magnetar merger product with the line of sight off the short GRB jet. We investigate possible light curves of these transients by considering different observer's viewing angles. We perform Monte Carlo simulations to calculate the peak luminosity function (LF) and event rate density of these X-ray transients. By considering that a fraction of massive neutron stars may be supra-massive and later collapse into black holes after spinning down, we investigate how the predicted LF depends on the equation of state (EoS) of the central object and the geometry of the system. In general, the LF can be fit by two log-normal distributions peaking around $10^{46.4}$ and $10^{49.6}$ $\rm erg\,s^{-1}$, corresponding to the trapped and free zones, respectively. For the majority of the EoS models, the current non-detection is consistent with having a free zone solid angle at most a few times of the solid angle of the short GRB jet. The event rate density of these X-ray transients is around a few tens of $\rm Gpc^{-3}yr^{-1}$ for luminosity above $10^{45}$ $\rm erg\,s^{-1}$. We predict that future X-ray telescopes (such as Einstein Probe) with sensitivity $\sim 10^{-11}$ $\rm erg\,s^{-1}\,cm^{-2}$ would detect as many as several tens of such transients per year per steradian. Within 200 Mpc, the aLIGO average range for NS-NS mergers, the estimated event rate of these transients is about 1 transient per year all sky.

Einstein's weak equivalence principle (WEP) states that any freely falling, uncharged test particle follows the same identical trajectory independent of its internal structure and composition. Since the polarization of a photon is considered to be part of its internal structure, we propose that polarized photons from astrophysical transients, such as gamma-ray bursts (GRBs) and fast radio bursts (FRBs), can be used to constrain the accuracy of the WEP through the Shapiro time delay effect. Assuming that the arrival time delays of photons with different polarizations are mainly attributed to the gravitational potential of the Laniakea supercluster of galaxies, we show that a strict upper limit on the differences of the parametrized post-Newtonian parameter $\gamma$ value for the polarized optical emission of GRB 120308A is $\Delta\gamma<1.2\times10^{-10}$, for the polarized gamma-ray emission of GRB 100826A is $\Delta\gamma<1.2\times10^{-10}$, and for the polarized radio emission of FRB 150807 is $\Delta\gamma<2.2\times10^{-16}$. These are the first direct verifications of the WEP for multiband photons with different polarizations. In particular, the result from FRB 150807 provides the most stringent limit to date on a deviation from the WEP, improving by one order of magnitude the previous best result based on Crab pulsar photons with different energies.

We performed a global study on the timing and spectral properties of type-B quasi-periodic oscillations (QPOs) in the outbursts of black hole X-ray binaries. The sample is built based on the observations of {\it Rossi X-ray Timing Explorer}, via searching in the literature in RXTE era for all the identified type-B QPOs. To enlarge the sample, we also investigated some type-B QPOs that are reported but not yet fully identified. Regarding to the time lag and hard/soft flux ratio, we found that the sources with type-B QPOs behave in two subgroups. In one subgroup, type-B QPO shows a hard time lag that firstly decreases and then reverse into a soft time lag along with softening of the energy spectrum. In the other subgroup, type-B QPOs distribute only in a small region with hard time lag and relatively soft hardness. These findings may be understood with a diversity of the homogeneity showing up for the hot inner flow of different sources. We confirm the universality of a positive relation between the type-B QPO frequency and the hard component luminosity in different sources. We explain the results by considering that the type-B QPO photons are produced in the inner accretion flow around the central black hole, under a local Eddington limit. Using this relationship, we derived an mass estimation of $9.3-27.1$ $\rm{M_{\odot}}$ for the black hole in H 1743-322.

A $\gamma$-ray transient, Swift J0644.5-5111, has been claimed to be associated with FRB 131104. The $\gamma$-ray energy output is estimated as $E_\gamma \approx 5\times 10^{51}$\,erg at the nominal $z\approx 0.55$ redshift implied by the dispersion measure of FRB 131104. However, a long-term radio imaging follow-up observations only place an upper limit on the radio afterglow flux of Swift J0644.5-5111. Applying the external shock model, we make a detailed constraint on the afterglow parameters for the FRB 131104/Swift J0644.5-5111 system. We find that for the commonly used microphysics shock parameters (e.g., $\epsilon_e=0.1$, $\epsilon_B=0.01$ and $p=2.3$), if the redshift value inferred from the DM value is correct to order of magnitude (i.e., $z>0.1$), the ambient medium number density should be $\leq 10^{-3}~\rm{cm^{-3}}$, which is the typical value for a compact binary merger environment but disfavors a massive star origin. Assuming a typical ISM density, one would require that the redshift of the FRB much smaller than the value inferred from DM value ($z\ll0.1$), implying a non-cosmological origin of DM. The constraints are much looser if one adopts smaller $\epsilon_B$ and $\epsilon_e$ values, as observed in some GRB afterglows. The FRB 131104/Swift J0644.5-5111 association remains plausible. We critically discuss possible progenitor models for the system.

To take advantage of the astrophysical potential of Gamma-Ray Bursts (GRBs), Chinese and French astrophysicists have engaged the SVOM mission (Space-based multi-band astronomical Variable Objects Monitor). Major advances in GRB studies resulting from the synergy between space and ground observations, the SVOM mission implements space and ground instrumentation. The scientific objectives of the mission put a special emphasis on two categories of GRBs: very distant GRBs at z$>$5 which constitute exceptional cosmological probes, and faint/soft nearby GRBs which allow probing the nature of the progenitors and the physics at work in the explosion. These goals have a major impact on the design of the mission: the on-board hard X-ray imager is sensitive down to 4 keV and computes on line image and rate triggers, and the follow-up telescopes on the ground are sensitive in the NIR. At the beginning of the next decade, SVOM will be the main provider of GRB positions and spectral parameters on very short time scale. The SVOM instruments will operate simultaneously with a wide range of powerful astronomical devices. This rare instrumental conjunction, combined with the relevance of the scientific topics connected with GRB studies, warrants a remarkable scientific return for SVOM. In addition, the SVOM instrumentation, primarily designed for GRB studies, composes a unique multi-wavelength observatory with rapid slew capability that will find multiple applications for the whole astronomy community beyond the specific objectives linked to GRBs. This report lists the scientific themes that will benefit from observations made with SVOM, whether they are specific GRB topics, or more generally all the issues that can take advantage of the multi-wavelength capabilities of SVOM.

A supramassive, strongly-magnetized millisecond neutron star (NS) has been proposed to be the candidate central engine of at least some short gamma-ray bursts (SGRBs), based on the "internal plateau" commonly observed in the early X-ray afterglow. While a previous analysis shows a qualitative consistency between this suggestion and the Swift SGRB data, the distribution of observed break time $t_b$ is much narrower than the distribution of the collapse time of supramassive NSs for the several NS equations-of-state (EoSs) investigated. In this paper, we study four recently-constructed "unified" NS EoSs, as well as three developed strange quark star (QS) EoSs within the new confinement density-dependent mass model. All the EoSs chosen here satisfy the recent observational constraints of the two massive pulsars whose masses are precisely measured. We construct sequences of rigidly rotating NS/QS configurations with increasing spinning frequency $f$, from non-rotating ($f = 0$) to the Keplerian frequency ($f = f_{\rm K}$), and provide convenient analytical parametrizations of the results. Assuming that the cosmological NS-NS merger systems have the same mass distribution as the Galactic NS-NS systems, we demonstrate that all except the BCPM NS EoS can reproduce the current $22\%$ supramassive NS/QS fraction constraint as derived from the SGRB data. We simultaneously simulate the observed quantities (the break time $t_b$, the break time luminosity $L_b$ and the total energy in the electromagnetic channel $E_{\rm total}$) of SGRBs, and find that while equally well reproducing other observational constraints, QS EoSs predict a much narrower $t_b$ distribution than that of the NS EoSs, better matching the data. We therefore suggest that the post-merger product of NS-NS mergers might be fast-rotating supramassive QSs rather than NSs.

Five TeV neutrino events weakly correlated with five gamma-ray bursts (GRBs) were detected recently by IceCube. This work is an attempt to show that if the GRB identifications are verified, the observed time delays between the TeV neutrinos and gamma-ray photons from GRBs provide attractive candidates for testing fundamental physics with high accuracy. Based on the assumed associations between the TeV neutrinos and GRBs, we find that the limiting velocity of the neutrinos is equal to that of photons to an accuracy of $\sim1.9\times10^{-15}-2.5\times10^{-18}$, which is about $10^{4}-10^{7}$ times better than the constraint obtained with the neutrino possibly from a blazar flare. In addition, we set the most stringent limits up to date on the energy scale of quantum gravity for both the linear and quadratic violations of Lorentz invariance, namely $E_{\rm QG, 1}>6.3\times10^{18}-1.5\times10^{21}$ GeV and $E_{\rm QG, 2}>2.0\times10^{11}-4.2\times10^{12}$ GeV, which are essentially as good as or are an improvement of one order of magnitude over the results previously obtained by the GeV photons of GRB 090510 and the PeV neutrino from a blazar flare. Assuming that the Shapiro time delay is caused by the gravitational potential of the Laniakea supercluster of galaxies, we also place the tightest limits to date on Einstein's weak equivalence principle through the relative differential variations of the parameterized post-Newtonian parameter $\gamma$ values for two different species of particles (i.e., neutrinos and photons), yielding $\Delta\gamma \sim 10^{-11}-10^{-13}$. However, it should be emphasized again that these limits here obtained are at best forecast of what could be achieved if the GRB/neutrino correlations would be finally confirmed.

A conservative constraint on the Einstein Weak Equivalence Principle (WEP) can be obtained under the assumption that the observed time delay between correlated particles from astronomical sources is dominated by the gravitational fields through which they move. Current limits on the WEP are mainly based on the observed time delays of photons with different energies. It is highly desirable to develop more accurate tests that include the gravitational wave (GW) sector. The detection by the advanced LIGO/VIRGO systems of gravitational waves will provide attractive candidates for constraining the WEP, extending the tests to gravitational interactions, with potentially higher accuracy. Considering the capabilities of the advanced LIGO/VIRGO network and the source direction uncertainty, we show that the joint detection of GWs and electromagnetic signals could probe the WEP to an accuracy down to $10^{-10}$, which is one order of magnitude tighter than previous limits, and seven orders of magnitude tighter than the multi-messenger (photons and neutrinos) results by supernova 1987A.

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.

We investigate joint low-latency gravitational wave (GW) detection and prompt electromagnetic (EM) follow-up observations of coalescing binary neutron stars (BNSs). Assuming that BNS mergers are associated with short duration gamma ray bursts (SGRBs), we evaluate if rapid EM follow-ups can capture the prompt emission, early engine activity or reveal any potential by-products such as magnetars or fast radio bursts. To examine the expected performance of extreme low-latency search pipelines, we simulate a population of coalescing BNSs and use these to estimate the detectability and localisation efficiency at different times before merger. Using observational SGRB flux data corrected to the range of the advanced GW interferometric detectors, we determine what EM observations could be achieved from low-frequency radio up to high energy $\gamma$-ray. We show that while challenging, breakthrough multi-messenger science is possible through low latency pipelines.

We investigate the prospects for joint low-latency gravitational wave (GW) detection and prompt electromagnetic (EM) follow-up observations of coalescing binary neutron stars (BNSs). For BNS mergers associated with short duration gamma-ray bursts (SGRBs), we for the first time evaluate the feasibility of rapid EM follow-ups to capture the prompt emission, early engine activity or reveal any potential by-products such as magnetars or fast radio bursts. To achieve our goal, we first simulate a population of coalescing BNSs using realistic distributions of source parameters and estimate the detectability and localisation efficiency at different times before merger. We then use a selection of facilities with GW follow-up agreements in place, from low-frequency radio to high energy $\gamma$-ray to assess the prospects of prompt follow-up. We quantify our assessment using observational SGRB flux data extrapolated to be within the horizon distances of the advanced GW interferometric detectors LIGO and Virgo and to the prompt phase immediately following the binary merger. Our results illustrate that while challenging, breakthrough multi-messenger science is possible with EM follow-up facilities with fast responses and wide fields-of-view. We demonstrate that the opportunity to catch the prompt stage ($<$ 5s) of SGRBs, can be enhanced by speeding up the detection pipelines of both GW observatories and EM follow-up facilities. We further show that the addition of an Australian instrument to the optimal detector network could possibly improve the angular resolution by a factor of two and thereby contribute significantly to GW-EM multi-messenger astronomy.

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.

The observed time delays between different energy bands from TeV blazars provide a new interesting way of testing the Einstein Equivalence Principle (EEP). If the whole time delay is assumed to be dominated by the gravitational field of the Milky Way, the conservative upper limit on the EEP can be estimated. Here we show that the strict limits on the differences of the parameterized post-Newtonian parameter $\gamma$ values are $\gamma_{\rm TeV}-\gamma_{\rm keV}<3.86\times10^{-3}$ for Mrk 421 and $\gamma_{\rm TeV}-\gamma_{\rm keV}<4.43\times10^{-3}$ for Mrk 501, while expanding the scope of the tested EEP energy range out to the TeV--keV range for the first time. With the small time lag from the 0.2--0.8 TeV and $>0.8$ TeV light curves of PKS 2155-304, a much more severe constraint on $\gamma$ differences of $\sim10^{-6}$ can be achieved, although the energy difference is of order of $\sim$ TeV. Furthermore, we can combine these limits on the energy dependence of $\gamma$ with the bound on the absolute $\gamma$ value $\gamma-1\sim0.3\%$ from light deflection measurements at the optical (eV) bands, and conclude that this absolute bound on $\gamma$ can be extended from optical to TeV energies.

The accuracy of Einstein's Equivalence Principle (EEP) can be tested with the observed time delays between correlated particles or photons that are emitted from astronomical sources. Assuming as a lower limit that the time delays are caused mainly by the gravitational potential of the Milky Way, we prove that fast radio bursts (FRBs) of cosmological origin can be used to constrain the EEP with high accuracy. Taking FRB 110220 and two possible FRB/gamma-ray burst (GRB) association systems (FRB/GRB 101011A and FRB/GRB 100704A) as examples, we obtain a strict upper limit on the differences of the parametrized post-Newtonian parameter $\gamma$ values as low as $\left[\gamma(1.23\; \rm GHz)-\gamma(1.45\; \rm GHz)\right]<4.36\times10^{-9}$. This provides the most stringent limit up to date on the EEP through the relative differential variations of the $\gamma$ parameter at radio energies, improving by 1 to 2 orders of magnitude the previous results at other energies based on supernova 1987A and GRBs.

Recently the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) published its first data release (DR1), which is ranked as the largest stellar spectra dataset in the world so far. We combine the PASTEL catalog and SIMBAD radial velocities as a testing standard to validate the DR1 stellar parameters (effective temperature $T_{\mathrm{eff}}$, surface gravity $\log g$, metallicity $\mathrm{[Fe/H]}$ and radial velocity $V_{\mathrm{r}}$). Through cross-identification of the DR1 catalogs and the PASTEL catalog, we obtain a preliminary sample of 422 stars. After removal of stellar parameter measurements from problematic spectra and applying effective temperature constraints to the sample, we compare the stellar parameters from DR1 with those from PASTEL and SIMBAD to prove that the DR1 results are reliable in restricted $T_{\mathrm{eff}}$ ranges. We derive standard deviations of 110 K, 0.19 dex, 0.11 dex and 4.91 $\mathrm{km\,s^{-1}}$ , for $T_{\mathrm{eff}}$, $\log g$, $[\mathrm{Fe/H}]$ when $T_{\mathrm{eff}}<8000\,\mathrm{K}$, and for $V_{\mathrm{r}}$ when $T_{\mathrm{eff}}<10000\,\mathrm{K}$, respectively. Systematic errors are negligible except for that of $V_{\mathrm{r}}$. Besides, metallicities in DR1 are systematically higher than those in PASTEL, in the range of PASTEL $[\mathrm{Fe/H}]<-1.5$.

The Einstein Equivalence Principle (EEP) can be probed with astrophysical sources emitting simultaneously different types of neutral particles, or particles with varying energies, by testing their time of flight through the same gravitational field. Here we use the time delays between correlated photons from cosmological transients to constrain the accuracy of the EEP. We take data from two gamma-ray bursts as an example, and use, as a lower limit to the theoretical time delays between different energies, delays arising from only the gravitational field of our own galaxy. We then show that the parameterized post-Newtonian parameter $\gamma$ is the same for photons over energy ranges between eV and MeV and between MeV and GeV to a part in $10^{-7}$, which is at least one order of magnitude better than previous limits. Combining this bound on the wavelength dependence of $\gamma$ with the absolute bound $|\gamma-1|<0.3\%$ from light-deflection measurements at optical (eV) wavelengths, we thus extend this absolute bound on $\gamma$ to GeV energies.

A generic synchrotron external shock model is the widely preferred paradigm used to interpret the broad-band afterglow data of gamma-ray bursts (GRBs), including predicted observable signatures from a reverse shock which have been confirmed by observations. Investigations of the nature of the reverse shock emission can provide valuable insights into the intrinsic properties of the GRB ejecta. Here we briefly review the standard and the extended models of the reverse shock emission, discussing the connection between the theory and observations, including the implications of the latest observational advances.

Within the framework of the external shock model of gamma-ray bursts (GRBs) afterglows, we perform a morphological analysis of the early optical lightcurves to directly constrain model parameters. We define four morphological types, i.e. the reverse shock dominated cases with/without the emergence of the forward shock peak (Type I/ Type II), and the forward shock dominated cases without/with $\nu_m$ crossing the band (Type III/IV). We systematically investigate all the Swift GRBs that have optical detection earlier than 500 s and find 3/63 Type I bursts (4.8%), 12/63 Type II bursts (19.0%), 30/63 Type III bursts (47.6%), 8/63 Type IV bursts (12.7%) and 10/63 Type III/IV bursts (15.9%). We perform Monte Carlo simulations to constrain model parameters in order to reproduce the observations. We find that the favored value of the magnetic equipartition parameter in the forward shock ($\epsilon_B^f$) ranges from $10^{-6}$ to $10^{-2}$, and the reverse-to-forward ratio of $\epsilon_B$ ($R_B$) is about 100. The preferred electron equipartition parameter $\epsilon_e^{r,f}$ value is 0.01, which is smaller than the commonly assumed value, e.g., 0.1. This could mitigate the so- called "efficiency problem" for the internal shock model, if $\epsilon_e$ during the prompt emission phase (in the internal shocks) is large (say, $\sim 0.1$). The preferred $R_B$ value is in agreement with the results in previous works that indicates a moderately magnetized baryonic jet for GRBs.

GRB 080503 is a short gamma-ray burst (GRB) detected by \emph{Swift} and has been classified as a compact-star-merger-origin GRB. The soft extended emission and the simultaneous late re-brightening in both the X-ray and optical afterglow lightcurves raise interesting questions regarding its physical origin. We show that the broad-band data of GRB 080503 can be well explained within the framework of the double neutron star merger model, provided that the merger remnant is a rapidly-rotating massive neutron star with an extremely high magnetic field (i.e. a millisecond magnetar). We show that the late optical re-brightening is consistent with the emission from a magnetar-powered "merger-nova". This adds one more case to the growing sample of merger-novae associated with short GRBs. The soft extended emission and the late X-ray excess emission are well connected through a magnetar dipole spin-down luminosity evolution function, suggesting that direct magnetic dissipation is the mechanism to produce these X-rays. The X-ray emission initially leaks from a hole in the merger ejecta pierced by the short GRB jet. The hole subsequently closes after the magnetar spins down and the magnetic pressure drops below ram pressure. The X-ray photons are then trapped behind the mergernova ejecta until the ejecta becomes optically thin at a later time. This explains the essentially simultaneous re-brightening in both the optical and X-ray lightcurves. Within this model, future gravitational wave sources could be associated with a bright X-ray counterpart along with the mergernova, even if the short GRB jet beams away from Earth.

Gamma-Ray Bursts are the most energetic explosions in the Universe, and are among the most promising for detecting multiple non-electromagnetic signals, including cosmic rays, high energy neutrinos and gravitational waves. The multi-GeV to TeV gamma-ray range of GRB could have significant contributions from hadronic interactions, mixed with more conventional leptonic contributions. This energy range is important for probing the source physics, including overall energetics, the shock parameters and the Lorentz factor. We discuss some of the latest observational and theoretical developments in the field.

The external forward shock (EFS) models have been the standard paradigm to interpret the broad-band afterglow data of gamma-ray bursts (GRBs). One prediction of the models is that some afterglow temporal breaks at different energy bands should be achromatic. Observations in the Swift era have revealed chromatic afterglow behaviors at least in some GRBs, casting doubts on the EFS origin of GRB afterglows. In this paper, we perform a systematic study to address the question: how bad/good are the external forward shock models? Our sample includes 85 GRBs well-monitored X-ray and optical lightcurves. Based on how well the data abide by the EFS models, we categorize them as: Gold sample: (Grade I and II) include 45/85 GRBs. They show evidence of, or are consistent with having, an achromatic break. The temporal/spectral behaviors in each afterglow segment are consistent with the predictions (closure relations) of the EFS models. Silver sample: (Grade III and IV) include 37/85 GRBs. They are also consistent with having an achromatic break, even though one or more afterglow segments do not comply with the closure relations. Bad sample: (Grade V), 3/85 shows direct evidence of chromatic behaviors, suggesting that the EFS models are inconsistent with the data. These are included in the Bad sample. We further perform statistical analyses of various observational properties ($\alpha$, $\beta$, $t_b$ and model parameters (energy injection index q, p, $\theta_j$, $\eta_\gamma$, etc) of the GRBs in the Gold Sample, and derive constraints on the magnetization parameter $\epsilon_B$ in the EFS. Overall, we conclude that the simplest EFS models can account for the multi-wavelength afterglow data of at least half of the GRBs. When more advanced modeling (e.g., long-lasting reverse shock, structured jets) is invoked, up to $>90 \%$ of the afterglows may be interpreted within the framework of the ESF models.

In quantum theory of gravity, we expect the Lorentz Invariance Violation (LIV) and the modification of the dispersion relation between energy and momentum for photons. The effect of the energy-dependent velocity due to the modified dispersion relation for photons was studied in the standard cosmological context by using a sample of Gamma Ray Bursts (GRBs). In this paper we mainly discuss the possible LIV effect by using different cosmological models for the accelerating universe. Due to the degeneracies among model parameters, the GRBs' time delay data are combined with the cosmic microwave background data from the Planck first year release, the baryon acoustic oscillation data at six different redshifts, as well as Union2 type Ia supernovae data, to constrain both the model parameters and the LIV effect. We find no evidence of LIV.

The GRB central engine intrinsic activity time $T_{\rm ce}$ is usually described through either the $\gamma$-ray duration $T_{90}$ or through a generalized burst duration $t_{\rm burst}$ which includes both the $\gamma$-ray emission and (when present) an extended flaring X-ray plateau. Here, we define a more specific operational description of $T_{\rm ce}$, and within the framework of the internal-external shock model, we develop a numerical code to study the relationship between $T_{90}$ and $T_{\rm ce}$, as well as between $t_{\rm burst}$ and $T_{\rm ce}$, for different initial conditions. We find that when $T_{\rm ce}\lesssim 10^4$ s, late internal collisions or refreshed external collisions result in values of $T_{\rm 90}$ and $t_{\rm burst}$ larger than $T_{\rm ce}$, usually by factors of $2-3$. For $T_{\rm ce}\gtrsim 10^4$ s, the $t_{\rm burst}$ is always a good estimator for $T_{\rm ce}$, while $T_{90}$ can underpredict $T_{\rm ce}$ when the late central engine activity is moderate. We find a clear bimodal distribution for $T_{\rm ce}$, based on our simulations as well as on the observational data for $T_{90}$ and $t_{\rm burst}$. We suggest that $t_{\rm burst}$ is a reliable measure for defining "ultra-long" GRBs. Bursts with $T_{90}$ of order $10^3$ s need not belong to a special population, while bursts with $t_{\rm burst} > 10^4$ s, where the late central engine activity is more moderate and shows up in X-rays, may represent a new population. These conclusions are insensitive to the initial conditions assumed in the models.

We present a wide dataset of gamma-ray, X-ray, UVOIR, and radio observations of the Swift GRB100814A. At the end of the slow decline phase of the X-ray and optical afterglow, this burst shows a sudden and prominent rebrightening in the optical band only, followed by a fast decay in both bands. The optical rebrightening also shows chromatic evolution. Such a puzzling behaviour cannot be explained by a single component model. We discuss other possible interpretations, and we find that a model that incorporates a long-lived reverse shock and forward shock fits the temporal and spectral properties of GRB100814A the best.

In view of the recent Fermi observations of GRB prompt emission spectra, we develop a theory of photosphere emission of a hybrid relativistic outflow with a hot fireball component (defined by dimensionless entropy $\eta$) and a cold Poynting-flux component (defined by magnetization $\sigma_0$ at the central engine). We consider the scenarios both without and with sub-photospheric magnetic dissipations. Based on a simplified toy model of jet dynamics, we develop two approaches: a "bottom-up" approach to predict the temperature (for a non-dissipative photosphere) and luminosity of the photosphere emission and its relative brightness for a given pair of $(\eta,\sigma_0)$; and a "top-down" approach to diagnose central engine parameters ($\eta$ and $\sigma_0$) based on the observed quasi-thermal photosphere emission properties. We show that a variety of observed GRB prompt emission spectra with different degrees of photosphere thermal emission can be reproduced by varying $\eta$ and $\sigma_0$ within the non-dissipative photosphere scenario. In order to reproduce the observed spectra, the outflows of most GRBs need to have a significant $\sigma$, both at the central engine, and at the photosphere. The $\sigma$ value at $10^{15}$ cm from the central engine (a possible non-thermal emission site) is usually also greater than unity, so that internal-collision-induced magnetic reconnection and turbulence (ICMART) may be the mechanism to power the non-thermal emission. We apply our top-down approach to GRB 110721A, and find that the temporal evolution behavior of its blackbody component can be well interpreted with a time-varying $(\eta,\sigma_0)$ at the central engine, instead of invoking a varying engine base size $r_0$ as proposed by previous authors.

The physical nature of fast radio bursts (FRBs) is not identified. Detecting electromagnetic counterparts in other wavelengths is essential to measure their distances and to settle down their physical nature. Assuming that at least some of them are of a cosmological origin, we calculate their afterglow lightcurves in multi-wavelengths (X-rays, optical and radio) by assuming a range of their total kinetic energies and redshifts. We focus on forward shock emission, but also consider the possibility that some of them might have bright reverse shock emission. In general, the FRB afterglows are too faint to be detected by current detectors. Only if an FRB has a very low radiative efficiency in radio (hence, a very large kinetic energy), and when it is close enough, can its afterglow be detected in the optical and radio bands. We discuss observational strategies to detect these faint afterglows using future telescopes such as LSST and EVLA.

Recently, both theoretical arguments and observational evidence suggested that a small fraction of fast radio bursts (FRBs) could be associated with gamma-ray bursts (GRBs). If such FRB/GRB association systems are commonly detected in the future, the combination of dispersion measure (DM) derived from FRBs and redshifts derived from GRBs makes these systems a plausible tool to conduct cosmography. We quantify uncertainties in deriving the redshift-dependent DM_{IGM} as a function of z, and test how well dark energy models can be constrained with Monte Carlo simulations. We show that with potentially several 10s of FRB/GRB systems detected in a decade or so, one may reach reasonable constraints on wCDM models. When combined with SN Ia data, unprecedented constraints on dark energy equation of state may be achieved, thanks to the prospects of detecting FRB/GRB systems at relatively high redshifts. The ratio between the mean value < DM_{IGM} (z)> and luminosity distance (D_{L} (z)) is insensitive to dark energy models. This gives the prospects of applying SN Ia data to calibrate < DM_{IGM} (z)> using a relatively small sample of FRB/GRB systems, allowing a reliable constraint on the baryon inhomogeneity distribution as a function of redshift. The methodology developed in this paper can also be applied, if the FRB redshifts can be measured by other means. Some caveats of putting this method into practice are also discussed.

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.

The Palomar Transient Factory (PTF) team recently reported the discovery of a rapidly fading optical transient source, PTF11agg. A long-lived scintillating radio counterpart was identified, but the search for a high energy counterpart showed negative results. The PTF team speculated that PTF11agg may represent a new class of relativistic outbursts. Here we suggest that a neutron star (NS)-NS merger system with a supra-massive magnetar central engine could be a possible source to power such a transient, if our line of sight is not on the jet axis direction of the system. These systems are also top candidates for gravitational wave sources to be detected in the advanced LIGO/Virgo era. We find that the PTF11agg data could be explained well with such a model, suggesting that at least some gravitational wave bursts due to NS-NS mergers may be associated with such a bright electromagnetic counterpart without a \gamma-ray trigger.

A massive millisecond magnetar may survive a merger of a neutron star (NS) binary, which would continuously power the merger ejecta. We develop a generic dynamic model for the merger ejecta with energy injection from the central magnetar. The ejecta emission (the "merger-nova") powered by the magnetar peaks in the UV band and the peak of lightcurve progressively shifts to an earlier epoch with increasing frequency. A magnetar-powered mergernova could have an optical peak brightness comparable to a supernova, which is a few tens or hundreds times brighter than the radioactive-powered merger-novae (the so-called macro-nova or kilo-nova). On the other hand, such a merger-nova would peak earlier and have a significantly shorter duration than that of a supernova. An early collapse of the magnetar could suppress the brightness of the optical emission and shorten its duration. Such millisecond-magnetar-powered merger-novae may be detected from NS-NS merger events without an observed short gamma-ray burst, and could be a bright electromagnetic counterpart for gravitational wave bursts due to NS-NS mergers. If detected, it suggests that the merger leaves behind a massive NS, which has important implications for the equation-of-state of nuclear matter.

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

As the technology of gravitational-wave and neutrino detectors becomes increasingly mature, a multi-messenger era of astronomy is ushered in. Advanced gravitational wave detectors are close to making a ground-breaking discovery of gravitational wave bursts (GWBs) associated with mergers of double neutron stars (NS-NS). It is essential to study the possible electromagnetic (EM) and neutrino emission counterparts of these GWBs. Recent observations and numerical simulations suggest that at least a fraction of NS-NS mergers may leave behind a massive millisecond magnetar as the merger product. Here we show that protons accelerated in the forward shock powered by a magnetar wind pushing the ejecta launched during the merger process would interact with photons generated in the dissipating magnetar wind and emit high energy neutrinos and photons. We estimate the typical energy and fluence of the neutrinos from such a scenario. We find that $\sim$PeV neutrinos could be emitted from the shock front as long as the ejecta could be accelerated to a relativistic speed. The diffuse neutrino flux from these events, even under the most optimistic scenarios, is too low to account for the two events announced by the IceCube Collaboration, but it is only slightly lower than the diffuse flux of GRBs, making it an important candidate for the diffuse background of $\sim$PeV neutrinos. The neutron-pion decay of these events make them a moderate contributor to the sub-TeV gamma-ray diffuse background.

We continue our systematic statistical study of various components in gamma-ray burst (GRB) optical lightcurves. We decompose the early onset bump and the late re-brightening bump with empirical fits and analyze their statistical properties. Among the 146 GRBs that have well-sampled optical lightcurves, the onset and re-brightening bumps are observed in 38 and 26 GRBs, respectively. It is found that the typical rising and decaying slopes for both the onset and re-brightening bumps are ~1.5 and -1.15, respectively. No early onset bumps in the X-ray band are detected to be associated with the optical onset bumps, while an X-ray re-brightening bump is detected for half of the re-brightening optical bumps. The peak luminosity is anti-correlated with the peak time, L_p\propto t_{p}^{-1.81+/-0.32} and L_p\propto t_{p}^{-0.83+/-0.17} for the onset and re-brightening bumps, respectively. Both L_p and the isotropic energy release of the onset bumps are correlated with E_{\gamma, iso}, whereas no similar correlation is found for the re-brightening bumps. Taking the onset bumps as probes for the properties of the fireballs and their ambient medium, we find that the typical power-law index of the relativistic electrons is 2.5 and the medium density profile behaves as n\propto r^{-1} within the framework of the synchrotron external shock models. With the medium density profile obtained from our analysis, we also confirm the correlation between initial Lorentz factor (\Gamma_0) and E_{\gamma, iso} in our previous work. The jet component that produces the re-brightening bump seems to be on-axis and independent of the prompt emission jet component. Its typical kinetic energy budget would be about one order of magnitude larger than the prompt emission component, but with a lower \Gamma_0, typically several tens.

If double neutron star mergers leave behind a massive magnetar rather than a black hole, a bright early afterglow can follow the gravitational wave burst (GWB) even if there is no short gamma-ray burst (SGRB) - GWB association or there is an association but the SGRB does not beam towards earth. Besides directly dissipating the proto-magnetar wind as suggested by Zhang, we here suggest that the magnetar wind could push the ejecta launched during the merger process, and under certain conditions, would reach a relativistic speed. Such a magnetar-powered ejecta, when interacting with the ambient medium, would develop a bright broad-band afterglow due to synchrotron radiation. We study this physical scenario in detail, and present the predicted X-ray, optical and radio light curves for a range of magnetar and ejecta parameters. We show that the X-ray and optical lightcurves usually peak around the magnetar spindown time scale (10^3-10^5s), reaching brightness readily detectable by wide-field X-ray and optical telescopes, and remain detectable for an extended period. The radio afterglow peaks later, but is much brighter than the case without a magnetar energy injection. Therefore, such bright broad-band afterglows, if detected and combined with GWBs in the future, would be a probe of massive millisecond magnetars and stiff equation-of-state for nuclear matter.

After 2010, three outbursts of H 1743-322 were detected by RXTE. We have carried out timing and spectral analysis of the data, emphasizing on the two with relatively complete evolution history presented in the RXTE/PCA observations. We then constitute an enlarged outburst sample for H 1743-322 which allows to investigate the spectral transitions in more details. We find that the spectral transitions to high-soft state constrain a region for four outbursts in hardness-intensity diagram. An extension of the region locates in the vicinity of the failed outburst in 2008, and excludes it from a successful group. We therefore suggest the failed outburst in 2008 may present the first almost successful outburst sample, which is important for modelling the outburst, especially upon the conditions required for transition to high-soft state.

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.

Well-sampled optical lightcurves of 146 gamma-ray bursts (GRBs) are compiled from the literature. By empirical fitting we identify eight possible emission components and summarize the results in a "synthetic" lightcurve. Both optical flare and early shallow-decay components are likely related to long-term central engine activities. We focus on their statistical properties in this paper. Twenty-four optical flares are obtained from 19 GRBs. The isotropic R-band energy is smaller than 1% of $E_{\gamma, \rm iso}$. The relation between isotropic luminosities of the flares and gamma-rays follows $L^{\rm F}_{\rm R, iso}\propto L_{{\gamma}, \rm iso}^{1.11\pm 0.27}$. Later flares tend to be wider and dimmer, i.e., $w^{\rm F}\sim t^{\rm F}_{\rm p}/2$ and $L^{\rm F}_{\rm R, iso}\propto [t^{\rm F}_{\rm p}/(1+z)]^{-1.15\pm0.15}$. The detection probability of the optical flares is much smaller than that of X-ray flares. An optical shallow decay segment is observed in 39 GRBs. The relation between the break time and break luminosity is a power-law, with an index of $-0.78\pm 0.08$, similar to that derived from X-ray flares. The X-ray and optical breaks are usually chromatic, but a tentative correlation is found. We suggest that similar to the prompt optical emission that tracks $\gamma$-rays, the optical flares are also related to the erratic behavior of the central engine. The shallow decay component is likely related to a long-lasting spinning-down central engine or piling up of flare materials onto the blastwave. Mixing of different emission components may be the reason of the diverse chromatic afterglow behaviors.

Gamma-ray bursts (GRBs) have variable lightcurves. Although most models attribute the observed variability to one physical origin (e.g. central engine activity, clumpy circumburst medium, relativistic turbulence), some models invoke two physically distinct variability components. We develop a method, namely, the stepwise filter correlation (SFC) method, to decompose the variability components in a GRB lightcurve. Based on a low-pass filter technique, we progressively filter the high frequency signals from the lightcurve, and then perform a correlation analysis between each adjunct pair of filtered lightcurves. Our simulations suggest that if a mock lightcurve contains a slow variability component superposed on a rapidly varying time sequence, the correlation coefficient as a function of the filter frequency would display a prominent dip feature around the frequency of the slow component. Through simulations, we demonstrate that this method can identify significant clustering structures of a lightcurve in the frequency domain, and proved that it can catch superposed signals that are otherwise not easy to retrieve based on other methods (e.g. the power density spectrum analysis method). We apply this method to 266 BATSE bright GRBs. We find that the majority of the bursts have clear evidence of such a superposition effect. We perform a statistical analysis of the identified variability components, and discuss the implications for GRB physics.

For the use of Gamma-Ray Bursts (GRBs) to probe cosmology in a cosmology-independent way, a new method has been proposed to obtain luminosity distances of GRBs by interpolating directly from the Hubble diagram of SNe Ia, and then calibrating GRB relations at high redshift. In this paper, following the basic assumption in the interpolation method that objects at the same redshift should have the same luminosity distance, we propose another approach to calibrate GRB luminosity relations with cosmographic fitting directly from SN Ia data. In cosmography, there is a well-known fitting formula which can reflect the Hubble relation between luminosity distance and redshift with cosmographic parameters which can be fitted from observation data. Using the Cosmographic fitting results from the Union set of SNe Ia, we calibrate five GRB relations using GRB sample at $z\leq1.4$ and deduce distance moduli of GRBs at $1.4< z \leq 6.6$ by generalizing above calibrated relations at high redshift. Finally, we constrain the dark energy parameterization models of the Chevallier-Polarski-Linder (CPL) model, the Jassal-Bagla-Padmanabhan (JBP) model and the Alam model with GRB data at high redshift, as well as with the Cosmic Microwave Background radiation (CMB) and the baryonic acoustic oscillation (BAO) observations, and we find the $\Lambda$CDM model is consistent with the current data in 1-$\sigma$ confidence region.

We present a systematic analysis of the spectral and temporal properties of 17 GRBs co-detected by GBM and LAT on board the Fermi satellite by May 2010. We performed a time-resolved spectral analysis of all the bursts with the finest temporal resolution allowed by statistics, in order to avoid temporal smearing of different spectral components. We found that the time-resolved spectra of 14 out of 17 GRBs are best modeled with the Band function over the entire Fermi spectral range, which may suggest a common origin for emissions detected by LAT and GBM. GRB 090902B and GRB 090510 require the superposition between an MeV component and an extra power law component, with the former having a sharp cutoff above E_p. For GRB 090902B, this MeV component becomes progressively narrower as the time bin gets smaller, and can be fit with a Planck function as the time bin becomes small enough. In general, we speculate that phenomenologically there may be three elemental spectral components : (I) a Band-function component (e.g. in GRB 080916C) that extends in a wide energy range and does not narrow with reducing time bins, which may be of the non-thermal origin; (II) a quasi-thermal component (e.g. in GRB 090902B) with the spectra progressively narrowing with reducing time bins; and (III) another non-thermal power law component extending to high energies. The spectra of different bursts may be decomposed into one or more of these elemental components. We compare this sample with the BATSE sample and investigate some correlations among spectral parameters. We discuss the physical implications of the data analysis results for GRB prompt emission, including jet compositions (matter-dominated vs. Poynting-flux-dominated outflow), emission sites (internal shock, external shock or photosphere), as well as radiation mechanisms (synchrotron, synchrotron self-Compton, or thermal Compton upscattering).

Self-gravitating accretion disks collapse to star-forming(SF) regions extending to the inner edge of the dusty torus in active galactic nuclei (AGNs). A full set of equations including feedback of star formation is given to describe the dynamics of the regions. We explore the role of supernovae explosion (SNexp), acting to excite turbulent viscosity, in the transportation of angular momentum in the regions within 1pc scale. We find that accretion disks with typical rates in AGNs can be driven by SNexp in the regions and metals are produced spontaneously. The present model predicts a metallicity--luminosity relationship consistent with that observed in AGNs. As relics of SF regions, a ring (or belt) consisting of old stars remains for every episode of supermassive black hole activity. We suggest that multiple stellar rings with random directions interact and form a nuclear star cluster after episodes driven by star formation.

It has been argued that the long gamma-ray burst (GRB) of GRB 060614 without an associated supernova (SN) has challenged the current classification and fuel model for long GRBs, and thus a tidal disruption model has been proposed to account for such an event. Since it is difficult to detect SNe for long GRBs at high redshift, the absence of an SN association cannot be regarded as the solid criterion for a new classification of long GRBs similar to GRB 060614, called GRB 060614-type bursts. Fortunately, we now know that there is an obvious periodic substructure observed in the prompt light curve of GRB 060614. We thus use such periodic substructure as a potential criterion to categorize some long GRBs into a new class of bursts, which might have been fueled by an intermediate-mass black hole (IMBH) gulping a star, rather than a massive star collapsing to form a black hole. Therefore, the second criterion to recognize for this new class of bursts is whether they fit the tidal disruption model. From a total of 328 Swift GRBs with accurately measured durations and without SN association, we find 25 GRBs satisfying the criteria for GRB 060614-type bursts: seven of them are with known redshifts and 18 with unknown redshifts. These new bursts are ~6% of the total Swift GRBs, which are clustered into two subclasses: Type I and Type II with considerably different viscous parameters of accretion disks formed by tidally disrupting their different progenitor stars. We suggest that the two different kinds of progenitors are solar-type stars and white dwarfs: the progenitors for four Type I bursts with viscous parameter of around 0.1 are solar-type stars, and the progenitors for 21 Type II bursts with viscous parameter of around 0.3 are white dwarfs. The potential applications of this new class of GRBs as cosmic standard candles are discussed briefly

Using 64 ms count data for long gamma-ray bursts ($T_{90}>$2.6s), we analyze the quantity, relative spectral lag (RSL), which is defined as $\tau_{31}/FWHM_{(1)}$, where $\tau_{31}$ is the spectral lag between energy bands 1 and 3, and $FWHM_{(1)}$ denotes full width at half maximum of the pulse in energy channel 1. To get insights into features of the RSLs, we investigate in detail all the correlations between them and other parameters with a sub-sample including nine long bursts. The general cross-correlation technique is adopted to measure the lags between two different energy bands. We can derive the below conclusions. Firstly, the distribution of RSLs is normal and oncentrates on around the value of 0.1. Secondly, the RSLs are weakly correlated with $FWHM$, asymmetry, peak flux ($F_{p}$), peak energy ($E_{p}$) and spectral indexes ($\alpha$ and $\beta$), while they are uncorrelated with $\tau_{31}$, hardness-ratio ($HR_{31}$) and peak time ($t_m$). The final but important discovery is that redshift ($z$) and peak luminosity ($L_{p}$) are strongly correlated with the RSLs which can be measured easily and directly. We find that the RSL is a good redshift and peak luminosity estimator.