The interstellar oxygen isotopic ratio of 18O/17O can reflect the relative amount of the secular enrichment by ejecta from high-mass versus intermediate-mass stars. Previous observations found a Galactic gradient of 18O/17O, i.e., low ratios in the Galactic center and large values in the Galactic disk, which supports the insideout formation scenario of our Galaxy. However, the observed objects are not many and, in particular, not so many at large galactocentric distances. Thus we started a systematic study on Galactic interstellar 18O/17O, through observations of C18O and C17O multi-transition lines toward a large sample of 286 sources (at least one order of magnitude larger than previous ones), from the Galactic center region to the far outer Galaxy (~22 kpc). Here we present our observations of J=1-0 lines of C18O and C17O, with the ARO12m and the IRAM 30m telescope. We detected successfully both C18O and C17O 1-0 lines for 34 sources among our IRAM30m sample of 50 targets and for 166 sources among our ARO12m sample of 260 targets. The C18O optical depth effect on our ratio results, evaluated by fitting results of C17O spectra with hyperfine components and our RADEX non-LTE model calculation for the strongest source, was found to be insignificant. Beam dilution does not seem to be a problem either, which was supported by the fact of no systematic variation between the isotopic ratio and the heliocentric distance, and consistent measured ratios from two telescopes for most of those detected sources. Our results, though there are still very few detections made for sources in the outer Galaxy, confirm the apparent 18O/17O gradient of 18O/17O = (0.10+-0.03)R_GC+(2.95+-0.30), with a Pearson's rank correlation coefficient R = 0.69. This is supported by the newest Galactic chemical evolution model including the impact of massive stellar rotators and novae.

Distance measurement is crucial to astronomy. Here we suggest a new conceptual method to measure the distance by using a local instrument. By engaging the double-slit interference and by considering the phase information of the light, the position of the intensity maximum is related to the distance of the source. Consequently, the precise measurement of the position can be used to measure the distance of the remote source.

The morphology of bipolar planetary nebulae (PNe) can be attributed to interactions between a fast wind from the central engine and dense toroidal shaped ejecta left over from common envelope (CE) evolution. Here we use the 3-D hydrodynamic AMR code AstroBEAR to study the possibility that bipolar PN outflows can emerge collimated even from an uncollimated spherical wind in the aftermath of a CE event. The output of a single CE simulation via the SPH code PHANTOM serves as the initial conditions. Four cases of winds, all with high enough momenta to account for observed high momenta preplanetary nebula outflows, are injected spherically from the region of the CE binary remnant into the ejecta. We compare cases with two different momenta and cases with no radiative cooling versus application of optically thin emission via a cooling curve to the outflow. Our simulations show that in all cases highly collimated bipolar outflows result from deflection of the spherical wind via the interaction with the CE ejecta. Significant asymmetries between the top and bottom lobes are seen in all cases. The asymmetry is strongest for the lower momentum case with radiative cooling. While real post CE winds may be aspherical, our models show that collimation via "inertial confinement" will be strong enough to create jet-like outflows even beginning with maximally uncollimated drivers. Our simulations reveal detailed shock structures in the shock focused inertial confinement (SFIC) model and develop a lens-shaped inner shock that is a new feature of SFIC driven bipolar lobes.

The recent discovery of a fast radio burst (FRB 200428) from the Galactic magnetar SGR 1935+2154 robustly indicated that FRB phenomena can sometimes be produced by magnetars, although it is uncertain whether the cosmological FRBs can share the same origin with this Galactic event. The association of FRB 200428 with an X-ray burst (XRB) further offers important implications for the physical processes responsible for the FRB phenomena. By assuming that the XRB emission is produced in the magnetosphere of the magnetar, we investigate the possibility of that the FRB emission is produced by the synchrotron maser (SM) mechanism, which is powered by a shock due to the collision of an $e^{\pm}$ ejecta with a baryonic cloud. It is found that this shock-powered SM model can in principle account for the FRB 200428 observations, if the collision just occurred on the line of sight and the ejecta lunched by magnetar bursts can have appropriate ingredients and structures. To be specific, a burst ejecta should consist of an ultra-relativistic and extremely highly collimated $e^{\pm}$ component and a sub-relativistic and wide-spreading baryonic component. The cloud blocking the $e^{\pm}$ ejecta is just a remnant of a previous baryonic ejecta. Meanwhile, as a result of the synchrotron emission of the shocked material, an intense millisecond X-ray pulse is predicted to overlap the magnetosphere XRB emission, which in principle provides a way to test the model. Additionally, the peak frequency of the SM radiation is constrained to be about a few hundred MHz and the radiation efficiency is around $10^{-4}$.

FRB 180916.J0158+65 has been found to repeatedly emit fast radio bursts with the period in roughly 16 days. We propose that such periodicity comes from orbit-induced, spin precession of the emitter, which is possibly a neutron star. Depending on the mass of the companion, the binary period ranges from several hundred to thousands of seconds. Such tight binaries have a relatively short lifetime, which does not likely come from the gravitational decay of a wide binary. We comment on the relation to GW190425 and the possibility in LISA and LIGO detections.

The reverse shock (RS) model is generally introduced to interpret the optical afterglows with the rapid rising and decaying, such as the early optical afterglow of GRB 990123 (which is also called optical flash). In this paper, we collected 11 gamma-ray bursts (GRBs) early optical afterglows, which have such signatures of dominant RS emission. Since the temporal slopes of the optical flashes are determined by both the medium density distribution index $k$ and the electron spectral index $p$, we apply the RS model of the thin shell case to the optical flashes and determine the ambient medium of the progenitors. We find that the $k$ value is in the range of 0 - 1.5. The $k$ value in this paper is consistent with the result in Yi et al. (2013), where the forward shock (FS) model was applied to some onset bumps. However, the method adopted in this paper is only applicable to GRB afterglows with significant sharp rising and decaying RS emission. Our results indicate that the RS model can also be applied to confirm the circumburst medium, further implying that GRBs may have diverse circumburst media.

The nuclide $^{14}$C can be produced in the atmosphere by high energy particles and $\gamma$-rays from high-energy phenomena. Through the carbon cycle, some of $^{14}$CO$_2$ produced in the atmosphere can be retained in annual tree rings. Four events of rapid increase of the $^{14}$C content occurred in AD 775, AD 994, BC 660 and BC 3371 were found. Recently, the data of Jull et al. (2020) was inconsistent with our records around BC 3371. We measured our sample again and found the $^{14}$C records are consistent with the value in Wang et al. (2017). Therefore, our $^{14}$C records are robust. The inconsistency may be caused by the difference of calendar ages for the wood samples, or the physical origin of the event. First, crossdating on ring width can be performed only between trees whose growth has the same environmental conditions. Because the master tree-ring for dendrochronology is lack for Chinese trees. The master tree-ring from California has to be used. Therefore, the calendar ages derived from dendrochronology may be not precise. Second, the $^{14}$C even may be not global. One evidence is the variation of $^{14}$C content around AD 1006. The $^{14}$C contents of Californian trees increase 12\textperthousand~ in two years, while Japanese trees show no $^{14}$C increase.

In this paper, we present a novel method to test the Einstein's Equivalence Principle (EEP) using (simultaneous) multi-wavelength radio observations of polarized blazars. We analyze simultaneous multi-wavelength polarization observations of 3C 279 at 22, 43, and 86 GHz obtained by two antennas of the Korean VLBI Network. We obtained 15 groups of polarization data, and applied the Metropolis-Hastings Markov Chain (MHMC) to simulate the parameters when considering the EEP effect and the simplest form of Faraday rotation (single external Faraday screen). The final results show the constraint of the parameterized post-Newtonian (PPN) parameter $\gamma$ discrepancy as $\Delta \gamma_{p} = (1.91\pm0.34)\times10^{-20}$. However, the single external Faraday screen is an oversimplification for blazars because there are numerous observations show complex Faraday rotation behavior for blazars due to internal/external Faraday dispersion, beam depolarization, etc. The value $\Delta \gamma_{p}$ results of this paper can only be considered as upper limits. Only if all other effects are revealed and considered, should the result be taken as a direct measurement of the violation of the EEP.

The extreme mass ratio inspiral (EMRI), defined as a stellar-mass compact object inspiraling into a supermassive black hole (SMBH), has been widely argued to be a low-frequency gravitational wave (GW) source. EMRIs providing accurate measurements of black hole mass and spin, are one of the primary interests for Laser Interferometer Space Antenna (LISA). However, it is usually believed that there are no electromagnetic (EM) counterparts to EMRIs. Here we show a new formation channel of EMRIs with tidal disruption flares as EM counterparts. In this scenario, flares can be produced from the tidal stripping of the helium (He) envelope of a massive star by an SMBH. The left compact core of the massive star will evolve into an EMRI. We find that, under certain initial eccentricity and semimajor axis, the GW frequency of the inspiral can enter LISA band within 10 $\sim$ 20 years, which makes the tidal disruption flare an EM precursor to EMRI. Although the event rate is just $2\times 10^{-4}~\rm Gpc^{-3}yr^{-1}$, this association can not only improve the localization accuracy of LISA and help to find the host galaxy of EMRI, but also serve as a new GW standard siren for cosmology.

We compute the forces, torque and rate of work on the companion-core binary due to drag in global simulations of common envelope (CE) evolution for three different companion masses. Our simulations help to delineate regimes when conventional analytic drag force approximations are applicable. During and just prior to the first periastron passage of the in-spiral phase, the drag force is reasonably approximated by conventional analytic theory and peaks at values proportional to the companion mass. Good agreement between global and local 3D "wind tunnel" simulations, including similar net drag force and flow pattern, is obtained for comparable regions of parameter space. However, subsequent to the first periastron passage, the drag force is up to an order of magnitude smaller than theoretical predictions, quasi-steady, and depends only weakly on companion mass. The discrepancy is exacerbated for larger companion mass and when the inter-particle separation reduces to the Bondi-Hoyle-Lyttleton accretion radius, creating a turbulent thermalized region. Greater flow symmetry during this phase leads to near balance of opposing gravitational forces in front of and behind the companion, hence a small net drag. The reduced drag force at late times helps explain why companion-core separations necessary for envelope ejection are not reached by the end of limited duration CE simulations.

Tremaine and Weinberg (TW) proposed a conceptually simple procedure relying on long-slit spectroscopy to measure the pattern speeds of bars ($\Omega_{\rm p}$) in disk galaxies. Using a simulated galaxy, we investigate the potential biases and uncertainties of TW measurements using increasingly popular integral-field spectrographs (IFSs), for which multiple pseudo-slits (and thus independent measurements) can be constructed with a single observation. Most important, to establish the spatial coverage required and ensure the validity of the measurements, the inferred $\Omega_{\rm p}$ must asymptotically converge as the (half-)length of each pseudo-slit used is increased. The requirement for our simulation is to reach $\approx1.3$ times the half-light radius, but this may vary from galaxy to galaxy. Only those slits located within the bar region yield accurate measurements. We confirm that the position angle of the disk is the dominant source of systematic error in TW $\Omega_{\rm p}$ measurements, leading to under/overestimates of tens of percents for inaccuracies of even a few degrees. Recasting the data so that the data grid aligns with the disk major axis leads to slightly reduced uncertainties. Accurate measurements are obtained only for well-defined ranges of the bar angle (relative to the galaxy major axis) $\phi_{\rm bar}$ and the inclination angle $i$, here $10\lesssim\phi_{\rm bar}\lesssim75^{\circ}$ and $105\lesssim\phi_{\rm bar}\lesssim170^{\circ}$ and $15\lesssim i\lesssim70^{\circ}$. The adopted (pseudo-)slit widths, spatial resolution, and (unless extremely aggressive) spatial binning of IFS data have no significant impact on the measurements. Our results thus provide useful guidelines for reliable and accurate direct $\Omega_{\rm p}$ measurements with IFS observations.

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.

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 spreading of the X-line out of the reconnection plane under a strong guide field is investigated using large-scale three-dimensional (3D) particle-in-cell (PIC) simulations in asymmetric magnetic reconnection. A simulation with a thick, ion-scale equilibrium current sheet (CS) reveals that the X-line spreads at the ambient ion/electron drift speeds, significantly slower than the Alfv\'en speed based on the guide field $V_{Ag}$. Additional simulations with a thinner, sub-ion-scale CS show that the X-line spreads at $V_{Ag}$ (Alfv\'enic spreading), much higher than the ambient species drifts. An Alfv\'enic signal consistent with kinetic Alfv\'en waves develops and propagates, leading to CS thinning and extending, which then results in reconnection onset. The continuous onset of reconnection in the signal propagation direction manifests as Alfv\'enic X-line spreading. The strong dependence on the CS thickness of the spreading speeds, and the X-line orientation are consistent with the collisionless tearing instability. Our simulations indicate that when the collisionless tearing growth is sufficiently strong in a thinner CS such that $\gamma/\Omega_{ci}\gtrsim\mathcal{O}(1)$, Alfv\'enic X-line spreading can take place. Our results compare favorably with a number of numerical simulations and recent magnetopause observations. A key implications is that the magnetopause CS is typically too thick for Alfv\'enic X-line spreading to effectively take place.

We collected the information of 304 gamma-ray bursts (GRBs) from the literature, and analyzed the correlations among the host galaxy offsets (the distance from the site of the GRB to the center of its host galaxy), $T_{\rm 90,i}$ (the duration $T_{\rm 90}$ in rest-frame), $T_{\rm R45,i}$ (the duration $T_{\rm R45}$ in rest-frame), $E_{\rm \gamma,iso}$ (the isotropic equivalent energy), $L_{\rm \gamma,iso}$ ($=E_{\rm \gamma,iso}/T_{\rm 90,i}$, the isotropic equivalent luminosity) and $L_{\rm pk}$ (peak luminosity). We found that $ T_{\rm 90,i}$, $ T_{\rm R45,i}$, $ E_{\rm \gamma,iso}$, $L_{\rm pk}$ have negative correlation with $ {\rm offset}$, which is consistent with origin of short GRBs (SGRBs) and long GRBs (LGRBs). On separate analysis, we found similar results for $\log E_{\rm \gamma,iso}$ - $\log {\rm (offset)}$ and $\log L_{\rm pk}$ - $\log {\rm (offset)}$ relations in case of SGRBs only, while no obvious relation for LGRBs. There is no correlations between offset and $L_{\rm \gamma,iso}$. We also put the special GRB 170817A {and GRB 060218A} on the plots. {The two GRBs both have low luminosity and small offset.} In the $ \log ({\rm offset})- \log T_{\rm 90,i}$ plot, we found GRB 170817A locates in between the two regions of SGRBs and LGRBs and it is the outlier in the $ {\rm offset}-E_{\rm \gamma, iso}$, $ {\rm offset}-L_{\rm \gamma, iso}$ and $ {\rm offset}-L_{\rm pk}$ plots. Together with GRB 060218A being an outlier in all plots, it indicates the speciality of GRBs 170817A and 060218A, and might imply more subgroups of the GRB samples.

In order to obtain an overview of the gamma-ray bursts (GRBs), we need a full sample. In this paper, we collected 6289 GRBs (from GRB 910421 to GRB 160509A) from the literature, including prompt emission, afterglow and host galaxy properties. We hope to use this large sample to reveal the intrinsic properties of GRB. We have listed all the data in machine readable tables, including the properties of the GRBs, correlation coefficients and linear regression results of two arbitrary parameters, and linear regression results of any three parameters. These machine readable tables could be used as a data reservoir for further studies on the classifications or correlations. One may find some intrinsic properties from these statistical results. With this comprehensive table, it is possible to find relations between different parameters, and to classify the GRBs into different kinds of sub-groups. With the completion, it may reveal the nature of GRBs and may be used as tools like pseudo-redshift indicators, standard candles, etc. All the machine readable data and statistical results are available on the website of the journal.

Peculiar velocities are a precious tool to study the large-scale distribution of matter in the local universe and test cosmological models. However, present measurements of peculiar velocities are based on empirical distance indicators, which introduce large error bars. Here we present a new method to measure the peculiar velocities, by directly estimating luminosity distances through waveform signals from inspiralling compact binaries and measuring redshifts from electromagnetic (EM) counterparts. In the future, with the distance uncertainty of GW events reducing to $0.1$ per cent by future GW detectors, the uncertainty of the peculiar velocity can be reduced to $10$ km/s at 100 mega parsecs. We find that dozens of GW events with EM counterparts can provide a Hubble constant $H_0$ uncertainty of $0.5\%$ and the growth rate of structure with a $0.6\%$ precision in the third-generation ground-base GW detectors, which can reconcile the $H_0$ tension and determine the origins for cosmic accelerated expansion.

Well-sampled optical light curves of 50 gamma-ray bursts (GRBs) with plateau features are compiled from the literature. By empirical fitting, we obtained the parameters of the optical plateaus, such as the decay slopes ($\alpha_{\rm 1}$ and $\alpha_{\rm 2}$), the break times ($T_{\rm b}$), and the corresponding optical fluxes ($F_{\rm b}$) at the break times. The break time of optical plateaus ranges from tens of seconds to $10^6$ seconds, with a typical value about $10^4$ seconds. We have calculated the break luminosity, and it mainly ranges from $10^{44}$ erg $s^{-1}$ to $10^{47}$ erg $s^{-1}$, which is generally two or three orders of magnitude less than the corresponding break luminosity of the X-ray afterglow plateaus. We reanalyzed the optical plateaus and also found that a significantly tighter correlation exists when we added the isotropic equivalent energy of GRBs $E_{\rm \gamma,iso}$ into the $L_{\rm b,z}-T_{\rm b,z}$ relation. The best fit correlation is obtained to be $L_{\rm b,z}\propto T_{\rm b,z}^{-0.9}E_{\rm \gamma,iso}^{0.4}$. We next explored the possible correlations among $L_{\rm b,z}$, $T_{\rm b,z}$ and $E_{\rm p,i}$, and found there is also a tight correlation between them, which takes the form of $L_{\rm b,z}\propto T_{\rm b,z}^{-0.9}E_{\rm p,i}^{0.5}$. We argue that these two tight $L_{\rm b,z}-T_{\rm b,z}-E_{\rm \gamma,iso}$ and $L_{\rm b,z}-T_{\rm b,z}-E_{\rm p,i}$ correlations are more physical, and it may be directly related to radiation physics of GRBs. The tight correlations are possible to be used as standard candles.

Type Ic supernovae (SNe) can be classified as: normal SNe Ic, type Ic super-luminous SNe (SLSNe Ic), X-ray flash (XRF) connected SNe, broad-line SNe Ic (SN Ic-BL), and gamma-ray burst (GRB) connected SNe. Here we suggest an inner connection for all kinds of SNe Ic based on a pair of jets being successfully launched: a normal SN Ic is a normal core collapsar without jets launched; a GRB associated SN Ic is a core collapsar with relativistic jets launched and successfully breaking out the envelope of the progenitor; an XRF associated SN Ic is a core collapsar with jets launched but can only develop a relativistic shock breakout; an SN Ic-BL is an off-axis GRB or XRF associated supernova; and an SLSN Ic is close to the XRF-SN Ic but the shock breakout is not relativistic and most of the jet energy is deposited into the supernova component. Based on the luminosity-distance diagram, we derived the luminosity function of all different types of SNe Ic as a whole. We also show that the normal SNe Ic and GRB connected SNe Ic have similar accumulative distributions.

The jet breaks in the afterglow lightcurves of short gamma-ray bursts (SGRBs), rarely detected so far, are crucial for estimating the half-opening angles of the ejecta ($\theta_{\rm j}$) and hence the neutron star merger rate. In this work we report the detection of jet decline behaviors in GRB 150424A and GRB 160821B and find $\theta_{\rm j}\sim 0.1$ rad. Together with five events reported before 2015 and other three "identified" recently (GRB 050709, GRB 060614 and GRB 140903A), we have a sample consisting of nine SGRBs and one long-short GRB with reasonably estimated $\theta_{\rm j}$. In particular, three {\it Swift} bursts in the sample have redshifts $z\leq 0.2$, with which we estimate the local neutron star merger rate density {to be $\sim 1109^{+1432}_{-657}~{\rm Gpc^{-3}~yr^{-1}}$ or $162^{+140}_{-83} {\rm Gpc^{-3}yr^{-1}}$ if the narrowly-beamed GRB 061201 is excluded}. Inspired by the typical $\theta_{\rm j}\sim 0.1$ rad found currently, we further investigate whether the off-beam GRBs (in the uniform jet model) or the off-axis events (in the structured jet model) can significantly enhance the GRB/GW association or not. For the former the enhancement is at most moderate, while for the latter the enhancement can be much greater and a high GRB/GW association probability of $\sim 10\%$ is possible. We also show that the data of GRB 160821B may contain a macronova/kilonova emission component with a temperature of $\sim 3100$ K at $\sim 3.6$ days after the burst and more data are needed to ultimately clarify.

GRB 170817A is the first short gamma-ray burst (GRB) with direct detection of the gravitational-wave radiation and also the spectroscopically identified macronova emission (i.e., AT 2017gfo). The prompt emission of this burst, however, is underluminous in comparison with the other short GRBs with known redshift. In this work, we examine whether GRB 170817A is indeed unique. We firstly show that GRB 130603B/macronova may be the on-axis "analogs" of GRB 170817A/AT 2017gfo, and the extremely dim { but long-lasting} afterglow emission of GRB 170817A may suggest a low number density ($\sim 10^{-5}~{\rm cm^{-3}}$) of its circumburst medium { and a structured outflow}. We then discuss whether GRB 070923, GRB 080121, GRB 090417A, GRB 111005A, and GRB 170817A form a new group of very nearby underluminous GRBs originated from neutron star mergers. If the short events GRB 070923, GRB 080121, and GRB 090417A are indeed at a redshift of $\sim 0.076,~0.046,~0.088$, respectively, their isotropic energies of the prompt emission are $\sim 10^{47}$ erg and thus comparable to the other two events. The non-detection of optical counterparts of GRB 070923, GRB 080121, GRB 090417A, and GRB 111005A, however, strongly suggests that the macronovae from neutron star mergers are significantly diverse in luminosities or, alternatively, there is the other origin channel (for instance, the white dwarf and black hole mergers). We finally suggest that GW170817/GRB 170817A are likely not alone and similar events will be detected by the upgraded/upcoming gravitational-wave detectors and the electromagnetic monitors.

GRB 111005A, one long duration gamma-ray burst (GRB) occurred within a metal-rich environment that lacks massive stars with $M_{\rm ZAMS}\geq 15M_\odot$, is not coincident with supernova emission down to stringent limit and thus should be classified as a "long-short" GRB (lsGRB; also known as SN-less long GRB or hybrid GRB), like GRB 060505 and GRB 060614. In this work we show that in the neutron star merger model, the non-detection of the optical/infrared emission of GRB 111005A requires a sub-relativistic neutron-rich ejecta with the mass of $\leq 0.01~M_\odot$, (significantly) less massive than that of GRB 130603B, GRB 060614 and GRB 050709. The lsGRBs are found to have a high rate density and the neutron star merger origin model can be unambiguously tested by the joint observations of the second generation gravitational wave (GW) detectors and the full-sky gamma-ray monitors such as Fermi-GBM and the proposing GECAM. If no lsGRB/GW association is observed in 2020s, alternative scenarios have to be systematically investigated. With the detailed environmental information achievable for the very-nearby events, a novel kind of merger or explosion origin may be identified.

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

On 17 August 2017, a gravitational wave event (GW170817) and an associated short gamma-ray burst (GRB 170817A) from a binary neutron star merger had been detected. The followup optical/infrared observations also identified the macronova/kilonova emission (AT2017gfo). In this work we discuss some implications of the remarkable GW170817/GRB 170817A/AT2017gfo association. We show that the $\sim 1.7$s time delay between the gravitational wave (GW) and GRB signals imposes very tight constraint on the superluminal movement of gravitational waves (i.e., the relative departure of GW velocity from the speed of light is $\leq 4.3\times 10^{-16}$) or the possible violation of weak equivalence principle (i.e., the difference of the gamma-ray and GW trajectories in the gravitational field of the galaxy and the local universe should be within a factor of $\sim 3.4\times 10^{-9}$). The so-called Dark Matter Emulators and a class of contender models for cosmic acceleration ("Covariant Galileon") are ruled out, too. The successful identification of Lanthanide elements in the macronova/kilonova spectrum also excludes the possibility that the progenitors of GRB 170817A are a binary strange star system. The high neutron star merger rate (inferred from both the local sGRB data and the gravitational wave data) together with the significant ejected mass strongly suggest that such mergers are the prime sites of heavy r-process nucleosynthesis.

Cosmic rays interact with the Earth's atmosphere to produce $^{14}$C, which can be absorbed by trees. Therefore, rapid increases of $^{14}$C in tree rings can be used to probe previous cosmic-ray events. By this method, three $^{14}$C rapidly increasing events have been found. Plausible causes of these events include large solar proton events, supernovae or short gamma-ray bursts. However, due to the lack of measurements of $^{14}$C by year, the occurrence frequency of such $^{14}$C rapidly increasing events is poorly known. In addition, rapid increases may be hidden in the IntCal13 data with five-year resolution. Here we report the result of $^{14}$C measurements using an ancient buried tree during the period between BC 3388 and 3358. We find a rapid increase of about 9\textperthousand~ in the $^{14}$C content from BC 3372 to BC 3371. We suggest that this event could originate from a large solar proton event.

Assuming that the early optical emission is dominated by the external reverse shock (RS) in the standard model of gamma-ray bursts (GRBs), we intend to constrain RS models with the initial Lorentz factor $\Gamma_0$ of the outflows based on the ROTSE-III observations. We consider two cases of the RS behavior: the relativistic shock and the non-relativistic shock. For homogeneous interstellar medium (ISM) and wind circum-burst environment, the constraints can be achieved by the fact that the peak flux $F_{\rm \nu}$ at the RS crossing time should be lower than the observed upper limit $F_{\rm \nu,limit}$. We consider the different spectral regimes that the observed optical frequency $\nu_{\rm opt}$ may locate in, which are divided by the orders for the minimum synchrotron frequency $\nu_{\rm m}$ and the cooling frequency $\nu_{\rm c}$. {\bf Considering the homogeneous and wind environment around GRBs, we find that the relativistic RS case can be constrained by the (upper and lower) limits of $\Gamma_0$ in a large range from about hundreds to thousands for 36 GRBs reported by ROTSE-III. The constraints on the non-relativistic RS case are achieved with limits of $\Gamma_0$ for 26 bursts ranging from $\sim 30$ to $\sim 350$.} The lower limits of $\Gamma_0$ achieved for the relativistic RS model is disfavored based on the previously discovered correlation between the initial Lorentz factor $\Gamma_0$ and the isotropic gamma-ray energy $E_{\rm \gamma, iso}$ released in prompt phase.

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

By fitting the bolometric light curves of 31 super-luminous supernovae (SLSNe) with the magnetar engine model, we derive the ejecta masses and magnetar parameters for these SLSNe. The lower boundary of magnetic field strengths of SLSN magnetars can be set just around the critical field strength $B_{\rm c}$ of electron Landau quantization. In more details, SLSN magnetars can further be divided into two subclasses of magnetic fields of $\sim(1-5)B_{\rm c}$ and $\sim(5-10)B_{\rm c}$, respectively. It is revealed that these two subclasses of magnetars are just associated with the slow-evolving and fast-evolving bolometric light curves of SLSNe. In comparison, the magnetars harbored in gamma-ray bursts (GRBs) and associated hypernovae are usually inferred to have much higher magnetic fields with a lower boundary about $\sim10B_{\rm c}$. This robustly suggests that it is the magnetic fields that play the crucial role in distinguishing SLSNe from GRBs/hypernovae. The rotational energy of SLSN magnetars are found to be correlated with the masses of supernova ejecta, which provides a clue to explore the nature of their progenitors. Moreover, the distribution of ejecta masses of SLSNe is basically intermediate between those of normal core-collapse supernovae and hypernovae. This could indicate an intrinsic connection among these different stellar explosions.

Based on the luminosity-distance diagram, we propose a method to quickly estimate the luminosity function for any certain astrophysical objects. Giving the mean distance between any two objects at a given luminosity range, we can find the relation between the mean distance and the luminosity, and consequently can obtain the luminosity function. Not like the straightforward counting method, this method does not need a complete sample. The only requirement is that the object distributes uniformly in space. We apply this method to a simulated sample, and find it can produce the luminosity function properly. This method can also be used for energy function.

One of the cornerstones of modern physics is Einstein's special relativity, with its constant speed of light and zero photon mass assumptions. Constraint on the rest mass m_{\gamma} of photons is a fundamental way to test Einstein's theory, as well as other essential electromagnetic and particle theories. Since non-zero photon mass can give rise to frequency-(or energy-) dependent dispersions, measuring the time delay of photons with different frequencies emitted from explosive astrophysical events is an important and model-independent method to put such a constraint. The cosmological gamma-ray bursts (GRBs), with short time scales, high redshifts as well as broadband prompt and afterglow emissions, provide an ideal testbed for m_{\gamma} constraints. In this paper we calculate the upper limits of the photon mass with GRB early time radio afterglow observations as well as multi-band radio peaks, thus improve the results of Schaefer (1999) by nearly half an order of magnitude.

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.

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

The Active Particle-induced X-ray Spectrometer (APXS) is one of the payloads on board the Yutu rover of Chang'E-3 mission. In order to assess the instrumental performance of APXS, a ground verification test was done for two unknown samples (basaltic rock, mixed powder sample). In this paper, the details of the experiment configurations and data analysis method are presented. The results show that the elemental abundance of major elements can be well determined by the APXS with relative deviations < 15 wt. % (detection distance = 30 mm, acquisition time = 30 min). The derived detection limit of each major element is inversely proportional to acquisition time and directly proportional to detection distance, suggesting that the appropriate distance should be < 50mm.

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

We propose a direct and model-independent method to constrain the Lorentz factor of a relativistically expanding object, like gamma-ray bursts. Only the measurements, such as thermal component of the emission, the distance and the variable time scale of the light curve, are used. If the uncertainties are considered, we will obtain lower limits of the Lorentz factor instead. We apply this method to GRB 090618 and get a lower limit of the Lorentz factor to be 22. The method can be used to any relativistically moving object, such as gamma-ray bursts, blazars, and soft gamma-ray repeaters, providing the thermal component of the emission being observed.

In observations, it is important to deproject the two fundamental quantities characterizing a bar, i.e., its length ($a$) and ellipticity ($e$), to face-on values before any careful analyses. However, systematic estimation on the uncertainties of the commonly used deprojection methods is still lacking. Simulated galaxies are well suited in this study. We project two simulated barred galaxies onto a 2D plane with different bar orientations and disk inclination angles ($i$). Bar properties are measured and deprojected with the popular deprojection methods in the literature. Generally speaking, deprojection uncertainties increase with increasing $i$. All the deprojection methods behave badly when $i$ is larger than $60^\circ$, due to vertical thickness of the bar. Thus, future statistical studies of barred galaxies should exclude galaxies more inclined than $60^\circ$. At moderate inclination angles ($i\leq60^\circ$), 2D deprojection methods (analytical and image stretching) and Fourier-based methods (Fourier decomposition and bar-interbar contrast) perform reasonably well with uncertainties $\sim10\%$ in both the bar length and ellipticity. Whereas the uncertainties of the 1D analytical deprojection can be as high as $100\%$ in certain extreme case. We find that different bar measurement methods show systematic differences in the deprojection uncertainties. We further discuss the deprojection uncertainty factors with the emphasis on the most important one, i.e., the 3D structure of the bar itself. We construct two triaxial toy bar models that can qualitatively reproduce the results of the 1D and 2D analytical deprojections; they confirm that the vertical thickness of the bar is the main source of uncertainties.

Anomalies of the cosmic microwave background (CMB) maps have been widely acquainted nowadays from the Wilkinson Microwave Anisotropy Probe (WMAP) satellite to the Planck satellite. One of the anomalies is a multipole alignment from $l=2$ to $l=5$. In our work, we investigate the angular distribution of gamma-ray bursts (GRBs) to find whether there is the same anomaly of GRB as CMB. We perform spherical harmonics expansion on GRB samples to derive coefficients of a few first terms of the expansion terms and find that there is rough multipole alignment from $l=2$ to $l=4$ while the dipole and $l=5$ multipole is in a distant direction, and that the quadrupole is obviously planar and the other ones are normal.

According to the analysis of the $^{14}$C content of two Japanese trees over a period of approximately 3000 years at high time resolution, Miyake (2012) found a rapid increase at AD 774-775 and later on at AD 992-993 (Miyake 2013). This corresponds to a high-energy event happened within one year that input $\gamma$-ray energy about 7$\times{}$10$^{24}$erg to the Earth, leaving the origin a mystery. Such strong event should have an unusual optical counterpart, and have been recorded in historical literature. We searched Chinese historical materials around AD 744-775 and AD 992-993, but no remarkable event was found except a violent thunderstorm in AD 775. However, the possibility of a thunderstorm containing so much energy is still unlikely. We conclude the event caused the $^{14}$C increase is still unclear. This event most probably has no optical counterpart, and short gamma-ray burst, giant flare of a soft gamma-ray repeater and terrestrial $\gamma$-ray flash may all be the candidates.

The three Mn I lines at 17325, 17339 and 17349 \AA~are among the 25 strongest lines ($\log(gf)>0.5$) in the H-band. They are all heavily broadened due to hyperfine structure and the profiles of these lines have so far not been understood. Earlier studies of these lines even suggested that they were blended. In this work, the profiles of these three infra-red (IR) lines have been studied theoretically and compared to experimental spectra to assist in the complete understanding of the solar spectrum in the IR. It is shown that the structure of these lines can not be described in the conventional way by the diagonal $A$ and $B$ hyperfine interaction constants. The off-diagonal hyperfine interaction not only has large impact on the energies of the hyperfine levels, but also introduces a large intensity redistribution among the hyperfine lines, changing the line profiles dramatically. By performing large-scale calculations of the diagonal and off-diagonal hyperfine interaction and $gf$-values between the upper and lower hyperfine levels and using a semi-empirical fitting procedure, agreement between our synthetic and experimental spectra was achieved. Furthermore, we compare our results with observations of stellar spectra. The spectra of the Sun and the K1.5 III red giant star Arcturus were modelled in the relevant region, $1.73-1.74$ $\mu$m using our theoretically predicted $gf$-values and energies for each individual hyperfine line. Satisfactory fits were obtained and clear improvements were found using our new data compared with the old available Mn I data. A complete list of energies and $gf$-values for all the $3d^54s({^7S})4d$ e$^{6}$D $-$ $3d^54s({^7S})4f$ w$^{6}$F hyperfine lines are available as supplementary material online, whereas only the stronger lines are presented and discussed in detail in this paper.

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

We present a new representation of light curves, which is quite different from the binning method. Instead of choosing uniform bins, the reciprocal of interval between two successive photons is adopted to represent the counting rate. A primary application of this light curve is to compute the power density spectrum by means of Lomb Periodogram and to find possible periods. To verify this new representation, we apply this method to artificial periodic data and some known periodic celestial objects, and the periods are all correctly found. Compared with the traditional fast Fourier transform method, our method does not rely on the bin size and has a spontaneously high time resolution, guaranteeing a wide frequency range in power density spectrum, and is especially useful when the photons are rare for its little information losses. Some other applications of the new light curve, like pulse identification, variability and spectral time lag, are also discussed.

On 13 December 2012, Chang'e-2 conducted a successful flyby of the near-Earth asteroid 4179 Toutatis at a closest distance of 770 $\pm$ 120 meters from the asteroid's surface. The highest-resolution image, with a resolution of better than 3 meters, reveals new discoveries on the asteroid, e.g., a giant basin at the big end, a sharply perpendicular silhouette near the neck region, and direct evidence of boulders and regolith, which suggests that Toutatis may bear a rubble-pile structure. Toutatis' maximum physical length and width are (4.75 $\times$ 1.95 km) $\pm$10$\%$, respectively, and the direction of the +$z$ axis is estimated to be (250$\pm$5$^\circ$, 63$\pm$5$^\circ$) with respect to the J2000 ecliptic coordinate system. The bifurcated configuration is indicative of a contact binary origin for Toutatis, which is composed of two lobes (head and body). Chang'e-2 observations have significantly improved our understanding of the characteristics, formation, and evolution of asteroids in general.

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.

Long-duration gamma-ray bursts (GRBs) at z < 1 are in most cases found to be accompanied by bright, broad-lined Type Ic supernovae (SNe Ic-BL). The highest-energy GRBs are mostly located at higher redshifts, where the associated SNe are hard to detect observationally. Here we present early and late observations of the optical counterpart of the very energetic GRB 130427A. Despite its moderate redshift z = 0.3399+/-0.0002, GRB 130427A is at the high end of the GRB energy distribution, with an isotropic-equivalent energy release of Eiso ~ 9.6x10^53 erg, more than an order of magnitude more energetic than other GRBs with spectroscopically confirmed SNe. In our dense photometric monitoring, we detect excess flux in the host-subtracted r-band light curve, consistent with what expected from an emerging SN, ~0.2 mag fainter than the prototypical SN 1998bw. A spectrum obtained around the time of the SN peak (16.7 days after the GRB) reveals broad undulations typical of SNe Ic-BL, confirming the presence of a SN, designated SN 2013cq. The spectral shape and early peak time are similar to those of the high expansion velocity SN 2010bh associated with GRB 100316D. Our findings demonstrate that high-energy long-duration GRBs, commonly detected at high redshift, can also be associated with SNe Ic-BL, pointing to a common progenitor mechanism.

Complex network approaches have been recently developed as an alternative framework to study the statistical features of time-series data. We perform a visibility-graph analysis on both the daily and monthly sunspot series. Based on the data, we propose two ways to construct the network: one is from the original observable measurements and the other is from a negative-inverse-transformed series. The degree distribution of the derived networks for the strong maxima has clear non-Gaussian properties, while the degree distribution for minima is bimodal. The long-term variation of the cycles is reflected by hubs in the network which span relatively large time intervals. Based on standard network structural measures, we propose to characterize the long-term correlations by waiting times between two subsequent events. The persistence range of the solar cycles has been identified over 15\,--\,1000 days by a power-law regime with scaling exponent $\gamma = 2.04$ of the occurrence time of the two subsequent and successive strong minima. In contrast, a persistent trend is not present in the maximal numbers, although maxima do have significant deviations from an exponential form. Our results suggest some new insights for evaluating existing models. The power-law regime suggested by the waiting times does indicate that there are some level of predictable patterns in the minima.

The X-ray emission from Swift J1644+57 is not steadily decreasing instead it shows multiple pulses with declining amplitudes. We model the pulses as reverse shocks from collisions between the late ejected shells and the externally shocked material, which is decelerated while sweeping the ambient medium. The peak of each pulse is taken as the maximum emission of each reverse shock. With a proper set of parameters, the envelope of peaks in the light curve as well as the spectrum can be modelled nicely.

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

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

The shallow decay phase of the early X-ray afterglow in gamma-ray bursts discovered by {\it Swift} is a widely discussed topic. As the spectral index does not change at the transferring of the shallow decay phase to a normal phase, it implies this transferring should be a dynamical change rather than a spectral evolution. We suggest both the shallow decay phase and the normal phase are from the external shock in a wind environment, while the transferring time is the deceleration time. We apply this model to GRBs 050319 and 081008, and find them can be well explained by choosing a proper set of parameters.

We compile 23 Gamma-ray Bursts (GRBs) and 21 blazars with estimated Doppler factors, and the Doppler factors of GRBs are estimated from their Lorentz factors by assuming their jet viewing angles \theta->0 degree. Using the conventional assumption that the prompt emission of GRBs is dominated by the synchrotron radiation, we calculate the synchrotron luminosity of GRBs from their total isotropic energy and burst duration. Intriguingly, we discover a uniform correlation between the synchrotron luminosity and Doppler factor, L_syn \propto D^3.1, for GRBs and blazars, which suggests that they may share some similar jet physics. One possible reason is that GRBs and blazars have, more or less, similar intrinsic synchrotron luminosities and both of them are strongly enhanced by the beaming effect. After Doppler and redshift-correction, we find that the intrinsic peak energy of the GRBs ranges from 0.1 to 3 keV with a typical value of 1 keV. We further correct the beaming effect for the observed luminosity of GRBs and find that there exists a positive correlation between the intrinsic synchrotron luminosity and peak energy for GRBs, which is similar to that of blazars. Our results suggest that both the intrinsic positive correlation and the beaming effect may be responsible for the observed tight correlation between the isotropic energy and the peak energy in GRBs (so called "Amati" relation).

Fermi observations of GeV emission from GRBs have suggested that the Lorentz factor of some GRBs is around a thousand or even higher. At the same time the same Fermi observations have shown an extended GeV emission indicating that this higher energy emission might be a part of the afterglow and it does not come from the same region as the lower energy prompt emission. If this interpretation is correct than we should reconsider the opacity limits on the Loretnz factor of the emitting regions which are based on a one-zone model. We describe here a two-zone model in which the GeV photons are emitted in a larger radius than the MeV photons and we calculate the optical depth for pair creation of a GeV photon passing the lower energy photons shell. We find that, as expected, the new two-zone limits on the Lorentz factor are significantly lower. The corresponding limits for the Fermi bursts are lower by a factor of five compared to the one-zone model and it is possible that both the MeV and GeV regions have relatively modest Lorentz factors (~200 - 400).

In the fireball model, it is more physically realistic that gamma-ray burst (GRB) ejecta have a range of bulk Lorentz factors (assuming $M\propto \Gamma^{-s}$). The low Lorentz factor part of the ejecta will catch up with the high Lorentz factor part when the latter is decelerated by the surrounding medium to a comparable Lorentz factor. Such a process will develop a long-lasting weak reverse shock until the whole ejecta are shocked. Meanwhile, the forward shocked materials are gradually supplied with energy from the ejecta that are catching-up, and thus the temporal decay of the forward shock emission will be slower than that without an energy supply. However, the reverse shock may be strong. Here, we extend the standard reverse-forward shock model to the case of radially nonuniform ejecta. We show that this process can be classified into two cases: the thick shell case and the thin shell case. In the thin shell case, the reverse shock is weak and the temporal scaling law of the afterglow is the same as that in Sari & Mesz (2000). However, in the thick shell case, the reverse shock is strong and thus its emission dominates the afterglow in the high energy band. Our results also show slower decaying behavior of the afterglow due to the energy supply by low Lorentz factor materials, which may help the understanding of the plateau observed in the early optical and X-ray afterglows.

While it is generally agreed that the emitting regions in Gamma-Ray Bursts (GRBs) move ultra relativistically towards the observer, different estimates of the initial Lorentz factors, $\Gamma_0$, lead to different, at times conflicting estimates. We show here that the quiet periods in which the signals goes down below the instrumental thresholds, put strong upper limits on the values of $\Gamma_0$. According to the standard internal-external shocks model an external shock should develop during the prompt stage. This external shock radiates in the hard X-rays to soft gamma-rays bands and this emission should be seen as a smooth background signal. The observed deep minima indicate that this contribution is negligible. This limits, in turn, $\Gamma_0$. We obtain upper limits on $\Gamma_0$ for several bursts with typical values around hundreds. We compare these values with those obtained by the other methods, which typically yield lower limits. The results are marginally consistent leaving only a narrow range of allowed values for $\Gamma_0$.

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

We calculate the high energy (sub-GeV to TeV) prompt and afterglow emission of GRB 080319B that was distinguished by a naked-eye optical flash and by an unusual strong early X-ray afterglow. There are three possible sources for high energy emission: the prompt optical and $\gamma$-ray photons IC scattered by the accelerated electrons, the prompt photons IC scattered by the early external reverse-forward shock electrons, and the higher band of the synchrotron and the synchrotron self-Compton emission of the external shock. There should have been in total {hundreds} high energy photons detectable for the Large Area Telescope (LAT) onboard the Fermi satellite, and {tens} photons of those with energy $> 10$ GeV. The $> 10$ GeV emission had a duration about twice that of the soft $\gamma$-rays. AGILE could have observed these energetic signals if it was not occulted by the Earth at that moment. The physical origins of the high energy emission detected in GRB 080514B, GRB 080916C and GRB 081024B are also discussed. These observations seem to be consistent with the current high energy emission models.

Inverse Compton (IC) scattering is one of two viable mechanisms that can produce the prompt gamma-ray emission in Gamma-Ray Bursts. IC requires low energy seed photons and a population of relativistic electrons that upscatter them. The same electrons can upscatter the gamma-ray photons to even higher energies in the TeV range. Using the current upper limits on the prompt optical emission we show here that under general conservative assumption the IC mechanism suffers from an "energy crisis". Namely, IC will over-produce a very high energy component that would carry much more energy than the observed prompt gamma-rays. Our analysis is general and it makes no assumptions on the specific mechanism that produces the relativistic electrons population.

The extremely bright optical flash that accompanied GRB 080319B suggested, at first glance, that the prompt $\gamma$-rays in this burst were produced by Synchrotron self Compton (SSC). We analyze here the observed optical and $\gamma$ spectrum. We find that the very strong optical emission poses, due to self absorption, very strong constraints on the emission processes and put the origin of the optical emission at a very large radius, almost inconsistent with internal shock. Alternatively it requires a very large random Lorentz factor for the electrons. We find that SSC could not have produced the prompt $\gamma$-rays. We also show that the optical emission and the $\gamma$ rays could not have been produced by synchrotron emission from two populations of electron within the same emitting region. Thus we must conclude that the optical and the $\gamma$-rays were produced in different physical regions. A possible interpretation of the observations is that the $\gamma$-rays arose from internal shocks but the optical flash resulted from external shock emission. This would have been consistent with the few seconds delay observed between the optical and $\gamma$-rays signals.

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

We analyze the hitherto available space-based X-ray data as well as ground-based optical data of the X-ray transient 080109/SN2008D. From the data we suggest that (i) The initial transient (\lesssim 800 sec) is attributed to the reverse shock emission of a mildly relativistic (\Gamma \sim a few) outflow stalled by the dense stellar wind. (ii) The subsequent X-ray afterglow (\lesssim 2\times 10^4 sec) can be ascribed to the forward shock emission of the outflow, with a kinetic energy \sim 10^{46} erg, when sweeping up the stellar wind medium. (iii) The late X-ray flattening (\gtrsim 2\times 10^4$ sec) is powered by the fastest non-decelerated component of SN2008D's ejecta. (iv) The local event rate of X-ray transient has a lower limit of \sim 1.6\times 10^4 yr^{-1} Gpc^{-3}, indicating a vast majority of X-ray transients have a wide opening angle of \gtrsim 100 degree. The off-axis viewing model is less likely. (v) Transient 080109/SN2008D may lead to a continuum from GRB-SN to under-luminous GRB-/XRF-SN to X-ray transient-SN and to ordinary Ibc SN (if not every Ibc SN has a relativistic jet), as shown in Figure 2 of this Letter.

By considering synchrotron radiative process in the internal shock model and assuming that all internal shocks are nearly equally energetic, we analyze the gamma-ray burst (GRB) emission at different radii corresponding to different observed times. We apply this model to GRB 050904 and find that our analytical results can provide a natural explanation for the multi-band observations of GRB 050904. This suggests that the X-ray flare emission and the optical emission of this burst could have originated from internal shocks being due to collisions among nearly-equally-energetic shells ejected from the central engine. Thus GRB 050904 appears to be a burst with super-long central engine activity.

Using the conservation of energy and momentum during collisions of any two shells, we consider the efficiency of gamma-ray bursts by assuming that the ejecta from the central engine are equally massive and have the same Lorentz factors. We calculate the efficiency and the final Lorentz factor of the merged whole shell for different initial diversities of Lorentz factors and for different microscopic radiative efficiency. As a result, a common high efficiency in the range of 0.1 to 0.9 is considerable, and a very high value near 100% is also reachable if the diversity of the Lorentz factors is large enough.

By neglecting sideways expansion of gamma-ray burst (GRB) jets and assuming their half-opening angle distribution, we estimate the detectability of orphan optical afterglows. This estimation is carried out by calculating the durations of off-axis optical afterglows whose flux density exceeds a certain observational limit. We show that the former assumption leads to more detectable orphans, while the latter suppresses the detectability strongly compared with the model with half-opening angle $\theta_j=0.1$. We also considered the effects of other parameters, and find that the effects of the ejecta energy $E_j$ and post-jet-break temporal index $-\alpha_2$ are important but that the effects of the electron-energy distribution index $p$, electron energy equipartition factor $\epsilon_e$, and environment density $n$ are insignificant. If $E_j$ and $\alpha_2$ are determined by other methods, one can constrain the half-opening angle distribution of jets by observing orphan afterglows. Adopting a set of "standard" parameters, the detectable rate of orphan afterglows is about $1.3\times 10^{-2} {deg}^{-2} {yr}^{-1}$, if the observed limiting magnitude is 20 in R-band.

Similar to the pulsar, the magnetic axis and the spin axis of the gamma-ray burst source may not lie on the same line. This may cause a ring-like jet due to collimation of the precessing magnetic axis. We analyze the tail emission from such a jet, and find that it has a shallow decay phase with temporal index equal to -1/2 if the Lorentz factor of the ejecta is not very high. This phase is consistent with the shallow decay phase of some early X-ray afterglow detected by {\it{swift}}. The ring-like jet has a tail cusp with sharp rising and very sharp decay. This effect can provide an explanation for the re-brightening and sharp decay of the X-ray afterglow of GRB 050709.

Using the fluid-magnetofluid formalism, we obtain axisymmetric stability criteria for a composite disc system consisting of stellar and gaseous magnetized singular isothermal discs (MSIDs). Lou & Zou recently constructed exact global stationary configurations for both axisymmetric and nonaxisymmetric coplanar MHD perturbations in such a composite MSID system and proposed the MHD $D_s-$criteria for axisymmetric stability by the hydrodynamic analogy. In a different perspective, we derive and analyze here the time-dependent WKBJ dispersion relation in the low-frequency and tight-winding regime to examine axisymmetric stability properties. By introducing a rotational Mach number $D_s$ for the ratio of the stellar rotation speed $V_s$ to the stellar velocity dispersion $a_s$, one readily determines the stable range of $D_s^2$ numerically to establish the $D_s-$criteria for axisymmetric MSID stability. Those MSID systems rotating either too fast (ring ragmentation) or too slow (Jeans collapse) are unstable. Our WKBJ results of (M)SID instability provide physically compelling explanations for the stationary analysis of Lou & Zou. We further introduce an effective MHD $Q$ parameter for a composite MSID system and compare with the earlier work of Elmegreen, Jog and Shen & Lou. As expected, an axisymmetric dark matter halo enhances the stability against axisymmetric disturbances in a composite partial MSID system.

When a cold shell sweeps up the ambient medium, a forward shock and a reverse shock will form. We analyze the reverse-forward shocks in a wind environment, including their dynamics and emission. An early afterglow is emitted from the shocked shell, e.g., an optical flash may emerge. The reverse shock behaves differently in two approximations: relativistic and Newtonian cases, which depend on the parameters, e.g., the initial Lorentz factor of the ejecta. If the initial Lorentz factor is much less than $114 E_{53}^{1/4} \Delta_{0,12}^{-1/4} A_{*,-1}^{-1/4}$, the early reverse shock is Newtonian. This may take place for the wider of a two-component jet, an orphan afterglow caused by a low initial Lorentz factor, and so on. The synchrotron self absorption effect is significant especially for the Newtonian reverse shock case, since the absorption frequency $\nu_a$ is larger than the cooling frequency $\nu_c$ and the minimum synchrotron frequency $\nu_m$ for typical parameters. For the optical to X-ray band, the flux is nearly unchanged with time during the early period, which may be a diagnostic for the low initial Lorentz factor of the ejecta in a wind environment. We also investigate the early light curves with different wind densities, and compare them with these in the ISM model.

We construct aligned and unaligned stationary perturbation configurations in a composite system of stellar and coplanarly magnetized gaseous singular isothermal discs (SIDs) coupled by gravity. In comparison with SID problems studied earlier, there exist three possible classes of stationary solutions allowed by more dynamic freedoms. Our exact global perturbation solutions and critical points are valuable for testing numerical magnetohydrodynamic codes. For galactic applications, our model analysis contains more realistic elements and offer useful insights for structures and dynamics of disc galaxies consisting of stars and magnetized gas.