It is generally believed that the variability of photospheric emission in gamma-ray bursts (GRBs) traces that of the jet power. This work further investigates the variability of photospheric emission in a variable jet. By setting a constant $\eta$ (dimensionless entropy of the jet), we find that the light curve of the photospheric emission shows a ``tracking'' pattern on the time profile of jet power. However, the relative variability is significantly low in the photospheric emission compared with that in the jet power. If the $\eta$ is genetic variable, the variability of the photospheric emission is not only limited by the jet power but also affected by $\eta$ strongly. It becomes complex and is generally different from that of the jet power. Moreover, the opposite phase may stand in the variabilities of the photospheric emission at different photon energies. We also find that the relative variability does not remain constant over the photon energies with an obvious reduction at a certain energy. This is consistent with the analysis of GRB 090902B in which an appreciable thermal component has been detected in a wide energy range. For several other GRBs coupling with the thermal component, we conservatively evaluate the variability of the thermal and non-thermal emission, respectively. Our results show that the relative variability of the thermal emission is likely comparable to that of the non-thermal emission for these bursts. In addition, the analysis of GRB~120323A reveals that the variability of the photospheric emission may be of the opposite phase from that of the non-thermal emission.

By assuming the formation of a black hole soon after the merger event of GW170817, Shibata et al. updated the constraints on the maximum mass ($M_\textrm{max}$) of a stable neutron star within $\lesssim$ 2.3 $M_{\odot}$, but there is no solid evidence to rule out $M_\textrm{max}>2.3~M_{\odot}$ from the point of both microphysical and astrophysical views. In order to explain massive pulsars, it is naturally expected that the equation of state (EOS) would become stiffer beyond a specific density. In this paper, we consider the possibility of EOSs with $M_\textrm{max}>2.3~M_{\odot}$, investigating the stiffness and the transition density in a polytropic model. Two kinds of neutron stars are considered, i.e., normal neutron stars (the density vanishes on gravity-bound surface) and strange stars (a sharp density discontinuity on self-bound surface). The polytropic model has only two parameter inputs in both cases: ($\rho_{\rm t}$, $\gamma$) for gravity-bound objects, while ($\rho_{\rm s}$, $\gamma$) for self-bound ones, with $\rho_{\rm t}$ the transition density, $\rho_{\rm s}$ the surface density and $\gamma$ the polytropic exponent. In the matter of $M_\textrm{max}>2.3~M_{\odot}$, it is found that the smallest $\rho_{\rm t}$ and $\gamma$ should be $\sim 0.50~\rho_0$ and $\sim 2.65$ for normal neutron stars, respectively, whereas for strange star, we have $\gamma > 1.40$ if $\rho_{\rm s} > 1.0~\rho_0$ and $\rho_{\rm s} < 1.58~\rho_0$ if $\gamma <2.0$ ($\rho_0$ is the nuclear saturation density). These parametric results could guide further research of the real EOS with any foundation of microphysics if a pulsar mass higher than $2.3~M_{\odot}$ is measured in the future. We also derive rough results of common neutron star radius range, which is $9.8~\rm{km} < R_{1.4} < 13.8~\rm{km}$ for normal neutron stars and $10.5~\rm{km} < R_{1.4} < 12.5~\rm{km}$ for strange stars.

The origin of gamma-ray burst (GRB) X-ray plateau, especially the internal plateau is still unclear, but could be related to GRB's central engine of magnetar. It is generally believed that the spin-down power of the magnetar is injected into forward external shock, however we propose here that most of the power will be dissipated behind the GRB jet through larger amplitude electromagnetic wave (LAEMW). The relevant physical conditions and observational implications are analyzed and discussed, and various kinds of X-ray light curves could be reproduced. Although it is still a matter of debate about the chromatic multi-band afterglow in the standard external afterglow fireball model, we can explain naturally this feature under this scenario. Furthermore, we predict that the X-ray emission of spin-down wind could possibly precede the prompt emission of GRB jet if the energy of LAEMW is dissipated first but shock-induced radiation in the jet is produced later. It is emphasized that both the GRB jet and the spin-down wind should have significant observational consequences in the magnetar scenario, and would be focused equally in GRB physics.

We propose a geometrical explanation for the periodically and nonperiodicly repeating fast radio bursts (FRBs) under neutron star (NS)-companion system. We suggest a critical binary separation, $r_{\rm c}$, within which the interaction between the NS and its companion can trigger FRBs. On one hand, for an elliptic orbit with the minimum and maximum binary separations, $r_{\rm min}$ and $r_{\rm max}$ (satisfying the condition $r_{\rm min}<r_{\rm c}<r_{\rm max}$), this binary system could reproduce periodical FRBs. A small separation $r$ ($<r_{\rm c}$) of such an elliptic orbit would result in an active period of periodical FRBs. On the other hand, if $r_{\rm max}<r_{\rm c}$, the modulation of orbital motion will not work due to persistent interaction, and this kind of repeating FRBs should be nonperiodic. It is estimated that the critical separation $r_{\rm c}\sim 10^{12}$ cm for stellar binary system.

The unknown equation of state (EoS) of neutron stars (NSs) is puzzling because of rich non-perturbative effects of strong interaction there. A method to constrain the EoS by using the detected X-ray plateaus of gamma-ray bursts (GRBs) is proposed in this paper. Observations show some GRB X-ray plateaus may be powered by strongly magnetized millisecond NSs. The properties of these NSs should then satisfy: (i) the spin-down luminosity of these NSs should be brighter than the observed luminosity of the X-ray plateaus; (ii) the total rotational energy of these NSs should be larger than the total energy of the X-ray plateaus. Through the case study of GRB 170714A, the moment of inertia of NSs is constrained as $I>1.0\times 10^{45}\left ( \frac{P_{\rm cri}}{1\;\rm ms} \right )^{2} \;\rm g\cdot cm^{2}$, where $P_{\rm cri}$ is the critical rotational period that an NS can achieve. The constraint of the radii of NSs according to GRB 080607 is shown in Table 1.

The fast radio burst, FRB 171019, was relatively bright when discovered first by ASKAP, but was identified as a repeater with three faint bursts detected later by GBT and CHIME. These observations lead to the discussion of whether the first bright burst shares the same mechanism with the following repeating bursts. A model of binary neutron star merger is proposed for FRB 171019, in which the first bright burst occurred during the merger event, while the subsequent repeating bursts are starquake-induced, and generally fainter, as the energy release rate for the starquakes can hardly exceed that of the catastrophic merger event. This scenario is consistent with the observation that no burst detected is as bright as the first one.

It is well known that collisionless systems are dissipation free from the perspective of particle collision and thus conserve entropy. On the other hand, processes such as magnetic reconnection and turbulence appear to convert large-scale magnetic energy into heat. In this paper, we investigate the energization and heating of collisionless plasma. The dissipation process is discussed in terms of fluid entropy in both isotropic and gyrotropic forms. Evolution equations for the entropy are derived and they reveal mechanisms that lead to changes in fluid entropy. These equations are verified by a collisionless particle-in-cell simulation of multiple reconnecting current sheets. In addition to previous findings regarding the pressure tensor, we emphasize the role of heat flux in the dissipation process.

Recent research shows that the spectral lag is closely related to the spectral evolution in GRBs. In this paper, we study the spectral lag for a radiating jet shell with a high-energy cut-off radiation spectrum. For the jet shell with a cut-off power-law spectrum, the spectral lag monotonically increases with the photon energy and levels off at a certain photon energy. It is the same for the jet shell with a Band cut-off spectrum (Bandcut). However, a turn-over from the positive lags to negative lags appears in the high-energy range for the jet shell with Bandcut, which is very similar to that observed in GRB~160625B. The dependence of the spectral lags on the spectral shape/evolution are studied in details. In addition, the spectral lag behavior observed in GRB~160625B is naturally reproduced based on our theoretical outcome.

The Lorentz factor of a relativistic jet and its evolution during the jet expansion are difficult to estimate, especially for the jets in gamma-ray bursts (GRBs). However, it is related to the understanding of jet physics. Owing to the absorption of two-photon pair production ($\gamma\gamma{\leftrightarrow}e^+e^-$), a high-energy spectral cutoff may appear in the radiation spectrum of GRBs. We search such kind of high-energy cutoff in GRB 160625B, which is one of the brightest bursts in recent years. It is found that the high-energy spectral cutoff is obvious for the first pulse in the second emission episode of GRB 160625B (i.e., $\sim$186-192 s after the burst first trigger), which is smooth and well-shaped. Then, we estimate the Lorentz factor and radiation location of the jet shell associated with the first pulse in the second emission episode of GRB 160625B. It is found that the radiation location increases with time. In addition, the Lorentz factor remains almost constant during the expansion of the jet shell. This reveals that the magnetization of the jet is low or intermediate in the emission region, event though the jet could be still Poynting flux dominated at smaller radii to avoid a bright thermal component in the emission episode.

Strict measurement of the speed of gravitational wave (GW) is very important for fundamental physics. In this paper, taking cosmological effect into account, we derive a more precise formula for calculating the speed of GW based on GW 170817-like events. We find that in the case of high redshift, the usual luminosity distance needs to be replaced by co-moving distance. We consider some possible electromagnetic signals which can be significantly reduce the intrinsic uncertainty between the emission times of the gravitational and electromagnetic signals. Hopefully, the measurement accuracy will be improved to $\sim 10^{-18}\;\rm m\; s^{-1}$.

Neutrino annihilation process on a hyperaccreting disk is one of the leading models to explain the generation of relativistic jets of gamma-ray bursts (GRBs). However, there is still uncertainty regarding the NAE of neutron star-accretion disc (NS-disc) system because of complicated microphysics processes and effects of strong magnetic field. In this paper, we investigate the NAE by assuming that the prompt jet of GRB 070110 is driven by neutrino pair annihilation in the NS-disc system. Our calculation shows $\eta_{\nu\bar{\nu}}> 1.2 \times 10^{-3}$ under the estimated accretion rate $\dot{M}\simeq 0.04\rm M_{\odot}\cdot s^{-1}$. Independent of the detailed accretion disc models, our result shows that the magnetic field may play an important role in the neutrino annihilation process on the hyperaccreting magnetized accretion disc. Compared with the theoretical value of $\eta_{\rm \nu \bar{\nu}}$ of the nonmagnetized BH-disc system, the NAE should increase significantly in the case of the NS-disc system if the GRB is powered by magnetar-disc system.

The measurement of the speed of gravitational waves (GWs) is useful to distinguish general relativity from massive gravity. We propose a new model-independent strategy to measure the speed of GWs with the distorted pulsars. Theoretically, when the the GW frequencies from a distorted pulsar are twice the frequencies of EM pulses, they should be emitted at the same time. By measuring the arrival times of these two signals emitted at the same time, the speed of GWs can be calculated with the time difference. Specifically, when the glitches of pulsars are taken into consideration, some pulsars in our Galaxy and nearby galaxies are potential to test our new strategy at high accuracy in the foreseeable future. On the other hand, the new method is meaningful as a motivation for future design of instruments.

We study the scalar stochastic gravitational-wave background (SGWB) from astrophysical sources, including compact binary mergers and stellar collapses, in the Bras-Dicke theory of gravity. By contrast to tensor waves, we found the scalar SGWB to predominantly arise from stellar collapses. These collapses not only take place at higher astrophysical rates, but emit more energy. This is because, unlike tensor radiation, which mainly starts from quadrupole order, the scalar perturbation can be excited by changes in the monopole moment. In particular, in the case of stellar collapse into a neutron star or a black hole, the monopole radiation, at frequencies below 100\,Hz, is dominated by the memory effect. At low frequencies, the scalar SGWB spectrum follows a power law of $\Omega_\text{S}\propto f^\alpha$, with $\alpha = 1$. We predict that $\Omega_\text{S}$ is inversely proportional to the square of $\omega_{\rm BD}+2$, with $\quad(\omega_{\rm BD}+2)^2\Omega_S(f=25\,{\rm Hz}) = 2.8\times 10^{-6}$. We also estimate the detectability of the scalar SGWB for current and third-generation detector networks, and the bound on $\omega_{\rm BD}$ that can be imposed from these observations.

Magnetic flux ropes are commonly observed throughout the heliosphere, and recent studies suggest that interacting flux ropes are associated with some energetic particle events. In this work, we carry out 2D particle-in-cell (PIC) simulations to study the coalescence of two magnetic flux ropes (or magnetic islands), and the subsequent plasma energization processes. The simulations are initialized with two magnetic islands embedded in a reconnecting current sheet. The two islands collide and eventually merge into a single island. Particles are accelerated during this process as the magnetic energy is released and converted to the plasma energy, including bulk kinetic energy increase by the ideal electric field, and thermal energy increase by the fluid compression and the non-ideal electric field. We find that contributions from these different energization mechanisms are all important and comparable with each other. Fluid shear and a non-gyrotropic pressure tensor also contribute to the energy conversion process. For simulations with different box sizes ranging from $L_x \sim 25$--$100 d_i$ and ion-to-electron mass ratios $m_i / m_e = 25$, 100 and 400, we find that the general evolution is qualitatively the same for all runs, and the energization depends only weakly on either the system size or the mass ratio. The results may help us understand plasma energization in solar and heliospheric environments.

It has been speculated that quantum gravity corrections may lead to modifications to space-time geometry near black hole horizons. Such structures may cause reflections to gravitational waves, causing {\it echoes} that follow the main gravitational waves from binary black hole coalescence. We show that such echoes, if exist, will give rise to a stochastic gravitational-wave background, which is very substantial if the near-horizon structure has a near unity reflectivity for gravitational waves, readily detectable by Advanced LIGO. In case reflectivity is much less than unity, the background will mainly be arising from the first echo, with a level proportional to the power reflectivity of the near-horizon structure, but robust against uncertainties in the location of the structure --- as long as it is very close to the horizon. Sensitivity of third-generation detectors allows the detection of a background that corresponds to power reflectivity $\sim 10^{-3}$, if the uncertainties in the binary black-hole merger rate can be removed. We note that the echoes do alter the $f^{2/3}$ power law of the background spectra at low frequencies, which is rather robust against the uncertainties.

Two-episode emission components separated by quiescent gaps in the prompt emission of gamma-ray bursts (GRBs) have been observed in the {\em Swift} era, but there is a lack of spectral information due to the narrow energy band of the {\em Swift}/Burst Alert Telescope. In this paper, a systematic analysis of the spectral and temporal properties of the prompt emission of 101 {\em Fermi}/Gamma-ray Burst Monitor detected long GRBs show the existence of two-episode emission components in the light curves, with quiescent times of up to hundreds of seconds. We focus on investigating the differences of those two emission episodes. We find that the light curves of the two emission components exhibit different behavior, e.g., a soft emission component that either precedes or follows the main prompt emission or that the intensity of the two emission episodes are comparable with each other. No statistically significant correlation in the duration of the two emission episodes can be claimed. We define a new parameter $\varepsilon$ as the ratio of the peak flux of the first and second emission episodes and find that a higher $\varepsilon$ corresponds to a larger fluence. The preferred spectral model in our analysis is a cutoff power-law model for most GRBs. The distribution of $E_p$ for episodes I and II range from tens of keV to 1000 keV with a lognormal fit and there are no significant differences between them. Moreover, we do not find significant relationships between $\varepsilon$ and $E_p$ for the two emission episodes. Those results suggest that these two-episode emission components likely share the same physical origin.

The favored progenitor model for Gamma-ray Bursts (GRBs) with Supernova (SN) association is the core collapse of massive stars. One possible outcome of such a collapse is a rapidly spinning, strongly magnetized neutron star ("magnetar"). We systematically analyze the multi-wavelength data of GRB/SN associations detected by several instruments before 2017 June. Twenty GRB/SN systems have been confirmed via direct spectroscopic evidence or a clear light curve bump, as well as some spectroscopic evidence resembling a GRB-SN. We derive/collect the basic physical parameters of the GRBs and the SNe, and look for correlations among these parameters. We find that the peak brightness, $\rm ^{56}Ni$ mass, and explosion energy of SNe associated with GRBs are statistically higher than other Type Ib/c SNe. A statistically significant relation between the peak energy of GRBs and the peak brightness of their associated SNe is confirmed. No significant correlations are found between the GRB energies (either isotropic or beaming-corrected) and the supernova energy. We investigate the energy partition within these systems and find that the beaming-corrected GRB energy of most systems is smaller than the SN energy, with less than 30\% of the total energy distributed in the relativistic jet. The total energy of the systems is typically smaller than the maximum available energy of a millisecond magnetar ($2\times 10^{52}$ erg), especially if aspherical SN explosions are considered. The data are consistent with-though not a proof of-the hypothesis that most, but not all, GRB/SN systems are powered by millisecond magnetars.

We study the gravitational wave (GW) production induced by the asymmetric jets of gamma-ray bursts (GRBs). The asymmetric jets result in a recoil force acted on the central compact object, whose motion leads to emission of GW. Under reasonable assumptions and simplifications, we derive the analytic form of the produce GWs. The amplitude of emitted GWs is estimated to be relatively low, but possibility exists that they can be detected by future experiments such as the Einstein Telescope. We find the dynamical properties of the central object, which is difficult to be studied via the electromagnetic (EW) channel, can be inferred by measuring the emitted GWs. Moreover, we find the emitted GWs can be used determine whether the relativistic jets is launched by the neutrino annihilation process or the Blandford-Znajek process, which cannot be clearly distinguished by the current GRB observations. Our work manifests the importance of the GW channel in multi-messenger astronomy. The physical information encoded in the GW and EW emissions of an astrophysical object is complementary to each other; in case some physics can not be effectively investigated using the EW channel alone, including the GW channel can be very helpful.

We compile the radio-optical-X-ray spectral energy distributions (SEDs) of 65 knots and 29 hotspots in 41 active galactic nucleus jets to examine their high energy radiation mechanisms. Their SEDs can be fitted with the single-zone leptonic models, except for the hotspot of Pictor A and six knots of 3C 273. The X-ray emission of one hotspot and 22 knots is well explained as synchrotron radiations under the equipartition condition; they usually have lower X-ray and radio luminosities than the others, which may be due to a lower beaming factor. An inverse Compton (IC) process is involved for explaining the X-ray emission of the other SEDs. Without considering the equipartition condition, their X-ray emission can be attributed to the synchrotron-self-Compton (SSC) process, but the derived jet power (P_jet) are not correlated with L_k and most of them are larger than L_k with more than three orders of magnitude, where L_k is the jet kinetic power estimated with their radio emission. Under the equipartition condition, the X-ray emission is well interpreted with the IC process to the cosmic microwave background photons (IC/CMB). In this scenario, the derived P_jet of knots and hotspots are correlated with and comparable to L_k. These results suggest that the IC/CMB model may be the promising interpretation of their X-ray emission. In addition, a tentative knot-hotspot sequence in the synchrotron peak-energy--peak-luminosity plane is observed, similar to the blazar sequence, which may be attributed to their different cooling mechanisms of electrons.

Motivating by the discovery of association between GW 170817 and sGRB 170817A, we present a comprehensive analysis for sGRBs observed with Fermi/GBM in 9 operation years and study the properties of sGRB 170817A -like events. We derive a catalog of 275 typical sGRBs and 48 sGRB 170817A-like weak events from the GBM data of 2217 GRBs. We visibly identify two patterns of their light curve, single episode (Pattern I, 61\% of the SGRBs) and multiple episodes (Pattern II, 39\% of the SGRBs). Their duration distribution shows a tentative bimodal feature. Their spectra can be fitted with a cutoff power-law model, except for 4 sGRBs, and the spectral indices normally distribute at $\Gamma=0.69\pm 0.40$. Their $E_p$ values show a tentative bimodal distribution with peaks at 145 keV and 500 keV. No correlation among $T_{90}$, $E_p$, and $\Gamma$ is found. GRB 170817A is a soft, weak sGRB with $ E_{p}=124\pm 106$ keV, $L_{\rm iso}=(5.67\pm4.65)\times10^{46}\rm ~erg~s^{-1}$, and $E_{\rm iso}=(3.23\pm2.65)\times10^{46}\rm ~erg$. It follows the $E_{\rm iso}-E_{\rm p}$ relation of typical short GRBs. Its lightcurve is of Pattern II. Two lightcurve patterns, together with the potential two components in the $E_{\rm p}$ and $T_{90}$ distributions, we suspect that the current sample may include two distinct types of sGRBs from different progenitors. sGRB 170817A-like events may be from NS-NS mergers and those sGRBs with a Pattern I lightcurve may be from another distinct type of compact binary.

A rapidly spinning, strongly magnetized neutron star is invoked as the central engine for some Gamma-ray bursts (GRBs), especially, the $"$internal plateau$"$ feature of X-ray afterglow. However, for these $"$internal plateau$"$ GRBs, how to produce their prompt emission remains an open question. Two different physical process have been proposed in the literature, (1) a new-born neutron star is surrounded by a hyper-accreting and neutrino cooling disk, the GRB jet can be powered by neutrino annihilation aligning the spin axis; (2) a differentially rotating millisecond pulsar was formed due to different angular velocity between the interior core and outer shell parts of the neutron star, which can power an episodic GRB jet. In this paper, by analyzing the data of one peculiar GRB 070110 (with internal plateau), we try to test which model being favored. By deriving the physical parameters of magnetar with observational data, the parameter regime for initial period ($P_{0\rm }$) and surface polar cap magnetic field ($B_{\rm p}$) of the central NS are $(0.96\sim 1.2 )~\rm ms$ and $(2.4\sim 3.7)\times 10^{14}~\rm G$, respectively. The radiative efficiency of prompt emission is about $\eta_{\gamma} \sim 6\%$. However, the radiative efficiency of internal plateau ($\eta_{\rm X}$) is larger than $31\%$ assuming the $M_{\rm NS}\sim1.4 M_{\odot}$ and $P_{0\rm }\sim1.2 ~\rm ms$. The clear difference between the radiation efficiencies of prompt emission and internal plateau implies that they maybe originated from different components (e.g. prompt emission from the relativistic jet powered by neutrino annihilation, while the internal plateau from the magnetic outflow wind).