A fast radio burst was recently detected to be associated with a hard X-ray burst from the Galactic magnetar SGR 1935+2154. Scenarios involving magnetars for FRBs are hence highly favored. In this work, we suggest that the impact between an asteroid and a magnetar could explain such a detection. According to our calculations, an asteroid of mass $10^{20}$ g will be disrupted at $7 \times 10^9$ cm when approaching the magnetar, the accreted material will flow along the magnetic field lines at Alfv\'en radius $\sim 10^7$ cm. After falling onto the magnetar surface, an instant accretion column will be formed, producing a Comptonized X-ray burst and an FRB in the magnetosphere. We show that all the observational features could be interpreted self-consistently in our scenario. Quasi-periodic oscillations in this specific X-ray burst may help to verify our scenario.

The Tianlai Cylinder Pathfinder is a radio interferometer array designed to test techniques for 21 cm intensity mapping in the post-reionization Universe, with the ultimate aim of mapping the large scale structure and measuring cosmological parameters such as the dark energy equation of state. Each of its three parallel cylinder reflectors is oriented in the north-south direction, and the array has a large field of view. As the Earth rotates, the northern sky is observed by drift scanning. The array is located in Hongliuxia, a radio-quiet site in Xinjiang, and saw its first light in September 2016. In this first data analysis paper for the Tianlai cylinder array, we discuss the sub-system qualification tests, and present basic system performance obtained from preliminary analysis of the commissioning observations during 2016-2018. We show typical interferometric visibility data, from which we derive the actual beam profile in the east-west direction and the frequency band-pass response. We describe also the calibration process to determine the complex gains for the array elements, either using bright astronomical point sources, or an artificial on site calibrator source, and discuss the instrument response stability, crucial for transit interferometry. Based on this analysis, we find a system temperature of about 90 K, and we also estimate the sensitivity of the array.

The characteristics of the spectral evolution of the prompt emission of gamma-ray bursts (GRBs), which are closely related to the radiation mechanism (synchrotron or photosphere), are still an unsolved subject. Here, by performing the detailed time-resolved spectral fitting of GRB 131231A, which has a very bright and well-defined single pulse, some interesting spectral evolution features have been found. (i) Both the low-energy spectral index $\alpha$ and the peak energy $E_{\rm p}$ exhibit the "flux-tracking" pattern ("double-tracking" characteristics). (ii) The parameter relations, i.e., $F$ (the energy flux)-$\alpha$, $F$-$E_{\rm p}$, and $E_{\rm p}$-$\alpha$, along with the analogous Yonetoku $E_{\rm p}$-$L_{\gamma,\rm iso}$ relation for the different time-resolved spectra, show strong monotonous (positive) correlations, both in the rising and the decaying phases. (iii) The values of $\alpha$ do not exceed the synchrotron limit ($\alpha$= -2/3) in all slices across the pulse, favoring the synchrotron origin. We argue that the one-zone synchrotron emission model with the emitter streaming away at a large distance from the central engine can explain all of these special spectral evolution characteristics.

Recently, GRB 170817A was confirmed to be associated with GW 170817, which was produced by a neutron star - neutron star merger. It indicates that at least some short gamma-ray bursts come from binary neutron star mergers. Theoretically, it is widely accepted that short gamma-ray bursts can be produced by two distinctly different mechanisms, binary neutron star mergers and neutron star - black hole mergers. These two kinds of bursts should be different observationally due to their different trigger mechanisms. Motivated by this idea, we collect a universal data set constituted of 51 short gamma-ray bursts observed by $Swift$/BAT, among which 14 events have extended emission component. We study the observational features of these 51 events statistically. It is found that our samples are consisted of two distinct groups. They clearly show a bimodal distribution when their peak photon fluxes at 15-150 keV band are plotted against the corresponding fluences. Most interestingly, all the 14 short bursts with extended emission lie in a particular region. When the fluences are plotted against the burst durations, short bursts with extended emission again tend to concentrate in the long duration segment. These features strongly indicate that short gamma-ray bursts really may come from two distinct types of progenitors. We argue that those short gamma-ray bursts with extended emission come from the coalescence of neutron stars, while the short gamma-ray bursts without extended emission come from neutron star - black hole mergers.

A plateau phase in the X-ray afterglow is observed in a significant fraction of gamma-ray bursts (GRBs). Previously, it has been found that there exists a correlation among three key parameters concerning the plateau phase, i.e., the end time of the plateau phase in the GRB rest frame ($T_{a}$), the corresponding X-ray luminosity at the end time ($L_{X}$) and the isotropic energy of the prompt GRB ($E_{\gamma,\rm{iso}}$). In this study, we systematically search through all the \emph{Swift} GRBs with a plateau phase that occurred between 2005 May and 2018 August. We collect 174 GRBs, with redshifts available for all of them. For the whole sample, the correlation between $L_{X}$, $T_{a}$ and $E_{\gamma,\rm{iso}}$ is confirmed, with the best fit relation being $L_{X}\propto T_{a}^{-1.01}E_{\gamma,\rm{iso}}^{0.84}$. Such an updated three-parameter correlation still supports that the central leftover after GRBs is probably a millisecond magnetar. It is interesting to note that short GRBs with duration less than 2 s in our sample also follow the same correlation, which hints that the merger production of two neutron stars could be a high mass magnetar, but not necessarily a black hole. Moreover, GRBs having an "internal" plateau (i.e., with a following decay index being generally smaller than -3) also obey this correlation. It further strengthens the idea that the internal plateau is due to the delayed collapse of a high mass neutron star into a black hole. The updated three-parameter correlation indicates that GRBs with a plateau phase may act as a standard candle for cosmology study.

Since the true ground state of the hadrons may be strange quark matter (SQM), pulsars may actually be strange stars rather than neutron stars. According to this SQM hypothesis, strange planets can also stably exist. The density of normal matter planets can hardly be higher than 30 g cm$^{-3}$. As a result, they will be tidally disrupted when its orbital radius is less than $\sim 5.6\times10^{10} \rm \, cm $, or when the orbital period ($P_{\rm orb}$) is less than $ \sim \rm 6100\, s $. On the contrary, a strange planet can safely survive even when it is very close to the host, due to its high density. The feature can help us identify SQM objects. In this study, we have tried to search for SQM objects among close-in exoplanets orbiting around pulsars. Encouragingly, it is found that four pulsar planets (XTE J1807-294 b, XTE J1751-305 b, PSR 0636 b, PSR J1807-2459A b) completely meet the criteria of $P_{\rm orb} < \rm 6100\, s $, and are thus good candidates for SQM planets. The orbital periods of two other planets (PSR J1719+14 b and PSR J2051-0827 b) are only slightly higher than the criteria. They could be regarded as potential candidates. Additionally, we find that the periods of five white dwarf planets (GP Com b, V396 Hya b, J1433 b, WD 0137-349 b, and SDSS J1411+2009 b) are less than 0.1 days. We argue that they might also be SQM planets. It is further found that the persistent gravitational wave emissions from at least three of these close-in planetary systems are detectable to LISA. More encouragingly, the advanced LIGO and Einstein Telescope are able to detect the gravitational wave bursts produced by the merger events of such SQM planetary systems, which will provide a unique test for the SQM hypothesis.

The broadband afterglow of GW170817/GRB 170817A is believed to be from an off-axis structured jet. The central engine of a gamma-ray burst usually launches a pair of outflows that move oppositely, it is reasonable to consider the emission from a double-sided structured jets, with a near-jet moving toward us and a counter-jet moving away from us. Assuming that the two branches of the jet have the same physical parameters, we have calculated their radio emission. It is found that the counter-jet component will emerge in the radio light curves $\sim$ 2500 days post-merger.It typically leads to a plateau in the light curve, thus could be hopefully revealed by high accuracy radio observations. We have also considered the possibilities that both branches have different parameters, and found that if some of the physical parameters of the counter-jet are evaluated more favorably than those of the near-jet, then the counter-jet emission will be enhanced and will clearly show up as an obvious plateau or even a rebrightening. For example, when the circum-burst medium encountered by the counter-jet is assumed to be 100 times denser than that of the near-jet, a remarkable radio plateau will appear at $\sim$ 600 days. However, in X-ray bands, the counter-jet component is generally too faint to be discerned. It is argued that the late radio observations of GW170817/GRB 170817A can help determine the key parameters and diagnose the environment of the event.

We perform a series of relativistic magnetohydrodynamics simulations to investigate how a hot magnetic jet propagates within the dynamical ejecta of a binary neutron star merger, with the focus on how the jet structure depends on the delay time of jet launching with respect to the merger time, $\Delta t_{\rm jet}$. We find that regardless of the jet launching delay time, a structured jet with an angle-dependent luminosity and Lorentz factor is always formed after the jet breaks out the ejecta. On the other hand, the jet launching delay time has an impact on the jet structure. If the jet launching delay time is relatively long, e.g., $\ge 0.5$ s, the line-of-sight material has a dominant contribution from the cocoon. On the other hand, for a relatively short jet launching delay time, the jet penetrates through the ejecta early on and develops an angular structure afterward. The line-of-sight ejecta is dominated by the structured jet itself. We discuss the case of GW170817/GRB 170817A within the framework of both long and short $\Delta t_{\rm jet}$. Future more observations of GW/GRB associations can help to differentiate between these two scenarios.

Although the physical origin of gamma-ray burst (GRB) prompt emission is still controversial, synchrotron radiation from accelerated electrons is a promising mechanism. It is believed that electrons are accelerated continuously by ultra-relativistic shocks or magnetic reconnections. At the same time, these electrons will be cooled via several processes (mainly adiabatic expansion, synchrotron radiation and inverse Compton scattering (ICS)), which regulate the distribution of electrons. An extra high-energy spectrum component is expected to be induced by ICS. However, the gamma-gamma annihilation effect may eliminate the high-energy photons and prevent the observers from distinguishing the extra component. We perform numerical calculations by taking these effects into account and discuss whether the extra ICS component could be observed. By exploring the plausible parameter space for relevant quantities of the GRB ejecta, we present the electron distributions and the corresponding spectra of synchrotron radiation and ICS. It is found that the extra component is observable only for ejecta with a rather large bulk Lorentz factor. A large Lorentz factor means the adiabatic expansion is the leading process in the electron cooling procedure, which makes the low-energy electron distribution spectrum to be relatively hard. Therefore the ICS component is more likely to be detected for GRBs that have a hard low-energy photon spectrum.

Currently, FRB 121102 is the only fast radio burst source that was observed to give out bursts repeatedly. It shows a high repeating rate, with more than one hundred bursts being spotted, but with no obvious periodicity in the activities. Thanks to its repetition, the source was well localized with a subarcsecond accuracy, leading to a redshift measurement of about 0.2. FRB 121102 is a unique source that can help us understand the enigmatic nature of fast radio bursts. In this study, we analyze the characteristics of the waiting times between bursts from FRB 121102. It is found that there is a clear bimodal distribution for the waiting times. While most waiting times cluster at several hundred seconds, a small portion of the waiting times are strikingly in the range of 2--40 millisecond. More interestingly, it is found that the waiting time does not correlate with the burst intensity, either for the preceding burst or for the subsequent burst. It strongly indicates that the repeating bursts should be generated by some external mechanisms, but not internal mechanisms. As a result, the models involving collisions between small bodies and neutron stars could be competitive mechanisms for this interesting source.

A large-scale ordered magnetic field plays a very important role in the formation and acceleration of a gamma-ray burst (GRB) jet. During the GRB prompt phase, some dissipation processes may happen and disturb the magnetic field, making the field become random to some extent. Therefore, a mixed magnetic field consisting of an ordered component and a random component is plausible for the early afterglow phase. Here we investigate the polarization evolution and light curve of an afterglow under such a 3-dimensional mixed magnetic field. Three kinds of ordered component (i.e., aligned, toroidal and radial) are discussed. We find that the three cases are distinguishable through jet polarization evolution. The polarization angle for a 3D mixed magnetic field with an aligned ordered component can evolve gradually but only changes abruptly by $90^\circ$ in the toroidal and radial cases. Furthermore, during the reverse shock crossing time, the polarization degree (PD) can be non-zero for the toroidal case but roughly zero for the radial case. Since an aligned component in a jet corresponds to a magnetar central engine and a toroidal component corresponds to a black hole, GRB central engines are distinguishable through polarization observations even if the magnetic field is mixed in a jet. In addition, our polarization calculation can be applied to GRB 170817A associated with GW170817. Using the recently-observed PD upper limit $12\%$ of GRB 170817A at $t=244$ days, the magnetic field strength ratio of the ordered to random components in this afterglow is constrained to be $\lesssim0.9$.

We constrain cosmological parameters using only Hubble parameter data and quantify the impact of future Hubble parameter measurements on parameter estimation for the most typical dark energy models. We first constrain cosmological parameters using 52 current Hubble parameter data including the Hubble constant measurement from the Hubble Space Telescope. Then we simulate the baryon acoustic oscillation signals from WFIRST (Wide-Field Infrared Survey Telescope) covering the redshift range of $z\in [0.5, 2]$ and the redshift drift data from E-ELT (European Extremely Large Telescope) in the redshift range of $z\in [2,5]$. It is shown that solely using the current Hubble parameter data could give fairly good constraints on cosmological parameters. Compared to the current Hubble parameter data, with the WFIRST observation the $H(z)$ constraints on dark energy would be improved slightly, while with the E-ELT observation the $H(z)$ constraints on dark energy is enormously improved.

The Tianlai Pathfinder is designed to demonstrate the feasibility of using a wide field of view radio interferometers to map the density of neutral hydrogen in the Universe after the Epoch of Reionizaton. This approach, called 21~cm intensity-mapping, promises an inexpensive means for surveying the large-scale structure of the cosmos. The Tianlai Pathfinder presently consists of an array of three, 15~m $\times$ 40~m cylinder telescopes and an array of sixteen, 6~m diameter dish antennas located in a radio-quiet part of western China. The two types of arrays were chosen to determine the advantages and disadvantages of each approach. The primary goal of the Pathfinder is to make 3D maps by surveying neutral hydrogen over large areas of the sky %$20,000 {\rm deg}^2$ in two different redshift ranges: first at $1.03 > z > 0.78$ ($700 - 800$~MHz) and later at $0.21 > z > 0.12$ ($1170 - 1270$~MHz). The most significant challenge to $21$~cm intensity-mapping is the removal of strong foreground radiation that dwarfs the cosmological signal. It requires exquisite knowledge of the instrumental response, i.e. calibration. In this paper, we provide an overview of the status of the Pathfinder and discuss the details of some of the analysis that we have carried out to measure the beam function of both arrays. We compare electromagnetic simulations of the arrays to measurements, discuss measurements of the gain and phase stability of the instrument, and provide a brief overview of the data processing pipeline.

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.

X-ray flares, lasting for $\sim 100 - 1000$ s in the X-ray band, are often observed following gamma-ray bursts (GRBs). The physical origin of X-ray flares is still unknown merely with the temporal/spectral information. On the other hand, some polarimeters are expected to be launched within several years thanks to the increasing interest on astronomical X-ray polarimetry. Here, by assuming that X-ray flares are synchrotron radiation from relativistic spherical shells, we show that the linear polarization degree during the rising phase of an X-ray flare is much higher for the emitting region with toroidal magnetic fields than that with random magnetic fields. In the decay phase of the flare, the evolution of the polarization degree is determined by the curvature effect of the emitting shell, which is a natural feature of jet scenarios for flares. Therefore, the measurement of the polarization of X-ray flares would provide a useful tool to probe the configuration of magnetic fields in the emission region, and may even help to test the curvature effect. The information on the magnetic configuration can further help us to understand the properties of GRB jets.

Fast radio bursts (FRBs) are bright radio pulses from the sky with millisecond durations and Jansky-level flux densities. Their origins are still largely uncertain. Here we suggest a new model for FRBs. We argue that the collision of a white dwarf with a black hole can give birth to a transient accretion disk, from which powerful episodic magnetic blobs will be launched. The collision between two consecutive magnetic blobs can result in a catastrophic magnetic reconnection, which releases a large amount of free magnetic energy and forms a forward shock. The shock propagates through the cold magnetized plasma within the blob in the collision region, radiating through synchrotron maser mechanism, which is responsible for a non-repeating FRB signal. Our calculations show that the theoretical energetics, radiation frequency, duration timescale, and event rate can be well consistent with the observational characteristic of FRBs.

Fast radio bursts (FRBs) are transient radio sources at cosmological distances. No counterparts in other bands have been observed for { non-repeating FRBs}. Here we suggest the collapse of strange star crusts as a possible origin for FRBs. Strange stars, which are composed of almost equal numbers of u, d, and s quarks, may be encapsulated by a thin crust of normal hadronic matter. When a strange star accretes matter from its environment, the crust becomes heavier and heavier. It may finally collapse, leading to the release of a large amount of magnetic energy and plenty of electron/positron pairs on a very short timescale. Electron/positron pairs in the polar cap region of the strange star can be accelerated to relativistic velocities, streaming along the magnetic field lines to form a thin shell. FRBs are produced by coherent emission from these electrons when the shell is expanding. Basic characteristics of observed FRBs can be explained in our model.

The brightness of the multi-wavelength afterglow of GRB 170817A is increasing unexpectedly even $\sim$ 160 days after the associated gravitational burst. Here we suggest that the brightening can be caused by a late-time energy injection process. We use an empirical expression to mimic the evolution of the injection luminosity, which is consisted of a power-law rising phase and a power-law decreasing phase. It is found that the power-law indices of the two phases are $0.92$ and $-2.8$, respectively, with the peak time of the injection being $\sim 110$ days. The energy injection could be due to some kind of accretion, with the total accreted mass being $\sim 0.006 M_\odot$. However, normal fall-back accretion, which usually lasts for a much shorter period, cannot provide a natural explanation. Our best-fit decaying index of $-2.8$ is also at odds with the expected value of $-5/3$ for normal fall-back accretion. Noting that the expansion velocities of the kilonova components associated with GW170817 is $0.1-0.3\,c$, we argue that there should also be some ejecta with correspondingly lower velocities during the coalescence of the double neutron star system. They are bound by the gravitational well of the remnant central compact object and might be accreted at a timescale of about 100 days, providing a reasonable explanation for the energy injection. Detailed studies on the long-lasting brightening of GRB 170817A thus may provide useful information on the matter ejection during the merger process of binary neutron stars.

Recent follow-up observations of the binary neutron star (NS) merging event GW170817/SGRB 170817A reveal that its X-ray/optical/radio emissions are brightening continuously up to $\sim 100$ days post-merger. This late-time brightening is unexpected from the kilonova model or the off-axis top-hat jet model for gamma-ray burst afterglows. In this paper, by assuming that the merger remnant is a long-lived NS, we propose that the interaction between an electron-positron-pair ($e^+e^-$) wind from the central NS and the jet could produce a long-lived reverse shock, from which a new emission component would rise and can interpret current observations well. The magnetic-field-induced ellipticity of the NS is taken to be $4 \times 10^{-5}$ in our modeling, so that the braking of the NS is mainly through the gravitational wave (GW) radiation rather than the magnetic dipole radiation, and the emission luminosity at early times would not exceed the observational limits. In our scenario, since the peak time of the brightening is roughly equal to the spin-down time scale of the NS, the accurate peak time may help constrain the ellipticity of the remnant NS. We suggest that radio polarization observations of the brightening would help to distinguish our scenario from other scenarios. Future observations on a large sample of short gamma-ray burst afterglows or detections of GW signals from merger remnants would test our scenario.

We analytically investigate the Schwinger pair production in the de Sitter background by using {\em the uniform asymptotic approximation method}, and show that the equation of motion in general has two turning points, and the nature of these points could be single, double, real or complex, depending on the choice of the free parameters involved in the theory. Different natures of these points lead to different electric currents. In particular, when $\beta \equiv m^2/H^2-9/4$ is positive, both turning points are complex, and the electric current due to the Schwinger process is highly suppressed, where $m$ and $H$ denote, respectively, the mass of the particle and the Hubble parameter. For the turning points to be real, it is necessary to have $\beta < 0$, and the more negative of $\beta$, the easier to produce particles. In addition, when $\beta < 0$, we also study the particle production when the electric field $E$ is very weak. We find that the electric current in this case is proportional to $E^{1/2 - \sqrt{|\beta|}}$, which is strongly enhanced in the weak electric field limit when $m < \sqrt{2} H$.

The low-energy spectra of gamma-ray bursts' (GRBs) prompt emission are closely related to the energy distribution of electrons, which is further regulated by their cooling processes. We develop a numerical code to calculate the evolution of the electron distribution with given initial parameters, in which three cooling processes (i.e., adiabatic, synchrotron and inverse Compton cooling) and the effect of decaying magnetic field are coherently considered. A sequence of results are presented by exploring the plausible parameter space for both the fireball and the Poynting-flux-dominated regime. Different cooling patterns for the electrons can be identified and they are featured by a specific dominant cooling mechanism. Our results show that the hardening of the low-energy spectra can be attributed to the dominance of synchrotron self-Compton cooling within the internal shock model, or to decaying synchrotron cooling within the Poynting-flux-dominated jet scenario. These two mechanisms can be distinguished by observing the hard low-energy spectra of isolated short pulses in some GRBs. The dominance of adiabatic cooling can also lead to hard low-energy spectra when the ejecta is moving at an extreme relativistic speed. The information from the time-resolved low-energy spectra can help to probe the physical characteristics of the GRB ejecta via our numerical results.

When an emitting spherical shell with a constant Lorentz factor turns off emission abruptly at some radius, its high-latitude emission would obey the relation of $\hat{\alpha}$ (the temporal index) = $2 + \hat{\beta}$ (the spectral index). However, this relation is violated by the X-ray fares in some gamma-ray bursts (GRBs), whose $\hat{\alpha}$ is much more steeper. We show that the synchrotron radiation should be anisotropic when the angular distribution of accelerated electrons has preferable orientation, and this anisotropy would naturally lead to the steeper decay for the high-latitude emission if the the intrinsic emission is limb-brightened. We use our simple toy model to reproduce the temporal and the spectral evolution of the X-ray flares. We show that our model can well interpret the steep decay of the X-ray flares in the three GRBs selected as example. Recent simulations on particle acceleration may support the specific anisotropic distribution of the electrons adopted in our work. Reversely, confirmation of the anisotropy in the radiation would give meaningful clues to the details of electron acceleration in the emitting region.

A wide bandwidth, dual polarized, modified four-square antenna is presented as a feed antenna for radio astronomical measurements. A linear array of these antennas is used as a line-feed for cylindrical reflectors for Tianlai, a radio interferometer designed for 21~cm intensity mapping. Simulations of the feed antenna beam patterns and scattering parameters are compared to experimental results at multiple frequencies across the 650 - 1420 MHz range. Simulations of the beam patterns of the combined feed array/reflector are presented as well.

A 3.56-hour white dwarf (WD) - M dwarf (MD) close binary system, AR Scorpii, was recently reported to show pulsating emission in radio, IR, optical, and UV, with a 1.97-minute period, which suggests the existence of a WD with a rotation period of 1.95 minutes. We propose a model to explain the temporal and spectral characteristics of the system. The WD is a nearly perpendicular rotator, with both open field line beams sweeping the MD stellar wind periodically. A bow shock propagating into the stellar wind accelerates electrons in the wind. Synchrotron radiation of these shocked electrons can naturally account for the broad-band (from radio to X-rays) spectral energy distribution of the system.

The prompt emission of gamma-ray bursts (GRBs) is characterized by rapid variabilities, which may be a direct reflection of the unsteady central engine. We perform a series of axisymmetric 2.5-dimensional simulations to study the propagation of relativistic, hydrodynamic, intermittent jets through the envelope of a GRB progenitor star. A realistic rapidly rotating star is incorporated as the background of jet propagation, and the star is allowed to collapse due to the gravity of the central black hole. By modeling the intermittent jets with constant-luminosity pulses with equal on and off durations, we investigate how the half-period, $T$, affects the jet dynamics. For relatively small $T$ values (e.g. 0.2 s), the jet breakout time $t_{\rm bo}$ depends on the opening angle of the jet, with narrower jets more penetrating and reaching the surface at shorter times. For $T \leq 1$ s, the reverse shock crosses each pulse before the jet penetrates through the stellar envelope. As a result, after the breakout of the first group of pulses at $t_{\rm bo}$, several subsequent pulses vanish before penetrating the star, causing a quiescent gap. For larger half-periods ($T=2.0, 4.0$ s), all the pulses can successfully penetrate through the envelope, since each pulse can propagate through the star before the reverse shock crosses the shell. Our results may interpret the existence of a weak precursor in many long GRBs, given that the GRB central engine injects intermittent pulses with half-period $T \leq 1$ s. The observational data seem to be consistent with such a possibility.

Optical re-brightenings in the afterglows of some gamma-ray bursts (GRBs) are unexpected within the framework of the simple external shock model. While it has been suggested that the central engines of some GRBs are newly born magnetars, we aim to relate the behaviors of magnetars to the optical re-brightenings. A newly born magnetar will lose its rotational energy in the form of Poynting-flux, which may be converted into a wind of electron-positron pairs through some magnetic dissipation processes. As proposed by Dai (2004), this wind will catch up with the GRB outflow and a long-lasting reverse shock would form. By applying this scenario to GRB afterglows, we find that the reverse shock propagating back into the electron-positron wind can lead to an observable optical re-brightening and a simultaneous X-ray plateau (or X-ray shallow decay). In our study, we select four GRBs, i.e., GRB 080413B, GRB 090426, GRB 091029, and GRB 100814A, of which the optical afterglows are well observed and show clear re-brightenings. We find that they can be well interpreted. In our scenario, the spin-down timescale of the magnetar should be slightly smaller than the peak time of the re-brightening, which can provide a clue to the characteristics of the magnetar.

It was found by Amati et al. in 2002 that for a small sample of 9 gamma-ray bursts, more distant events appear to be systematically harder in the soft gamma-ray band. Here, we have collected a larger sample of 65 gamma-ray bursts, whose time integrated spectra are well established and can be well fitted with the so called Band function. It is confirmed that a correlation between the redshifts ($z$) and the low-energy indices ($\alpha$) of the Band function does exist, though it is a bit more scattered than the result of Amati et al. This correlation can not be simply attributed to the effect of photon reddening. Furthermore, correlations between $\alpha$ and $E_{\rm peak}$ (the peak energy in the $\nu F_{\nu}$ spectrum in the rest frame), $\alpha$ and $E_{\rm iso}$ (the isotropic energy release), $\alpha$ and $L_{\rm iso}$ (the isotropic luminosity) are also found, which indicate that these parameters are somehow connected. The results may give useful constraints on the physics of gamma-ray bursts.

The detection of optical re-brightenings and X-ray plateaus in the afterglows of gamma-ray bursts (GRBs) challenges the generic external shock model. Recently, we have developed a numerical method to calculate the dynamic of the system consisting of a forward shock and a reverse shock. Here, we briefly review the applications of this method in the afterglow theory. By relating these diverse features to the central engines of GRBs, we find that the steep optical re-brightenings would be caused by the fall-back accretion of black holes, while the shallow optical re-brightenings are the consequence of the injection of the electron-positron-pair wind from the central magnetar. These studies provide useful ways to probe the characteristics of GRB central engines.

The shallow decay phase and flares in the afterglows of gamma-ray bursts (GRBs) is widely believed to be associated with the later activation of central engine. Some models of energy injection involve with a continuous energy flow since the GRB trigger time, such as the magnetic dipole radiation from a magnetar. However, in the scenario involving with a black hole accretion system, the energy flow from the fall-back accretion may be delayed for a fall-back time $\sim t_{\rm fb}$. Thus we propose a delayed energy injection model, the delayed energy would cause a notable rise to the Lorentz factor of the external shock, which will "generate" a bump in the multiple band afterglows. If the delayed time is very short, our model degenerates to the previous models. Our model can well explain the significant re-brightening in the optical and infrared light curves of GRB 081029 and GRB 100621A. A considerable fall-back mass is needed to provide the later energy, this indicates GRBs accompanied with fall-back material may be associated with a low energy supernova so that fraction of the envelope can be survived during eruption. The fall-back time can give meaningful information of the properties of GRB progenitor stars.

As a new kind of radio transient sources detected at $\sim 1.4$ GHz, fast radio bursts are specially characterized by their short durations and high intensities. Although only ten events are detected so far, fast radio bursts may actually frequently happen at a rate of $\sim 10^{3}$ --- $10^4~\rm{sky}^{-1}~\rm{day}^{-1}$. We suggest that fast radio bursts can be produced by the collisions between neutron stars and asteroids. This model can naturally explain the millisecond duration of fast radio bursts. The energetics and event rate can also be safely accounted for. Fast radio bursts thus may be one side of the multifaces of the neutron star-small body collision events, which are previously expected to lead to X-ray/gamma-ray bursts or glitch/anti-glitches.

We explore the impact of the Sandage-Loeb (SL) test on the precision of cosmological constraints for $f(T)$ gravity theories. The SL test is an important supplement to current cosmological observations because it measures the redshift drift in the Lyman-$\alpha$ forest in the spectra of distant quasars, covering the "redshift desert" of $2 \lesssim z \lesssim5$. To avoid data inconsistency, we use the best-fit models based on current combined observational data as fiducial models to simulate 30 mock SL test data. We quantify the impact of these SL test data on parameter estimation for $f(T)$ gravity theories. Two typical $f(T)$ models are considered, the power-law model $f(T)_{PL}$ and the exponential-form model $f(T)_{EXP}$. The results show that the SL test can effectively break the existing strong degeneracy between the present-day matter density $\Omega_m$ and the Hubble constant $H_0$ in other cosmological observations. For the considered $f(T)$ models, a 30-year observation of the SL test can improve the constraint precision of $\Omega_m$ and $H_0$ enormously but cannot effectively improve the constraint precision of the model parameters.

By detecting redshift drift in the spectra of Lyman-$\alpha$ forest of distant quasars, Sandage-Loeb (SL) test directly measures the expansion of the universe, covering the "redshift desert" of $2 \lesssim z \lesssim5$. Thus this method is definitely an important supplement to the other geometric measurements and will play a crucial role in cosmological constraints. In this paper, we quantify the ability of SL test signal by a CODEX-like spectrograph for constraining interacting dark energy. Four typical interacting dark energy models are considered: (i) $Q=\gamma H\rho_c$, (ii) $Q=\gamma H\rho_{de}$, (iii) $Q=\gamma H_0\rho_c$, and (iv) $Q=\gamma H_0\rho_{de}$. The results show that for all the considered interacting dark energy models, relative to the current joint SN+BAO+CMB+$H_0$ observations, the constraints on $\Omega_m$ and $H_0$ would be improved by about 60\% and 30--40\%, while the constraints on $w$ and $\gamma$ would be slightly improved, with a 30-yr observation of SL test. We also explore the impact of SL test on future joint geometric observations. In this analysis, we take the model with $Q=\gamma H\rho_c$ as an example, and simulate future SN and BAO data based on the space-based project WFIRST. We find that in the future geometric constraints, the redshift drift observations would help break the geometric degeneracies in a meaningful way, thus the measurement precisions of $\Omega_m$, $H_0$, $w$, and $\gamma$ could be substantially improved using future probes.

Fast radio bursts (FRBs) are newly discovered radio transient sources. Their high dispersion measures indicate an extragalactic origin. But due to the lack of observational data in other wavelengths, their progenitors still remain unclear. Here we suggest the collisions between neutron stars and asteroids/comets as a promising mechanism for FRBs. During the impact process, a hot plasma fireball will form after the material of the small body penetrates into the neutron star surface. The ionized matter inside the fireball will then expand along the magnetic field lines. Coherent radiation from the thin shell at the top of the fireball will account for the observed FRBs. Our scenario can reasonably explain the main features of FRBs, such as their durations, luminosities, and the event rate. We argue that for a single neutron star, FRBs are not likely to happen repeatedly in a forseeable time span since such impacts are of low probability. We predict that faint remnant X-ray emissions should be associated with FRBs, but it may be too faint to be detected by detectors at work.

We suggest that the collision of a small solid body with a pulsar can lead to an observable glitch/anti-glitch. The glitch amplitude depends on the mass of the small body and the impact parameter as well. In the collision, a considerable amount of potential energy will be released either in the form of a short hard X-ray burst or as a relatively long-lasting soft X-ray afterglow. The connection between the glitch amplitude and the X-ray energetics can help to diagnose the nature of these timing anomalies.

In recent years, more and more gamma-ray bursts with late rebrightenings in multi-band afterglows unveil the late-time activities of the central engines. GRB 100814A is a special one among the well-sampled events, with complex temporal and spectral evolution. The single power-law shallow decay index of the optical light curve observed by GROND between 640 s and 10 ks is $\alpha_{\rm opt} = 0.57 \pm 0.02$, which apparently conflicts with the simple external shock model expectation. Especially, there is a remarkable rebrightening in the optical to near infrared bands at late time, challenging the external shock model with synchrotron emission coming from the interaction of the blast wave with the surrounding interstellar medium. In this paper, we invoke a magnetar with spin evolution to explain the complex multi-band afterglow emission of GRB 100814A. The initial shallow decay phase in optical bands and the plateau in X-ray can be explained as due to energy injection from a spin-down magnetar. At late time, with the falling of materials from the fall-back disk onto the central object of the burster, angular momentum of the accreted materials is transferred to the magnetar, which leads to a spin-up process. As a result, the magnetic dipole radiation luminosity will increase, resulting in the significant rebrightening of the optical afterglow. It is shown that the observed multi-band afterglow emission can be well reproduced by the model.

Strange quark matter (SQM) may be the true ground state of hadronic matter, indicating that the observed pulsars may actually be strange stars, but not neutron stars. According to this SQM hypothesis, the existence of a hydrostatically stable sequence of strange quark matter stars has been predicted, ranging from 1 --- 2 solar mass strange stars, to smaller strange dwarfs and even strange planets. While gravitational wave (GW) astronomy is expected to open a new window to the universe, it will shed light on the searching for SQM stars. Here we show that due to their extreme compactness, strange planets can spiral very close to their host strange stars, without being tidally disrupted. Like inspiraling neutron stars or black holes, these systems would serve as a new kind of sources for GW bursts, producing strong gravitational waves at the final stage. The events occurring in our local Universe can be detected by the upcoming gravitational wave detectors, such as Advanced LIGO and the Einstein Telescope. This effect provides a unique probe to SQM objects and is hopefully a powerful tool for testing the SQM hypothesis.

We study an extensive sample of 87 GRBs for which there are well sampled and simultaneous optical and X-ray light-curves. We extract the cleanest possible signal of the afterglow component, and compare the temporal behaviors of the X-ray light-curve, observed by Swift XRT, and optical data, observed by UVOT and ground-based telescopes for each individual burst. Overall we find 62\% GRBs that are consistent with the standard afterglow model. When more advanced modeling is invoked, up to 91\% of the bursts in our sample may be consistent with the external shock model. A large fraction of these bursts are consistent with occurring in a constant interstellar density medium (ISM) (61\%) while only 39\% of them occur in a wind-like medium. Only 9 cases have afterglow light-curves that exactly match the standard fireball model prediction, having a single power law decay in both energy bands which are observed during their entire duration. In particular, for the bursts with chromatic behavior additional model assumptions must be made over limited segments of the light-curves in order for these bursts to fully agree with the external shock model. Interestingly, for 54\% of the X-ray and 40\% of the optical band observations the end of the shallow decay ($t^{\sim-0.5}$) period coincides with the jet break ($t^{\sim-p}$) time, causing an abrupt change in decay slope. The fraction of the burst that consistent with the external shock model is independent of the observational epochs in the rest frame of GRBs. Moreover, no cases can be explained by the cooling frequency crossing the X-ray or optical band.

The modified holographic Ricci dark energy coupled to interacting relativistic and non-relativistic dark matter is considered in the nonflat Friedmann-Robertson-Walker universe. Through examining the deceleration parameter, one can find that the transition time of the Universe from decelerating to accelerating phase in the interacting holographic Ricci dark energy model is close to that in the $\Lambda$ cold dark matter model. The evolution of modified holographic Ricci dark energy's state parameter and the evolution of dark matter and dark energy's densities shows that the dark energy holds the dominant position from the near past to the future. By studying the statefinder diagnostic and the evolution of the total pressure, one can find that this model could explain the Universe's transition from the radiation to accelerating expansion stage through the dust stage. According to the $Om$ diagnostic, it is easy to find that when the interaction is weak and the proportion of relativistic dark matter in total dark matter is small, this model is phantom-like. Through our studying, we find the interaction and the relativistic dark matter's proportion all have great influence on the evolution of the Universe.

Previous studies have shown that for the Supernova Legacy Survey three-year (SNLS3) data there is strong evidence for the redshift-evolution of color-luminosity parameter $\beta$ of type Ia supernovae (SN Ia). In this paper, we explore the effects of varying $\beta$ on the cosmological constraints of holographic dark energy (HDE) model. In addition to the SNLS3 data, we also use Planck distance prior data of cosmic microwave background (CMB), as well as galaxy clustering (GC) data extracted from Sloan Digital Sky Survey (SDSS) data release 7 and Baryon Oscillation Spectroscopic Survey (BOSS). We find that, for the both cases of using SN data alone and using SN+CMB+GC data, involving an additional parameter of $\beta$ can reduce $\chi^2$ by $\sim$ 36; this shows that $\beta$ deviates from a constant at 6$\sigma$ confidence levels. Adopting SN+CMB+GC data, we find that compared to the constant $\beta$ case, varying $\beta$ yields a larger fractional matter density $\Omega_{m0}$ and a smaller reduced Hubble constant $h$; moreover, varying $\beta$ significantly increases the value of HDE model parameter $c$, leading to $c\approx 0.8$, consistent with the constraint results obtained before Planck. These results indicate that the evolution of $\beta$ should be taken into account seriously in the cosmological fits. In addition, we find that relative to the differences between the constant $\beta$ and varying $\beta(z)$ cases, the effects of different light-curve fitters on parameter estimation are very small.

The Sandage-Loeb (SL) test directly measures the expansion rate of the universe in the redshift range of $2\lesssim z\lesssim 5$ by detecting redshift drift in the spectra of Lyman-$\alpha$ forest of distant quasars. We discuss the impact of the future SL test data on parameter estimation for the $\Lambda$CDM, the $w$CDM, and the $w_0w_a$CDM models. To avoid the potential inconsistency with other observational data, we take the best-fitting dark energy model constrained by the current observations as the fiducial model to produce 30 mock SL test data. The SL test data provide an important supplement to the other dark energy probes, since they are extremely helpful in breaking the existing parameter degeneracies. We show that the strong degeneracy between $\Omega_m$ and $H_0$ in all the three dark energy models is well broken by the SL test. Compared to the current combined data of type Ia supernovae, baryon acoustic oscillation, cosmic microwave background, and Hubble constant, the 30-yr observation of SL test could improve the constraints on $\Omega_m$ and $H_0$ by more than 60\% for all the three models. But the SL test can only moderately improve the constraint on the equation of state of dark energy. We show that a 30-yr observation of SL test could help improve the constraint on constant $w$ by about 25\%, and improve the constraints on $w_0$ and $w_a$ by about 20\% and 15\%, respectively. We also quantify the constraining power of the SL test in the future high-precision joint geometric constraints on dark energy. The mock future supernova and baryon acoustic oscillation data are simulated based on the space-based project JDEM. We find that the 30-yr observation of SL test would help improve the measurement precision of $\Omega_m$, $H_0$, and $w_a$ by more than 70\%, 20\%, and 60\%, respectively, for the $w_0w_a$CDM model.

It has been found that, for the Supernova Legacy Survey three-year (SNLS3) data, there is strong evidence for the redshift-evolution of color-luminosity parameter $\beta$. In this paper, adopting the $w$-cold-dark-matter ($w$CDM) model and considering its interacting extensions (with three kinds of interaction between dark sectors), we explore the evolution of $\beta$ and its effects on parameter estimation. In addition to the SNLS3 data, we also take into account the Planck distance priors data of the cosmic microwave background (CMB), the galaxy clustering (GC) data extracted from SDSS DR7 and BOSS, as well as the direct measurement of Hubble constant from the Hubble Space Telescope (HST) observation. We find that, for all the interacting dark energy (IDE) models, adding a parameter of $\beta$ can reduce $\chi^2$ by $\sim$ 34, indicating that $\beta_1 = 0$ is ruled out at 5.8$\sigma$ confidence level (CL). Furthermore, it is found that varying $\beta$ can significantly change the fitting results of various cosmological parameters: for all the dark energy models considered in this paper, varying $\beta$ yields a larger $\Omega_{c0}$ and a larger $w$; on the other side, varying $\beta$ yields a smaller $h$ for the $w$CDM model, but has no impact on $h$ for the three IDE models. This implies that there is a degeneracy between $h$ and $\gamma$. Our work shows that the evolution of $\beta$ is insensitive to the interaction between dark sectors, and then highlights the importance of considering $\beta$'s evolution in the cosmology fits.

The detection of the B-mode polarization of the cosmic microwave background (CMB) by the BICEP2 experiment implies that the tensor-to-scalar ratio $r$ should be involved in the base standard cosmology. In this paper, we extend the $\Lambda$CDM+$r$+neutrino/dark radiation models by replacing the cosmological constant with the dynamical dark energy with constant $w$. Four neutrino plus dark energy models are considered, i.e., the $w$CDM+$r+\sum m_\nu$, $w$CDM+r + $N_{\rm eff}$, $w$CDM+r + $\sum m_\nu$ + $N_{\rm eff}$, and $w$CDM+r + $N_{\rm eff}$ + $m_{\nu,{\rm sterile}}^{\rm eff}$ models. The current observational data considered in this paper include the Planck temperature data, the WMAP 9-year polarization data, the baryon acoustic oscillation data, the Hubble constant direct measurement data, the Planck Sunyaev-Zeldovich cluster counts data, the Planck CMB lensing data, the cosmic shear data, and the BICEP2 polarization data. We test the data consistency in the four cosmological models, and then combine the consistent data sets to perform joint constraints on the models. We focus on the constraints on the parameters $w$, $\sum m_\nu$, $N_{\rm eff}$, and $m_{\nu,{\rm sterile}}^{\rm eff}$.

Re-brightening bumps are frequently observed in gamma-ray burst (GRB) afterglows. Many scenarios have been proposed to interpret the origin of these bumps, of which a blast wave encountering a density-jump in the circumburst environment has been questioned by recent works. We develop a set of differential equations to calculate the relativistic outflow encountering the density-jump by extending the work of Huang et al. (1999). This approach is a semi-analytic method and is very convenient. Our results show that late high-amplitude bumps can not be produced under common conditions, only short plateau may emerge even when the encounter occurs at early time ($< 10^4$ s). In general, our results disfavor the density-jump origin for those observed bumps, which is consistent with the conclusion drawn from full hydrodynamics studies. The bumps thus should be due to other scenarios.

The Sandage-Loeb (SL) test is a unique method to probe dark energy in the "redshift desert" of $2\lesssim z\lesssim 5$, and thus it provides an important supplement to the other dark energy probes. Therefore, it is of great importance to quantify how the future SL test data impact on the dark energy constraints. To avoid the potential inconsistency in data, we use the best-fitting model based on the other geometric measurements as the fiducial model to produce 30 mock SL test data. The 10-yr, 20-yr, and 30-yr observations of SL test are analyzed and compared in detail. We show that compared to the current combined data of type Ia supernovae, baryon acoustic oscillation, cosmic microwave background, and Hubble constant, the 30-yr observation of SL test could improve the constraint on $\Omega_m$ by about $80%$ and the constraint on $w$ by about $25%$. Furthermore, the SL test can also improve the measurement of the possible direct interaction between dark energy and dark matter. We show that the SL test 30-yr data could improve the constraint on $\gamma$ by about $30%$ and $10%$ for the $Q=\gamma H\rho_c$ and $Q=\gamma H\rho_{de}$ models, respectively.

GRB 120326A is an unusual gamma-ray burst (GRB) which has a quite long plateau and a very late rebrightening both in X-ray and optical bands. The similar behavior of the optical and X-ray light curves suggests that they maybe have a common origin. The long plateau starts from several hundred seconds and ends at tens of thousands seconds. The peak time of the late rebrightening is about 30000 s. We analyze the energy injection model by means of numerical and analytical solutions, considering both the wind environment and ISM environment for GRB afterglows. We especially study the influence of the injection starting time, ending time, stellar wind density (or density of the circumburst environment), and injection luminosity on the shape of the afterglow light curves, respectively. We find that the light curve is largely affected by the parameters in the wind model. There is a "bump" at the late time only in the wind model too. In the wind case, it is interesting that the longer the energy injected, the more obvious the rebrightening will be. We also find the peak time of bump is determined by the stellar wind density. We use the late continuous injection model to interpret the unusual afterglow of GRB 120326A. The model can well fit the observational data, however, we find that the time scale of the injection must be larger than ten thousands seconds. This implies that the time scale of the central engine activity must be more than ten thousands seconds. This can give useful constraints on the central engine of GRBs. We consider a new born millisecond pulsar with strong magnetic field as the central engine. On the other hand, our results suggest that the circumburst environment of GRB 120326A is very likely a stellar wind.

Glitches have been frequently observed in neutron stars. Previously these glitches unexceptionally manifest as sudden spin-ups that can be explained as due to impulsive transfer of angular momentum from the interior superfluid component to the outer solid crust. Alternatively, such spin-up glitches may also be due to large-scale crust-cracking events. However, an unprecedented anti-glitch was recently reported for the magnetar 1E 2259+586. In this case, the magnetar clearly exhibited a sudden spin-down, strongly challenging previous glitch theories. Here we show that the anti-glitch can be well explained by the collision of a small solid body with the magnetar. The intruder has a mass of about $1.1 \times 10^{21}$ g. Its orbital angular momentum is assumed to be antiparallel to that of the spinning magnetar, so that the sudden spin-down can be naturally accounted for. The observed hard X-ray burst and decaying softer X-ray emission associated with the anti-glitch can also be reasonably explained. Our study indicates that a completely different type of glitches as due to collisions between small bodies and neutron stars should exist and may have already been observed previously. It also hints a new way for studying the capture events by neutron stars: through accurate timing observations of pulsars.