Owing to the detection of an extremely bright fast radio burst (FRB) 200428 associated with a hard X-ray counterpart from the magnetar soft gamma-ray repeater (SGR) 1935+2154, the distance of SGR 1935+2154 potentially hosted in the supernova remnant (SNR) G57.2+0.8 can be revisited. Under the assumption that the SGR and the SNR are physically related, in this Letter, by investigating the dispersion measure (DM) of the FRB contributed by the foreground medium of our Galaxy and the local environments and combining with other observational constraints, we find that the distance of SGR 1935+2154 turns out to be $9.0\pm2.5\,$kpc and the SNR radius falls into $10$ to $18\,$pc since the local DM contribution is as low as $0-18\,$pc cm$^{-3}$. These results are basically consistent with the previous studies. In addition, an estimate for the Faraday rotation measure of the SGR and SNR is also carried out.

Very recently, an extremely bright fast radio burst (FRB) 200428 with two sub-millisecond pulses was discovered to come from the direction of the Galactic magnetar SGR 1935+2154, and an X-ray burst (XRB) counterpart was detected simultaneously. These observations favor magnetar-based interior-driven models. In this Letter, we propose a different model for FRB 200428 associated with an XRB from SGR 1935+2154, in which a magnetar with high proper velocity encounters an asteroid of mass $\sim10^{20}\,$g. This infalling asteroid in the stellar gravitational field is first possibly disrupted tidally into a great number of fragments at radius $\sim {\rm a\,\,few}$ times $10^{10}\,$cm, and then slowed around the Alfv$\acute{\rm e}$n radius by an ultra-strong magnetic field and in the meantime two major fragments of mass $\sim 10^{17}\,$g that cross magnetic field lines produce two pulses of FRB 200428. The whole asteroid is eventually accreted onto the poles along magnetic field lines, impacting the stellar surface, creating a photon-e$^\pm$ pair fireball trapped initially in the stellar magnetosphere, and further leading to an XRB. We show that this gravitationally-powered model can interpret all of the observed features self-consistently.

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

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

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

The periodic activity of repeating fast radio burst (FRB) 180916.J0158+65 was recently reported by the CHIME/FRB Collaboration team. 28 bursts from this source not only show a $\sim16$-day period with an active phase of $\sim 4.0$ days but also exhibit a broken power law in differential energy distribution. In this paper, we suggest that FRB 180916.J0158+65-like periodic FRBs would provide a unique probe of extragalactic asteroid belts (EABs), based on our previously-proposed pulsar-EAB impact model, in which repeating FRBs arise from an old-aged, slowly-spinning, moderately-magnetized pulsar traveling through an EAB around another stellar-mass object. These two objects form a binary and thus the observed period is in fact the orbital period. We show that this model can be used to well interpret all of the observed data of FRB 180916.J0158+65. Furthermore, we constrain the EAB's physical properties and find that (1) the outer radius of the EAB is at least an order of magnitude smaller than that of its analogue in the solar system, (2) the differential size distribution of the EAB's asteroids at small diameters (large diameters) is shallower (steeper) than that of solar-system small objects, and (3) the two belts have a comparable mass.

In Dec. 2011 PSR B0540-69 experienced a spin-down rate transition (SRT), after which the spin-down power of the pulsar increased by ~36%. About 1000 days after the SRT, the X-ray luminosity of the associated pulsar wind nebula (PWN) was found to brighten by 32+/-8%. After the SRT, the braking index n of PSR B0540-69 changes from n=2.12 to n=0.03 and then keeps this value for about five years before rising to n=0.9 in the following years. We find that most of the current models have difficulties in explaining the measured braking index. One exceptive model of the braking index evolution is the increasing dipole magnetic field of PSR B0540-69. We suggest that the field increase may result from some instabilities within the pulsar core that enhance the poloidal component at the price of toroidal component of the magnetic field. The increasing dipole magnetic field will result in the X-ray brightening of the PWN. We fit the PWN X-ray light curve by two models: one assumes a constant magnetic field within the PWN during the brightening and the other assumes an enhanced magnetic field proportional to the energy density of the PWN. It appears that the two models fit the data well, though the later model seems to fit the data a bit better. This provides marginal observational evidence that magnetic field in the PWN is generated by the termination shock. Future high-quality and high-cadence data are required to draw a solid conclusion.

The first repeating fast radio burst (FRB), FRB 121102, was found to be associated with a spatially coincident, persistent nonthermal radio source, but the origin of the persistent emission remains unknown. In this paper, we propose that the persistent emission is produced via synchrotron-heating process by multiple bursts of FRB 121102 in a self-absorbed synchrotron nebula. As a population of bursts of the repeating FRB absorbed by the synchrotron nebula, the energy distribution of electrons in the nebula will change significantly. As a result, the spectrum of the nebula will show a hump steadily. For the persistent emission of FRB 121102, the total energy of bursts injecting into the nebula is required to be about $3.3\times10^{49}\,\unit{erg}$, the burst injection age is over $6.7\times 10^4\,\unit{yr}$, the nebula size is $\sim0.02\,\unit{pc}$, and the electron number is about $3.2\times10^{55}$. We predict that as more bursts inject, the brightness of the nebula would be brighter than the current observation, and meanwhile, the peak frequency would become higher. Due to the synchrotron absorption of the nebula, some low-frequency bursts would be absorbed, which may explain why most bursts were detected above $\sim1~\unit{GHz}$.

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.

It is widely believed that magnetars could be born in core-collapse supernovae (SNe), binary neutron star (BNS) or binary white dwarf (BWD) mergers, or accretion-induced collapse (AIC) of white dwarfs. In this paper, we investigate whether magnetars could also be produced from neutron star--white dwarf (NSWD) mergers, motivated by FRB 180924-like fast radio bursts (FRBs) possibly from magnetars born in BNS/BWD/AIC channels suggested by \cite{mar19}. By a preliminary calculation, we find that NSWD mergers with unstable mass transfer could result in the NS acquiring an ultra-strong magnetic field via the dynamo mechanism due to differential rotation and convection or possibly via the magnetic flux conservation scenario of a fossil field. If NSWD mergers can indeed create magnetars, then such objects could produce at least a subset of FRB 180924-like FRBs within the framework of flaring magnetars, since the ejecta, local environments, and host galaxies of the final remnants from NSWD mergers resemble those of BNS/BWD/AIC channels. This NSWD channel is also able to well explain both the observational properties of FRB 180924-like and FRB 180916.J0158+65-like FRBs within a large range in local environments and host galaxies.

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.

We present spectroscopy of stars in the immediate vicinity of the dwarf nova (DN) KZ Gem to confirm its identification, which had been ambiguous in the literature. Analysis of 73 radial velocities spanning from 2014 to 2019 provides a high-precision orbital period of 0.2224628(2)\,d ($\sim5.34$\,hr) and shows KZ\,Gem to be a double-lined DN. Time series photometry taken from 2016 to 2018 shows a variable double-hump modulation with a full amplitude of $\sim0.3$\,mag, along with five Gaussian-like transient events lasting $\sim30$\,min or more. Using the light curve code XRBinary and nonlinear fitting code NMfit, we obtain an optimized binary model of the dwarf nova (DN) KZ Gem, from time series photometry, consisting of a Roche-lobe-filling K type dwarf with a mass transfer rate of $2.7\,-\,7.9\times10^{-10}\,{\rm M}_{\odot}\,{\rm yr}^{-1}$ to a large, cool and thick disk surrounding a white dwarf, in an orbit with an inclination of $51^{\circ}.6(\pm1^{\circ}.4)$. Two hotspots on the disk are demonstrated to cause the observed variations in the ellipsoidal modulations from the secondary star. This physical model is compatible with the Gaia distance of KZ\,Gem.

Besides light curves and spectra, polarization provides a different powerful tool of studying the $\gamma-$ray burst (GRB) prompt phase. Compared with the time-integrated and energy-integrated polarization, time-resolved and energy-resolved polarization can deliver more physical information about the emitting region. Here we investigate time-resolved and energy-resolved polarization of GRB prompt emission using the synchrotron models. We find that the equal arrival time surface effect is very important in shaping the PD curves when the physical conditions of emitting region changes violently with radius. Polarization properties are neither correlated with the spectral lag nor the peak energy evolution patterns. Polarization properties with a mixed magnetic field are very similar to those for a corresponding ordered magnetic field but the former has a smaller polarization degree. The emission at the MeV peak can be highly polarized for a synchrotron model while it is unpolarized as predicted by a dissipative photosphere model. Future energy-resolved polarization observations can distinguish between these two models.

Young neutron stars (NSs) born in core-collapse explosions are promising candidates for the central engines of fast radio bursts (FRBs), since the first localized repeating burst FRB 121102 happens in a star forming dwarf galaxy, which is similar to the host galaxies of superluminous supernovae (SLSNe) and long gamma-ray bursts (LGRBs). However, FRB 180924 and FRB 190523 are localized to massive galaxies with low rates of star formation, compared with the host of FRB 121102. Meanwhile, the offsets between the bursts and host centers are about 4 kpc and 29 kpc for FRB 180924 and FRB 190523, respectively. These properties of hosts are similar to short gamma-ray bursts \textbf{(SGRBs)}, which are produced by mergers of binary neutron star (BNS) or neutron star-black hole (NS-BH). Therefore, the NSs powering FRBs may be formed in BNS mergers. In this paper, we study the BNS merger rates, merger times, and predict their most likely merger locations for different types of host galaxies using population synthesis method. We find that the BNS merger channel is consistent with the recently reported offsets of FRB 180924 and FRB 190523. The offset distribution of short GRBs is well reproduced by population synthesis using galaxy model which is similar to GRB hosts. The event rate of FRBs (including non-repeating and repeating), is larger than those of BNS merger and short GRBs, which requires a large fraction of observed FRBs emitting several bursts. Using curvature radiation by bunches in NS magnetospheres, we also predict the observational properties of FRBs from BNS mergers, including the dispersion measure, and rotation measure. At late times ($t\geq1$yr), the contribution to dispersion measure and rotation measure from BNS merger ejecta could be neglected.

In this paper, we present the study for the energy reservoir powering the light curves (LCs) of PS1-12cil and SN 2012aa which are superluminous and luminous supernovae (SNe), respectively. The multi-band and bolometric LCs of these two SNe show unusual secondary bumps after the main peaks. The two-peaked LCs cannot be explained by any simple energy-source models (e.g., the $^{56}$Ni cascade decay model, the magnetar spin-down model, and the ejecta-circumstellar medium interaction model). Therefore, we employ the $^{56}$Ni plus ejecta-circumstellar medium (CSM) interaction (CSI) model, the magnetar plus CSI model, and the double CSI model to fit their bolometric LCs, and find that both these two SNe can be explained by the double CSI model and the magnetar plus CSI model. Based on the modeling, we calculate the the time when the shells were expelled by the progenitors: provided that they were powered by double ejecta-shell CSI, the inner and outer shells might be expelled $\sim 0.2-3.6$ and $\sim 2-25$ years before the explosions of the SNe, respectively; the shells were expelled $\sim 2-20$ years before the explosions of the SNe if they were powered by magnetars plus CSI.

The isotopes $^{60}$Fe and $^{26}$Al originate from massive stars and their supernovae, reflecting ongoing nucleosynthesis in the Galaxy. We studied the gamma-ray emission from these isotopes at characteristic energies 1173, 1332, and 1809 keV with over 15 years of SPI data, finding a line flux in $^{60}$Fe combined lines of $(0.31\pm 0.06) \times 10^{-3}$ ph cm$^{-2}$ s$^{-1}$ and the $^{26}$Al line flux of $(16.8\pm 0.7) \times 10^{-4}$ ph cm$^{-2}$ s$^{-1}$ above the background and continuum emission for the whole sky. Based on the exponential-disk grid maps, we characterise the emission extent of $^{26}$Al to find scale parameters $R_0 =7.0^{+1.5}_{-1.0}$ kpc and $z_0=0.8^{+0.3}_{-0.2}$ kpc, however the $^{60}$Fe lines are too weak to spatially constrain the emission. Based on a point source model test across the Galactic plane, the $^{60}$Fe emission would not be consistent with a single strong point source in the Galactic center or somewhere else, providing a hint for a diffuse nature. We carried out comparisons of emission morphology maps using different candidate-source tracers for both $^{26}$Al and $^{60}$Fe emissions, and suggests that the $^{60}$Fe emission is more likely to be concentrated towards the Galactic plane. We determine the $^{60}$Fe/$^{26}$Al $\gamma$-ray flux ratio at $(18.4\pm4.2)\,\%$ , when using a parameterized spatial morphology model. Across the range of plausible morphologies, it appears possible that $^{26}$Al and $^{60}$Fe are distributed differently in the Galaxy. Using the best fitting maps for each of the elements, we constrain flux ratios in the range 0.2--0.4. We discuss its implications for massive star models and their nucleosynthesis.

The properties of fast radio bursts (FRBs) indicate that the physical origin of this type of astrophysical phenomenon is related to neutron stars. The first detected repeating source, FRB 121102, is associated with a persistent radio counterpart. In this paper, we propose that this radio counterpart could arise from a pulsar wind nebula powered by a magnetar without surrounding supernova ejecta. Its medium is a stratified structure produced by a progenitor wind. The model parameters are constrained by the spectrum of the counterpart emission, the size of the nebula, and the large but decreasing rotation measure (RM) of the repeating bursts. In addition, the observed dispersion measure is consistent with the assumption that all of the RM comes from the shocked medium.

Two accretion columns have been argued to form over the surface of a newborn millisecond magnetar for an extremely high accretion rate $\gtrsim1.8\times10^{-2}M_\odot\ {\rm s^{-1}}$ that may occur in the core-collapse of a massive star. In this paper, we investigate the characteristics of these accretion columns and their gravitational wave (GW) radiation. For a typical millisecond magnetar (surface magnetic field strength $B\sim10^{15}$ G and initial spin period $P\sim1$ ms), we find (1) its accretion columns are cooled via neutrinos and can reach a height $\sim1$ km over the stellar surface; (2) its column-induced characteristic GW strain is comparable to the sensitivities of the next generation ground-based GW detectors within a horizon $\sim1$ Mpc; (3) the magnetar can survive only a few tens of seconds; (4) during the survival timescale, the height of the accretion columns increases rapidly to the peak and subsequently decreases slowly; (5) the column mass, characteristic GW strain, and maximum GW luminosity have simultaneous peaks in a similar rise-fall evolution. In addition, we find that the magnetar's spin evolution is dominated by the column accretion torque. A possible association with failed supernova is also discussed.

We report our photometric and spectroscopic follow-up observations of the superoutbursts of three dwarf novae (GWAC\,180415A, GWAC\,181017A and GWAC\,181211A) identified independently by the Ground Wide-angle Cameras system, one of the ground-based instruments of the China-France SVOM mission. Based on a combination of our photometry and that taken from the AAVSO, our period analysis of the superhumps enables us to determine the mass ratios to be 0.0967-0.1163, 0.1879-0.1883 and 0.0981-0.1173 for GWAC\,180415A, GWAC\,181017A and GWAC\,181211A, respectively. GWAC\,180415A can be firmly identified as a WZ sge-type dwarf novae due to its long duration ($\sim2$ weeks) multiple rebrightenings with amplitudes of 3-4 magnitudes, the early superhump associated with a double-wave modulation and the low mass ratio. The inferred low mass ratio and location in the $\varepsilon-P_{\mathrm{orb}}$ diagram suggest that GWAC\,181211A is a WZ sge-type dwarf novae candidate. The measured Balmer decrements suggest the Balmer line emission is produced from an optical thick region in GWAC\,180415A and GWAC\,181017A, and from an optical thin region in GWAC\,181211A.

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.

It was recently proposed that a significant fraction of ultraluminous X-ray sources (ULXs) actually host a neutron star (NS) accretor. We have performed a systematic study on the NS ULX population in Milky Way-like galaxies, by combining binary population synthesis and detailed stellar evolution calculations. Besides a normal star, the ULX donor can be a helium star (the hydrogen envelope of its progenitor star was stripped during previous common envelope evolution) if the NS is accreting at a super-Eddington rate via Roche lobe overflow. We find that the NS$-$helium star binaries can significantly contribute the ULX population, with the overall number of about several in a Milky Way-like galaxy. Our calculations show that such ULXs are generally close systems with orbital period distribution peaked at $ \sim 0.1 $ day (with a tail up to $ \sim100 $ days), and the helium stars have relatively low masses distributing with a maximum probability at $ \sim 1M_{\odot} $.

It is well-known that ordinary supernovae (SNe) are powered by 56Ni cascade decay. Broad-lined type Ic SNe (SNe Ic-BL) are a subclass of SNe that are not all exclusively powered by 56Ni decay. It was suggested that some SNe Ic-BL are powered by magnetar spin-down. iPTF16asu is a peculiar broad-lined type Ic supernova discovered by the intermediate Palomar Transient Factory. With a rest-frame rise time of only 4 days, iPTF16asu challenges the existing popular models, for example, the radioactive heating (56Ni-only) and the magnetar+56Ni models. Here we show that this rapid rise could be attributed to interaction between the SN ejecta and a pre-existing circumstellar medium ejected by the progenitor during its final stages of evolution, while the late-time light curve can be better explained by energy input from a rapidly spinning magnetar. This model is a natural extension to the previous magnetar model. The mass-loss rate of the progenitor and ejecta mass are consistent with a progenitor that experienced a common envelope evolution in a binary. An alternative model for the early rapid rise of the light curve is the cooling of a shock propagating into an extended envelope of the progenitor. It is difficult at this stage to tell which model (interaction+magnetar+56Ni or cooling+magnetar+56Ni) is better for iPTF16asu. However, it is worth noting that the inferred envelope mass in the cooling+magnetar+56Ni is very high.

We extract 18 candidate short gamma-ray bursts (SGRBs) with precursors from 660 SGRBs observed by {\em Fermi} and {\em Swift} satellites, and carry out a comprehensive analysis on their temporal and spectral features. We obtain the following results: (1) For a large fraction of candidates, the main burst durations are longer than their precursor durations, comparable to their quiescent times from the end of precursors to the beginning of their main bursts. (2) The average flux of precursors tends to increase as their main bursts brighten. (3) As seen from the distributions of hardness ratio and spectral fitting, the precursors are slightly spectrally softer with respect to the main bursts. Moreover, a significant portion of precursors and all main bursts favor a non-thermal spectrum. (4) The precursors might be a probe of the progenitor properties of SGRBs such as the magnetic field strength and the crustal equation of state if they arise from some processes before mergers of binary compact objects rather than post-merger processes.

For a binary composed of a spinning black hole (BH) (with mass $\gtrsim 7M_\odot$) and a strongly magnetized neutron star (NS) (with surface magnetic field strength $B_{\rm S,NS}\gtrsim10^{12}$\,G and mass $\sim 1.4M_\odot$), the NS as a whole will possibly eventually plunge into the BH. During the inspiral phase, the spinning BH could be charged to the Wald charge quantity $Q_{\rm W}$ until merger in an electro-vacuum approximation. During the merger, if the spinning charged BH creates its own magnetosphere due to an electric field strong enough for pair cascades to spark, the charged BH would transit from electro-vacuum to force-free cases and could discharge in a time $\gtrsim1~{\rm ms}$. As the force-free magnetosphere is full of a highly conducting plasma, the magnetic flux over the NS's caps would be retained outside the BH's event horizon under the frozen-in condition. Based on this scenario, we here investigate three possible energy-dissipation mechanisms that could produce electromagnetic (EM) counterparts in a time interval of the BH's discharge post a BH-NS merger-induced gravitational wave event: (1) magnetic reconnection at the BH's poles would occur, leading to a millisecond bright EM signal, (2) a magnetic shock in the zone of closed magnetic field lines due to the detachment and reconnection of the entire BH magnetic field would probably produce a bright radio emission, e.g., a fast radio burst, and (3) the Blandford-Znajek mechanism would extract the BH's rotational energy, giving rise to a millisecond-duration luminous high-energy burst. We also calculate the luminosities due to these mechanisms as a function of BH's spin for different values of $B_{\rm S,NS}$.

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.

Very recently \citet{XueYQ2019} reported an important detection of the X-ray transient, CDF-S XT2, whose light curve is analogous to X-ray plateau features of gamma-ray burst afterglows. They suggested that this transient is powered by a remnant stable magnetar from a binary neutron star merger since several pieces of evidence (host galaxy, location, and event rate) all point toward such an assumption. In this paper, we revisit this scenario and confirm that this X-ray emission can be well explained by the internal gradual magnetic dissipation process in an ultra-relativistic wind of the newborn magnetar. We show that both the light curve and spectral evolution of CDF-S XT2 can be well fitted by such a model. Furthermore, we can probe some key properties of the central magnetar, such as its initial spin period, surface magnetic field strength and wind saturation Lorentz factor.

Plateaus are common in X-ray afterglows of gamma-ray bursts. Among a few scenarios for the origin of them, the leading one is that there exists a magnetar inside and persistently injects its spin-down energy into an afterglow. In previous studies, the radiation efficiency of this process is assumed to be a constant $\gtrsim0.1$, which is quite simple and strong. In this work we obtain the efficiency from a physical point of view and find that this efficiency strongly depends on the injected luminosity. One implication of this result is that those X-ray afterglow light curves which show steeper temporal decay than $t^{-2}$ after the plateau phase can be naturally understood now. Also, the braking indexes deduced from afterglow fitting are found to be larger than those in previous studies, which are more reasonable for newborn magnetars.

A central compact object (CCO, e.g. a black hole) with an accretion disk has been suggested as the common central engine of various astrophysical phenomena, such as gamma-ray bursts (GRBs), tidal disruption events (TDEs) and active galactic nuclei (AGNs). A jet powered by such a system might precess due to the misalignment of the angular momenta of the CCO and accretion disk. Some quasi-periodic behaviors observed in the light curves of these phenomena can be well interpreted within the framework of a precessing jet model. In this paper, we study the emission polarization of precessing jets in the three kinds of phenomena. The polarization angle also shows a gradual change for the synchrotron emission in both the random and toroidal magnetic field configurations with the precessing jet, while it can only change abruptly by $90^\circ$ for the non-precessing top-hat jet. Polarization properties are periodic due to the assumptions made in our model. The polarization observations are crucial to confirm the precession nature of jets in GRBs, TDEs and AGNs.

The Large High Altitude Air Shower Observatory (LHAASO) project is a new generation multi-component instrument, to be built at 4410 meters of altitude in the Sichuan province of China, with the aim to study with unprecedented sensitivity the spec trum, the composition and the anisotropy of cosmic rays in the energy range between 10$^{12}$ and 10$^{18}$ eV, as well as to act simultaneously as a wide aperture (one stereoradiant), continuously-operated gamma ray telescope in the energy range between 10$^{11}$ and $10^{15}$ eV. The experiment will be able of continuously surveying the TeV sky for steady and transient sources from 100 GeV to 1 PeV, t hus opening for the first time the 100-1000 TeV range to the direct observations of the high energy cosmic ray sources. In addition, the different observables (electronic, muonic and Cherenkov/fluorescence components) that will be measured in LHAASO will allow to investigate origin, acceleration and propagation of the radiation through a measurement of energy spec trum, elemental composition and anisotropy with unprecedented resolution. The remarkable sensitivity of LHAASO in cosmic rays physics and gamma astronomy would play a key-role in the comprehensive general program to explore the High Energy Universe. LHAASO will allow important studies of fundamental physics (such as indirect dark matter search, Lorentz invariance violation, quantum gravity) and solar and heliospheric physics. In this document we introduce the concept of LHAASO and the main science goals, providing an overview of the project.

Over the past decade, many rapidly evolving optical transients (REOTs), whose rise and decline timescales are significantly shorter than those of canonical supernovae (SNe), have been discovered and studied. Some REOTs have high peak luminosity ($\gtrsim 10^{43}$ erg s$^{-1}$), disfavoring the radioactivity-powered-SN model that has been widely adopted to explain normal SNe. In this paper, we study three luminous REOTs (PS1-10bjp, PS1-11bbq, and PS1-13ess) and use a model involving magnetar energy input to fit their bolometric light curves and temperature evolution. We find that core-collapse SNe (CCSNe) powered by magnetars with $P_0 \approx 18$--34 ms and $B_p \approx$ (2.5--5.8) $\times10^{15}$ G can reproduce their bolometric light curves as well as the temperature evolution. The inferred values of ejecta mass are $\sim 0.40$--0.46 M$_\odot$, slightly smaller than that of the Type Ic SN 1994I, indicating that they can be explained by the magnetar-powered stripped CCSN model.

SN~2007D is a nearby (redshift $z = 0.023146$), luminous Type Ic supernova (SN) having a narrow light curve (LC) and high peak luminosity. Previous research based on the assumption that it was powered by the $^{56}$Ni cascade decay suggested that the inferred $^{56}$Ni mass and the ejecta mass are $\sim 1.5$M$_{\odot}$ and $\sim 3.5$M$_{\odot}$, respectively. In this paper, we employ some multiband LC models to model the $R$-band LC and the color ($V-R$) evolution of SN~2007D to investigate the possible energy sources powering them. We find that the pure $^{56}$Ni model is disfavored; the multiband LCs of SN~2007D can be reproduced by a magnetar whose initial rotational period $P_{0}$ and magnetic field strength $B_p$ are $7.28_{-0.21}^{+0.21}$ (or $9.00_{-0.42}^{+0.32}$) ms and $3.10_{-0.35}^{+0.36}\times 10^{14}$ (or $2.81_{-0.44}^{+0.43}\times 10^{14}$) G, respectively. By comparing the spectrum of SN~2007D with that of some superluminous SNe (SLSNe), we find that it might be a luminous SN like several luminous ``gap-filler" optical transients that bridge ordinary and SLSNe, rather than a genuine SLSN.

The formation of neutron stars (NSs), both from collapses of massive stars and mergers of compact objects, can be usually indicated by bright transients emitted from explosively-ejected material. In particular, if the newborn NSs can rotate at a millisecond period and have a sufficiently high magnetic field, then the spin-down of the NSs would provide a remarkable amount of energy to the emitting material. As a result, super-luminous supernovae could be produced in the massive stellar collapse cases, while some unusual fast evolving and luminous optical transients could arise from the cases of NS mergers and accretion-induced collapses of white dwarfs. In all cases, if the dipolar magnetic fields of the newborn NSs can be amplified to be as high as $10^{15}$ G, a relativistic jet could be launched and then a gamma-ray burst can be produced as the jet successfully breaks out from the surrounding nearly-isotropic ejected material.

Fast radio bursts (FRBs) are bright milliseconds radio transients with large dispersion measures (DMs), and can be used as potential astrophysical and cosmological tools. Many models for FRB progenitors have been proposed. However, none of them can explain all the observational properties, meaning that FRBs remain one of the most intriguing mysteries in astronomy. The statistical properties of radio bursts can unveil the underlying physics. Here, we report statistical results of the repeating FRB 121102, and show that FRB 121102 and solar type III radio bursts share four statistical properties: power-law frequency distributions for energies, fluxes, durations and waiting times. All of the distributions can be explained by avalanche models of self-organized criticality (SOC) systems. It is well known that solar type III radio burst arises from the nonlinear conversion of Langmuir waves generated by two-stream instability of electron beams, which are accelerated by bursty magnetic reconnections. The similarities support that repeating FRBs are coherent emissions from magnetic-reconnection-driven beams in the magnetospheres of magnetars.

GRB 170817A, detected by Fermi-GBM 1.7\,s after the merger of a neutron star (NS) binary, provides the first direct evidence for a link between such a merger and a short-duration gamma-ray burst. The X-ray observations after GRB 170817A indicate a possible X-ray flare with a peak luminosity $L_{\rm peak} \sim 2\times 10^{39}\,{\rm erg\,s}^{-1}$ near day 156. Here we show that this X-ray flare may be understood based on a slim disc around a compact object. On the one hand, there exists the maximal accretion rate $\dot M_{\rm max}$ for the slim disc, above which an optically thick outflow is significant and radiation from the disc is obscured. Based on the energy balance analysis, we find that $\dot M_{\rm max}$ is in the range of $\sim 4\dot M_{\rm Edd}$ and $\sim 21\dot M_{\rm Edd}$ when the angular velocity of the slim disc is between $\rm (1/5)^{1/2}\Omega_K$ and $\rm \Omega_K$ (where $\dot M_{\rm Edd}$ is the Eddington accretion rate and $\Omega_K$ is the Keplerian angular velocity). With $\dot M_{\rm max}$, the slim disc can provide a luminosity $\sim L_{\rm peak}$ for a compact object of $2.5 M_{\sun}$. On the other hand, if the merger of two NSs forms a typical neutrino-dominated accretion disc whose accretion rate $\dot M$ follows a power-law decline with an index $-1.8$ , then the system must pass through the outflow regime and enter the slim disc in $\sim 11-355$ days. These results imply that a post-merger slim accretion disc could account for the observed late-time $L_{\rm peak}$.

The mergers of black hole (BH)--neutron star (NS) binaries have been one of the most interesting topics in astrophysics, because such events have been thought to possibly produce multimessenger signals including gravitational waves and broadband electromagnetic (EM) waves. In this paper, we investigate EM emission from the inspiral of a binary composed of a spinning BH and a magnetized NS. Observationally, the BH is usually more massive than $\sim7M_\odot$ and the NS has a mass $\simeq 1.4M_\odot$. During the inspiral of such a binary, the BH will accumulate more and more charges based on the charging scenario of Wald, even though the BH will eventually swallow the NS whole inevitably. We calculate the emission luminosities and energies through three energy dissipation mechanisms: magnetic dipole radiation, electric dipole radiation, and magnetic reconnection. We show that magnetic dipole radiation due to the spin of the increasingly charged BH and magnetic reconnection in between the BH and the NS could be most significant at the final inspiral stage. We find that if the BH is rapidly rotating and the NS is strongly magnetized, these mechanisms would lead to a detectable EM signal (e.g., a short-duration X-ray transient).

Observations of short-duration gamma-ray bursts and their afterglows show that a good fraction (perhaps $\gtrsim50\%$) of binary neutron star mergers lead to strongly magnetized, rapidly rotating pulsars (including millisecond magnetars), no matter whether the pulsar remnants are short- or long-lived. Such compact objects are very likely to have significant radial oscillations and high interior temperatures, as indicated in recent numerical simulations. In this paper, we have investigated rotation-induced gravitational radiation from possibly existing, radially oscillating pulsars after binary neutron star mergers, and find that this mechanism can efficiently damp the radial oscillations. The resulting gravitational waves (GWs) could have a non-negligible contribution to the high-frequency spectrum. We provide an order-of-magnitude estimate of the event rate and suggest that such GW events would be detectable with the advanced LIGO/Virgo or next-generation detectors. Our discussion can also be applied to newborn, radially oscillating, millisecond pulsars formed through the other astrophysical processes.

Supernovae (SNe) are the most brilliant optical stellar-class explosions. Over the past two decades, several optical transient survey projects discovered more than $\sim 100$ so-called superluminous supernovae (SLSNe) whose peak luminosities and radiated energy are $\gtrsim 7\times 10^{43}$ erg s$^{-1}$ and $\gtrsim 10^{51}$ erg, at least an order of magnitude larger than that of normal SNe. According to their optical spectra features, SLSNe have been split into two broad categories of type I that are hydrogen-deficient and type II that are hydrogen-rich. Investigating and determining the energy sources of SLSNe would be of outstanding importance for understanding the stellar evolution and explosion mechanisms. The energy sources of SLSNe can be determined by analyzing their light curves (LCs) and spectra. The most prevailing models accounting for the SLSN LCs are the $^{56}$Ni cascade decay model, the magnetar spin-down model, the ejecta-CSM interaction model, and the jet-ejecta interaction model. In this \textit{review}, we present several energy-source models and their different combinations.

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

Gamma-ray bursts are the strongest explosions in the Universe since the Big Bang, believed to be produced either in forming black holes at the end of massive star evolution or merging of compact objects. Spectral and timing properties of gamma-ray bursts suggest that the observed bright gamma-rays are produced in the most relativistic jets in the Universe; however, the physical properties, especially the structure and magnetic topologies in the jets are still not well known, despite several decades of studies. It is widely believed that precise measurements of the polarization properties of gamma-ray bursts should provide crucial information on the highly relativistic jets. As a result there have been many reports of gamma-ray burst polarization measurements with diverse results, see, however many such measurements suffered from substantial uncertainties, mostly systematic. After the first successful measurements by the GAP and COSI instruments, here we report a statistically meaningful sample of precise polarization measurements, obtained with the dedicated gamma-ray burst polarimeter, POLAR onboard China's Tiangong-2 spacelab. Our results suggest that the gamma-ray emission is at most polarized at a level lower than some popular models have predicted; although our results also show intrapulse evolution of the polarization angle. This indicates that the low polarization degrees could be due to an evolving polarization angle during a gamma-ray burst.

The central engine of gamma-ray bursts (GRBs) is poorly constrained. There exist two main candidates: a fast-rotating black hole and a rapidly spinning magnetar. Furthermore, the X-ray plateaus are widely accepted by the energy injection into the external shock. In this paper, we systematically analyze the \emph{Swift}/XRT light curves of 101 GRBs having plateau phases and known redshifts (before May 2017). Since a maximum energy budget ($\sim2\times10^{52}$ erg) exists for magnetars but not for a black hole, this provides a good clue to identify the type of GRB central engine. We calculate the isotropic kinetic energy $E_{\rm K,iso}$ and the isotropic X-ray energy release $E_{\rm X,iso}$ for individual GRB. We identify three categories based on how likely a black hole harbor at central engine: 'Gold' (9 out of 101, both $E_{\rm X,iso}$ and $E_{\rm K,iso}$ exceed the energy budget), 'Silver' (69 out of 101, $E_{\rm X,iso}$ less than the limit but $E_{\rm K,iso}$ is greater than the limit), and 'Bronze' (23 out of 101, the energies are not above the limit). We then derive and test the black hole parameters with the Blandford Znajek mechanism, and find that the observations of the black hole candidate ('Gold'+'Silver') samples are consistent with the expectations of the black hole model. Furthermore, we also test the magnetar candidate ('Bronze') sample with the magnetar model, and find that the magnetar surface magnetic led ($B_{p}$) and initial spin period ($P_{0}$) are consistent with the expectations of the magnetar model. Our analysis indicates that, if the magnetar wind is isotropic, a magnetar central engine is possible for 20\% of the analyzed GRBs. For most GRBs a black hole is most likely operating.

In this White Paper we present the potential of the enhanced X-ray Timing and Polarimetry (eXTP) mission for studies related to Observatory Science targets. These include flaring stars, supernova remnants, accreting white dwarfs, low and high mass X-ray binaries, radio quiet and radio loud active galactic nuclei, tidal disruption events, and gamma-ray bursts. eXTP will be excellently suited to study one common aspect of these objects: their often transient nature. Developed by an international Consortium led by the Institute of High Energy Physics of the Chinese Academy of Science, the eXTP mission is expected to be launched in the mid 2020s.

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

We investigate emission signatures of binary compact star gravitational wave sources consisting of strongly magnetized neutron stars (NSs) and/or white dwarfs (WDs) in their late-time inspiral phase. Because of electromagnetic interactions between the magnetospheres of the two compact stars, a substantial amount of energy will be extracted, and the resultant power is expected to be $\sim 10^{38} - 10^{44}$ erg/s in the last few seconds before the two stars merge, when the binary system contains a NS with a surface magnetic field $10^{12}$ G. The induced electric field in the process can accelerate charged particles up to the EeV energy range. Synchrotron radiation is emitted from energetic electrons, with radiative energies reaching the GeV energy for binary NSs and the MeV energy for NS - WD or double WD binaries. In addition, a blackbody component is also presented and it peaks at several to hundreds keV for binary NSs and at several keV for NS - WD or double WD binaries. The strong angular dependence of the synchrotron radiation and the isotropic nature of the blackbody radiation lead to distinguishable modulation patterns between the two emission components. If coherent curvature radiation is presented, fast radio bursts could be produced. These components provide unique simultaneous electromagnetic signatures as precursors of gravitational wave events associated with magnetized compact star mergers and short gamma ray bursts (e.g., GRB 100717).

Recent wide-field surveys discovered new types of peculiar optical transients that showed diverse behaviors of the evolution of photospheric properties. We develop a general theory of homologous explosions with constant opacity, paying special attention on the evolution of the photospheric radius $R_{\rm ph}$. We find that regardless of the density distribution profile, $R_{\rm ph}$ always increases early on and decreases at late times. This result does not depend on the radiation and cooling processes inside the ejecta.The general rising/falling behavior of $R_{\rm ph}$ can be used to quickly diagnose whether the source originates from a supernova-like explosion. The shape of the $R_{\rm ph}$ evolution curve depends on the density profile, so the observations may be used to directly diagnose the density profile as well as the temperature profile of the ejecta. All the well-monitored supernovae show such a $R_{\rm ph}$ rising/falling behavior, which is consistent with our theory. The recently discovered peculiar transient AT2018cow showed a continuous decay of $R_{\rm ph}$, which is disfavored to be of a supernova-like explosion origin. Our result therefore supports the interpretation of this transient as a tidal disruption event.

The Intermediate Palomar Transient Factory reported the discovery of an unusual type II-P supernova iPTF14hls. Instead of a ~100-day plateau as observed for ordinary type II-P supernovae, the light curve of iPTF14hls has at least five distinct peaks, followed by a steep decline at ~1000 days since discovery. Until 500 days since discovery, the effective temperature of iPTF14hls is roughly constant at 5000-6000K . In this paper we propose that iPTF14hls is likely powered by intermittent fallback accretion. It is found that the light curve of iPTF14hls can be well fit by the usual t^{-5/3} accretion law until ~1000 days post discovery when the light curve transitions to a steep decline. To account for this steep decline, we suggest a power-law density profile for the late accreted material, rather than the constant profile as appropriated for the t^{-5/3} accretion law. Detailed modeling indicates that the total fallback mass is ~0.2M_{sun}, with an ejecta mass M_{ej}~21M_{sun}. We find the third peak of the light curve cannot be well fit by the fallback model, indicating that there could be some extra rapid energy injection. We suggest that this extra energy injection may be a result of a magnetic outburst if the central object is a neutron star. These results indicate that the progenitor of iPTF14hls could be a massive red supergiant.

In this paper, we revisit the scenario that an internal gradual magnetic dissipation takes place within the wind from a newborn millisecond magnetar can be responsible for gamma-ray burst production. We show that a combination of two emission components in this model, i.e., the photospheric emission from the wind and the synchrotron radiation within the magnetic reconnection region, can give a reasonable fit to the observed spectrum of the prompt emission phase of GRB 160804A. We obtain the physical parameters through a Monte Carlo procedure and deduce the initial spin period and magnetic field of the central magnetar. Furthermore, the independent afterglow fitting analysis gives a consistent result, adding great credibility to this scenario. In addition, we predict a subclass of GRBs called bursts from such a Magnetar wind Internal Gradual MAgnetic Dissipation (abbreviated as "MIGMAD bursts") that have several distinctive properties.

The gravitational wave event GW170817 associated with the short gamma-ray burst (GRB) 170817A confirms that binary neutron star (BNS) mergers are one of the origins of short GRBs. The associated kilonova emission, radioactively powered by nucleosynthesized heavy elements, was also detected. Black hole-neutron star (BH-NS) mergers have been argued to be another promising origin candidate of short GRBs and kilonovae. Numerical simulations show that the ejecta in BH-NS mergers is geometrically much more anisotropic than the BNS merger case. In this paper, we investigate observational signatures of kilonova emission from the anisotropic ejecta in BH-NS mergers. We find that a bump appears on the bolometric luminosity light curve due to the inhomogeneous mass distribution in the latitudinal direction. The decay slope of the single-band light curve becomes flatter and the spectrum also deviates from a single-temperature blackbody radiation spectrum due to the gradient in the velocity distribution of the ejecta. Future detection or non-detection of such signatures would be useful to test the mass ejection geometry in BH-NS mergers.

Recently, a high-energy neutrino event IceCube-170922A in the spatial and temporal coincidence with the flaring gamma-ray blazar TXS 0506+056 was reported. A neutrino outburst between September 2014 and March 2015 was discovered in the same direction by a further investigation of $9.5$ years of IceCube data, while the blazar is in a quiescent state during the outburst with a gamma-ray flux only about one-fifth of the neutrino flux. In this letter, we propose the neutrino outburst originates from the interaction between a relativistic jet and a dense gas cloud which may be formed via the tidally disrupted envelope of a red giant being blown by the impact of the jet. Gamma-ray photons and electron/positron pairs that are produced correspondingly will induce electromagnetic cascades. Comptonization of the cascade emission inside the cloud forms an X-ray photon field with Wien distribution. GeV flux is suppressed due to the absorption by the Comptonized photon field and, as a result, a hard spectrum above 10 GeV is formed. The gamma-ray spectrum predicted in our model is consistent with the Fermi-LAT data of TXS 0506+056.

We used the light curve code XRBinary to model the quiescent K2 light curves of three low-inclination cataclysmic variables (CVs): 1RXS\,J0632+2536 (J0632+2536), RZ\,Leo, TW\,Vir and the pre-CV WD\,1144+011. Optimized light curve models were obtained using a nonlinear fitting code NMfit and visualized by Phoebe 2.0. The disk model of J0632+2536 shows that one hotspot at the edge of the disk is enough to describe its light curve, while the other two dwarf nova (DN): RZ\,Leo and TW\,Vir require two hotspots. A typical pre-CV model with a weak irradiation effect for WD\,1144+011 can explain its single-hump modulation, and the newly observed spectrum confirms its previous classification. The synthetic analyses for the DN clearly indicate that phase zero of the double-hump modulations occurs around the secondary minimum and the primary hump is mainly caused by the hotspot at the edge of the disk. The quiescent disk has a flat temperature distribution with a power index of $\sim0.11$. The disk model of RZ\,Leo implies a truncated disk, supporting its previously speculated classification as an intermediate polar (IP). Except for the IP model of RZ\,Leo, which lacks a component related to the inferred accretion curtain, the models of J0632+2536, TW\,Vir and WD\,1144+011 are consistent with results from the Gaia mission. The derived masses and radii of the secondaries of the three DN are consistent with the semi-empirical relations for CV donor stars, while their effective temperatures are higher than the predictions. Irradiation of the donor stars is investigated to explain this discrepancy.

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

Multi-band observations of the fast radio burst (FRB) 131104 show that this burst may be associated with a gamma-ray transient entitled Swift J0644.5-5111. Follow-up observations for potential X-ray and radio counterparts of FRB 131104/Swift J0644.5-5111 got null results and provided the upper limits of the emission flux at 5.5 GHz, 7.5 GHz, $U$-band, and X-ray band. By assuming this association and using these upper limits, environmental properties (the fraction of energyin a magnetic field $\varepsilon_{\rm B}$ and the number density $n$) of the progenitor system of FRB 131104/Swift J0644.5-5111 were constrained in the context of the standard afterglow model that neglects the non-relativistic effect and jet effect by several groups. In this paper, we adopt a unified afterglow model that takes into account the non-relativistic effect and jet effect and use the upper limits of four bands (5.5 GHz, 7.5 GHz, $U$-band, and X-ray) to obtain more stringent constraints on the parameter space spanned by $\varepsilon_{\rm B}$ and $n$. We thus suggest that FRB 131104/Swift J0644.5-5111 might originate from a black hole-neutron star merger event. Moreover, we calculate multi-band emissions from a kilonova powered by the radioactivity of $r$-process elements synthesized in the ejected neutron-rich material and find that the $U$-band emission from the putative kilonova is significantly lower than the upper limit of the observations.

In a strong magnetic field, a high-energy photon can be absorbed and then produce an electron-positron pair. The produced electron/positron will in turn radiate a high-energy photon via synchrotron radiation, which then initiates a cascade. We built a one-dimensional Monte-Carlo code to study the development of the cascade especially after it reaches the saturated status, when almost all the energy of the primary particles transfers to the photons. The photon spectrum in this status has a cut-off due to the absorption by magnetic fields, which is much sharper than the exponential one. Below the cut-off, the spectral energy distribution (SED) manifest itself as a broken power-law with a spectral index of $0.5$ and $0.125$, respectively, below and above the broken energy. The SED can be fitted by a simple analytical function, which is solely determined by the product of the cascade scale $R$ and the magnetic field perpendicular to the motion of the particle B_{\perp}, with an accuracy better than 96\%. The similarity of the spectrum to that from the cascade in an isotropic black-body photon field is also studied.

The successful joint observation of the gravitational wave event GW170817 and its multi-wavelength electromagnetic counterparts first enables human to witness a definite merger event of two neutron stars (NSs). This historical event confirms the origin of short-duration gamma-ray bursts (GRBs), and in particular, identifies the theoretically-predicted kilonova phenomenon that is powered by radioactive decays of $r$-process heavy elements. However, whether a long-lived remnant NS could be formed during this merger event remains unknown, although such a central engine has been suggested by afterglow observations of some short-duration GRBs. By invoking this long-lived remnant NS, we here propose a model of hybrid energy sources for the kilonova AT2017gfo associated with GW 170817. While the early emission of AT2017gfo is still powered radioactively as usually suggested, its late emission is primarily caused by delayed energy injection from the remnant NS. In our model, only one single opacity is required and an intermediate value of $\kappa\simeq0.97\,\rm cm^2g^{-1}$ is revealed, which could be naturally provided by lanthanide-rich ejecta that is deeply ionized by the emission from a wind of the NS. These self-consistent results indicate that a long-lived remnant NS, which must own a very stiff equation of state, had been formed during the merger event of GW170817. This provides a very stringent constraint on the strong interaction in nuclear-quark matter. It is further implied that such GW events could provide a probe of the early spin and magnetic evolutions of NSs, e.g., the burying of surface magnetic fields.

Gamma-ray bursts (GRBs) have been widely argued to originate from binary compact object mergers or core collapses of massive stars. Jets from these systems may have two components, an inner, narrow sub-jet and an outer, wider sub-jet. Such a jet subsequently interacts with its ambient gas, leading to a reverse shock (RS) and a forward shock (FS). The magnetic field in the narrow sub-jet is very likely to be mixed by an ordered component and a random component during the afterglow phase. In this paper, we calculate light curves and polarization evolution of optical afterglows with this mixed magnetic field in the RS region of the narrow sub-jet in a two-component jet model. The resultant light curve has two peaks: an early peak arises from the narrow sub-jet and a late-time rebrightening is due to the wider sub-jet. We find the polarization degree (PD) evolution under such a mixed magnetic field confined in the shock plane is very similar to that under the purely ordered magnetic field condition. The two-dimensional `mixed' magnetic fields confined in the shock plane are essentially the ordered magnetic fields only with different configurations. The position angle (PA) of the two-component jet can change gradually or abruptly by $90^\circ$. In particular, an abrupt $90^\circ$ change of the PA occurs when the PD changes from its decline phase to rise phase.

CRTS J035905.9+175034 is the first eclipsing SU UMa system for which a superoutburst has been observed by Kepler in the short-cadence mode. The light curve contains one superoutburst, eight normal outbursts (including a precursor to the superoutburst), and several minioutbursts that are present before -- but not after -- the superoutburst. The superoutburst began with a precursor normal outburst, and shortly after the peak of the precursor, the system developed large-amplitude superhumps that achieved their maximum amplitude after just three superhump cycles. The period excess of the initial superhump period relative to the orbital period implies a mass ratio of 0.281 +/- 0.015, placing it marginally above most theoretical predictions of the highest-possible mass ratio for superhump formation. In addition, our analysis of the variations in eclipse width and depth, as well as the hot spot amplitudes, generally provides substantiation of the thermal-tidal instability model. The K2 data, in conjunction with our ground-based time-resolved spectroscopy and photometry from 2014-2016, allows us to determine many of the fundamental parameters of this system.

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

The high-energy (> 100MeV) emission observed by Fermi-LAT during the prompt phase of some luminous gamma-ray bursts (GRBs) could arise from the cascade induced by interactions between accelerated protons and the radiation field of GRBs. The photomeson process, which is usually suggested to operate in such a hadronic explanation, requires a rather high proton energy (> 10^17eV) for an efficient interaction. However, whether GRBs can accelerate protons to such a high energy is far from guaranteed, although they have been suggested as the candidate source for ultrahigh-energy cosmic rays. In this work, we revisit the hadronic model for the prompt high-energy emission of GRBs with a smaller maximum proton energy than the usually adopted value estimated from the Bohm condition. In this case, the Bethe-Heitler pair production process becomes comparably important or even dominates over the photomeson process. We show that with a relatively low maximum proton energy with a Lorentz factor of 10^5 in the comoving frame, the cascade emission can still reproduce various types of high-energy spectrum of GRBs. For most GRBs without high-energy emission detected, the maximum proton energy could be even lower and relax the constraints on the parameters of GRB jet resulting from the fact of non-detection of GRB neutrinos by IceCube.

A subset of short gamma-ray bursts (sGRBs) have been found to be characterized by near-infrared/optical bumps at $\sim 1$ days, some of which exhibit almost concurrent X-ray flares. Although the near-infrared/optical bumps may be a signature of kilonovae, the X-ray flares are not consistent with kilonovae. It is widely believed that sGRBs are produced by the mergers of double compact objects, during which sub-relativistic ejecta are launched. In this paper we propose that the above optical/X-ray features are indicative of the formation of long-lived magnetars following the mergers of double neutron stars. Observations and theoretical works imply that the spin-down power of the magnetars is injected into the ejecta as ultra-relativistic electron-positron pairs, i.e. pulsar wind nebulae (PWNe). Here we suggest such a PWN model and find that the optical bump and X-ray flare observed in GRB 050724 can be well understood in this PWN model. We show that the optical bump and X-ray flare may have different origins. Our results strengthen the evidence for the formation of magnetars in double neutron star mergers and justify the validity of the PWN model.

In this paper, we investigate two hydrogen-poor superluminous supernovae (SLSNe) iPTF15esb and iPTF13dcc whose light curves (LCs) show significant deviation from the smooth rise and fall. The LC of iPTF15esb exhibits two peaks and a post-peak plateau, and furthermore the late-time spectrum of iPTF15esb shows a strong, broad H$\alpha$ emission line. The early-time LC of iPTF13dcc shows a long duration bump followed by the second peak. Here we propose an ejecta-circumstellar medium (CSM) interaction model involving multiple shells/winds and use it to explain the LCs of iPTF15esb and iPTF13dcc. We find that the theoretical LCs reproduced by this model can well match the observations of iPTF15esb and iPTF13dcc. Based on our results, we infer that the progenitors have undergone multiple violent mass-loss processes before the SN explosion. In addition, we find that the variation trend of our inferred densities of the shells is consistent with that predicted by the stellar mass-loss history before an SN explosion. Further investigations for other bumpy SLSNe/SNe would shed light on their nature and provide a probe for the mass-loss history of their progenitors.

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

Broad-lined type Ic supernovae (SNe Ic-BL) are a subclass of rare core collapse SNe whose energy source is debated in the literature. Recently a series of investigations on SNe Ic-BL with the magnetar (plus 56Ni) model were carried out. Evidence for magnetar formation was found for the well-observed SNe Ic-BL 1998bw and 2002ap. In this paper we systematically study a large sample of SNe Ic-BL not associated with gamma-ray bursts. We use photospheric velocity data determined in a homogeneous way. We find that the magnetar+56Ni model provides a good description of the light curves and velocity evolution of our sample of SNe Ic-BL, although some SNe (not all) can also be described by the pure-magnetar model or by the two-component pure-56Ni model (3 out of 12 are unlikely explained by two-component model). In the magnetar+56Ni model, the amount of 56Ni required to explain their luminosity is significantly reduced, and the derived initial explosion energy is, in general, in accordance with neutrino heating. Some correlations between different physical parameters are evaluated and their implications regarding magnetic field amplification and the total energy reservoir are discussed.

Spectra of 38 candidate or known cataclysmic variables are presented. Most are candidate dwarf novae or systems containing possible highly magnetic white dwarfs, while a few (KR Aur, LS Peg, V380 Oph and V694 Mon) are previously known objects caught in unusual states. Individual spectra are used to confirm a dwarf nova nature or other classification while radial velocities of 15 systems provide orbital periods and velocity amplitudes that aid in determining the nature of the objects. Our results substantiate a polar nature for four objects, find an eclipsing SW Sex star below the period gap, another as a likely intermediate polar, as well as two dwarf novae with periods in the middle of the gap.

Very recently, the gravitational-wave (GW) event GW170817 was discovered to be associated with the short gamma-ray burst (GRB) 170817A. Multi-wavelength follow-up observations were carried out, and X-ray, optical, and radio counterparts to GW170817 were detected. The observations undoubtedly indicate that GRB 170817A originates from a binary neutron star merger. However, the GRB falls into the low-luminosity class that could have a higher statistical occurrence rate and detection probability than the normal (high-luminosity) class. This implies the possibility that GRB 170817A is intrinsically powerful, but we are off-axis and only observe its side emission. In this Letter, we provide a timely modeling of the multi-wavelength afterglow emission from this GRB and the associated kilonova signal from the merger ejecta, under the assumption of a structured jet, a two-component jet, and an intrinsically less-energetic quasi-isotropic fireball, respectively. Comparing the afterglow properties with the multi-wavelength follow-up observations, we can distinguish between these three models. Furthermore, a few model parameters (e.g., the ejecta mass and velocity) can be constrained.

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

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.

We study the most luminous known supernova (SN) associated with a gamma-ray burst (GRB), SN 2011kl. The photospheric velocity of SN 2011kl around peak brightness is $21,000\pm7,000$ km s$^{-1}$. Owing to different assumptions related to the light-curve (LC) evolution (broken or unbroken power-law function) of the optical afterglow of GRB 111209A, different techniques for the LC decomposition, and different methods (with or without a near-infrared contribution), three groups derived three different bolometric LCs for SN 2011kl. Previous studies have shown that the LCs without an early-time excess preferred a magnetar model, a magnetar+$^{56}$Ni model, or a white dwarf tidal disruption event model rather than the radioactive heating model. On the other hand, the LC shows an early-time excess and dip that cannot be reproduced by the aforementioned models, and hence the blue-supergiant model was proposed to explain it. Here we reinvestigate the energy sources powering SN 2011kl. We find that the two LCs without the early-time excess of SN 2011kl can be explained by the magnetar+$^{56}$Ni model, and the LC showing the early excess can be explained by the magnetar+$^{56}$Ni model taking into account the cooling emission from the shock-heated envelope of the SN progenitor, demonstrating that this SN might primarily be powered by a nascent magnetar.

The LIGO detection of gravitational waves (GW) from merging black holes in 2015 marked the beginning of a new era in observational astronomy. The detection of an electromagnetic signal from a GW source is the critical next step to explore in detail the physics involved. The Antarctic Survey Telescopes (AST3), located at Dome A, Antarctica, is uniquely situated for rapid response time-domain astronomy with its continuous night-time coverage during the austral winter. We report optical observations of the GW source (GW~170817) in the nearby galaxy NGC 4993 using AST3. The data show a rapidly fading transient at around 1 day after the GW trigger, with the $i$-band magnitude declining from $17.23\pm0.13$ magnitude to $17.72\pm0.09$ magnitude in $\sim 1.8$ hour. The brightness and time evolution of the optical transient associated with GW~170817 are broadly consistent with the predictions of models involving merging binary neutron stars. We infer from our data that the merging process ejected about $\sim 10^{-2}$ solar mass of radioactive material at a speed of up to $30\%$ the speed of light.

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

A binary neutron star (BNS) merger has been widely argued to be one of the progenitors of a short gamma-ray burst (SGRB). This central engine can be verified if its gravitational-wave (GW) event is detected simultaneously. Once confirmed, this kind of association will be a landmark in multi-messenger astronomy and will greatly enhance our understanding of the BNS merger processes. Due to the limited detection horizon of BNS mergers for the advanced LIGO/Virgo GW observatories, we are inclined to local SGRBs within few hundreds of mega-parsecs. Since normal SGRBs rarely fall into such a close range, to make it more observationally valuable, we have to focus on low-luminosity SGRBs which have a higher statistical occurrence rate and detection probability. However, there is a possibility that an observed low-luminosity SGRB is intrinsically powerful but we are off-axis and only observe its side emission. In this paper, we provide some theoretical predictions of both the off-axis afterglow emission from a nearby SGRB under the assumption of a structured jet and the macronova signal from the ejecta of this GW-detectable BNS merger. From the properties of the afterglow emission, we could distinguish an off-axis normal SGRB from an intrinsically low-energy quasi-isotropic class. Furthermore, with follow-up multi-wavelength observations, a few parameters for BNS mergers (e.g. the medium density and the ejecta mass and velocity) would be constrained.

Recently, a precise (sub-arcsecond) localization of the repeating fast radio burst (FRB) 121102 has led to the discovery of persistent radio and optical counterparts, the identification of a host dwarf galaxy at a redshift of $z=0.193$, and several campaigns of searches for higher-frequency counterparts, which gave only upper limits on the emission flux. Although the origin of FRBs remains unknown, most of the existing theoretical models are associated with pulsars, or more specifically, magnetars. In this paper, we explore persistent high-energy emission from a rapidly rotating highly magnetized pulsar associated with FRB 121102 if internal gradual magnetic dissipation occurs in the pulsar wind. We find that the efficiency of converting the spin-down luminosity to the high-energy (e.g., X-ray) luminosity is generally much smaller than unity, even for a millisecond magnetar. This provides an explanation for the non-detection of high-energy counterparts to FRB 121102. We further constrain the spin period and surface magnetic field strength of the pulsar with the current high-energy observations. In addition, we compare our results with the constraints given by the other methods in previous works and would expect to apply our new method to some other open issues in the future.

The quiescent periodic photometric modulations of two low-inclination cataclysmic variables observed in Kepler K2 Campaigns 0 and 1, KZ Gem and TW Vir, are investigated. A phase-correcting method was successfully used to detect the orbital modulations of KZ Gem and TW Vir and improve their orbital periods. The light curve morphologies of both CVs were further analyzed by defining flux ratios and creating colormaps. KZ Gem shows ellipsoidal modulations with an orbital period of 0.22242(1) day, twice the period listed in the updated RK catalogue (Edition 7.24). With this newly determined period, KZ Gem is no longer a CV in the period gap, but a long-period CV. A part of the quiescent light curve of TW Vir that had the highest stability was used to deduce its improved orbital period of 0.182682(3) day. The flat patterns shown in the colormaps of the flux ratios for KZ Gem demonstrate the stability of their orbital modulations, while TW Vir show variable orbital modulations during the K2 datasets. In TW Vir, the single versus double-peaked nature of the quiescent orbital variations before and after superoutburst may be related to the effect of the superoutburst on the accretion disk.

Type I X-ray bursts in a low-mass X-ray binary (LMXB) are caused by unstable nuclear burning of accreted materials. Semi-analytical and numerical studies of unstable nuclear burning have successfully reproduced partial properties of this kind of burst. However, some other properties (e.g. the waiting time) are not well explained. In this paper, we find that the probability distributions of fluence, peak count, rise time, duration and waiting time can be described as power-law-like distributions. This indicates that type I X-ray bursts may be governed by a self-organized criticality (SOC) process. The power-law index of waiting time distribution (WTD) is around $-1$, which is not predicted by any current waiting time model. We propose a physical burst rate model, in which the mean occurrence rate is inversely proportional to time $\lambda\propto t^{-1}$. In this case, the WTD is well explained by a non-stationary Poisson process within the SOC theory. In this theory, the burst size is also predicted to follow a power-law distribution, which requires that the emission area possesses only part of the neutron star surface. Furthermore, we find that the WTDs of some astrophysical phenomena can also be described by similar occurrence rate models.

We statistically study gamma-ray burst (GRB) optical flares from the Swift/UVOT catalog. We compile 119 optical flares, including 77 flares with redshift measurements. Some tight correlations among the time scales of optical flares are found. For example, the rise time is correlated with the decay time, and the duration time is correlated with the peak time of optical flares. These two tight correlations indicate that longer rise times are associated with longer decay times of optical flares, and also suggest that broader optical flares peak at later times, which are consistent with the corresponding correlations of X-ray flares. We also study the frequency distributions of optical flare parameters, including the duration time, rise time, decay time, peak time and waiting time. Similar power-law distributions for optical and X-ray flares are found. Our statistic results imply that GRB optical flares and X-ray flares may share the similar physical origin and both of them are possibly related to central engine activities.

Recently, a hydrogen-poor superluminous supernova (SLSN) iPTF15esb at redshift $z=0.224$ was reported, whose light curve (LC) and spectrum show several unusual characteristics. Its late-time spectrum shows a strong, broad H$\alpha $ emission line and the bolometric LC exhibits two peaks and a post-peak plateau. Here we propose an ejecta-circumstellar medium (CSM) interaction model involving multiple shells and winds to explain this double-peak SN. We find that the theoretical LC reproduced by this model can well match the observations of iPTF15esb. Based on this result, we infer that the progenitor has undergone at least three violent mass-loss processes before the SN explosion. Furthermore, we find that the masses of the CSM wind and shells from the outermost shell (the first eruption) to the innermost wind (the final eruption) decrease sequentially but their densities increase. The variation trend of the inferred densities of the shells and wind is consistent with the stellar structure before an SN explosion. Further investigations for similar SLSNe would provide a probe for the mass-loss history of their progenitors.

High-energy neutrinos are expected to originate from different stages in a gamma-ray burst (GRB) event. In this work we revisit the dissipative photospheric scenario, in which the GRB prompt emission is produced around the photospheric radius. Meanwhile, possible dissipation mechanisms (e.g., internal shocks or magnetic reconnection) could accelerate cosmic-rays (CRs) to ultra-high energies and then produce neutrinos via hadronuclear and photohadronic processes, which are referred to as prompt neutrinos. In this paper, we obtain the prompt neutrino spectrum of a single GRB within a self-consistent analytical framework, in which the jet-cocoon structure and possible collimation effects are included. We investigate a possible neutrino signal from the cocoon, which has been ignored in the previous studies. We show that if a GRB event happens at a distance of the order of Mpc, there is a great chance to observe the neutrino emission from the cocoon by IceCube, which is even more promising than jet neutrinos since the opening angle of the cocoon is much larger. We also figure out the diffuse neutrino flux of GRB cocoons and find that it could be comparable with that of the jets. Our results are consistent with the latest result reported by the IceCube collaboration that no significant correlation between neutrino events and observed GRBs is seen in the new data (Aartsen et al. 2017).

Previous works show that a correlation among radio luminosity, X-ray luminosity, and black hole (BH) mass from stellar-mass BHs in X-ray binaries to supermassive BHs in active galactic nuclei (AGNs), which leads to the so-called fundamental plane of BH activity. However, there are two competing explanations for this fundamental plane, including the jet-dominated model and the disk-jet model. Thus, the physical origin of this fundamental plane remains unknown. In this paper, we show that the X-ray luminosities, radio luminosities and BH masses of gamma-ray bursts (GRBs) and M82 X-1 also show a similar distribution. The universal scaling law among stellar-mass, intermediate and supermassive BH systems, together with the fact that radio and X-ray emission of GRBs originates from relativistic jets, reveals that the fundamental plane of BH activity is controlled by a jet, i.e., the radio and X-ray emission is mainly from the jet. Our work also suggests that the jets are scale-invariant with respect to the BH mass.

In this paper we collect 19 hydrogen-deficient superluminous supernovae (SLSNe) and fit their light curves, temperature evolution, and velocity evolution based on the magnetar-powered model. To obtain the best-fitting parameters, we incorporate the Markov Chain Monte Carlo approach. We get rather good fits for 7 events ($\chi^2$/d.o.f = $0.24-0.96$) and good fits for other 7 events ($\chi^2$/d.o.f = $1.37-3.13$). We find that the initial periods ($P_{0}$) and magnetic strength ($B_{p}$) of the magnetars supposed to power these SLSNe are in the range of $\sim 1.2-8.3$ ms and $\sim (0.2-8.8)\times 10^{14}\,$G, respectively; the inferred masses of the ejecta of these SLSNe are between 1 and 27.6 $M_{\odot }$, and the values of the gamma-ray opacity $\kappa_{\gamma }$ are between 0.01 and 0.82 cm$^2$ ~g$^{-1}$. We also calculate the fraction of the initial rotational energy of the magnetars harbored in the centers of the remnants of these SLSNe that is converted to the kinetic energy of the ejecta and find that the fraction is $\sim 19-97$ for different values of $P_{0}$ and $B_{p}$, indicating that the acceleration effect cannot be neglected. Moreover, we find that the initial kinetic energies of most of these SLSNe are so small ($\lesssim 2\times 10^{51}$ erg) that they can be easily explained by the neutrino-driven mechanism. These results can help clarify some important issues related to the energy-source mechanisms and explosion mechanisms and reveal the nature of SLSNe.

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

In this paper we investigate effects of the $r$-mode instability on a newborn rapidly-rotating magnetar with fall-back accretion. Such a magnetar could usually occur in core-collapse supernovae and gamma-ray bursts. We find that the magnetar's spin and $r$-mode evolution are influenced by accretion. If the magnetar is sufficiently spun up to a few milliseconds, gravitational radiation leads to the growth of the $r$-mode amplitude significantly. The maximum $r$-mode amplitude reaches an order $\sim 0.001$ when the damping due to the growth of a toroidal magnetic field balances the growth of the $r$-mode amplitude. If such a sufficiently spun-up magnetar was located at a distance less than 1\,Mpc, then gravitational waves would be detectable by the Einstein Telescope but would have an extremely low event rate. However, if the spin-up is insufficient, the growth of the $r$-mode amplitude is mainly due to the accretion torque. In this case, the maximum $r$-mode amplitude is of the order of $\sim 10^{-6}-10^{-5}$.

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

The similarity of the host galaxy of FRB 121102 with those of long gamma-ray bursts and Type I super-luminous supernovae suggests that this FRB could be associated with a young magnetar. By assuming the FRB emission to be produced within the magnetosphere, we derive a lower limit on the age of the magnetar, after which GHz emission is able to escape freely from the dense relativistic wind of the magnetar. Another lower limit is obtained by requiring the dispersion measure contributed by the electron/positron pair wind to be consistent with the observations of the host galaxy. Furthermore, we also derive some upper limits on the magnetar age with discussions on possible energy sources of the FRB emission and the recently-discovered persistent radio counterpart. As a result, some constraints on model parameters are addressed by reconciling the lower limits with the possible upper limits that are derived with an assumption of rotational energy source.

In this paper, we propose a new scenario in which a rapidly-rotating strongly-magnetized pulsar without any surrounding supernova ejecta produces fast radio bursts (FRBs) repeatedly via some mechanisms, and meanwhile, an ultra-relativistic electron/positron pair wind from the pulsar sweeps up its ambient dense interstellar medium, giving rise to a non-relativistic pulsar wind nebula (PWN). We show that the synchrotron radio emission from such a PWN is bright enough to account for the recently-discovered persistent radio source associated with the repeating FRB 121102 in reasonable ranges of the model parameters. In addition, our PWN scenario is consistent with the non-evolution of the dispersion measure inferred from all the repeating bursts observed in four years.

Muon pairs can be produced in the annihilation of ultra-high energy (UHE, $E \gtrsim 10^{18} \,\mathrm{eV}$) photons with low energy cosmic background radiation in the intergalactic space, giving birth to neutrinos. Although the branching ratio of muon pair production is low, products of other channels, which are mainly electron/positron pairs, will probably transfer most of their energies into the new generated UHE photon in the subsequent interaction with the cosmic background radiation via Compton scattering in deep Klein-Nishina regime. The regeneration of these new UHE photons then provides a second chance to produce the muon pairs, enhancing the neutrino flux. We investigate the neutrino production in the propagation of UHE photons in the intergalactic space at different redshifts, considering various competing processes such as pair production, double pair production for UHE photons, and triplet production and synchrotron radiation for UHE electrons. Regarding the least energetic outgoing particles as energy loss, we obtain the effective penetration length of the leading particle, as well as energy loss rate including the neutrino emission rate in the cascade process. Finally, we find that an extra component of UHE neutrinos will arise from the propagation of UHE cosmic rays due to the generated UHE photons and electron/positrons. However, the flux of this component is quite small, with a flux of at most $10\%$ of that of the conventional cosmogenic neutrino at a few EeV, in the absence of a strong intergalactic magnetic field and a strong cosmic radio background. The precise contribution of extra component depends on several factors, e.g., cosmic radio background, intergalactic magnetic field, the spectrum of proton, which will be discussed in this work.

Time-tag ultraviolet data obtained on the Hubble Space Telescope in 2013 reveal interesting variability related to the white dwarf spin in the two cataclysmic variables RZ Leo and CC Scl. RZ Leo shows a period at 220s and its harmonic at 110s, thus identifying it as a likely Intermediate Polar (IP). The spin signal is not visible in a short single night of ground based data in 2016, but the shorter exposures in that dataset indicate a possible partial eclipse. The much larger UV amplitude of the spin signal in the known IP CC Scl allows the spin of 389s, previously only seen at outburst, to be visible at quiescence. Spectra created from the peaks and troughs of the spin times indicate a hotter temperature of several thousand degrees during the peak phases, with multiple components contributing to the UV light.

Broad-lined type Ic supernovae (SNe Ic-BL) are peculiar stellar explosions that distinguish themselves from ordinary SNe. Some SNe Ic-BL are associated with long-duration (\gtrsim 2 s) gamma-ray bursts (GRBs). Black holes and magnetars are two types of compact objects that are hypothesized to be central engines of GRBs. In spite of decades of investigations, no direct evidence for the formation of black holes or magnetars has been found for GRBs so far. Here we report the finding that the early peak (t \lesssim 50 days) and late-time (t \gtrsim 300 days) slow decay displayed in the light curves of both SNe 1998bw (associated with GRB 980425) and 2002ap (not GRB-associated) can be attributed to magnetar spin-down with initial rotation period P0 \sim 20 ms, while the intermediate-time (50 \lesssim t \lesssim 300 days) linear decline is caused by radioactive decay of 56Ni. The connection between the early peak and late-time slow decline in the light curves is unexpected in alternative models. We thus suggest that GRB 980425 and SN 2002ap were powered by magnetars.

We work on a GRB sample whose initial Lorentz factors ($\Gamma_0$) are constrained with the afterglow onset method and the jet opening angles ($\theta_{\rm j}$) are determined by the jet break time. We confirm the $\Gamma_0$ - $E_{\gamma,\rm iso}$ correlation by Liang et al. (2010), and the $\Gamma_0$ - $L_{\gamma,\rm iso}$ correlation by L{\"u} et al. (2012). Furthermore, we find correlations between $\Gamma_0$ and the beaming corrected $\gamma$-ray energy ($E_{\gamma}$) and mean $\gamma$-ray luminosity ($L_{\gamma}$). By also including the kinetic energy of the afterglow, we find rough correlations (with larger scatter) between $\Gamma_0$ and the total ($\gamma$-ray plus kinetic) energy and the total mean luminosity, both for isotropic values and beaming corrected values: these correlations allow us to test the data with GRB central engine models. Limiting our sample to the GRBs that likely have a black hole central engine, we compare the data with theoretical predictions of two types of jet launching mechanisms from BHs, i.e. the non-magnetized $\nu \bar\nu$ - annihilation mechanism, and the strongly magnetized Blandford-Znajek (BZ) mechanism. We find that the data are more consistent with the latter mechanism, and discuss the implications of our findings for GRB jet composition.

Broad-lined type Ic supernovae (SNe Ic-BL) are of great importance because their association with long-duration gamma-ray bursts (LGRBs) holds the key to deciphering the central engine of LGRBs, which refrains from being unveiled despite decades of investigation. Among the two popularly hypothesized types of central engine, i.e., black holes and strongly magnetized neutron stars (magnetars), there is mounting evidence that the central engine of GRB-associated SNe (GRB-SNe) is rapidly rotating magnetars. Theoretical analysis also suggests that magnetars could be the central engine of SNe Ic-BL. What puzzled the researchers is the fact that light curve modeling indicates that as much as 0.2-0.5 solar mass of 56Ni was synthesized during the explosion of the SNe Ic-BL, which is unfortunately in direct conflict with current state-of-the-art understanding of magnetar-powered 56Ni synthesis. Here we propose a dynamic model of magnetar-powered SNe to take into account the acceleration of the ejecta by the magnetar, as well as the thermalization of the injected energy. Assuming that the SN kinetic energy comes exclusively from the magnetar acceleration, we find that although a major fraction of the rotational energy of the magnetar is to accelerate the SNe ejecta, a tiny fraction of this energy deposited as thermal energy of the ejecta is enough to reduce the needed 56Ni to 0.06 solar mass for both SN 1997ef and SN 2007ru. We therefore suggest that magnetars could power SNe Ic-BL both in aspects of energetics and of 56Ni synthesis.

Recently, DeLaunay et al. (2016) discovered a gamma-ray transient, Swift J0644.5-5111, associated with the fast radio burst (FRB) 131104. They also reported follow-up broadband observations beginning two days after the FRB and provided upper limits on a putative afterglow of this transient. In this paper, we show that if such a transient drives a relativistic shock as in a cosmological gamma-ray burst (GRB), these upper limits are consistent with an environment of which density is much less than that of an interstellar medium but typical for the outskirts' density of a galaxy when the typical values of three microphysical parameters of the shock are taken. This appears to be inconsistent with the catastrophic event models in which the central engine of Swift J0644.5-5111 is surrounded by an interstellar medium, but together with the properties of the gamma-ray transient, favors the binary neutron star merger origin. We further constrain the physical parameters of the postmerger object by assuming that Swift J0644.5-5111 results from internal dissipation of a spinning-down pulsar wind, and we find that the postmerger object is an ultra-strongly magnetized, very rapidly rotating pulsar. This merger event should have given birth to a gravitational wave burst, an FRB, and a short GRB or an extended X-ray/gamma-ray emission if a relativistic jet of the GRB is missed. Such "triplets" would be testable in the near future with the advanced LIGO and Virgo gravitational-wave observatories.

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

The merger of two white dwarfs is expected to result in a central fast rotating core surrounded by a debris disk, in which magnetorotational instabilities give rise to a hot magnetized corona and a magnetized outflow. The dissipation of magnetic energy via reconnection could lead to the acceleration of cosmic-rays in the expanding material, which would result in high energy neutrinos. We discuss the possibility of using these neutrino signals as probes of the outflow dynamics, magnetic energy dissipation rate and cosmic-ray acceleration efficiency. Importantly, the accompanying high-energy gamma-rays are absorbed within these sources because of the large optical depth, so these neutrino sources can be regarded as hidden cosmic-ray accelerators that are consistent with the non-detection of gamma-rays with Fermi-LAT. While the cosmic-ray generation rate is highly uncertain, if it reaches $\sim10^{45}\,\rm erg\,Mpc^{-3}\,yr^{-1}$, the diffuse neutrino flux could contribute a substantial fraction of the IceCube observations. We also evaluate the prospect of observing individual merger events, which provides a means for testing such sources in the future.

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

eXTP is a science mission designed to study the state of matter under extreme conditions of density, gravity and magnetism. Primary targets include isolated and binary neutron stars, strong magnetic field systems like magnetars, and stellar-mass and supermassive black holes. The mission carries a unique and unprecedented suite of state-of-the-art scientific instruments enabling for the first time ever the simultaneous spectral-timing-polarimetry studies of cosmic sources in the energy range from 0.5-30 keV (and beyond). Key elements of the payload are: the Spectroscopic Focusing Array (SFA) - a set of 11 X-ray optics for a total effective area of about 0.9 m^2 and 0.6 m^2 at 2 keV and 6 keV respectively, equipped with Silicon Drift Detectors offering <180 eV spectral resolution; the Large Area Detector (LAD) - a deployable set of 640 Silicon Drift Detectors, for a total effective area of about 3.4 m^2, between 6 and 10 keV, and spectral resolution <250 eV; the Polarimetry Focusing Array (PFA) - a set of 2 X-ray telescope, for a total effective area of 250 cm^2 at 2 keV, equipped with imaging gas pixel photoelectric polarimeters; the Wide Field Monitor (WFM) - a set of 3 coded mask wide field units, equipped with position-sensitive Silicon Drift Detectors, each covering a 90 degrees x 90 degrees FoV. The eXTP international consortium includes mostly major institutions of the Chinese Academy of Sciences and Universities in China, as well as major institutions in several European countries and the United States. The predecessor of eXTP, the XTP mission concept, has been selected and funded as one of the so-called background missions in the Strategic Priority Space Science Program of the Chinese Academy of Sciences since 2011. The strong European participation has significantly enhanced the scientific capabilities of eXTP. The planned launch date of the mission is earlier than 2025.

In this review we present a progress report of the connection between long-duration gamma-ray bursts (GRBs) and their accompanying supernovae (SNe). The analysis is from the point of view of an observer, with much of the emphasis placed on how observations, and the modelling of observations, have constrained what we known about GRB-SNe. We discuss their photometric and spectroscopic properties, their role as cosmological probes, including their measured luminosity$-$decline relationships, and how they can be used to measure the Hubble constant. We present a statistical analysis of their bolometric properties, and use this to determine the properties of the "average" GRB-SNe: which has a kinetic energy of $E_{\rm K} \approx 2.5\times10^{52}$ erg, an ejecta mass of $M_{\rm ej} \approx 6$ M$_{\odot}$, a nickel mass of $M_{\rm Ni} \approx 0.4$ M$_{\odot}$, a peak photospheric velocity of $v_{\rm ph} \approx 21,000$ km s$^{-1}$, a peak bolometric luminosity of $L_{\rm p} \approx 1\times10^{43}$ erg s$^{-1}$, and it reaches peak bolometric light in $t_{\rm p} \approx 13$ days. We discuss their geometry, consider the various physical processes that are thought to power the luminosity of GRB-SNe, and whether differences exist between GRB-SNe and the SNe associated with ultra-long duration GRBs. We discuss how observations of the environments of GRB-SNe further constrain the physical properties of their progenitor stars, and give an overview of the current theoretical paradigms of their suspected central engines. We also present an overview of the radioactively powered transients that have been photometrically associated with short-duration GRBs. We conclude the review by discussing what additional research is needed to further our understanding of GRB-SNe, in particular the role of binary-formation channels and the connection of GRB-SNe with superluminous SNe (abridged).

Almost all superluminous supernovae (SLSNe) whose peak magnitudes are $\lesssim -21$ mag can be explained by the $^{56}$Ni-powered model, magnetar-powered (highly magnetized pulsar) model or ejecta-circumstellar medium (CSM) interaction model. Recently, iPTF13ehe challenges these energy-source models, because the spectral analysis shows that $\sim 2.5M_\odot$ of $^{56}$Ni have been synthesized but are inadequate to power the peak bolometric emission of iPTF13ehe, while the rebrightening of the late-time light-curve (LC) and the H$\alpha$ emission lines indicate that the ejecta-CSM interaction must play a key role in powering the late-time LC. Here we propose a triple-energy-source model, in which a magnetar together with some amount ($\lesssim 2.5M_\odot$) of $^{56}$Ni may power the early LC of iPTF13ehe while the late-time rebrightening can be quantitatively explained by an ejecta-CSM interaction. Furthermore, we suggest that iPTF13ehe is a genuine core-collapse supernova rather than a pulsational pair-instability supernova candidate. Further studies on similar SLSNe in the future would eventually shed light on their explosion and energy-source mechanisms.

Very recently Spitler et al. (2016) and Scholz et al. (2016) reported their detections of sixteen additional bright bursts from the direction of the fast radio burst (FRB) 121102. This repeating FRB is inconsistent with all the catastrophic event models put forward previously for hypothetically non-repeating FRBs. Here we propose a different model, in which highly magnetized pulsars travel through asteroid belts of other stars. We show that a repeating FRB could originate from such a pulsar encountering lots of asteroids in the belt. During each pulsar-asteroid impact, an electric field induced outside the asteroid has such a large component parallel to the stellar magnetic field that electrons are torn off the asteroidal surface and accelerated to ultra-relativistic energies instantaneously. Subsequent movement of these electrons along magnetic field lines will cause coherent curvature radiation, which can account for all the properties of an FRB. In addition, this model can self-consistently explain the typical duration, luminosity, and repetitive rate of the seventeen bursts of FRB 121102. The predicted occurrence rate of repeating FRB sources may imply that our model would be testable in the next few years.

The X-ray afterglows of almost one half of gamma-ray bursts have been discovered to have a shallow decay phase by the {\em Swift} satellite, whose origin remains mysterious. Two main models have been proposed to explain this phase, relativistic wind bubbles (RWBs) and structured ejecta, which could originate from millisecond magnetars and rapidly-rotating black holes, respectively. Based on these models, we here investigate polarization evolution in the shallow decay phase of X-ray and optical afterglows. We find that in the RWB model, a significant bump of the polarization degree evolution curve appears during the shallow decay phase of both optical and X-ray afterglows, while the polarization position angle changes its direction by $90^\circ$ abruptly. In the structured ejecta model, however, the polarization degree does not evolve significantly during the shallow decay phase of afterglows, no matter whether the magnetic field configuration in the ejecta is random or globally large-scale. Therefore, we conclude that these two models for the shallow decay phase and relevant central engines would be testable with future polarization observations.

Star-forming and starburst galaxies are considered as one of the viable candidate sources of the high-energy cosmic neutrino background detected in IceCube. We revisit contributions of supernova remnants (SNRs) and hypernova remnants (HNRs) in such galaxies to the diffuse high-energy neutrino and gamma-ray backgrounds, in light of the latest Fermi data above 50GeV. We also take into account possible time dependent effects of the cosmic-ray (CR) acceleration during the SNR evolution. CRs accelerated by the SNR shocks can produce high-energy neutrinos up to $\sim100$ TeV energies, but CRs from HNRs can extend the spectrum up to PeV energies. We show that, only if HNRs are dominant over SNRs, the diffuse neutrino background above 100 TeV can be explained without contradicting the gamma-ray data. However, the neutrino data around 30 TeV remain unexplained, which might suggest a different population of gamma-ray dark CR sources. Alternatively, we consider possible contributions of Pop-III HNRs up to $z\lesssim10$, and show that they are not constrained by the gamma-ray data, and thus could contribute to the diffuse high-energy backgrounds if their explosion energy reaches ${\mathcal E}_{\rm POP-III}\sim{\rm a~few}\times10^{53}$erg. More conservatively, our results suggest that the explosion energy of POP-III HNRs is ${\mathcal E}_{\rm POP-III}\lesssim7\times{10}^{53}$erg.

Context: The LIGO consortium announced the first direct detection of gravitation wave event GW150914 from two merging black holes; however the nature of the black holes are still not clear. Aims: We study whether electromagnetic radiation can be detected from merging stellar mass black binaries like GW150914. Methods: We briefly investigate the possible growth and merging processes of the two stellar mass black holes in the merging event of GW150914 detected by aLIGO, as clocked by a distant external observer. Our main results are: (1) The description of the black hole growth using stationary metric of a pre-existing black hole predicts strong electromagnetic radiation from merging black holes, which is inconsistent with GW150914; (2) Only gravitational wave radiation can be produced in the coalescence of two black holes such as that in the GW150914 event, if the black hole growth is described using time-dependent metric considering the influence of the in-falling matter onto a pre-existing black hole, as clocked by a distant external observer. Conclusions: Future high sensitivity detections of gravitational waves from merging black holes might be used to probe matter distribution and space-time geometry in the vicinity of the horizon. Perhaps the GW150914-like events can be identified with traditional astronomy observations only if the black holes are embedded in extremely dense medium before their final merge, when very strong electromagnetic radiation is produced and can escape from the system.

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

In this paper we propose that a fast radio burst (FRB) could originate from the magnetic interaction between double neutron stars (NSs) during their final inspiral within the framework of a unipolar inductor model. In this model, an electromotive force is induced on one NS to accelerate electrons to an ultra-relativistic speed instantaneously. We show that coherent curvature radiation from these electrons moving along magnetic field lines in the magnetosphere of the other NS is responsible for the observed FRB signal, that is, the characteristic emission frequency, luminosity, duration and event rate of FRBs can be well understood. In addition, we discuss several implications of this model, including double-peaked FRBs and possible associations of FRBs with short-duration gamma-ray bursts and gravitational wave events.

The study of short-duration gamma-ray bursts provides growing evidence that a good fraction of double neutron star mergers lead to the formation of stable millisecond magnetars. The launch of Poynting flux by the millisecond magnetars could leave distinct electromagnetic signatures that reveal the energy dissipation processes in the magnetar wind. In previous studies (Wang & Dai 2013b; Wang et al. 2015), we assume that the magnetar wind becomes completely lepton-dominated so that electrons/positrons in the magnetar wind are accelerated by a diffusive shock. However, theoretical modeling of pulsar wind nebulae shows that in many cases the magnetic field energy in the pulsar wind may be strong enough to suppress diffusive shock acceleration. In this paper, we investigate the reverse shock emission as well as the forward shock emission with an arbitrary magnetization parameter $\sigma$ of a magnetar wind. We find that the reverse shock emission strongly depends on $\sigma$, and in particular, $\sigma \sim 0.3$ leads to the strongest reverse shock emission. Future observations would be helpful to diagnose the composition of the magnetar wind.

Fast radio bursts (FRBs) are radio bursts characterized by millisecond durations, high Galactic latitude positions, and high dispersion measures. Very recently, the cosmological origin of FRB 150418 has been confirmed by \cite{kea16}, and FRBs are now strong competitors as cosmological probes. The simple sharp feature of the FRB signal is ideal for them to probe some of the fundamental laws of physics. Here we show that by analyzing the delay time between different frequencies, the FRB data can place stringent upper limits on the rest mass of the photon. For FRB 150418 at $z=0.492$, one can potentially reach $m_{\gamma} \leq 5.2 \times 10^{-47}$ g, which is $10^{20}$ times smaller than the rest mass of electron, and is about $10^3$ times smaller than that obtained using other astrophysical sources with the same method.

We present photometry of the intermediate polar FO Aquarii obtained as part of the K2 mission using the Kepler space telescope. The amplitude spectrum of the data confirms the orbital period of 4.8508(4) h, and the shape of the light curve is consistent with the outer edge of the accretion disk being eclipsed when folded on this period. The average flux of FO Aquarii changed during the observations, suggesting a change in the mass accretion rate. There is no evidence in the amplitude spectrum of a longer period that would suggest disk precession. The amplitude spectrum also shows the white dwarf spin period of 1254.3401(4) s, the beat period of 1351.329(2) s, and 31 other spin and orbital harmonics. The detected period is longer than the last reported period of 1254.284(16) s, suggesting that FO Aqr is now spinning down, and has a positive $\dot{P}$. There is no detectable variation in the spin period over the course of the K2 observations, but the phase of the spin cycle is correlated with the system brightness. We also find the amplitude of the beat signal is correlated with the system brightness.

Millisecond magnetars can be formed via several channels: core-collapse of massive stars, accretion-induced collapse of white dwarfs (WDs), double WD mergers, double neutron star (NS) mergers, and WD-NS mergers. Because the mass of ejecta from these channels could be quite different, their light curves are also expected to be diverse. We evaluate the dynamic evolution of optical transients powered by millisecond magnetars. We find that the magnetar with short spin-down timescale converts its rotational energy mostly into the kinetic energy of the transient, while the energy of a magnetar with long spin-down timescale goes into radiation of the transient. This leads us to speculate that hypernovae could be powered by magnetars with short spin-down timescales. At late times the optical transients will gradually evolve into a nebular phase because of the photospheric recession. We treat the photosphere and nebula separately because their radiation mechanisms are different. In some cases the ejecta could be light enough that the magnetar can accelerate it to a relativistic speed. It is well known that the peak luminosity of a supernova (SN) occurs when the luminosity is equal to the instantaneous energy input rate, as shown by Arnett (1979). We show that photospheric recession and relativistic motion can modify this law. The photospheric recession always leads to a delay of the peak time $t_{\mathrm{pk}}$ relative to the time $t_{\times }$ at which the SN luminosity equals the instantaneous energy input rate. Relativistic motion, however, may change this result significantly.

There is growing evidence that a stable magnetar could be formed from the coalescence of double neutron stars. In previous papers, we investigated the signature of formation of stable millisecond magnetars in radio and optical/ultraviolet bands by assuming that the central rapidly rotating magnetar deposits its rotational energy in the form of a relativistic leptonized wind. We found that the optical transient PTF11agg could be the first evidence for the formation of post-merger millisecond magnetars. To enhance the probability of finding more evidence for the post-merger magnetar formation, it is better to extend the observational channel to other photon energy bands. In this paper we propose to search the signature of post-merger magnetar formation in X-ray and especially gamma-ray bands. We calculate the SSC emission of the reverse shock powered by post-merger millisecond magnetars. We find that the SSC component peaks at $1\,{\rm GeV}$ in the spectral energy distribution and extends to $\gtrsim 10\,{\rm TeV}$ for typical parameters. These energy bands are quite suitable for Fermi/LAT and CTA, which, with their current observational sensitivities, can detect the SSC emission powered by post-merger magnetars up to $1\,{\rm Gpc}$. NuSTAR, sensible in X-ray bands, can detect the formation of post-merger millisecond magnetars at redshift $z\sim 1$. Future improvement in sensitivity of CTA can also probe the birth of post-merger millisecond magnetars at redshift $z\sim 1$. However, because of the $\gamma$-$\gamma$ collisions, strong high-energy emission is clearly predicted only for ejecta masses lower than $10^{-3}M_\odot$.

There are 15 cataclysmic variables (CVs) observed in the first two campaigns of the K2 mission. In this paper, the eight CVs showing distinct features are analyzed in detail. Among these eight, modulations during quiescence are evident at the known orbital periods in the SU UMa stars QZ Vir and RZ Leo, and at our newly determined orbital periods of 1RXS J0632+2536 and WD 1144+011. The periodogram analysis for the quiescent light curve of QZ Vir reveals multi-period modulations and the coexistence of orbital and superhump periods. The phased orbital light curves for the other 3 CVs in quiescence display wide (about half cycle) and shallow (< 0.5 mag) eclipse features. Besides these modulations, their quiescent light curves reveal several transient events: a sudden decrease of system light in 1RXS J0632+2536, a low level flare-like event in QZ Vir, a short brightening event in RZ Leo and a temporary disappearance of the orbital modulation in WD 1144+011. The two known dwarf novae UV Gem and TW Vir and the CVs USNO-B1.01144-00115322 and CSS130516:111236.7+002807 show outbursts, including 1 complete and 3 incomplete normal outbursts and 2 complete superoutbursts. An incomplete but typical normal outburst confirms the dwarf nova identification of the USNO-B1.01144-00115322. The one complete normal outburst in UV Gem possibly provides the orbital period, since its modulations are shorter than the previously observed superhump period. The superoutbursts of TW Vir and CSS130516:111236.7+002807, along with their corresponding superhump periods, indicate that both objects are SU UMa stars. The derived superhump period of CSS130516:111236.7+002807 is 1.44 hr, implying that this new SU UMa star is close to the period minimum.

Fast radio bursts (FRBs) are mysterious transient sources. If extragalactic, as suggested by their relative large dispersion measures, their brightness temperatures must be extremely high. Some FRB models (e.g. young pulsar model, magnetar giant flare model, or supra-massive neutron star collapse model) suggest that they may be associated with a synchrotron nebula. Here we study a synchrotron-heating process by an FRB in a self-absorbed synchrotron nebula. If the FRB frequency is below the synchrotron self-absorption frequency of the nebula, electrons in the nebula would absorb FRB photons, leading to a harder electron spectrum and enhanced self-absorbed synchrotron emission. In the meantime, the FRB flux is absorbed by the nebula electrons. We calculate the spectra of FRB-heated synchrotron nebulae, and show that the nebula spectra would show a significant hump in several decades near the self-absorption frequency. Identifying such a spectral feature would reveal an embedded FRB in a synchrotron nebula.

In this paper we show that the most luminous supernova discovered very recently, ASASSN-15lh, could have been powered by a newborn ultra-strongly-magnetized pulsar, which initially rotates near the Kepler limit. We find that if this pulsar is a neutron star, its rotational energy could be quickly lost as a result of gravitational-radiation-driven r-mode instability; if it is a strange quark star, however, this instability is highly suppressed due to a large bulk viscosity associated with the nonleptonic weak interaction among quarks and thus most of its rotational energy could be extracted to drive ASASSN-15lh. Therefore, we conclude that such an ultra-energetic supernova provides a possible signature for the birth of a strange quark star.

The central engine and jet composition of gamma-ray bursts (GRBs) remain mysterious. Here we suggest that observations on polarization evolution of early optical afterglows may shed light on these questions. We first study the dynamics of a reverse shock and a forward shock that are generated during the interaction of a relativistic jet and its ambient medium. The jet is likely magnetized with a globally large-scale magnetic field from the central engine. The existence of the reverse shock requires that the magnetization degree of the jet should not be high ($\sigma\leq 1$), so that the jet is mainly composed of baryons and leptons. We then calculate the light curve and polarization evolution of an early optical afterglow, and find that when the polarization position angle changes by $90^\circ$ during the early afterglow, the polarization degree is zero for a toroidal magnetic field but is very likely to be non-zero for an aligned magnetic field. This result would be expected to provide a probe for the central engine of GRBs, because an aligned field configuration could originate from a magnetar central engine and a toroidal field configuration could be produced from a black hole via the Blandford-Znajek mechanism. Finally, for such two kinds of magnetic field configurations, we fit the observed data of the early optical afterglow of GRB 120308A equally well.

In this paper, we use two model-independent methods to standardize long gamma-ray bursts (GRBs) using the $E_{\rm iso}-E_{\rm p}$ correlation, where $E_{\rm iso}$ is the isotropic-equivalent gamma-ray energy and $E_{\rm p}$ is the spectral peak energy. We update 42 long GRBs and try to make constraint on cosmological parameters. The full sample contains 151 long GRBs with redshifts from 0.0331 to 8.2. The first method is the simultaneous fitting method. The extrinsic scatter $\sigma_{\rm ext}$ is taken into account and assigned to the parameter $E_{\rm iso}$. The best-fitting values are $a=49.15\pm0.26$, $b=1.42\pm0.11$, $\sigma_{\rm ext}=0.34\pm0.03$ and $\Omega_m=0.79$ in the flat $\Lambda$CDM model. The constraint on $\Omega_m$ is $0.55<\Omega_m<1$ at the 1$\sigma$ confidence level. If reduced $\chi^2$ method is used, the best-fit results are $a=48.96\pm0.18$, $b=1.52\pm0.08$ and $\Omega_m=0.50\pm0.12$. The second method is using type Ia supernovae (SNe Ia) to calibrate the $E_{\rm iso}-E_{\rm p}$ correlation. We calibrate 90 high-redshift GRBs in the redshift range from 1.44 to 8.1. The cosmological constraints from these 90 GRBs are $\Omega_m=0.23^{+0.06}_{-0.04}$ for flat $\Lambda$CDM, and $\Omega_m=0.18\pm0.11$ and $\Omega_{\Lambda}=0.46\pm0.51$ for non-flat $\Lambda$CDM. For the combination of GRB and SNe Ia sample, we obtain $\Omega_m=0.271\pm0.019$ and $h=0.701\pm0.002$ for the flat $\Lambda$CDM, and for the non-flat $\Lambda$CDM, the results are $\Omega_m=0.225\pm0.044$, $\Omega_{\Lambda}=0.640\pm0.082$ and $h=0.698\pm0.004$. These results from calibrated GRBs are consistent with that of SNe Ia. Meanwhile, the combined data can improve cosmological constraints significantly, comparing to SNe Ia alone. Our results show that the $E_{\rm iso}-E_{\rm p}$ correlation is promising to probe the high-redshift universe.

In this paper, we obtain the NIRB and SBGWs from the early stars, which are constrained by the observation of reionization and star formation rate. We study the transition from Pop III to Pop II stars via the star formation model of different population, which takes into account the reionization and the metal enrichment evolution. We calculate the two main metal pollution channels arising from the supernova-driven protogalactic outflows and "genetic channel". We obtain the SFRs of Pop III and Pop II and their NIRB and SBGWs radiation. We predict that the upper limit of metallicity in metal-enriched IGM (the galaxies whose polluted via "genetic channel") reaches $Z_{\rm crit}=10^{-3.5}Z_{\odot}$ at $z\sim13$ ($z\sim11$), which is consistent with our star formation model. We constrain on the SFR of Pop III stars from the observation of reionization. The peak intensity of NIRB is about $0.03-0.2~nW m^{-2}{sr}^{-1}$ at $\sim 1 \mu m$ for $z>6$. The prediction of NIRB signal is consistent with the metallicity evolution. We also obtain the gravitational wave background from the black holes formed by these early stars. The predicted gravitational wave background has a peak amplitude of $\Omega_{GW}\simeq8\times10^{-9}$ at $\nu=158$ Hz for Pop II star remnants. However, the background generated by Pop III.2 stars is much less than Pop II stars, with a peak amplitude of $\Omega_{GW}\simeq1.2\times10^{-11}$ at $\nu=28~Hz$. The background of Pop III.1 shifts to lower frequencies, and the amplitude of $\Omega _{GW}$ for Pop III.1 stars shows a minimum value at $\nu\simeq 10$ Hz, due to the lack of gravitational wave signals from the stars with $140~M_{\odot}<M_\ast<260~M_{\odot}$.

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

X-ray flares were discovered in the afterglow phase of gamma-ray bursts (GRBs) by the {\em Swift} satellite a decade ago and known as a canonical component in GRB X-ray afterglows. In this paper, we constrain the Lorentz factors of GRB X-ray flares using two different methods. For the first method, we estimate the lower limit on the bulk Lorentz factor with the flare duration and jet break time. In the second method, the upper limit on the Lorentz factor is derived by assuming that the X-ray flare jet has undergone saturated acceleration. We also re-estimate the initial Lorentz factor with GRB afterglow onsets, and find the coefficient of the theoretical Lorentz factor is 1.67 rather than the commonly used 2 for interstellar medium (ISM) and 1.44 for the wind case. We find that the correlation between the limited Lorentz factor and the isotropic radiation energy of X-ray flares in the ISM case is more consistent with that of prompt emission than the wind case in a statistical sense. For a comparison, the lower limit on Lorentz factor is statistically larger than the extrapolation from prompt bursts in the wind case. Our results indicate that X-ray flares and prompt bursts are produced by the same physical mechanism.

We provide a general analysis on the properties of emitting material of some rapidly evolving and luminous transients discovered recently with the Pan-STARRS1 Medium Deep Survey. It is found that these transients are probably produced by a low-mass non-relativistic outflow that is continuously powered by a newly born, rapidly spinning, and highly magnetized neutron star. Such a system could originate from an accretion-induced collapse of a white dwarf or a merger of a neutron star-neutron star binary. Therefore, observations to these transients would be helpful for constraining white dwarf and neutron star physics and/or for searching and identifying gravitational wave signals from the mergers.

Most type-Ic core-collapse supernovae (CCSNe) produce $^{56}$Ni and neutron stars (NSs) or black holes (BHs). The dipole radiation of nascent NSs has usually been neglected in explaining supernovae (SNe) with peak absolute magnitude $M_{\rm peak}$ in any band are $\gtrsim -19.5$~mag, while the $^{56}$Ni can be neglected in fitting most type-Ic superluminous supernovae (SLSNe Ic) whose $M_{\rm peak}$ in any band are $\lesssim -21$~mag, since the luminosity from a magnetar (highly magnetized NS) can outshine that from a moderate amount of $^{56}$Ni. For luminous SNe Ic with $-21 \lesssim M_{\rm peak}\lesssim -19.5$~mag, however, both contributions from $^{56}$Ni and NSs cannot be neglected without serious modeling, since they are not SLSNe and the $^{56}$Ni mass could be up to $\sim 0.5 M_{\odot}$. In this paper we propose a unified model that contain contributions from both $^{56}$Ni and a nascent NS. We select three luminous SNe Ic-BL, SN~2010ay, SN~2006nx, and SN~14475, and show that, if these SNe are powered by $^{56}$Ni, the ratio of $M_{\rm Ni}$ to $M_{\rm ej}$ are unrealistic. Alternatively, we invoke the magnetar model and the hybrid ($^{56}$Ni + NS) model and find that they can fit the observations, indicating that our models are valid and necessary for luminous SNe Ic. Owing to the lack of late-time photometric data, we cannot break the parameter degeneracy and thus distinguish among the model parameters, but we can expect that future multi-epoch observations of luminous SNe can provide stringent constraints on $^{56}$Ni yields and the parameters of putative magnetars.

Gamma-ray bursts (GRBs) are the most violent explosions in the Universe and can be used to explore the properties of high-redshift universe. It is believed that the long GRBs are associated with the deaths of massive stars. So it is possible to use GRBs to investigate the star formation rate (SFR). In this paper, we use Lynden-Bell's $c^-$ method to study the luminosity function and rate of \emph{Swift} long GRBs without any assumptions. We find that the luminosity of GRBs evolves with redshift as $L(z)\propto g(z)=(1+z)^k$ with $k=2.43_{-0.38}^{+0.41}$. After correcting the redshift evolution through $L_0(z)=L(z)/g(z)$, the luminosity function can be expressed as $\psi(L_0)\propto L_0^{-0.14\pm0.02}$ for dim GRBs and $\psi(L_0)\propto L_0^{-0.70\pm0.03}$ for bright GRBs, with the break point $L_{0}^{b}=1.43\times10^{51}~{\rm erg~s^{-1}}$. We also find that the formation rate of GRBs is almost constant at $z<1.0$ for the first time, which is remarkably different from the SFR. At $z>1.0$, the formation rate of GRB is consistent with the SFR. Our results are dramatically different from previous studies. Some possible reasons for this low-redshift excess are discussed. We also test the robustness of our results with Monte Carlo simulations. The distributions of mock data (i.e., luminosity-redshift distribution, luminosity function, cumulative distribution and $\log N-\log S$ distribution) are in good agreement with the observations. Besides, we also find that there are remarkable difference between the mock data and the observations if long GRB are unbiased tracers of SFR at $z<1.0$.

There is a sign that long-duration gamma-ray bursts (GRBs) originate from the core collapse of massive stars. During a jet puncturing through the progenitor envelope, high energy neutrinos can be produced by the reverse shock formed at the jet head. It is suggested that low-luminosity GRBs (LL-GRBs) are possible candidates of this high energy neutrino precursor up to $\sim {\rm PeV}$. Before leaving the progenitor, these high energy neutrinos must oscillate from one flavor to another with matter effect in the envelope. Under the assumption of a power-law stellar envelope density profile $\rho \propto r^{-\alpha}$ with an index $\alpha$, we study the properties of ${\rm TeV-PeV}$ neutrino oscillation. We find that adiabatic conversion is violated for these neutrinos so we do certain calibration of level crossing effect. The resonance condition is reached for different energies at different radii. We notice that the effective mixing angles in matter for ${\rm PeV}$ neutrinos are close to zero so the transition probabilities from one flavor to another are almost invariant for ${\rm PeV}$ neutrinos. We plot all the transition probabilities versus energy of ${\rm TeV-PeV}$ neutrinos from the birth place to the surface of the progenitor. With an initial flavor ratio $\phi_{\nu_e}^0:\phi_{\nu_\mu}^0:\phi_{\nu_\tau}^0=1:2:0$, we plot how the flavor ratio evolves with energy and distance when neutrinos are still in the envelope, and further get the ratio when they reach the Earth. For ${\rm PeV}$ neutrinos, the ratio is always $\phi_{\nu_e}:\phi_{\nu_\mu}:\phi_{\nu_\tau}\simeq0.30:0.37:0.33$ on Earth. In addition, we discuss the dependence of the flavor ratio on energy and $\alpha$ and get a pretty good result. This dependence may provide a promising probe of the progenitor structure.

Gamma-ray bursts (GRBs) are the most luminous electromagnetic explosions in the Universe, which emit up to $8.8\times10^{54}$ erg isotropic equivalent energy in the hard X-ray band. The high luminosity makes them detectable out to the largest distances yet explored in the Universe. GRBs, as bright beacons in the deep Universe, would be the ideal tool to probe the properties of high-redshift universe: including the cosmic expansion and dark energy, star formation rate, the reionization epoch and the metal enrichment history of the Universe. In this article, we review the luminosity correlations of GRBs, and implications for constraining the cosmological parameters and dark energy. Observations show that the progenitors of long GRBs are massive stars. So it is expected that long GRBs are tracers of star formation rate. We also review the high-redshift star formation rate derived from GRBs, and implications for the cosmic reionization history. The afterglows of GRBs generally have broken power-law spectra, so it is possible to extract intergalactic medium (IGM) absorption features. We also present the capability of high-redshift GRBs to probe the pre-galactic metal enrichment and the first stars.

There is accumulating evidence that at least a fraction of binary neutron star mergers result in rapidly spinning magnetars, with subrelativistic neutron-rich ejecta as massive as a small fraction of solar mass. The ejecta could be heated continuously by the Poynting flux emanated from the central magnetars. Such Poynting flux could become lepton-dominated so that a reverse shock develops. It was demonstrated that such a picture is capable of accounting for the optical transient PTF11agg (Wang & Dai 2013b). In this paper we investigate the X-ray and ultraviolet (UV) radiation as well as the optical and radio radiation studied by Wang & Dai (2013b). UV emission is particularly important because it has the right energy to ionize the hot ejecta at times $t\lesssim 600$ s. It is thought that the ejecta of binary neutron star mergers are a remarkably pure sample of r-process material, about which our understanding is still incomplete. In this paper we evaluate the possibility of observationally determining the bound-bound and bound-free opacities of the r-process material by timing the X-ray, UV, and optical radiation. It is found that these timings depend on the opacities weakly and therefore only loose constraints on the opacities can be obtained.

Recently, researches performed by two groups have revealed that the magnetar spin-down energy injection model with full energy trapping can explain the early-time light curves of SN 2010gx, SN 2013dg, LSQ12dlf, SSS120810 and CSS121015, but fails to fit the late-time light curves of these Superluminous Supernovae (SLSNe). These results imply that the original magnetar-powered model is challenged in explaining these SLSNe. Our paper aims to simultaneously explain both the early- and late-time data/upper limits by considering the leakage of hard emissions. We incorporate quantitatively the leakage effect into the original magnetar-powered model and derive a new semi-analytical equation. Comparing the light curves reproduced by our revised magnetar-powered model to the observed data and/or upper limits of these five SLSNe, we found that the late-time light curves reproduced by our semi-analytical equation are in good agreement with the late-time observed data and/or upper limits of SN 2010gx, CSS121015, SN 2013dg and LSQ12dlf and the late-time excess of SSS120810, indicating that the magnetar-powered model might be responsible for these SLSNe and that the gamma ray and X-ray leakage are unavoidable when the hard photons were down-Comptonized to softer photons. To determine the details of the leakage effect and unveil the nature of SLSNe, more high quality bolometric light curves and spectra of SLSNe are required.

X-ray flares have been discovered in black hole systems, such as gamma-ray bursts, the tidal disruption event Swift J1644+57, the supermassive black hole Sagittarius A$^*$ at the center of our Galaxy, and some active galactic nuclei. Their occurrences are always companied by relativistic jets. However, it is still unknown whether there is a physical analogy among such X-ray flares produced in black hole systems spanning nine orders of magnitude in mass. Here we report the observed data of X-ray flares, and show that they have three statistical properties similar to solar flares, including power-law distributions of energies, durations, and waiting times, which both can be explained by a fractal-diffusive self-organized criticality model. These statistical similarities, together with the fact that solar flares are triggered by a magnetic reconnection process, suggest that all of the X-ray flares are consistent with magnetic reconnection events, implying that their concomitant relativistic jets may be magnetically dominated.

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

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.

Relativistic jets are universal in long-duration gamma-ray burst (GRB) models. Before breaking out, they must propagate in the progenitor envelope along with a forward shock and a reverse shock forming at the jet head. Both electrons and protons will be accelerated by the shocks. High energy neutrinos could be produced by these protons interacting with stellar materials and electron-radiating photons. The jet will probably be collimated, which may have a strong effect on the final neutrino flux. Under the assumption of a power-law stellar-envelope density profile $\rho \propto r^{-\alpha}$ with an index $\alpha$, we calculate this neutrino emission flux by these shocks for low-luminosity GRBs (LL-GRBs) and ultra-long GRBs (UL-GRBs) in different collimation regimes, using the jet propagation framework developed by \citet{bro11}. We find that LL-GRBs and UL-GRBs are capable for detectable high energy neutrinos up to $\sim {\rm PeV}$, and obtain the final neutrino spectrum. Besides, we conclude that larger $\alpha$ corresponds to greater neutrino flux at high energy end ($\sim {\rm PeV}$) and higher maximum neutrino energy as well. However, such differences are so small that it is not promising for us to distinguish from observations, given the energy resolution we have now.

We investigate the mass distribution of long gamma-ray burst (GRB) host galaxies and the redshift distribution of long GRBs by considering that long GRBs occur in low-metallicity environments. We calculate the upper limit on the stellar mass of a galaxy which can produce long GRBs by utilizing the mass-metallicity (M-Z) relation of galaxies. After comparing with the observed GRB host galaxies masses, we find that the observed GRB host galaxy masses can fit the predicted masses well if GRBs occur in low-metallicity $12+\log\rm(O/H)_{\rm KK04}<8.7$. GRB host galaxies have low metallicity, low mass, and high star formation rate compared with galaxies of seventh data release of the Sloan Digital Sky Survey. We also study the cumulative redshift distribution of the latest \emph{Swift} long GRBs by adding dark GRBs and 10 new GRBs redshifts from TOUGH survey. The observed discrepancy between the GRB rate and the star formation history can be reconciled by considering that GRBs tend to occur in low-metallicity galaxies with $12+\log\rm(O/H)_{\rm KK04}<8.7$. We conclude that the metallicity cutoff that can produce long GRBs is about $12+\log\rm(O/H)_{\rm KK04}<8.7$ from the host mass distribution and redshift distribution.

Active galactic nuclei (AGNs) and gamma-ray bursts (GRBs) are powerful astrophysical events with relativistic jets. In this Letter the broadband spectral properties are compared between GRBs and the well-observed blazars. The distribution of GRBs are consistent with the well-known blazar sequence including the $\nu L_\nu(5\rm GHz)-\alpha_{\rm RX}$ and $\nu L_\nu(5\rm GHz)-\nu_{\rm peak}$ correlations, where $\alpha_{\rm RX}$ is defined as the broadband spectral slope in radio-to-X-ray bands, and $\nu_{\rm peak}$ is defined as the spectral peak frequency. Moreover, GRBs occupy the low radio luminosity end of these sequences. These two correlations suggest that GRBs could have a similar radiation process with blazars both in the prompt emission and afterglow phases, i.e., synchrotron radiation.

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

In this paper we give a brief review of our recent studies on the long and short gamma-ray bursts (GRBs) detected Swift, in an effort to understand the puzzle of classifying GRBs. We consider that it is still an appealing conjecture that both long and short GRBs are drawn from the same parent sample by observational biases.

Many long-duration gamma-ray bursts (GRBs) were observed by {\it Swift}/XRT to have plateaus in their X-ray afterglow light curves. This plateau phase has been argued to be evidence for long-lasting activity of magnetar (ultra-strongly magnetized neutron stars) central engines. However, the emission efficiency of such magnetars in X-rays is still unknown. Here we collect 24 long GRB X-ray afterglows showing plateaus followed by steep decays. We extend the well-known relationship between the X-ray luminosity ${L_{\mathrm{X}}}$ and spin-down luminosity $L_{\mathrm{sd}}$ of pulsars to magnetar central engines, and find that the initial rotation period $P_{0}$ ranges from 1 ms to 10 ms and that the dipole magnetic field $B$ is centered around $10^{15}$ G. These constraints not only favor the suggestion that the central engines of some long GRBs are very likely to be rapidly rotating magnetars but also indicate that the magnetar plateau emission efficiency in X-rays is close to 100%.

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

The equation of state (EOS), $w(z)$, is the most important parameter of dark energy. We reconstruct the evolution of this EOS in a model-independent way using the latest cosmic microwave background (CMB) data from Planck and other observations, such as type Ia supernovae (SNe Ia), the baryonic acoustic oscillation measurements (SDSS, 6dF, BOSS, and WiggleZ), and the Hubble parameter value $H(z)$. The results show that the EOS is consistent with the cosmological constant at the $2\sigma$ confidence level, not preferring a dynamical dark energy. The uncorrelated EOS of dark energy constraints from Planck CMB data are much tighter than those from the WMAP 9-year CMB data.

A long-duration gamma-ray burst (GRB) has been widely thought to arise from the collapse of a massive star, and it has been suggested that its ambient medium is a homogenous interstellar medium (ISM) or a stellar wind. There are two shocks when an ultra-relativistic fireball that has been ejected during the prompt gamma-ray emission phase sweeps up the circumburst medium: a reverse shock that propagates into the fireball, and a forward shock that propagates into the ambient medium. In this paper, we investigate the temporal evolution of the dynamics and emission of these two shocks in an environment with a general density distribution of $n\propto R^{-k}$ (where $R$ is the radius) by considering thick-shell and thin-shell cases. A GRB afterglow with one smooth onset peak at early times is understood to result from such external shocks. Thus, we can determine the medium density distribution by fitting the onset peak appearing in the light curve of an early optical afterglow. We apply our model to 19 GRBs, and find that their $k$ values are in the range of 0.4 - 1.4, with a typical value of $k\sim1$, implying that this environment is neither a homogenous interstellar medium with $k=0$ nor a typical stellar wind with $k=2$. This shows that the progenitors of these GRBs might have undergone a new mass-loss evolution.

Twenty-one new optical light curves, including five curves obtained in 2009 and sixteen curves detected from the AAVSO International Database spanning from 1977 to 2011, demonstrate 16 new primary minimum light times in the high state. Furthermore, seven newly found low-state transient events from 2006 to 2009 were discovered, consisting of five Gaussian-shaped events and two events with an exponential form with decay timescales of $<$0.005 days; these timescales are one order of magnitude shorter than those of previous X-ray flare events. In the state transition, two special events were detected: a "disrupted event" with an amplitude of $\sim$ 2 mag and a duration of $\sim$ 72 minutes and continuing R-band twin events larger than all known R-band flares detected in M-type red dwarfs. All 45 available high-state data points spanning over 35 yr were used to construct an updated O-C diagram of \objectname{AM Herculis}, which clearly shows a significant sine-like variation with a period of 12-15 yr and an amplitude of 6-9 minutes. Using the inspected physical parameters of the donor star, the secular variation in the O-C diagram cannot be interpreted by any decided angular momentum loss mechanism, but can satisfy the condition $\tau_{\dot{\rm M}_{2}}\simeq\tau_{\rm KH}>>\tau_{\dot{\rm R}_{\rm 2}}$, which is required by numerical calculations of the secular evolution of cataclysmic variables. In order to explain the prominent periodic modulation, three plausible mechanisms - spot motion, the light travel-time effect, and magnetic active cycles - are discussed in detail.

Based on the stiff equations of state of neutron stars (NS) and the discovery of high-mass NSs, a NS-NS merger will leave behind, with high probabilities, a rapidly rotating massive magnetar. The central magnetar will dissipate its rotational energy to the outflow by injecting Poynting flux, which will become lepton-dominated so that a long-lasting reverse shock (RS) is developed. We calculate the emission of the RS as well as the emission of forward shock (FS) and find that in most cases the RS emission is stronger than FS emission. It is found that the recently discovered transient, PTF11agg, can be neatly accounted for by the RS emission powered by a millisecond magnetar. Other alternative models have been considered and cannot explain the observed light curves well. We therefore suggest that PTF11agg be the first evidence for RS emission from a postmerger millisecond magnetar.

X-ray flares detected in nearly half of gamma-ray burst (GRB) afterglows are one of the most intriguing phenomena in high-energy astrophysics. All the observations indicate that the central engines of bursts, after the gamma-ray emission has ended, still have long periods of activity, during which energetic explosions eject relativistic materials, leading to late-time X-ray emission. It is thus expected that X-ray flares provide important clues to the nature of the central engines of GRBs, and more importantly, unveil the physical mechanism of the flares themselves, which has so far remained mysterious. Here we report statistical results of X-ray flares of GRBs with known redshifts, and show that X-ray flares and solar flares share three statistical properties: power-law frequency distributions for energies, durations, and waiting times. All of the distributions can be well understood within the physical framework of a self-organized criticality (SOC) system. The statistical properties of X-ray flares of GRBs are similar to solar flares, and thus can both be attributed to a SOC process. Both types of flares may be driven by a magnetic reconnection process, but X-ray flares of GRBs are produced in ultra-strongly magnetized millisecond pulsars or long-term hyperaccreting disks around stellar-mass black holes.

Recently, the GeV radiation during the X-ray flare activity in GRB 100728A was detected by Feimi/LAT. Here we study the dynamics and emission properties of collision between two homogeneous shells based on the late internal shock model. The GeV photons can be produced from X-ray flare photons up-scattered by relativistic electrons accelerated by forward-reverse shocks, where involved radiative processes include synchrotron self-Compton and crossing inverse-Compton scattering. By analytical and numerical calculations, the observed spectral properties in GRB 100728A can be well explained.

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

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

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

The possibility that we live in a special place in the universe, close to the center of a large, radially inhomogeneous void, has attracted attention recently as an alternative to dark energy or modified gravity to explain the accelerating universe. We show that the distribution of orientations of galaxy pairs can be used to test the Copernican principle that we are not in a central or special region of Universe. The popular void models can not fit both the latest type Ia supernova, cosmic microwave background data and the distribution of orientations of galaxy pairs simultaneously. Our results rule out the void models at the $4\sigma$ confidence level as the origin of cosmic acceleration and favor the Copernican principle.

We explore high-redshift gamma-ray bursts (GRBs) as promising tools to probe pre-galactic metal enrichment. We utilize the bright afterglow of a Pop III GRB exploding in a primordial dwarf galaxy as a luminous background source, and calculate the strength of metal absorption lines that are imprinted by the first heavy elements in the intergalactic medium (IGM). To derive the GRB absorption line diagnostics, we use an existing highly-resolved simulation of the formation of a first galaxy which is characterized by the onset of atomic hydrogen cooling in a halo with virial temperature >10^4 K. We explore the unusual circumburst environment inside the systems that hosted Pop III stars, modeling the density evolution with the self-similar solution for a champagne flow. For minihalos close to the cooling threshold, the circumburst density is roughly proportional to (1+z) with values of about a few cm^{-3}. In more massive halos, corresponding to the first galaxies, the density may be larger, n>100 cm^{-3}. The resulting afterglow fluxes may be detectable with the JWST and VLA in the near-IR and radio wavebands, respectively, out to redshift z>20. We predict that the maximum of the afterglow emission shifts from near-IR to millimeter bands with peak fluxes from mJy to Jy at different observed times. GRBs are ideal tools for probing the metal enrichment in the early IGM, due to their high luminosities and featureless power-law spectra. The metals in the first galaxies produced by the first supernova (SN) explosions are likely to reside in low-ionization stages. We show that if the afterglow can be observed sufficiently early, analysis of the metal lines can distinguish whether the first heavy elements were produced in a PISN, or a core-collapse (Type II) SN, thus constraining the initial mass function of the first stars.

Soft extended emission (EE) following initial hard spikes up to 100 seconds was observed with {\em Swift}/BAT for about half of short-type gamma-ray bursts (SGRBs). This challenges the conversional central engine models of SGRBs, i.e., compact star merger models. In the framework of the black hole-neutron star merger models, we study the roles of the radial angular momentum transfer in the disk and the magnetic barrier around the black hole for the activity of SGRB central engines. We show that the radial angular momentum transfer may significantly prolong the lifetime of the accretion process and multiple episodes may be switched by the magnetic barrier. Our numerical calculations based on the models of the neutrino-dominated accretion flows suggest that the disk mass is critical for producing the observed EE. In case of the mass being $\sim 0.8M_{\odot}$, our model can reproduce the observed timescale and luminosity of both the main and EE episodes in a reasonable parameter set. The predicted luminosity of the EE component is lower than the observed EE with about one order of magnitude and the timescale is shorter than 20 seconds if the disk mass being $\sim 0.2M_{\odot}$. {\em Swift}/BAT-like instruments may be not sensitive enough to detect the EE component in this case. We argue that the EE component would be a probe for merger process and disk formation for compact star mergers.

The shallow decay phase or plateau phase of early afterglows of gamma-ray bursts (GRBs), discovered by Swift, is currently understood as being due to energy injection to a relativistic blast wave. One natural scenario for energy injection invokes a millisecond magnetar as the central engine of GRBs, because the conventional model of a pulsar predicts a nearly constant magnetic-dipole-radiation luminosity within the spin-down timescale. However, we note that significant brightening occurs in some early afterglows, which apparently conflicts with the above scenario. Here we propose a new model to explain this significant brightening phenomena by considering a hyperaccreting fallback disk around a newborn millisecond magnetar. We show that for typical values of the model parameters, sufficient angular momentum of the accreted matter is transferred to the magnetar and spins it up. It is this spin-up that leads to a dramatic increase of the magnetic dipole radiation luminosity with time and thus significant brightening of an early afterglow. Based on this model, we carry out numerical calculations and fit well early afterglows of 12 GRBs assuming sufficiently strong fallback accretion. If the accretion is very weak, our model turns out to be the conventional energy-injection scenario of a pulsar. Therefore, our model can provide a unified explanation for the shallow decay phase, plateaus, and significant brightening of early afterglows.

We report here the discovery of an optical flare observed in R band from the red-dwarf eclipsing binary CU Cnc whose component stars are at the upper boundary of full convection (M1=0.43 and M2=0.4M0, M0 is the solar mass). The amplitude of the flare is the largest among those detected in R band (~0.52mag) and the duration time is about 73 minutes. As those observed on the Sun, quasi-periodic oscillations were seen during and after the flare. Three more R-band flares were found by follow up monitoring. In total, this binary was monitored photometrically by using R filter for 79.9 hours, which reveals a R-band flare rate about 0.05 flares per hour. These detections together with other strong chromospheric and coronal activities, i.e., very strong H_alpha and H_beta emission features and an EUV and X-ray source, indicate that it has very strong magnetic activity. Therefore, the apparent faintness (~1.4 magnitude in V) of CU Cnc compared with other single red dwarfs of the same mass can be plausibly explained by the high coverage of the dark spots.

The increasingly deep limit on the neutrino emission from gamma-ray bursts (GRBs) with IceCube observations has reached the level that could put useful constraints on the fireball properties. We first present a revised analytic calculation of the neutrino flux, which predicts a flux an order of magnitude lower than that obtained by the IceCube collaboration. For benchmark model parameters (e.g. the bulk Lorentz factor is \Gamma=10^{2.5}, the observed variability time for long GRBs is t_v=0.01 s and the ratio between the energy in accelerated protons and in radiation is \eta_p=10 for every burst) in the standard internal shock scenario, the predicted neutrino flux from 215 bursts during the period of the 40-string and 59-string configurations is found to be a factor of ~3 below the IceCube sensitivity. However, if we accept the recently found inherent relation between the bulk Lorentz factor and burst energy, the expected neutrino flux increases significantly and the spectral peak shifts to lower energy. In this case, the non-detection then implies that the baryon loading ratio should be \eta_p<10 if the variability time of long GRBs is fixed to t_v=0.01 s. Instead, if we relax the standard internal shock scenario but keep to assume \eta_p=10, the non-detection constrains the dissipation radius to be R>4x10^{12} cm assuming the same dissipation radius for every burst and benchmark parameters for fireballs. We also calculate the diffuse neutrino flux from GRBs for different luminosity functions existing in the literature. The expected flux exceeds the current IceCube limit for some luminosity functions, and thus the non-detection constrains \eta_p<10 in such cases when the variability time of long GRBs is fixed to t_v=0.01 s.

By using six new determined mid-eclipse times together with those collected from the literature, we found that the Observed-Calculated (O-C) curve of RR Cae shows a cyclic change with a period of 11.9 years and an amplitude of 14.3s, while it undergoes an upward parabolic variation (revealing a long-term period increase at a rate of dP/dt =+4.18(+-0.20)x10^(-12). The cyclic change was analyzed for the light-travel time effect that arises from the gravitational influence of a third companion. The mass of the third body was determined to be M_3*sin i' = 4.2(+-0.4) M_{Jup} suggesting that it is a circumbinary giant planet when its orbital inclination is larger than 17.6 degree. The orbital separation of the circumbinary planet from the central eclipsing binary is about 5.3(+-0.6)AU. The period increase is opposite to the changes caused by angular momentum loss via magnetic braking or/and gravitational radiation, nor can it be explained by the mass transfer between both components because of its detached configuration. These indicate that the observed upward parabolic change is only a part of a long-period (longer than 26.3 years) cyclic variation, which may reveal the presence of another giant circumbinary planet in a wide orbit.

We report here the tentative discovery of a Jovian planet in orbit around the rapidly pulsating subdwarf B-type (sdB-type) eclipsing binary NY Vir. By using new determined eclipse times together with those collected from the literature, we detect that the observed-calculated (O-C) curve of NY Vir shows a small-amplitude cyclic variation with a period of 7.9\,years and a semiamplitude of 6.1\,s, while it undergoes a downward parabolic change (revealing a period decrease at a rate of $\dot{P}=-9.2\times{10^{-12}}$). The periodic variation was analyzed for the light-travel time effect via the presence of a third body. The mass of the tertiary companion was determined to be $M_3\sin{i^{\prime}}=2.3(\pm0.3)$\,$M_{Jupiter}$ when a total mass of 0.60\,$M_{\odot}$ for NY Vir is adopted. This suggests that it is most probably a giant circumbinary planet orbiting NY Vir at a distance of about 3.3 astronomical units (AU). Since the rate of period decrease can not be explained by true angular momentum loss caused by gravitational radiation or/and magnetic braking, the observed downward parabolic change in the O-C diagram may be only a part of a long-period (longer than 15 years) cyclic variation, which may reveal the presence of another Jovian planet ($\sim2.5$$M_{Jupiter}$) in the system.

Gamma-ray bursts (GRBs) are attractive objects for constraining the nature of dark energy in a way complementary to other cosmological probes, especially at high redshifts. However, the apparent magnitude of distant GRBs can be distorted by the gravitational lensing from the density fluctuations along the line of sight. We investigate the gravitational lensing effect on the cosmological parameters and dark energy equation of state from GRBs. We first calibrated the GRB luminosity relations without assuming any cosmological models. The luminosity distances of low-redshift GRBs were calibrated with the cosmography method using a latest type Ia supernova (SNe Ia) sample. We find that the gravitational lensing has non-negligible effects on the determination of cosmological parameters and dark energy. The gravitational lensing shifts the best-fit constraints on cosmological parameters and dark energy. Because high-redshift GRBs are more likely to be reduced, the most probable value of the observed matter density $\Omega_M$ is slightly lower than its actual value. In the $\Lambda$CDM model, we find that the matter density parameter $\Omega_M$ will shift from 0.30 to 0.26 after including the gravitational lensing effect. The gravitational lensing also affects the dark energy equation of state by shifting it to a more negative value. We constrain the dark energy equation of state out to redshift $z\sim 8$ using GRBs for the first time, and find that the equation of state deviates from $\Lambda$CDM at the $1\sigma$ confidence level, but agrees with $w=-1$ at the $2\sigma$ confidence level.

In an effort to understand the puzzle of classifying gamma-ray bursts (GRBs), we perform a systematic study of {\it Swift} GRBs and investigate several short GRB issues. Though short GRBs have a short ($\lesssim2$ s) prompt duration as monitored by the Burst Alert Telescope, the composite light curves including both the prompt and afterglow emission suggest that most of the short GRBs have a similar radiative feature to long GRBs. Further, some well-studied short GRBs might also have an intrinsically long prompt duration, which renders them as a type of short GRB imposters. Genuine short GRBs detected by {\it Swift} might be rare that discriminating the observed short GRBs is, not surprisingly, troublesome. In particular, the observational biases in the host identification and redshift measurement of GRBs should be taken with great caution. The redshift distribution which has been found to be different for long and short GRBs might have been strongly affected by the measurement methods. We find that the redshifts measured from the presumed host galaxies of long and short GRBs appear to have a similar distribution.

Low-luminosity gamma-ray bursts (GRBs) with luminosity . 10^49erg/s probably consititute a distinct population from the classic high-luminosity GRBs. They are the most luminous objects detected so far within ~ 100 Mpc, the horizon distance of ultra-high energy cosmic rays (UHECRs), so they are considered to be candidate sources of UHECRs. It was recently argued that the energy production rate in UHECRs is much larger than that in gamma-ray photons of long GRBs measured by the Fermi satellite, which, if true, would challenge the view that GRBs can be the sources of UHECRs. We here suggest that many of the low luminosity GRBs, due to their low luminosity, can not trigger the current GRB detectors and hence their contribution to the local gamma-ray energy production rate is missing. We find that the real local energy production rate by low-luminosity GRBs, taking into account the missing part, which constitutes a dominant fraction of the total amount, could be sufficient to account for the flux of UHECRs. Due to the low-luminosity, only intermediate-mass or heavy nuclei can be accelerated to ~ 10^20 eV. We discuss the acceleration and survival of these UHE nuclei in low-luminosity GRBs, especially in those missing low-luminosity GRBs. At last, the accompanying diffuse neutrino flux from the whole low-luminosity GRB population is calculated.

The recently discovered high-energy transient Swift J164449.3+573451 (Sw J1644+57) is thought to arise from the tidal disruption of a passing star by a dormant massive black hole. Modeling of the broadband emission suggests the presence of a powerful relativistic jet, which contributes dominantly to the observed X-ray emission. Here we suggest that protons can be accelerated to ultra-high energies by internal shocks occurring in the jets, but their flux is insufficient to account for the observed flux of ultra-high energy cosmic rays. High energy protons can produce ~0.1-10 PeV neutrinos through photomeson interactions with X-ray photons. The large X-ray fluence (7x10^-4 erg cm^-2) and high photopion efficiency, together with the insignificant cooling of secondary mesons, result in bright neutrino emission expected from Sw J1644+57 if the jet composition is matter-dominated. One to several neutrinos may be detected by a Km^3-scale detector from one tidal disruption event similar to Sw J1644+57, thereby providing a powerful probe of the composition of the jets.

We present broadband (radio, optical, X-ray and GeV) fits to the afterglow light curves and spectra of three long-duration gamma-ray bursts (GRBs 080916C, 090902B, and 090926A) detected by the Gamma-Ray Burst Monitor (GBM) and Large Area Telescope (LAT) instruments on the Fermi satellite. Using the observed broadband data, we study the origin of the high energy emission, and suggest that the early-time GeV emission and the late-time radio, optical, and X-ray afterglows can be understood as being due to synchrotron emission from an external forward shock caused by structured ejecta propagating in a wind bubble jumping to a homogeneous density medium. If the ceasing time for majority of the energy injection is assumed to be close to the deceleration time of the forward shock, the structured ejecta with continuous energy injection to the forward shock can well explain the early rising feature of the GeV mission from these burst, and the density-jump medium can account for some certain plateaus or flares in the late afterglows. From our fits, we find that, on one hand, the external shock origin of the GeV photons will make the optical depth have not significant contribution to the early LAT rising part, which will loosen strong constraint of lower limits of Lorentz factor. On the other hand, these Fermi-LAT events preferentially occur in a low-density circumburst environment, in which case the Klein-Nishina cutoff will significantly suppress the Self-Synchrotron Compton (SSC) radiation. Such an environment might result from superbubbles or low-metallicity progenitor stars (which have a low mass-loss rate at late times of stellar evolution) of type Ib/c supernovae.

Several interesting luminosity correlations among gamma-ray burst (GRB) variables have been recently discussed extensively. In this paper, we derive the six luminosity correlations ($\tlag-L$, $V-L$, $\epkk-L$, $\epkk-E_\gamma$, $\trt-L$, $\epkk-E_{\gamma, \mathrm{iso}}$) from the light curves and spectra of the latest 116 long GRBs, including the time lag ($\tlag$) between low and high photon energy light curves, the variability ($V$) of the light curve, the peak energy of the spectrum ($\epkk$), and the minimum rise time ($\trt$) of the peaks. We find that the intrinsic scatter of the $V-L$ correlation is too large and there seems no inherent correlation between the two parameters using the latest GRB data. The other five correlations indeed exist when the sample is enlarged. The $\epkk-E_\gamma$ correlation has a significantly lower intrinsic scatter compared to the other correlations. We divide the full data into four redshift bins when testing possible evolution of the correlations with redshift. We find no statistically significant evidence for the redshift evolution of the luminosity correlations. To avoid the circularity problem when constraining the cosmological parameters, we simultaneously minimize $\chi^2$ with respect to both correlation parameters $a$, $b$ and the cosmological parameters using the maximum likelihood method. For the flat $\Lambda$CDM, the best fit is $\omm=0.31^{+0.13}_{-0.10}$. We also constrain the possible evolution of the equation of state (EOS) of the dark energy using the GRBs together with the Union2 compilation of SNe Ia and the $H(z)$ data. The result is consistent with the cosmological constant at $2 \sigma$ confidence level and mainly due to the GRB data, the dark energy EOS shows slight deviation from -1 at $z \geq0.5$ as was persistently presented with many previous data sets.

Using the precise times of mid-egress of the eclipsing polar HU Aqr, we discovered that this polar is orbited by two or more giant planets. The two planets detected so far have masses of at least 5.9 and 4.5\,M_{Jup}. Their respective distances from the polar are 3.6 AU and 5.4 AU with periods of 6.54 and 11.96 years, respectively. The observed rate of period decrease derived from the downward parabolic change in O-C curve is a factor 15 larger than the value expected for gravitational radiation. This indicates that it may be only a part of a long-period cyclic variation, revealing the presence of one more planet. It is interesting to note that the two detected circumbinary planets follow the Titus-Bode law of solar planets with n=5 and 6. We estimate that another 10 years of observations will reveal the presence of the predicted third planet.

Recent studies suggest that Swift gamma-ray bursts (GRBs) may not trace an ordinary star formation history. Here we show that the GRB rate turns out to be consistent with the star formation history with an evolving stellar initial mass function (IMF). We first show that the latest Swift sample of GRBs reveals an increasing evolution in the GRB rate relative to the ordinary star formation rate at high redshifts. We then assume only massive stars with masses greater than the critical value to produce GRBs, and use an evolving stellar IMF suggested by Dav\'{e} (2010) to fit the latest GRB redshift distribution. This evolving IMF would increase the relative number of massive stars, which could lead to more GRB explosions at high redshifts. We find that the evolving IMF can well reproduce the observed redshift distribution of Swift GRBs.

In this paper we present the spectrum of synchro-curvature self-Compton (SCSC) radiation of relativistic electrons with a power-law distribution of Lorentz factors. We find that the resulting spectrum is significantly different from that of either synchrotron self-Compton or curvature self-Compton radiation if both the curvature radius of the magnetic field and the cyclotron radius of the electrons are within some proper ranges. The effects of electrons' cooling and drifting, the low-energy self absorption in seed spectra, and the Klein-Nishina cutoff are also discussed, in order to get an accurate picture. We take gamma-ray bursts (GRBs) as our example environment for discussions. The results would be considered as a universal approach of the self-Compton emission of relativistic electrons moving in curved magnetic fields, and thus could be applied to many astrophysical phenomena, including GRBs, active galactic nuclei (AGNs), and pulsars.

Hyperaccreting disks around neutron stars or magnetars cooled via neutrino emission can be the potential central engine of GRBs. The neutron-star disk can cool more efficiently, produce much higher neutrino luminosity and neutrino annihilation luminosity than its black hole counterpart with the same accretion rate. The neutron star surface boundary layer could increase the annihilation luminosity as well. An ultra relativistic jet via neutrino annihilation can be produced along the stellar poles. Moreover, we investigate the effects of strong fields on the disks around magnetars. In general, stronger fields give higher disk densities, pressures, temperatures and neutrino luminosity; the neutrino annihilation mechanism and the magnetically-driven pulsar wind which extracts the stellar rotational energy can work together to generate and feed an even stronger ultra-relativistic jet along the stellar magnetic poles.

Context: Swift observations suggest that the X-ray afterglow emission of some gamma-ray bursts (GRB) may have internal origins, and the conventional external shock (ES) cannot be the exclusive source of the afterglow emission. Aims: If the central compact objects of some GRBs are millisecond magentars, the magnetar winds could play an important role in the (internal) X-ray afterglow emission, which is our focus here. Methods: The dynamics and the synchrotron radiation of the termination shock (TS) of the magmnetar winds, as well as the simultaneous GRB ES, are investigated by considering the magnetization of the winds. Results: As a result of the competition between the emission of the wind TS and the GRB ES, two basic types of X-ray afterglows are predicted, i.e., the TS-dominated and the ES-dominated types. Moreover, our results also show that both of the two types of afterglows have a shallow-decay phase and a normal-decay one, as observed by the \textit{Swift} satellite. This indicates that some observed X-ray afterglows could be (internally) produced by the magnetar winds, but not necessarily GRB ESs.

Combining with our newest CCD times of light minimum of EM Cygni, all 45 available times of light minimum including 7 data with large scatters are compiled and the updated O-C analysis is made. The bestfit for the O-C diagram of EM Cygni is a quadratic-plus-sinusoidal fit. The secular orbital period decrease rate -2.5(\pm 0.3)x10^{-11} s s^{-1} means that magnetic braking effect with a rate of mass loss via stellar wind, 2.3x10^{-10}Msunyr^{-1}, is needed for explaining the observed orbital period decrease. Moreover, for explaining the significant cyclical period change with a period of \sim 17.74(\pm 0.01)yr shown in the O-C diagram, magnetic activity cycles and light travel-time effect are discussed in detail. The O-C diagram of EM Cygni cannot totally rule the possibility of multi-periodic modulation out due to the gaps presented after 25000 cycles. Based on the hypothesis of a K-type third star in literature, light trave-time effect may be a more plausible explanation. However, the low orbital inclination of the third body (\sim 7.4 degree) suggests that the hypothetic K-type third star may be captured by EM Cygni. But assuming the spectral contamination from a block of circumbinary material instead of a K-type third star, the third star may be a brown dwarf in case of the coplanar orbit with parent binary.

Four newest CCD eclipse timings of the white dwarf for polar UZ Fornacis and Six updated CCD mid-eclipse times for SW Sex type nova-like V348 Puppis are obtained. The detailed O-C analyses for both CVs inside period gap are made. Orbital period increases at a rate of $2.63(\pm0.58)\times10^{-11} s\;s^{-1}$ for UZ Fornacis and of $5.8(\pm1.9)\times10^{-12} s\;s^{-1}$ for V348 Puppis, respectively, are discovered in their new O-C diagrams. However, the conservative mass transfer from the secondary to the massive white dwarf cannot explain the observed orbital period increases for both CVs, which are regarded as part of modulations at longer periods. Moreover, the O-C diagram of UZ Fornacis shows a possible cyclical change with a period of $\sim23.4(\pm5.1)yr$. For explaining the observed cyclical period changes in UZ Fornacis, both mechanisms of magnetic activity cycles in the late-type secondary and the light travel-time effect are regarded as two probable causes. Not only does the modulation period 23.4yr obey the empirical correlation derived by \cite{lan99}, but also the estimated fractional period change $\Delta P/P\sim7.3\times10^{-7}$ displays a behavior similar to that of the CVs below the period gap. On the other hand, a calculation for the light travel-time effect implies that the tertiary component in UZ Fornacis may be a brown dwarf with a high confidence level, when the orbital inclination of the third body is larger than $16^{\circ}$.

The two CCD photometries of the intermediate polar TV Columbae are made for obtaining the two updated eclipse timings with high precision. There is an interval time \sim 17yr since the last mid-eclipse time observed in 1991. Thus, the new mid-eclipse times can offer an opportunity to check the previous orbital ephemerides. A calculation indicates that the orbital ephemeris derived by Augusteijn et al. (1994) should be corrected. Based on the proper linear ephemeris (Hellier, 1993), the new orbital period analysis suggests a cyclical period variation in the O-C diagram of TV Columbae. Using Applegate's mechanism to explain the periodic oscillation in O-C diagram, the required energy is larger than that a M0-type star can afford over a complete variation period \sim 31.0(\pm 3.0)yr. Thus, the light travel-time effect indicates that the tertiary component in TV Columbae may be a dwarf with a low mass, which is near the mass lower limit \sim 0.08Msun as long as the inclination of the third body high enough.

(Abridged) The hyperaccreting neutron star or magnetar disks cooled via neutrino emission can be a candidate of gamma-ray burst (GRB) central engines. The strong field $\geq10^{15}-10^{16}$ G of the magnetar can play a significant role in affecting the disk properties and even lead to the funnel accretion process. We investigate the effects of strong fields on the disks around magnetars, and discuss implications of such accreting magnetar systems for GRB and GRB-like events. We discuss quantum effects of the strong fields on the disk, and use the MHD conservation equations to describe the behavior of the disk flow coupled with a large scale field, which is generated by the star-disk interaction. In general, stronger fields give higher disk densities, pressures, temperatures and neutrino luminosity, and change the electron fraction and degeneracy state significantly. A magnetized disk is always viscously stable outside the Alfv\'{e}n radius, but will be thermally unstable near the Alfv\'{e}n radius where the magnetic field plays a more important role in transferring the angular momentum and heating the disk than the viscous stress. The funnel accretion process will be only important for an extremely strong field, which creates a magnetosphere inside the Alfv\'{e}n radius and truncates the plane disk. Because of higher temperature and more concentrated neutrino emission of the magnetar surface ring-like belt region covered by funnel accretion, the neutrino annihilation rate from the accreting magnetars can be much higher than that from accreting neutron stars without fields. Furthermore, the neutrino annihilation mechanism and the magnetically-driven pulsar wind from the magnetar surface can work together to generate and feed an ultra-relativistic jet along the stellar magnetic poles.

A model of jet precession driven by a neutrino-cooled disc around a spinning black hole is present in order to explain the temporal structure and spectral evolution of gamma-ray bursts (GRBs). The differential rotation of the outer part of a neutrino dominated accretion disc may result in precession of the inner part of the disc and the central black hole, hence drives a precessed jet via neutrino annihilation around the inner part of the disc. Both analytic and numeric results for our model are present. Our calculations show that a black hole-accretion disk system with black hole mass $M \simeq 3.66 M_\odot$, accretion rate $\dot{M} \simeq 0.54 M_\odot \rm s^{-1}$, spin parameter $a=0.9$ and viscosity parameter $\alpha=0.01$ may drive a precessed jet with period P=1 s and luminosity $L=10^{51}$ erg s$^{-1}$, corresponding to the scenario for long GRBs. A precessed jet with $P=0.1$s and $L=10^{50}$ erg s$^{-1}$ may be powered by a system with $M \simeq 5.59 M_\odot$, $\dot{M} \simeq 0.74 M_\odot \rm s^{-1}$, $a=0.1$, and $\alpha=0.01$, possibly being responsible for the short GRBs. Both the temporal and spectral evolution in GRB pulse may explained with our model. GRB central engines likely power a precessed jet driven by a neutrino-cooled disc. The global GRB lightcurves thus could be modulated by the jet precession during the accretion timescale of the GRB central engine. Both the temporal and spectral evolution in GRB pulse may be due to an viewing effect due to the jet precession.

In this paper we investigate the deceleration, jerk and snap parameters to distinguish between the dark energy and modified gravity models by using high redshift gamma-ray bursts (GRBs) and supernovae (SNe). We first derive the expressions of deceleration, jerk and snap parameters in dark energy and modified gravity models. In order to constrain the cosmographic parameters, we calibrate the GRB luminosity relations without assuming any cosmological models using SNe Ia. Then we constrain the models (including dark energy and modified gravity models) parameters using type Ia supernovae and gamma-ray bursts. Finally we calculate the cosmographic parameters. GRBs can extend the redshift - distance relation up to high redshifts, because they can be detected to high redshifts.We find that the statefinder pair (r,s) could not be used to distinguish between some dark energy and modified gravity models, but these models could be differentiated by the snap parameter. Using the model-independent constraints on cosmographic parameters, we conclude that the \LambdaCDM model is consistent with the current data.

The high-redshift star formation rate (SFR) is difficult to measure directly even by modern approaches. Long-duration gamma-ray bursts (GRBs) can be detected to the edge of the visible universe because of their high luminorsities. The collapsar model of long gamma-ray bursts indicates that they may trace the star formation history. So long gamma-ray bursts may be a useful tool of measuring the high-redshift SFR. Observations show that long gamma-ray bursts prefer to form in a low-metallicity environment. We study the high-redshift SFR up to z~8.3 considering the Swift GRBs tracing the star formation history and the cosmic metallicity evolution in different background cosmological models including $\Lambda$CDM, quintessence, quintessence with a time-varying equation of state, and brane-world model. We use latest Swift GRBs including two highest-$z$ GRBs, GRB 080913 at $z=6.7$ and GRB 090423 at $z=8.3$. We find that the SFR at $z>4$ shows a steep decay with a slope of $\sim -5.0$ in $\Lambda$CDM. In the other three models, the high-redshift SFR is slightly different from $\Lambda$CDM model, and also shows a steep decay.

We revisit the vertical structure of neutrino-dominated accretion flows in spherical coordinates. We stress that the flow should be geometrically thick when advection becomes dominant. In our calculation, the luminosity of neutrino annihilation is enhanced by one or two orders of magnitude. The empty funnel along the rotation axis can naturally explain the neutrino annihilable ejection.

Extended high-energy(>100MeV) gamma-ray emission that lasts much longer than the prompt sub-MeV emission has been detected from quite a few gamma-ray bursts (GRBs) by Fermi Large Area Telescope (LAT) recently. A plausible scenario is that this emission is the afterglow synchrotron emission produced by electrons accelerated in the forward shocks. In this scenario, the electrons that produce synchrotron high-energy emission also undergo inverse-Compton (IC) loss and the IC scattering with the synchrotron photons should be in the Klein-Nishina regime. Here we study effects of the Klein-Nishina scattering on the high-energy synchrotron afterglow emission. We find that, at early times the Klein-Nishina suppression effect on those electrons that produce the high-energy emission is usually strong and therefore their inverse-Compton loss is small with a Compton parameter Y < a few for a wide range of parameter space. This leads to a relatively bright synchrotron afterglow at high energies that can be detected by Fermi LAT. As the Klein-Nishina suppression effect weakens with time, the inverse-Compton loss increases and could dominate over the synchrotron loss in some parameter space. This will lead to a faster temporal decay of the high-energy synchrotron emission than what is predicted by the standard synchrotron model, which may explain the observed rapid decay of the early high-energy gamma-ray emission in GRB090510 and GRB090902B.

Existing theory and models suggest that a Type I (merger) GRB should have a larger jet beaming angle than a Type II (collapsar) GRB, but so far no statistical evidence is available to support this suggestion. In this paper, we obtain a sample of 37 beaming angles and calculate the probability that this is true. A correction is also devised to account for the scarcity of Type I GRBs in our sample. The probability is calculated to be 83% without the correction and 71% with it.

Newborn neutron stars surrounded by hyperaccreting and neutrino-cooled disks may exist in some gamma-ray bursts (GRBs) and/or supernovae (SNe). In this paper we further study the structure of such a neutron-star disk based on the two-region (i.e., inner & outer) disk scenario following our previous work, and calculate the neutrino annihilation luminosity from the disk in various cases. We investigate the effects of the viscosity parameter, energy parameter (measuring the neutrino cooling efficiency of the inner disk) and outflow strength on the structure of the entire disk as well as the effect of emission from the neutron star surface boundary emission on the total neutrino annihilation rate. The inner disk satisfies the entropy-conservation or the advection-dominated self-similar structure depending on the energy parameter. An outflow from the disk decreases the density and pressure but increases the thickness of the disk. Moreover, compared with the black-hole disk, the neutrino annihilation luminosity above the neutron-star disk is higher, and the neutrino emission from the boundary layer could increase the neutrino annihilation luminosity by about one order of magnitude higher than the disk without boundary emission. The neutron-star disk with the advection-dominated inner disk could produce the highest neutrino luminosity while the disk with an outflow has the lowest. Although a heavily mass-loaded outflow from the neutron star surface at early times of neutron star formation prevents the outflow material from being accelerated to a high bulk Lorentz factor, an energetic ultrarelativistic jet via neutrino annihilation can be produced above the stellar polar region at late times if the disk accretion rate and the neutrino emission luminosity from the surface boundary layer are sufficiently high.

Fermi observations of high-energy gamma-ray emission from GRB 080916C shows that its spectrum is consistent with the Band function from MeV to tens of GeV. Assuming one single emission mechanism dominates in the whole energy range, we show that this spectrum is consistent with synchrotron origin by shock-accelerated electrons. The simple electron inverse-Compton model and the hadronic model are found to be less viable. In the synchrotron scenario, the synchrotron self-Compton scattering is likely to be in the Klein-Nishina regime and therefore the resulting high-energy emission is subdominant, even though the magnetic field energy density is lower than that in relativistic electrons. The Klein-Nishina inverse-Compton cooling may also affect the low-energy electron number distribution and hence results in a low-energy synchrotron photon spectrum $n(\nu)\propto\nu^{-1}$ below the peak energy. Under the framework of the electron synchrotron interpretation, we constrain the shock microphysical parameters and derive a lower limit of the upstream magnetic fields. The detection of synchrotron emission extending to about 70 GeV in the source frame in GRB 080916C favors the Bohm diffusive shock acceleration if the bulk Lorentz factor of the relativistic outflow is not significantly greater than thousands.

Recently, it was suggested that a photospheric component that results from the internal dissipation occurring in the optically thick inner parts of relativistic outflows may be present in the prompt $\gamma$/X-ray emission of gamma-ray bursts or X-ray flashes. We explore high-energy neutrino emission in this dissipative photosphere model, assuming that the composition of the outflow is baryon-dominated. We find that neutrino emission from proton-proton collision process forms an interesting signature in the neutrino spectra. Under favorable conditions for the shock dissipation site, these low-energy neutrinos could be detected by ${\rm km^3}$ detectors, such as Icecube. Higher energies ($\ga10$ TeV) neutrino emission from proton-proton collision and photo-pion production processes could be significantly suppressed for dissipation at relatively small radii, due to efficient Bethe-Heitler cooling of protons and/or radiative cooling of the secondary mesons in the photosphere radiation. As the dissipation shocks continue further out, high energy neutrinos from photo-pion production process becomes dominant.

In the popular internal shock model for the prompt emission of gamma-ray bursts (GRBs), collisions between a series of relativistic shells generate lots of paired forward and reverse shocks. We show that the synchrotron emission produced by the forward and reverse shocks respectively could peak at two quite different energy bands if the Lorentz factors of these two types of shocks are significantly different with each other (e.g., one shock is relativistic and the other is Newtonian). We then investigate whether this scenario is applicable to the case of GRB 080319B and find that a bimodal distribution of the shell Lorentz factors, peaking at $\sim400$ and $\sim10^5$, is required. In addition, this scenario predicts an accompanying inverse-Compton (IC) GeV emission with a luminosity comparable to (not much higher than) that of the synchrotron MeV emission, which can be tested with future \textit{Fermi} observations.

We study afterglow flares of gamma-ray bursts (GRBs) in the framework of the late internal shock (LIS) model based on a careful description for the dynamics of a pair of shocks generated by a collision between two homogeneous shells,. First, by confronting the model with some fundamental observational features of X-ray flares, we find some constraints on the properties of the pre-collision shells that are directly determined by the central engine of GRBs. Second, high energy emission associated with X-ray flares, which arises from synchrotron self-Compton (SSC) emission of LISs, is investigated in a wide parameter space. The predicted flux of high energy flares may reach as high as $\sim 10^{-8}\rm erg cm^{-2}s^{-1}$, which is likely to be detectable with the Large Area Telescope (LAT) aboard \textit{the Fermi Space Telescope}

Lorentz invariant violation (LIV) test is very important to study in the new physics. All the known astrophysical constraints either have a very small examinable parameter space, or are only suitable for some special theoretical models. Here we suggest that it is possible to detect the time-delay of ultra-high-energy cosmic-rays (UHECRs) directly. We discuss some difficulties in our method, including the intergalactic magnetic fields. It seems that none of them are crucial, hence this method could give a larger examinable parameter space and a stronger constraint on LIV.

(Abridged) We study the effects of a global magnetic field on viscously-rotating and vertically-integrated accretion disks around compact objects using a self-similar treatment. We extend Akizuki & Fukue's work (2006) by discussing a general magnetic field with three components ($r, \phi, z$) in advection-dominated accretion flows (ADAFs). We also investigate the effects of a global magnetic field on flows with convection. For these purposes, we first adopt a simple form of the kinematic viscosity $\nu=\alpha c_{s}^{2}/\Omega_{K}$ to study magnetized ADAFs. Then we consider a more realistic model of the kinematic viscosity $\nu=\alpha c_{s}H$, which makes the infall velocity increase but the sound speed and toroidal velocity decrease. We next use two methods to study magnetized flows with convection, i.e., we take the convective coefficient $\alpha_{c}$ as a free parameter to discuss the effects of convection for simplicity. We establish the $\alpha_{c}-\alpha$ relation for magnetized flows using the mixing-length theory and compare this relation with the non-magnetized case. If $\alpha_{c}$ is set as a free parameter, then $|v_{r}|$ and $c_{s}$ increase for a large toroidal magnetic field, while $|v_{r}|$ decreases but $|v_{\phi}|$ increases (or decreases) for a strong and dominated radial (or vertical) magnetic field with increasing $\alpha_{c}$. In addition, the magnetic field makes the $\alpha_{c}-\alpha$ relation be distinct from that of non-magnetized flows, and allows the $\rho\propto r^{-1}$ or $\rho\propto r^{-2}$ structure for magnetized non-accreting convection-dominated accretion flows with $\alpha+g\alpha_{c}< 0$ (where $g$ is the parameter to determine the condition of convective angular momentum transport).

(Abridged) We here study the structure of a hyperaccretion disk around a neutron star. We consider a steady-state hyperaccretion disk around a neutron star, and as a reasonable approximation, divide the disk into two regions, which are called inner and outer disks. The outer disk is similar to that of a black hole and the inner disk has a self-similar structure. In order to study physical properties of the entire disk clearly, we first adopt a simple model, in which some microphysical processes in the disk are simplified, following Popham et al. and Narayan et al. Based on these simplifications, we analytically and numerically investigate the size of the inner disk, the efficiency of neutrino cooling, and the radial distributions of the disk density, temperature and pressure. We see that, compared with the black-hole disk, the neutron star disk can cool more efficiently and produce a much higher neutrino luminosity. Finally, we consider an elaborate model with more physical considerations about the thermodynamics and microphysics in the neutron star disk (as recently developed in studying the neutrino-cooled disk of a black hole), and compare this elaborate model with our simple model. We find that most of the results from these two models are basically consistent with each other.

The observations of a nearby low-luminosity gamma-ray burst (GRB) 060218 associated with supernova SN 2006aj may imply an interesting astronomical picture where a supernova shock breakout locates behind a relativistic GRB jet. Based on this picture, we study neutrino emission for early afterglows of GRB 060218-like GRBs, where neutrinos are expected to be produced from photopion interactions in a GRB blast wave that propagates into a dense wind. Relativistic protons for the interactions are accelerated by an external shock, while target photons are basically provided by the incoming thermal emission from the shock breakout and its inverse-Compton scattered component. Because of a high estimated event rate of low-luminosity GRBs, we would have more opportunities to detect afterglow neutrinos from a single nearby GRB event of this type by IceCube. Such a possible detection could provide evidence for the picture described above.

We investigate the effect of X-ray echo emission in gamma-ray bursts (GRBs). We find that the echo emission can provide an alternative way of understanding X-ray shallow decays and jet breaks. In particular, a shallow decay followed by a "normal" decay and a further rapid decay of X-ray afterglows can be together explained as being due to the echo from prompt X-ray emission scattered by dust grains in a massive wind bubble around a GRB progenitor. We also introduce an extra temporal break in the X-ray echo emission. By fitting the afterglow light curves, we can measure the locations of the massive wind bubbles, which will bring us closer to finding the mass loss rate, wind velocity, and the age of the progenitors prior to the GRB explosions.

In this paper we analyze the constraints on the property of dark energy from cosmological observations. Together with SNe Ia Gold sample, WMAP, SDSS and 2dFGRS data, we include 69 long Gamma-Ray Bursts (GRBs) data in our study and perform global fitting using Markov Chain Monte Carlo (MCMC) technique. Dark energy perturbations are explicitly considered. We pay particular attention to the time evolution of the equation of state of dark energy parameterized as $w_{DE}=w_0+w_a(1-a)$ with $a$ the scale factor of the universe, emphasizing the complementarity of high redshift GRBs to other cosmological probes. It is found that the constraints on dark energy become stringent by taking into account high redshift GRBs, especially for $w_a$, which delineates the evolution of dark energy.

The origin of the ultrahigh-energy (UHE) cosmic rays (CRs) from the second knee ($\sim6\times10^{17}$eV) above in the CR spectrum is still unknown. Recently, there has been growing evidence that a peculiar type of supernovae, called hypernovae, are associated with sub-energetic gamma-ray bursts (GRBs), such as SN1998bw/GRB980425 and SN2003lw/GRB031203. Such hypernovae appear to have high (up to mildly relativistic) velocity ejecta, which may be linked to the sub-energetic GRBs. Assuming a continuous distribution of the kinetic energy of the hypernova ejecta as a function of its velocity $E_k\propto (\Gamma\beta)^{-\alpha}$ with $\alpha\sim 2$, we find that 1) the external shock wave produced by the high velocity ejecta of a hypernova can accelerate protons up to energies as high as $10^{19} {\rm eV}$; 2) the cosmological hypernova rate is sufficient to account for the energy flux above the second knee; and 3) the steeper spectrum of CRs at these energies can arise in these sources. In addition, hypernovae would also give rise to a faint diffuse UHE neutrino flux, due to $p\gamma$ interactions of the UHE CRs with hypernova optical-UV photons.

The property of dark energy and the physical reason for acceleration of the present universe are two of the most difficult problems in modern cosmology. The dark energy contributes about two-thirds of the critical density of the present universe from the observations of type-Ia supernova (SNe Ia) and anisotropy of cosmic microwave background (CMB).The SN Ia observations also suggest that the universe expanded from a deceleration to an acceleration phase at some redshift, implying the existence of a nearly uniform component of dark energy with negative pressure. We use the ``gold'' sample containing 157 SNe Ia and two recent well-measured additions, SNe Ia 1994ae and 1998aq to explore the properties of dark energy and the transition redshift. For a flat universe with the cosmological constant, we measure $\Omega_{M}=0.28_{-0.05}^{+0.04}$, which is consistent with Riess et al. The transition redshift is $z_{T}=0.60_{-0.08}^{+0.06}$. We also discuss several dark energy models that define the $w(z)$ of the parameterized equation of state of dark energy including one parameter and two parameters ($w(z)$ being the ratio of the pressure to energy density). Our calculations show that the accurately calculated transition redshift varies from $z_{T}=0.29_{-0.06}^{+0.07}$ to $z_{T}=0.60_{-0.08}^{+0.06}$ across these models. We also calculate the minimum redshift $z_{c}$ at which the current observations need the universe to accelerate.

The widely existing shallow decay phase of the X-ray afterglows of gamma-ray bursts (GRBs) is generally accepted to be due to long-lasting energy injection. The outflows carrying the injecting energy, based on the component that is dominative in energy, fall into two possible types: baryon-dominated and lepton-dominated ones. The former type of outflow could be ejecta that is ejected during the prompt phase of a GRB and consists of a series of baryonic shells with a distribution of Lorentz factors, and the latter type could be an electron-positron-pair wind that is driven by the post-burst central engine. We here provide a unified description for the dynamics of fireballs based on these two types of energy injection, and calculate the corresponding high-energy photon emission by considering synchrotron radiation and inverse Compton scattering (including synchrotron self-Compton and combined inverse-Compton) of electrons. We find that, in the two energy-injection models, there is a plateau (even a hump) in high-energy light curves during the X-ray shallow decay phase. In particular, a considerable fraction of the injecting energy in the lepton-dominated model can be shared by the long-lasting reverse shock since it is relativistic. Furthermore, almost all of the energy of the reverse shock is carried by leptons, and thus the inverse-Compton emission is enhanced dramatically. Therefore, this model predicts more significant high-energy afterglow emission than the baryon-dominated model. We argue that these observational signatures would be used to discriminate between different energy-injection models in the upcoming {\em Gamma-ray Large Area Space Telescope} (GLAST) era.

Swift/BAT has detected ~ 200 long-duration GRBs, with redshift measurements for ~50 of them. We derive the luminosity function (Phi^{HL}) and the local event rate (rho_0^{HL}) of the conventional high luminosity (HL) GRBs by using the z-known Swift GRBs. Our results are generally consistent with that derived from the CGRO/BATSE data. However, the fact that Swift detected a low luminosity (LL) GRB, GRB 060218, at z=0.033 within ~ 2 year of operation, together with the previous detection of the nearby GRB 980425, suggests a much higher local rate for these LL-GRBs. We explore the possibility that LL-GRBs as a distinct GRB population from the HL-GRBs. We find that rho_0^{LL} is 325_{-177}^{+352} Gpc^{-3} yr^{-1}, which is much higher than rho_0^{HL}(1.12_{-0.20}^{+0.43} Gpc^{-3} yr^{-1}). This rate is ~ 0.7% of the local Type Ib/c SNe. Our results, together with the finding that less than 10% of Type Ib/c SNe are associated with off-beam GRBs, suggest that LL-GRBs have a beaming factor typically less than 14, or a jet angle typically wider than 31^{o}. The high local GRB rate, the small beaming factor, and low luminosity make the LL-GRBs distinct from the HL-GRBs. Although the current data could not fully rule out the possibility that both HL- and LL-GRBs are the same population, our results suggest that LL-GRBs are likely a unique GRB population and the observed low redshift GRB sample is dominated by the LL-GRBs.

It has been widely shown that the cosmological parameters and dark energy can be constrained by using data from type-Ia supernovae (SNe Ia), the cosmic microwave background (CMB) anisotropy, the baryon acoustic oscillation (BAO) peak from Sloan Digital Sky Survey (SDSS), the X-ray gas mass fraction in clusters, and the linear growth rate of perturbations at z=0.15 as obtained from the 2dF Galaxy Redshift Survey. Recently, gamma-ray bursts (GRBs) have also been argued to be promising standard candles for cosmography. In this paper, we present constraints on the cosmological parameters and dark energy by combining a recent GRB sample including 69 events with the other cosmological probes. First, we find that for the LambdaCDM cosmology this combination makes the constraints stringent and the best fit is close to the flat universe. Second, we fit the flat Cardassian expansion model and find that this model is consistent with the LambdaCDM cosmology. Third, we present constraints on several two-parameter dark energy models and find that these models are also consistent with the LambdaCDM cosmology. Finally, we reconstruct the dark energy equation-of-state parameter w(z) and the deceleration parameter q(z). We see that the acceleration could have started at a redshift from z_T=0.40_{-0.08}^{+0.14} to z_T=0.65_{-0.05}^{+0.10}. This difference in the transition redshift is due to different dark energy models that we adopt. The most stringent constraint on w(z) lies in the redshift range z\sim 0.3-0.6.

The postburst object of a GRB is likely to be a highly magnetized, rapidly rotating compact object (e.g., a millisecond magnetar), which could produce an ultrarelativistic electron-positron-pair wind. The interaction of such a wind with an outwardly expanding fireball ejected during the burst leads to a relativistic wind bubble (RWB). We numerically calculate the dynamics and radiative properties of RWBs and use this model to explain the shallow decay phase of the early X-ray afterglows observed by Swift. We find that RWBs can fall into two types: forward-shock-dominated and reverse-shock-dominated bubbles. Their radiation during a period of $\sim 10^{2}-10^{5}$ seconds is dominated by the shocked medium and the shocked wind, respectively, based on different magnetic energy fractions of the shocked materials. For both types, the resulting light curves always have a shallow decay phase. In addition, we provide an example fit to the X-ray afterglows of GRB 060813 and GRB 060814 and show that they could be produced by forward-shock-dominated and reverse-shock-dominated bubbles, respectively. This implies that, for some early afterglows (e.g., GRB 060814), the long-lasting reverse shock emission is strong enough to explain their shallow decay phase.

The afterglows of gamma-ray bursts (GRBs) have commonly been assumed to be due to shocks sweeping up the circum-stellar medium. However, most GRBs have been found in dense star-forming regions where a significant fraction of the prompt X-ray emission can be scattered by dust grains. Here we revisit the behavior of dust scattering of X-rays in GRBs. We find that the features of some X-ray afterglows from minutes to days after the gamma-ray triggers are consistent with the scattering of prompt X-ray emission from GRBs off host dust grains. This implies that some of the observed X-ray afterglows (especially those without sharp rising and decaying flares) could be understood with a dust-scattering--driven emission model.

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.

The known encounter velocity (6.1kms-1) and particle incidence angle (perpendicular) between the Stardust spacecraft and the dust emanating from the nucleus of comet Wild 2 fall within a range that allows simulation in laboratory light gas gun experiments designed to validate analytical methods for the interpretation of dust impacts on the aluminum foil components of the Stardust collector. Buckshot of a wide size, shape and density range of mineral, glass, polymer and metal grains, have been fired to impact perpendicularly upon samples of Stardust Al1100 foil, tightly wrapped onto aluminium alloy plate as an analogue of foil on the spacecraft collector. We have not yet been able to produce laboratory impacts by projectiles with weak and porous aggregate structure, as may occur in some cometary dust grains. In this report we present information on crater gross morphology and its dependence on particle size and density, the pre-existing major and trace element composition of the foil, geometrical issues for energy dispersive X-ray analysis of the impact residues in scanning electron microscopes, and the modification of dust chemical composition during creation of impact craters as revealed by analytical transmission electron microscopy. Together, these observations help to underpin the interpretation of size, density and composition for particles impacted upon the Stardust aluminum foils.

The spectral and temporal properties of the non-thermal emission ofthe nearby XRF 060218 in 0.3-150 keV band are studied. We show that both the spectral energy distribution and the light curve properties suggest the same origin of the non-thermal emission detected by {\em Swift} BAT and XRT. This event has the longest pulse duration and spectral lag observed to date among the known GRBs. The pulse structure and its energy dependence are analogous to typical GRBs. By extrapolating the observed spectral lag to the {\em CGRO/BATSE} bands we find that the hypothesis that this event complies with the same luminosity-lag relation with bright GRBs cannot be ruled out at $2\sigma$ significance level. These intriguing facts, along with its compliance with the Amati-relation, indicate that XRF 060218 shares the similar radiation physics as typical GRBs.

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.

It is long speculated that long duration gamma-ray bursts (GRBs) originate froma relativistic jet emerging from a collapsing massive star progenitor. Although associations of core-collapsing supernovae with long GRB afterglows have been identified in a number of systems, including the latest X-ray flash (XRF) 060218/SN 2006aj connection detected by Swift, direct evidence of a relativistic jet emerging from a collapsing star is still lacking. Here we report the detection of a thermal emission component (high-T component) accompanying the prompt X-ray emission of XRF 060218, with temperature kT_H=1.21+0.22/-0.24 keV and effective blackbody radius R_H~ 5\times 10^{9} cm. This high-T component co-exists with another low-T thermal component as reported by Campana et al. 2006 for at least 2700 seconds, but evolves independently with respect to the low-T component by tracing the lightcurve of the non-thermal component. We identify this high-T thermal component as the emission of a hot cocoon surrounding the relativistic jet, as expected from the theoretical models.

Although the link between long Gamma Ray Bursts (GRBs) and supernovae (SNe) has been established, hitherto there have been no observations of the beginning of a supernova explosion and its intimate link to a GRB. In particular, we do not know however how a GRB jet emerges from the star surface nor how a GRB progenitor explodes. Here we report on observations of the close GRB060218 and its connection to SN2006aj. In addition to the classical non-thermal emission, GRB060218 shows a thermal component in its X-ray spectrum, which cools and shifts into the optical/UV band as time passes. We interpret these features as arising from the break out of a shock driven by a mildly relativistic shell into the dense wind surrounding the progenitor. Our observations allow us for the first time to catch a SN in the act of exploding, to directly observe the shock break-out and to provide strong evidence that the GRB progenitor was a Wolf-Rayet star.

The gamma-ray burst (GRB) 060218/SN 2006aj is a peculiar event, with the second lowest redshift, low luminosity, long duration, chromatic lightcurve features, and in particular, the presence of a thermal component in the X-ray and UV-optical spectra. Thanks to detailed temporal and spectral coverage of the {\em Swift} observatory, the abundant data allow the GRB prompt emission to be modelled in great detail for the first time. The low flux of prompt UV/optical emission disfavors the conventional internal shock/synchrotron radiation models, which generally predict strong UV/optical emission. Here we show that the unusual prompt emission of GRB 060218 can be produced by inverse-Compton scattering of shock-accelerated relativistic electrons off the detected thermal photons. A pair of (forward plus reverse) shocks form when a relativistic outflow interacts with a preexisting slower shell. The observed gamma-ray emission and X-ray emission arise from the reverse-shocked and forward-shocked regions, respectively. A fit to the data requires an initially increasing outflow luminosity, which is consistent with the prediction of the popular collapsar model.

Recent observations support the suggestion that short-duration gamma-ray bursts are produced by compact star mergers. The X-ray flares discovered in two short gamma-ray bursts last much longer than the previously proposed postmerger energy release time scales. Here we show that they can be produced by differentially rotating, millisecond pulsars after the mergers of binary neutron stars. The differential rotation leads to windup of interior poloidal magnetic fields and the resulting toroidal fields are strong enough to float up and break through the stellar surface. Magnetic reconnection--driven explosive events then occur, leading to multiple X-ray flares minutes after the original gamma-ray burst.

A new method of measuring cosmology with gamma-ray bursts(GRBs) has been proposed by Liang and Zhang recently. In this method, only observable quantities including the rest frame peak energy of the \nu F_\nu spectrum (E'_p), the isotropic energy of GRB (E_{\gamma,iso}), and the break time of the optical afterglow light curves in the rest frame (t'_b) are used. By considering this method we constrain the cosmological parameters and the redshift at which the universe expanded from the deceleration to acceleration phase. We add five recently-detected GRBs to the sample and derive E_{\gamma, iso}/10^{52} ergs=(0.93\pm 0.25)\times (E'_p/100 keV)^{1.91\pm 0.32}\times (t'_b/day)^{-0.93\pm 0.38} for a flat cosmology with \Omega_M=0.28 and H_0=71.0 km s^{-1} Mpc^{-1}. This relation is independent of the medium density around bursts and the efficiency of conversion of the explosion energy to gamma-ray energy that are difficult to measure. We regard the E_{\gamma,iso}(E'_p, t'_b) relationship as a standard candle and find 0.05<\Omega_M<0.48 and \Omega_\Lambda<1.15 (at the 1\sigma confidence level). In a flat universe with the cosmological constant we obtain 0.25<\Omega_M <0.46 and 0.54<\Omega_\Lambda<0.78 at the 1\sigma confidence level. The transition redshift is z_T=0.69_{-0.12}^{+0.11}. Combining 20 GRBs and 157 type Ia supernovae, we find \Omega_M=0.29 \pm0.03 for a flat universe and the transition redshift is z_T=0.61_{-0.05}^{+0.06}, which is slightly larger than the value found by considering SNe Ia alone. In particular, We also discuss several dark-energy models in which the equation of state w(z) is arameterized and investigate constraints on the cosmological parameters in detail.

We analyze several recently detected gamma-ray bursts (GRBs) with late X-ray flares in the context of late internal shock and late external shock models. We find that the X-ray flares in GRB 050421 and GRB 050502B originate from late internal shocks, while the main X-ray flares in GRB 050406 and GRB 050607 may arise from late external shocks. Under the assumption that the central engine has two periods of activities, we get four basic types of X-ray light curves. The classification of these types depends on which period of activities produces the prompt gamma-ray emission (Type 1 and Type 2: the earlier period; Type 3 and Type 4: the late period), and on whether the late ejecta catching up with the early ejecta happens earlier than the deceleration of the early ejecta (Type 1 and Type 3) or not (Type 2 and Type 4). We find that the X-ray flare caused by a late external shock is a special case of Type 1. Our analysis reveals that the X-ray light curves of GRBs 050406, 050421, and 050607 can be classified as Type 1, while the X-ray light curve of GRB 050502B is classified as Type 2. However, the X-ray light curve of GRB 050406 is also likely to be Type 2. We also predict a long-lag short-lived X-ray flare caused by the inner external shock, which forms when a low baryon-loading long-lag late ejecta decelerates in the non-relativistic tail of an outer external shock driven by an early ejecta.

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 present a general method for testing gamma-ray burst (GRB) jet structure and carry out a comprehensive analysis about the prevalent jet structure models. According to the jet angular energy distribution, we can not only derive the expected distribution of the GRB isotropic-equivalent energy release for any possible jet structure, but also obtain a two-dimensional distribution including redshift z. By using the Kolmogorov-Smirnov test we compare the predicted distribution with the observed sample, and find that the power-law structured jet model is most consistent with the current sample and that the uniform jet model is also plausible. However, this conclusion is tentative because of the small size and the inhomogeneity of this sample. Future observations (e.g., Swift) will provide a larger and less biased sample for us to make a robust conclusion by using the procedure proposed in this paper.

The prompt localization of gamma-ray burst (GRB) 050525A by {\em Swift} allowed the rapid follow-up of the afterglow. The observations revealed that the optical afterglow had a major rebrightening starting at $\sim 0.01$ days and ending at $\sim 0.03$ days, which was followed by an initial power-law decay. Here we show that this early emission feature can be interpreted as the reverse shock emission superposed by the forward shock emission in an interstellar medium environment. By fitting the observed data, we further constrain some parameters of the standard fireball-shock model: the initial Lorentz factor of the ejecta $\gamma_0>120$, the magnetic energy fraction $\epsilon_B>4\times10^{-6}$, and the medium density $n<2 {\rm cm^{-3}}$. These limits are consistent with those from the other very-early optical afterglows observed so far. In principle, a wind environment for GRB 050525A is disfavored.

Three different methods of measuring cosmology with gamma-ray bursts (GRBs) have been proposed since a relation between the $\gamma$-ray energy $E_{\gamma} $ of a GRB jet and the peak energy $E_p$ of the $\nu F_{\nu}$ spectrum in the burst frame was reported by Ghirlanda and coauthors. In Method I, to calculate the probability for a favored cosmology, only the contribution of the $E_\gamma-E_p$ relation that is already best fitted for this cosmology is considered. We apply this method to a sample of 17 GRBs, and obtain the mass density $\Omega_M=0.15^{+0.45}_{-0.13}$ ($1\sigma$) for a flat $\Lambda$CDM universe. In Method II, to calculate the probability for some certain cosmology, contributions of all the possible $E_\gamma-E_p$ relations that are best fitted for their corresponding cosmologies are taken into account. With this method, we find a constraint on the mass density $0.14<\Omega_M<0.69$ ($1 \sigma$) for a flat universe. In Method III, to obtain the probability for some cosmology, contributions of all the possible $E_\gamma-E_p$ relations associated with their unequal weights are considered. With this method, we obtain an inspiring constraint on the mass density $0.16<\Omega_M<0.45$ ($1 \sigma$) for a flat universe, and a $\chi^2_{dof}=19.08/15=1.27$ for the concordance model of $\Omega_M=0.27$. (abridged)

The brightest giant flare from the soft $\gamma$-ray repeater (SGR) 1806-20 was detected on 2004 December 27. The isotropic-equivalent energy release of this burst is at least one order of magnitude more energetic than those of the two other SGR giant flares. Starting from about one week after the burst, a very bright ($\sim 80$ mJy), fading radio afterglow was detected. Follow-up observations revealed the multi-frequency light curves of the afterglow and the temporal evolution of the source size. Here we show that these observations can be understood in a two-component explosion model. In this model, one component is a relativistic collimated outflow responsible for the initial giant flare and the early afterglow, and another component is a subrelativistic wider outflow responsible for the late afterglow. We also discuss triggering mechanisms of these two components within the framework of the magnetar model.

Recent follow-up observations of the December 27 giant flare of SGR 1806-20 have detected a multiple-frequency radio afterglow from 240 MHz to 8.46 GHz, extending in time from a week to about a month after the flare. The angular size of the source was also measured for the first time. Here we show that this radio afterglow gives the first piece of clear evidence that an energetic blast wave sweeps up its surrounding medium and produces a synchrotron afterglow, the same mechanism as established for gamma-ray burst afterglows. The optical afterglow is expected to be intrinsically as bright as $m_R\simeq13$ at $t\la 0.1$ days after the flare, but very heavy extinction makes the detection difficult because of the low galactic latitude of the source. Rapid infrared follow-up observations to giant flares are therefore crucial for the low-latitude SGRs, while for high-latitude SGRs (e.g. SGR 0526-66), rapid follow-ups should result in identification of their possible optical afterglows. Rapid multi-wavelength follow-ups will also provide more detailed information of the early evolution of a fireball as well as its composition.

Soft gamma-ray repeaters are at determined distances and their positions are known accurately. If observed, afterglows from their soft gamma-ray bursts will provide important clues to the study of the so called "classical gamma-ray bursts". On applying the popular fireball/blastwave model of classical gamma-ray bursts to soft gamma-ray repeaters, it is found that their X-ray and optical afterglows are detectable. Monitoring of the three repeaters is solicited.

Polarization behaviors of optical afterglows from two-component gamma-ray burst jets are investigated, assuming various configurations for the two components. In most cases, the observed polarization is dominated by the inner narrow component for a long period. Interestingly, it is revealed that different assumptions about the lateral expansion of the jet can lead to different evolutions of the position angle of polarization. The observed afterglow light curve and polarization behaviors of GRB 020813 can be well explained by the two-component jet model. Particularly, the model is able to explain the constancy of the observed position angle in this event, given that the line of sight is slightly outside the narrow component.

Optical and radio afterglows arising from the shocks by relativistic conical ejecta running into pre-burst massive stellar winds are revisited. Under the homogeneous thin-shell approximation and a realistic treatment for the lateral expansion of jets, our results show that a notable break exists in the optical light curve in most cases we calculated in which the physical parameters are varied within reasonable ranges. For a relatively tenuous wind which cannot decelerate the relativistic jet to cause a light curve break within days, the wind termination shock due to the ram pressure of the surrounding medium occurs at a small radius, namely, a few times 10^{17} cm. In such a structured wind environment, the jet will pass through the wind within several hours and run into the outer uniform dense medium. The resulting optical light curve flattens with a shallower drop after the jet encounters the uniform medium, and then declines deeply, triggered by runaway lateral expansion.

In this paper we study the possibility of probing for the absolute neutrino mass and its variation with short Gamma Ray Burst (GRB). We have calculated the flight time difference between a massive neutrino and a photon in two different approaches to mass varying neutrinos. Firstly we parametrize the neutrino mass as a function of redshift in a model independent way, then we consider two specific models where the neutrino mass varies during the evolution of the Quintessence fields. Our calculations show in general the value of the time delay is changed substantially relative to a constant neutrino mass. Furthermore our numerical results show that the flight time delay in these models is expected to be larger than the duration time of the short GRB, which opens a possibility of testing the scenario of mass varying neutrino with the short GRB.

(abridged) Afterglow light curves are constructed analytically for realistic gamma-ray burst remnants decelerating in either a homogeneous interstellar medium or a stellar wind environment, taking into account the radiative loss of the blast wave, which affects the temporal behaviors significantly. Inverse Compton scattering is considered. The inverse Compton effect prolongs the fast-cooling phase markedly, during which the relativistic shock is semi-radiative and the radiation efficiency is approximately constant. It is further shown that the shock is still semi-radiative for quite a long time after it transits into the slow-cooling phase. The temporal decaying index of the X-ray afterglow light curve in this semi-radiative phase is more consistent with the observed $<\alpha_{X}> \sim 1.3$ than the commonly used adiabatic one. To manifest as a bump or even dominant in the X-ray afterglows during the relativistic stage, it is required that the density should be larger than about 1-10 cm$^{-3}$ in the interstellar medium case, or the wind parameter $A_{\ast}$ should be larger than unity in the stellar wind case.

We derive a relation between the flux density $F_{\nu,j}$ at the light-curve break of a gamma-ray burst (GRB) afterglow and the break time $t_{j}$. The break is due to the transition from the spherical-like to jet-like evolution of the afterglow, when the Lorentz factor of the jet equals the inverse of the initial half-opening angle, i.e., $\gamma=1/\theta_0$. We show that this relation indeed behaves as $F_{\nu,j}\propto t_{j}^{-p}$ among GRBs for the slow-cooling case, where $p$ is the power-law index of electron distribution. A statistical analysis of the optical jet breaks of nine GRBs gives $p=2.10\pm 0.21$, which is consistent with the shock acceleration theory. The value of $p$ derived in this way is different from the observed temporal index $\alpha_2$ ($F_{\nu}\propto t^{-\alpha_{2}}$) of the late-time light curve after $t_{j}$, which suffers several uncertainties from the unclear hydrodynamics of the sideways expansion and exhibits a large dispersion. Our results not only confirm that the remnants of GRBs are standard candles, but also provide the first evidence that the physical parameters of relativistic shocks are universal, with the favored values $\epsilon_{e}\sim 0.1$ and $\epsilon_{B}\sim 10^{-3}$.

An $E_{\gamma,{\rm jet}}\propto {E'_p}^{1.5}$ relationship with a small scatter for current $\gamma$-ray burst (GRB) data was recently reported, where $E_{\gamma,{\rm jet}}$ is the beaming-corrected $\gamma$-ray energy and $E'_p$ is the $\nu F_\nu$ peak energy in the local observer frame. By considering this relationship for a sample of 12 GRBs with known redshift, peak energy, and break time of afterglow light curves, we constrain the mass density of the universe and the nature of dark energy. We find that the mass density $\Omega_M=0.35\pm^{0.15}_{0.15}$ (at the $1\sigma$ confident level) for a flat universe with a cosmological constant, and the $w$ parameter of an assumed static dark-energy equation of state $w=-0.84\pm^{0.57}_{0.83}$ ($1\sigma$). Our results are consistent with those from type Ia supernovae. A larger sample established by the upcoming {\em Swift} satellite is expected to provide further constraints.

Based on the standard energy reservoir of gamma-ray burst (GRB) jets and the relationship between peak energy ($E_p$) of $\nu F_{\nu}$ spectrum and the equivalent-isotropic energy ($E_{\rm{iso}}$), we test the GRB probability distributions as a function of viewing angle ($\theta$) and of viewing angle together with redshift $z$, $P(\theta)$ and $P(\theta,z)$, predicted by the universal-structured jet (USJ) model with simulations on observational bases for a detection threshold of $S>4 \times 10^{-5}$ erg cm$^{-2}$, where $S$ is the total fluence in an energy band $>20$ keV. We simulate a sample including $10^{6}$ GRBs, and then derive $P(\theta)$ and $P(\theta,z)$. We show that both simulated and theoretical $P(\theta)$ and $P(\theta,z)$ are consistent with each other. This consistency seems to be regarded as support for the USJ model. A small difference between the theoretical and simulated results might be due to the observed $E_p$ distribution used in our simulations, which perhaps suffers a little of biases for $E_p<100$ keV.

We study the peak energy ($E_{\rm{p}}$) distribution of the $\nu F_{\nu}$ spectra of gamma-ray bursts (GRBs) and X-ray flashes (XRFs) with a sample of 57 bursts observed by {\em High Energy Transient Explorer 2} ({\em HETE-2}) French Gamma Telescope and discuss its implications for the jet structure models. Combining the observed $E_{\rm{p}}$ distribution of {\em HETE--2} GRBs/XRFs with that of BATSE GRBs, we find that the observed $E_{\rm{p}}$ distribution of GRBs/XRFs is a bimodal one with peaks of $\lesssim 30$ keV and $\sim 160-250$ keV. According to the recently-discovered equivalent-isotropic energy-$E_{\rm{p}}$ relationship, such a bimodal distribution implies a two-component structure of GRB/XRF jets. A simple simulation analysis shows that this structured jet model does roughly reproduce a bimodal distribution with peaks of $\sim 15$ and $\sim 200$ keV. We argue that future observations of the peak of $\sim 15$ keV in the $E_{\rm{p}}$ distribution would be evidence supporting this model. {\em Swift}, which covers an energy band of 0.2--150 keV, is expected to provide a key test for our results.

Using a sample of 2408 time-resolved spectra for 91 BATSE gamma-ray bursts (GRBs) presented by Preece et al., we show that the relation between the isotropic-equivalent luminosity (L_iso) and the spectral peak energy (E_p) in the cosmological rest frame, L_iso \propto E_p^2, not only holds within these bursts, but also holds among these GRBs, assuming that the burst rate as a function of redshift is proportional to the star formation rate. The possible implications of this relation for the emission models of GRBs are discussed. We suggest that both the kinetic-energy-dominated internal shock model and the magnetic-dissipation-dominated external shock model can well interpret this relation. We constrain the parameters for these two models, and find that they are in a good agreement with the parameters from the fittings to the afterglow data (abridged).

Highly magnetized, rapidly rotating compact objects are widely argued as central energy sources of $\gamma$-ray bursts (GRBs). After the GRB, such a magnetar-like object may directly lose its rotational energy through some magnetically-driven processes, which produce an ultrarelativistic wind dominated possibly by the energy flux of electron-positron pairs. The interaction of such a wind with an outward-expanding fireball leads to a relativistic wind bubble, being regarded as a relativistic version of the well-studied Crab Nebula. We here explore the dynamics of this wind bubble and its emission signatures. We find that when the injection energy significantly exceeds the initial energy of the fireball, the bulk Lorentz factor of the wind bubble decays more slowly than before, and more importantly, the reverse-shock emission could dominate the afterglow emission, which yields a bump in afterglow light curves. In addition, high polarization of the bump emission would be expected if a toroidal magnetic field in the shocked wind dominates over the random component.

Optical afterglows from two-component jets under various configurations are investigated numerically. Generally, the light curve is characterized by a rapid rebrightening when the observer is off-axis with respect to the narrow component, with the amplitude and peak time depending on detailed parameters. We further show that the optical afterglow of XRF 030723, especially its notable and rapid rebrightening, can be well explained by a typical two-component jet. This X-ray flash, together with GRB 030329, strongly hints the two-component jet model as a unified picture for X-ray flashes and gamma-ray bursts. With a narrow but ultra-relativistic inner outflow and a wide but less energetic outer ejecta, a two-component jet will be observed as a typical gamma-ray burst if our line of sight is within the angular scope of the narrow outflow. Otherwise, if the line of sight is within or slightly beyond the cone of the wide component, an X-ray flash will be detected.

It has been suggested that electromagnetic cascade of very high energy gamma-rays from gamma-ray bursts (GRBs) in the infrared/microwave background can produce delayed MeV-GeV photons. This delay could be caused by the angular spreading effect of the scattered microwave photons or deflection of the secondly pairs due to intergalactic magnetic field. Very high energy TeV photons of GRBs could be produced by a few mechanisms including the proton-synchrotron radiation and electron inverse Compton emission from GRB internal shocks as well as external shocks. We suggest that the information provided by the delayed emission could give constraints on models for TeV gamma-rays. A more accurate calculation of the delayed time caused by the angular spreading effect is presented by considering recent observations of the extragalactic infrared background and the theoretic high-redshift infrared background. We also suggest that the dependence of the maximum time delay of scattered photons on their energies, if determined by future GLAST detector, could differentiate the two mechanisms causing the time delay.

In many GRBs, afterglows have been observed for months or even years. It deserves noting that at such late stages, the remnants should have entered the deep Newtonian phase, during which the majority of shock-accelerated electrons will no longer be highly relativistic. However, a small portion of electrons are still ultra-relativistic and capable of emitting synchrotron radiation. Under the assumption that the electrons obey a power-law distribution according to their kinetic energy (not simply the Lorentz factor), we calculate optical afterglows from both isotropic fireballs and beamed ejecta, paying special attention to the late stages. In the beamed cases, it is found that the light curves are universally characterized by a flattening during the deep Newtonian phase. Implication of our results on orphan afterglows is also addressed.

For High frequency BL Lac objects (HBLs) like H1426+428, a significant fraction of their TeV gamma-rays emitted are likely to be absorbed in interactions with the diffuse IR background, yielding $e^\pm$ pairs. The resulting $e^\pm$ pairs generate one hitherto undiscovered GeV emission by inverse Compton scattering with the cosmic microwave background photons (CMBPs). We study such emission by taking the 1998-2000 CAT data, the reanalyzed 1999 & 2000 HEGRA data and the corresponding intrinsic spectra proposed by Aharonian et al. (2003a). We numerically calculate the scattered photon spectra for different intergalactic magnetic field (IGMF) strengths. If the IGMF is about $10^{-18}{\rm G}$ or weaker, there comes very strong GeV emission, whose flux is far above the detection sensitivity of the upcoming satellite GLAST! Considered its relatively high redshift ($z=0.129$), the detected GeV emission in turn provides us a valuable chance to calibrate the poor known spectral energy distribution of the intergalactic infrared background, or provides us some reliable constraints on the poorly known IGMF strength.

The shock model has successfully explained the observed behaviors of afterglows from long gamma-ray bursts (GRBs). Here we use it to investigate the so-called early afterglows from short GRBs, which arises from blast waves that are not decelerated considerably by their surrounding medium. We consider a nearby medium loaded with $e^{\pm}$ pairs (Beloborodov 2002). The temporal behaviors show first a soft-to-hard spectral evolution, from the optical to hard X-ray, and then a usual hard-to-soft evolution after the blast waves begin to decelerate. The light curves show variability, and consist of two peaks. The first peak, due to the pair effect, can be observed in the X-ray, though too faint and too short in the optical. The second peak will be easily detected by {\it Swift}. We show that detections of the double-peak structure in the light curves of early afterglows are very helpful to determine all the shock parameters of short GRBs, including both the parameters of the relativistic source and the surroundings. Besides, from the requirement that the forward-shock emission in short GRBs should be below the BATSE detection threshold, we give a strong constraint on the shock model parameters. In particular, the initial Lorentz factor of the source is limited to be no more than $\sim 10^3$, and the ambient medium density is inferred to be low, $n\la 10^{-1}$ cm$^{-3}$.

Gamma-ray bursts (GRBs) are believed to originate from ultra-relativistic winds/fireballs to avoid the "compactness problem". However, the most energetic photons in GRBs may still suffer from $\gamma-\gamma$ absorption leading to electron/positron pair production in the winds/fireballs. We show here that in a wide range of model parameters, the resulting pairs may dominate those electrons associated with baryons. Later on, the pairs would be carried into a reverse shock so that a shocked pair-rich fireball may produce a strong flash at lower frequencies, i.e. in the IR band, in contrast with optical/UV emission from a pair-poor fireball. The IR emission would show a 5/2 spectral index due to strong self-absorption. Rapid responses to GRB triggers in the IR band would detect such strong flashes. The future detections of many IR flashes will infer that the rarity of prompt optical/UV emissions is in fact due to dust obscuration in the star formation regions.

Magnetars can accelerate cosmic rays to high energies through the unipolar effect, and are also copious soft photon emitters. We show that young, fast-rotating magnetars whose spin and magnetic moment point in opposite directions emit high energy neutrinos from their polar caps through photomeson interactions. We identify a neutrino cut-off band in the magnetar period-magnetic field strength phase diagram, corresponding to the photomeson interaction threshold. Within uncertainties, we point out four possible neutrino emission candidates among the currently known magnetars, the brightest of which may be detectable for a chance on-beam alignment. Young magnetars in the universe would also contribute to a weak diffuse neutrino background, whose detectability is marginal, depending on the typical neutrino energy.

The optical flash accompanying GRB 990123 is believed to be powered by the reverse shock of a thin shell. With the best fitted physical parameters for GRB 990123 (Panaitescu & Kumar 2001) and the assumption that the parameters in the optical flash are the same as those in the afterglow, we show that: 1) the shell is thick but not thin, and we have provided the light curve for the thick shell case which coincides with the observation; 2) the theoretical peak flux of the optical flash accounts for only 3\times10^{-4} of the observed. In order to compensate this divergency, the physical parameters the electron energy ration and the magnetic ratio $\epsilon_e$, $\epsilon_B$ should be 0.61, 0.39 respectively, which are much different from those in the late afterglow.

The idea that gamma-ray bursts might be a kind of phenomena associated with neutron star kicks was first proposed by Dar & Plaga (1999). Here we study this mechanism in more detail and point out that the neutron star should be a high speed one (with proper motion larger than $\sim 1000$ km/s). It is shown that the model agrees well with observations in many aspects, such as the energetics, the event rate, the collimation, the bimodal distribution of durations, the narrowly clustered intrinsic energy, and the association of gamma-ray bursts with supernovae and star forming regions. We also discuss the implications of this model on the neutron star kick mechanism, and suggest that the high kick speed were probably acquired due to the electromagnetic rocket effect of a millisecond magnetar with an off-centered magnetic dipole.

We present an explanation for the unusual temporal feature of the GRB 030226 afterglow. The R-band afterglow of this burst faded as ~ t^{-1.2} in ~ 0.2 days after the burst, rebrightened during the period of ~ 0.2 - 0.5 days, and then declined with ~ t^{-2.0}. To fit such a light curve, we consider an ultrarelativistic jetted blast wave expanding in a density-jump medium. The interaction of the blast wave with a large density jump produces relativistic reverse and forward shocks. In this model, the observed rebrightening is due to emissions from these newly forming shocks, and the late-time afterglow is caused by sideways expansion of the jet. Our fitting implies that the progenitor star of GRB 030226 could have produced a stellar wind with a large density jump prior to the GRB onset.

The interaction of a relativistic fireball with its ambient medium is described through two shocks: a reverse shock that propagates into the fireball, and a forward shock that propagates into the medium. The observed optical flash of GRB 990123 has been considered to be the emission from such a reverse shock. The observational properties of afterglows suggest that the progenitors of some GRBs may be massive stars and their surrounding media may be stellar winds. We here study very early afterglows from the reverse and forward shocks in winds. An optical flash mainly arises from the relativistic reverse shock while a radio flare is produced by the forward shock. The peak flux densities of optical flashes are larger than 1 Jy for typical parameters, if we do not take into account some appropriate dust obscuration along the line of sight. The radio flare always has a long lasting constant flux, which will not be covered up by interstellar scintillation. The non-detections of optical flashes brighter than about 9th magnitude may constrain the GRBs isotropic energies to be no more than a few $10^{52}$ ergs and wind intensities to be relatively weak.

Large-scale, decelerating, relativistic X-ray jets due to material ejected from the black-hole candidate X-ray transient and microquasar XTE J1550-564 has been recently discovered with Chandra by Corbel et al. (2002). We find that the dynamical evolution of the eastern jet at the late time is consistent with the well-known Sedov evolutionary phase. A trans-relativistic external shock dynamic model by analogy with the evolution of gamma-ray burst remnants, is shown to be able to fit the observation data reasonably well. The inferred interstellar medium density around the source is well below the canonical value $n_ISM \sim 1 cm^{-3}$. We find that the emission from the continuously shocked interstellar medium (forward shock region) decays too slowly to be a viable mechanism for the eastern X-ray jet. However, the rapidly fading X-ray emission can be interpreted as synchrotron radiation from the non-thermal electrons in the adiabatically expanding ejecta. These electrons were accelerated by the reverse shock (moving back into the ejecta) which becomes important when the inertia of the swept external matter leads to an appreciable slowing down of the original ejecta. To ensure the dominance of the emission from the shocked ejecta over that from the forward shock region during the period of the observations, the magnetic field and electron energy fractions in the forward shock region must be far below equipartition. Future continuous, follow-up multi-wavelength observations of new ejection events from microquasars up to the significant deceleration phase should provide more valuable insight into the nature of the interaction between the jets and external medium.

The absorption feature detected in the prompt X-ray emission of GRB990705 has important consequences for its circum-burst environment and therefore on its afterglows. Here we investigate whether the circum-burst environment constrained by the absorption feature could be consistent with the observed H-band afterglow, which exhibits an earlier power law decay but a much faster decay about one day after the burst. Two possible geometries of the afterglow-emitting regions are suggested: 1) afterglow emission produced by the impact of the fireball on the surrounding torus, which serves as the absorbing material of the X-ray feature, as would be expected in the models involving that a supernova explosion precedes the gamma-ray burst by some time; 2)afterglow emission produced in the dense circum-burst medium inside the torus. In case 1), the faster decay at the later time is attributed to the disappearance of the shock due to the counter-pressure in the hot torus illuminated by the burst and afterglow photons. For case 2), the circum-burst medium density is found to be very high ($n\ga 10^4-10^5 {\rm cm^{-3}}$) if the emitting plasma is a jet or even higher if it is spherical. Future better observations of afterglows of GRBs that have absorption features might make it possible to make a more definite choice between these two scenarios.

Synchrotron self-Compton (SSC) process in the reverse shocks of gamma-ray bursts is suggested to be responsible for the observed prompt high-energy gamma-ray emissions from several gamma-ray bursts. We find that the SSC emission from the reverse shocks dominates over other emission processes in energy bands from tens of MeV to tens of GeV, for a wide range of shock parameters. This model is favorable for escape of energetic photons from the emitting regions due to a lower internal pair-production optical depth, as the characteristic size of the reverse shock region is much larger than that of internal shocks. We predict that, in this model, the prompt high-energy emissions are correlated with the prompt optical flashes, which can be test in the forthcoming GLAST era.

In previous works, a generic dynamical model has been suggested by Huang et al., which is shown to be correct for both adiabatic and radiative blastwaves, and in both ultra-relativistic and non-relativistic phases. In deriving their equations, Huang et al. have assumed that the radiative efficiency of the fireball is constant. They then applied their model directly to realistic cases where the radiative efficiency evolves with time. In this paper, we abandon the above assumption and re-derive a more accurate dynamical equation for gamma-ray burst remnants. Numerical results show that Huang et al.'s model is accurate enough in general cases.

Several high-frequency peaked BL Lac objects such as Mrk 501 are strong TeV emitters. However, a significant fraction of the TeV gamma rays emitted are likely to be absorbed in interactions with the diffuse IR background, yielding electron-positron pairs. Hence, the observed TeV spectrum must be steeper than the intrinsic one. Using the recently derived intrinsic $\gamma$-ray spectrum of Mrk 501 during its 1997 high state, we study the inverse-Compton scattering of cosmic microwave photons by the resulting electron-positron pairs, which implies the existence of a hitherto undiscovered GeV emission. The typical duration of the GeV emission is determined by the flaring activity time and the energy-dependent magnetic deflection time. We numerically calculate the scattered photon spectrum for different intergalactic magnetic field (IGMF) strengths, and find a spectral turnover and flare duration at GeV energies which are dependent on the field strength. We also estimate the scattered photon flux in the quiescent state of Mrk 501. The GeV flux levels predicted are consistent with existing EGRET upper limits, and should be detectable above the synchrotron -- self Compton (SSC) component with the {\em Gamma-Ray Large Area Space Telescope} ({\em GLAST}) for IGMFs $\lesssim 10^{-16}$ G, as expected in voids. Such detections would provide constraints on the strength of weak IGMFs.

Gamma-ray bursts may occur within pulsar wind bubbles (PWBs) under a number of scenarios, such as the supranova-like models in which the progenitor pulsar drives a powerful wind shocking against the ambient medium before it comes to death and produces a fireball. We here study the early afterglow emission from GRBs expanding into such a PWB environment. Different from the usual cold GRB external medium, the PWBs consist of a hot electron-positron medium with typical 'thermal' Lorentz factor of the order of gamma_w, the Lorentz factor of the pulsar particle wind. After GRB blast waves shock these hot electron-positron pairs, they will emit synchrotron radiation peaking at GeV bands. It is shown that GeV photons suffer negligible absorption by the soft photons radiation field in PWBs. Thus, strong GeV emissions in the early afterglow phases are expected, providing a plausible explanation for the long-duration GeV emission from GRB940217 detected by EGRET. Future GLAST may have the potential to test this GRB-PWB interaction model.

It is widely believed that multiwavelength afterglows of gamma-ray bursts (GRBs) originate from relativistic blast waves. We here show that in such blast waves, a significant fraction of the energy of shock-accelerated protons would be lost due to pion production by interactions with afterglow photons. This could lead to long-term production of $10^{16}$--$10^{18}$ eV neutrinos and sub-TeV $\gamma$-rays that accompany with usual afterglows, provided that the protons are accelerated to $10^{19}$ eV in the blast waves.

Based on a refined generic dynamical model, we investigate afterglows from jetted gamma-ray burst (GRB) remnants numerically. In the relativistic phase, the light curve break could marginally be seen. However, an obvious break does exist at the transition from the relativistic phase to the non-relativistic phase, which typically occurs at time 10 to 30 days. It is very interesting that the break is affected by many parameters, especially by the electron energy fraction (xi_e), and the magnetic energy fraction (xi_B^2). Implication of orphan afterglow surveys on GRB beaming is investigated. The possible existence of a kind of cylindrical jets is also discussed.

We generally analyze prompt high-energy emission above a few hundreds of GeV due to synchrotron self-Compton scattering in internal shocks. However, such photons cannot be detected because they may collide with cosmic infrared background photons, leading to electron/positron pair production. Inverse-Compton scattering of the resulting electron/positron pairs off cosmic microwave background photons will produce delayed MeV-GeV emission, which may be much stronger than a typical high-energy afterglow in the external shock model. We expand on the Cheng & Cheng model by deriving the emission spectrum and duration in the standard fireball shock model. A typical duration of the emission is ~ 10^3 seconds, and the time-integrated scattered photon spectrum is nu^{-(p+6)/4}, where p is the index of the electron energy distribution behind internal shocks. This is slightly harder than the synchrotron photon spectrum, nu^{-(p+2)/2}. The lower energy property of the scattered photon spectrum is dependent on the spectral energy distribution of the cosmic infrared background radiation. Therefore, future observations on such delayed MeV-GeV emission and the higher-energy spectral cutoff by the Gamma-Ray Large Area Space Telescope (GLAST) would provide a probe of the cosmic infrared background radiation.

It is believed that orphan afterglow searches can help to measure the beaming angle in gamma-ray bursts (GRBs). Great expectations have been put on this method. We point out that the method is in fact not as simple as we originally expected. Due to the baryon-rich environment that is common to almost all popular progenitor models, there should be many failed gamma-ray bursts, i.e., fireballs with Lorentz factor much less than 100 -- 1000, but still much larger than unity. In fact, the number of failed gamma-ray bursts may even be much larger than that of successful bursts. Owing to the existence of these failed gamma-ray bursts, there should be many orphan afterglows even if GRBs are due to isotropic fireballs, then the simple discovery of orphan afterglows never means that GRBs be collimated. Unfortunately, to distinguish a failed-GRB orphan and a jetted but off-axis GRB orphan is not an easy task. The major problem is that the trigger time is unknown. Some possible solutions to the problem are suggested.

We analyze in detail the hydrodynamics and afterglow emission of an ultrarelativistic blast wave when it expands in a density-jump medium. Such a medium is likely to appear in the vicinity of gamma-ray bursts (GRBs) associated with massive stars. The interaction of the blast wave with this medium is described through a reverse shock and a forward shock. We show that the reverse shock is initially relativistic if the factor of a density jump ($\alpha$) is much larger than 21, and Newtonian if $1<\alpha\ll 21$. We also calculate light curves of the afterglow emission during the interaction if the reverse shock is relativistic, and find that the optical flux density initially decays abruptly, then rises rapidly, and finally fades based on a power-law, which could be followed by an abrupt decay when the reverse shock has just crossed the originally swept-up matter. Therefore, one property of an afterglow occurring in a large-density-jump medium is an abrupt drop followed by a bump in the light curve and thus provides a probe of circumburst environments. In addition, this property could not only account for the optical afterglows of GRB 970508 and GRB 000301C but also explain the X-ray afterglow of GRB 981226.

We discuss the GRB afterglow emission from a relativistic jet that is initially in the radiative regime in which the accelerated electrons are fast cooling. We note that such a ``semiradiative'' jet decelerates faster than an adiabatic jet does. We also take into account the effect of strong inverse-Compton scattering on the cooling frequency in the synchrotron component and therefore on the light curve decay index. We find that there are two kinds of light-curve break for the jet effect. The first is an ``adiabatic break'' if the electrons become slow cooling before the jet enters a spreading phase, and the second is a ``radiative break'' which appears on the contrary case. We then show how a relativistic jet evolves dynamically and derive the overall temporal synchrotron emission in both cases, focusing on the change in light curve decay index around the break time. Finally, in view of our results, we rule out two cases for relativistic jets to account for the observed light curve breaks in a few afterglows: (i) an adiabatic jet with strong Compton cooling (Y>1) and with the cooling frequency nu_c locating in the observed energy range; (ii) a radiative jet with a significant fraction of total energy occupied by electrons (epsilon_e 1).

We derive light curves of the afterglow emission from highly collimated jets if the power-law index ($p$) of the electron energy distribution is above 1 but below 2. We find (1) below the characteristic synchrotron frequency, the light curve index depends generally on $p$. (2) As long as the jet expansion is spherical, the light curve index above the characteristic frequency increases slowly as the spectral index of the emission increases. (3) Once the jet enters the spreading phase, the high-frequency emission flux decays as $\propto t^{-(p+6)/4}$ rather than $\propto t^{-p}$. All these results differ from those in the case of $p>2$. We compare our analytical results with the observations on the GRB 010222 afterglow, and conclude that the jet model may be unable to explain the observed data. Thus, a more promising explanation for this afterglow seems to be the expansion of a relativistic fireball or a mildly collimated jet in a dense medium.

Assuming an adiabatic evolution of a gamma-ray burst (GRB) fireball interacting with an external medium, we calculate the hydrodynamics of the fireball with energy injection from a strongly magnetic millisecond pulsar through magnetic dipole radiation, and obtain the light curve of the optical afterglow from the fireball by synchrotron radiation. Results are given both for a homogeneous external medium and for a wind ejected by GRB progenitor. Our calculations are also available in both ultra-relativistic and non-relativistic phases. Furthermore, the observed R-band light curve of GRB{000301C} can be well fitted in our model, which might provide a probe of the properties of GRB progenitors.

We present a new model for gamma-ray bursts (GRBs) that are not only associated with supernovae but also have small baryon contamination. In this model, we assume a newborn neutron star to move outward at a kick velocity of $\sim 10^3 {\rm km} {\rm s}^{-1}$ in the supernova ejecta. We find that such a neutron star still hypercritically accretes its surrounding supernova matter. Once the stellar mass increases to some critical mass, the neutron star will undergo a phase transition to become a strange star, leading to an energy release of a few $10^{52}$ ergs. The phase transition, if possibly occuring just near the supernova front, will first result in an ultra-relativistic fireball and then a GRB. This provides a plausible explanation for the GRB-supernova association. We estimate the burst rate to be $\sim 10^{-6}$ per year per galaxy. Our model also predicts other possiblities. For example, if the resulting fireballs have a Lorentz factor of the order of a few, they will produce X-ray GRBs observed by BeppoSAX. We find the rate of such bursts to be $\sim 10^{-5}$ per year per galaxy.

We calculate the spectra of inverse Compton (IC) emissions in gamma-ray burst (GRB) shocks produced when relativistic ejecta encounters the external interstellar medium, assuming a broken power-law approximation to the synchrotron seed spectrum. Four IC processes, including the synchrotron self-Compton (SSC) processes in GRB forward and reverse shocks, and two combined-IC processes (i.e. scattering of reverse shock photons on the electrons in forward shocks and forward shock photons on the electrons in reverse shocks), are considered. We find that the SSC emission from reverse shocks dominates over other emission processes in energy bands from tens of MeV to tens of GeV, for a wide range of shock parameters. This mechanism may be responsible for the prompt high energy gamma-rays detected by the Energetic Gamma Ray Experiment Telescope (EGRET). At TeV energy bands, however, the combined-IC emissions and/or the SSC emission from the forward shocks become increasingly dominant for a moderately steep distribution of shocked electrons.

The afterglow emission from a spreading jet expanding in a circumstellar cloud is discussed. Prompt X-ray radiation and a strong UV flash from the reverse shock produced by the interaction of the jet with the cloud may destroy and clear the dust out to about 30 pc within the initial solid angle of the jet. As the sideways expansion of the jet becomes significant, most of the optical radiation from the high-latitude part of the jet may be absorbed by the dust outside the initial solid angle of the jet, but only the radiation from the part within the initial solid angle can be observed. We analytically show that the flux of the observational radiation decays as $\propto t^{-(p+1)}$ (where $p$ is the power-law index of the electron distribution) in the relativistic phase. This preliminary result motivates us to perform numerical calculations. Our results show that one break in the optical afterglow ligh curve extends over a factor of $\sim 3$ in time rather than one decade in time in the previous jet model. These results may provide a way to judge whether GRBs locate in dense clouds or not. Finally, we carry out a detailed modelling for the R-band afterglow of GRB 000926.

Currently popular models for progenitors of gamma-ray bursts (GRBs) are the mergers of compact objects and the explosions of massive stars. These two cases have distinctive environments for GRBs: compact object mergers occur in the interstellar medium (ISM) and the explosions of massive stars occur in the preburst stellar wind. We here discuss neutrino afterglows from reverse shocks as a result of the interaction of relativistic fireballs with their surrounding wind matter. After comparing with the analytical result of Waxman & Bahcall (2000) for the homogeneous ISM case, we find that the differential spectrum of neutrinos with energy from $\sim 3\times 10^{15}$ to $\sim 3\times 10^{17}$ eV in the wind case is softer by one power of the energy than in the ISM case. Furthermore, the expected flux of upward moving muons produced by neutrino interactions below a detector on the surface of the Earth in the wind case is $\sim 5$ events per year per km$^2$, which is about one order of magnitude larger than in the ISM case. In addition, these properties are independent of whether the fireballs are isotropic or beamed. Therefore, neutrino afterglows, if detected, may provide a way of distinguishing between GRB progenitor models based on the differential spectra of neutrinos and their event rates in a detector.

Some progenitor models of gamma-ray bursts (GRBs) (e.g., collapsars) may produce anisotropic jets in which the energy per unit solid angle is a power-law function of the angle ($\propto\theta^{-k}$). We calculate light curves and spectra for GRB afterglows when such jets expand either in the interstellar medium or in the wind medium. In particular, we take into account two kinds of wind: one ($n\propto r^{-3/2}$) possibly from a typical red supergiant star and another ($n\propto r^{-2}$) possibly from a Wolf-Rayet star. We find that in each type of medium, one break appears in the late-time afterglow light curve for small $k$ but becomes weaker and smoother as $k$ increases. When $k\ge 2$, the break seems to disappear but the afterglow decays rapidly. Thus, one expects that the emission from expanding, highly anisotropic jets provides a plausible explanation for some rapidly fading afteglows whose light curves have no break. We also present good fits to the optical afterglow light curve of GRB 991208. Finally, we argue that this burst might arise from a highly anisotropic jet expanding in the wind ($n\propto r^{-3/2}$) from a red supergiant to interpret the observed radio-to-optical-band afterglow data (spectrum and light curve).

In the standard scenario of the fireball model of gamma-ray bursts(GRBs), the huge initial energy release produces a relativistic blast wave expanding into the external medium and a reverse shock moving into and heating the fireball ejecta. We calculate the high energy gamma-ray emission due to inverse Compton scattering of the synchrotron photons from relativistic electrons in the reverse shock. Under the favorable values of the physical parameters of the GRBs and the interstellar medium, our result shows that during the prompt phase, this emission dominates over the component from the forward shock at high energy gamma-ray bands. This mechanism can excellently account for the observations of the prompt high energy gamma-rays detected by EGRET, such as from GRB930131.

Gamma-ray bursts with long durations are widely thought to arise from the collapse of massive stars, where the wind environment is unavoidable. It is also believed that $\gamma$-ray bursts come from jets. Considering these two points in this paper, we calculate the evolution of a highly collimated jet that expands in a stellar wind environment and the expected afterglow from such a jet. We use a set of refined dynamical equations and a realistic lateral speed of the jet, and find: (1) There is no observable break at the time when the Lorentz factor of the jet is equal to the inverse of its initial half-opening angle. (2) No obvious break appears at the time when the blast wave transits from the relativistic to the non-relativistic phase. (3) For the wind case, there is no flattening tendency even up to $10^9$ s. (4) Compared with the homogeneous medium case, our calculated flux is weaker in the stellar wind case. Finally, we find that two kinds of GRB models (neutron star mergers and massive star collapses) may be discriminated in our numerical results.

There is growing evidence that some long gamma-ray bursts (GRBs) arise from the core collapse of massive stars, and thus it is inevitable that the environments of these GRBs are preburst stellar winds or dense media. We studied, for the first time, the wind model for afterglows based on the Blandford-McKee self-similar solution of a relativistic shock, and suggested that GRB 970616 is an interactor with a stellar wind. We also proposed a dense medium model for some afterglows, e.g., the steepening in the light curve of the R-band afterglow of GRB 990123 may be caused by the adiabatic shock which has evolved from an ultrarelativistic phase to a nonrelativistic phase in a dense medium. We further discussed the dense medium model in more details, and investigated the effects of synchrotron self absorption and energy injection. A shock in a dense medium becomes nonrelativistic rapidly after a short relativistic phase. The afterglow from the shock at the nonrelativistic stage decays more rapidly than at the relativistic stage. Since some models for GRB energy sources predicted that a strongly magnetic millisecond pulsar may be born during GRB formation, we discussed the effect of such a pulsar on the evolution of the nonrelativistic shock through magnetic dipole radiation. We found that in the pulsar energy injection case, the dense medium model fits very well all the observational data of GRB 980519. Recently, we combined the dense medium model with the pulsar energy injection effect to provide a good fit to the optical afterglow data of GRB 000301C.

We have constrained the intrinsic parameters, such as the magnetic energy density fraction ($\epsilon_{B}$), the electron energy density fraction ($\epsilon_e$), the initial Lorentz factor ($\Gamma_0$) and the Lorentz factor of the reverse external shock ($\Gamma_{rs}$), of GRB990123, in terms of the afterglow information (forward shock model) and the optical flash information (reverse shock model). Our result shows: 1) the inferred values of $\epsilon_e$ and $\epsilon_B$ are consistent with the suggestion that they may be universal parameters, comparing to those inferred for GRB970508; 2) the reverse external shock may have become relativistic before it passed through the ejecta shell. Other instrinsic parameters of GRB990123, such as energy contained in the forward shock $E$ and the ambient density $n$ are also determined and discussed in this paper.

Conversion from neutron stars to strange stars as a possible mechanism of cosmological gamma-ray bursts (GRBs) has been discussed in previous works, although the existence of strange stars is still an open question. On the basis of this mechanism, we here outline an explanation of the connection between supernovae (SNe) and GRBs, which has got increasing evidence recently. An asymmetric but normal SN explosion leaves a massive ($\geq1.8{\rm M_\odot}$) and rapidly rotating neutron star, which then converts to a strange star few days later, due to its rapid spindown. The accompanied fireball, which can be accelerated to ultra-relativistic velocity ($\Gamma_0\sim 100$) due to the very low baryon contamination of the strange star, flows out along the direction of the high-velocity SN jet and subsequently produces a GRB and the following low energy afterglows by interacting with the surrounding stellar wind. We will also expect a very luminous supernova like SN1998bw, if a large fraction of the conversion energy finally turns into the kinetic energy of the supernova ejecta.

GRB{980519} is characterized by its rapidly declining optical and X-ray afterglows. Explanations of this behavior include models invoking a dense medium environment which makes the shock wave evolve quickly into the sub-relativistic phase, a jet-like outflow, and a wind-shaped circumburst medium environment. Recently, Frail {et al}. (1999a) found that the latter two cases are consistent with the radio afterglow of this burst. Here, by considering the trans-relativistic shock hydrodynamics, we show that the dense medium model can also account for the radio light curve quite well. The potential virtue of the dense medium model for GRB{980519} is that it implies a smaller angular size of the afterglow, which is essential for interpreting the strong modulation of the radio light curve. Optical extinction due to the dense medium is not important if the prompt optical-UV flash accompanying the $\gamma$-ray emission can destroy dust by sublimation out to an appreciable distance. Comparisons with some other radio afterglows are also discussed.

Due to some refinements in the dynamics, we can follow the overall evolution of a realistic jet numerically till its bulk velocity being as small as $\beta c \sim 10^{-3} c$. We find no obvious break in the optical light curve during the relativistic phase itself. However, an obvious break does exist at the transition from the relativistic phase to the non-relativistic phase, which typically occurs at time $t \sim 10^6$ --- $10^{6.5}$ s (i.e., 10 --- 30 d). The break is affected by many parameters, such as the electron energy fraction $\xi_{\rm e}$, the magnetic energy fraction $\xi_{\rm B}^2$, the initial half opening angle $\theta_0$, and the medium number density $n$. Increase of any of them to a large enough value will make the break disappear. Although the break itself is parameter-dependent, afterglows from jetted GRB remnants are uniformly characterized by a quick decay during the non-relativistic phase, with power law timing index $\alpha \geq 2.1$. This is quite different from that of isotropic fireballs, and may be of fundamental importance for determining the degree of beaming in $\gamma$-ray bursts observationally.

Whether gamma-ray bursts are highly beamed or not is a very difficult but important problem that we are confronted with. Some theorists suggest that beaming effect usually leads to a sharp break in the afterglow light curve during the ultra-relativistic phase, with the breaking point determined by $\gamma = 1 / \theta_0$, where $\gamma$ is the Lorentz factor of the blastwave and $\theta_0$ is the initial half opening angle of the ejecta, but numerical studies tend to reject the suggestion. We note that previous studies are uniformly based on dynamics that is not proper for non-relativistic blastwaves. Here we investigate the problem in more detail, paying special attention to the transition from the ultra-relativistic phase to the non-relativistic phase. Due to some crucial refinements in the dynamics, we can follow the overall evolution of a realistic jet till its velocity is as small as $\beta c \sim 10^{-3} c$. We find no obvious break in the optical light curve during the relativistic phase itself. However, an obvious break does appear at the transition from the relativistic phase to the Newtonian phase if the physical parameters involved are properly assumed. Generally speaking, the Newtonian phase is characterized by a sharp decay of optical afterglows, with the power law timing index $\alpha \sim 1.8$ - 2.1. This is due to the quick lateral expansion at this stage. The quick decay of optical afterglows from GRB 970228, 980326, and 980519, and the breaks in the optical light curves of GRB 990123 and 990510 may indicate the presence of highly collimated $\gamma$-ray burst ejecta.

Recent observations on the GRB 000301c afterglow show that three breaks appear in the R-band light curve, and in particular the decay slope at late times is as steep as -3.0. This unusual afterglow is clearly inconsistent with the standard afterglow shock model. Here we propose a non-standard model for the unusual R-band afterglow of GRB 000301c. In this model, an ultra-relativistic shock in a dense medium ("dirty environment")rapidly evolved to the non-relativistic phase in initial 1 day. During such a phase, the shock happened to be caught up with by two energetic shells ejected from the central engine at two different times, and the shock was refreshed, leading to two flattenings of the light curve. After each interaction between the shock and shell, the afterglow decayed as $\propto t^{-3.0}$ if the electron distribution index of the shocked medium, $p\approx 3.4$, derived from the optical spectrum. Therefore, this model can provide an excellent explanation for the flattening and steepening features of the GRB 000301c optical afterglow light curve. We further point out that the energy injection shells ejected from the central engine at later times may be material shells (e.g., in the massive star progenitor models related to black holes) or radiation shells (e.g., in the millisecond pulsar progenitor models).

Based on the refined dynamical model proposed by us earlier for beamed $\gamma$-ray burst ejecta, we carry out detailed numerical procedure to study those $\gamma$-ray bursts with rapidly fading afterglows (i.e., $\sim t^{-2}$). It is found that optical afterglows from GRB 970228, 980326, 980519, 990123, 990510 and 991208 can be satisfactorily fitted if the $\gamma$-ray burst ejecta are highly collimated, with a universal initial half opening angle $\theta_0 \sim 0.1$. The obvious light curve break observed in GRB 990123 is due to the relativistic-Newtonian transition of the beamed ejecta, and the rapidly fading optical afterglows come from synchrotron emissions during the mildly relativistic and non-relativistic phases. We strongly suggest that the rapid fading of afterglows currently observed in some $\gamma$-ray bursts is evidence for beaming in these cases.

In a recent paper (Dai & Lu 1999), we have proposed a simple model in which the steepening in the light curve of the R-band afterglow of the gamma-ray burst (GRB) 990123 is caused by the adiabatic shock which has evolved from an ultrarelativistic phase to a nonrelativistic phase in a dense medium. We find that such a model is quite consistent with observations if the medium density is about $3\times 10^6 {\rm cm}^{-3}$. Here we discuss this model in more details. In particular, we investigate the effects of synchrotron self absorption and energy injection. A shock in a dense medium becomes nonrelativistic rapidly after a short relativistic phase. The afterglow from the shock at the nonrelativistic stage decays more rapidly than at the relativistic stage. Since some models for GRB energy sources predict that a strongly magnetic millisecond pulsar may be born during the formation of GRB, we discuss the effect of such a pulsar on the evolution of the nonrelativistic shock through magnetic dipole radiation. We find that after the energy which the shock obtains from the pulsar is much more than the initial energy of the shock, the afterglow decay will flatten significantly. When the pulsar energy input effect disappears, the decay will steepen again. These features are in excellent agreement with the afterglows of GRB 980519, GRB 990510 and GRB 980326. Furthermore, our model fits very well all the observational data of GRB 980519 including the last two detections.

Quiescence and burst emission and relativistic particle winds of soft gamma-ray repeaters (SGRs) have been widely interpreted to result from ultrastrongly magnetized neutron stars. In this magnetar model, the magnetic energy and gravitational energy of the neutron stars are suggested as the energy sources of all the emission and winds. However, Harding, Contopoulos & Kazanas (1999) have shown that the magnetic field should be only ~ 3 x 10^{13} G in order to match the characteristic spin-down timescale of SGR 1806-20 and its SNR age. Here we argue that if the magnetic field is indeed so weak, the previously suggested energy sources seem problematic. We further propose a plausible model in which SGR pulsars are young strange stars with superconducting cores and with a poloidal magnetic field of ~ 3 x 10^{13} G. In this model, the movement of the flux tubes not only leads to crustal cracking of the stars, giving rise to deconfinement of crustal matter to strange matter, but also to movement of internal magnetic toruses with the flux tubes. The former process will result in burst and quiescence emission and the latter process will produce steady relativistic winds which will power the surrounding supernova remnants.

Both observations and theoretical models suggest that the emissions in gamma-ray bursts (GRBs) and the afterglows are beamed. We argue that the recent polarization measured in the afterglows gives further evidence of beaming in GRBs. In this approach, we adopted the polarization-magnitude relation of BL Lacertae objects to 4 GRBs with available polarization measurements and found that the data of the 4 GRBs are consistent with the relation of BL Lacertae objects. This result suggests that the emissions in GRBs are beamed and there are similarities between GRBs and BL Lacertae objects (or generally blazars). We suggest that the weak polarization of GRBs may result from low percentage of orderly magnetic field in the shock front.

We here propose a two-step model for gamma-ray bursts (GRBs) associated with supernovae. In the first step, the core collapse of a star with mass $\ge 19M_\odot$ leads to a massive neutron star and a normal supernova, and subsequently hypercritical accretion of the neutron star from the supernova ejecta may give rise to a jet through neutrino annihilation and/or Poynting flux along the stellar rotation axis. However, because of too much surrounding matter, this jet rapidly enters a nonrelativistic phase and evolves to a large bubble. In the second step, the neutron star promptly implodes to a rapidly rotating black hole surrounded by a torus when the mass of the star increases to the maximum mass and meanwhile its rotation frequency increases to the upper limit due to the accreted angular momentum. The gravitational binding energy of the torus may be dissipated by a magnetized relativistic wind, which may then be absorbed by the supernova ejecta, thus producing an energetic hypernova. The rotational energy of the black hole may be extracted by the Blandford-Znajek's mechanism, leading to another jet. This jet is relatively free of baryons and thus may be accelerated to an ultrarelativistic phase because the first jet has pushed out of its front matter and left a baryon-free exit. We expect that the second jet generates a GRB and its afterglow. Our two-step model may alleviate the baryon-contamination problem suffered possibly from in the hypernova models. Furthermore, this model not only accounts for association of several GRBs with supernovae but also explains well the features of the afterglows of these bursts.

It is suggested that the inner energetic engine of Gamma ray burst (GRB) may be the result of the transition of normal hadron to quark-gluon plasma (QGP) in rapidly-rotating and spin-down newborn neutron star. When such a nascent neutron star slows down through dipole electromagnetic and quadruple gravitational radiation, the increasing center density may reach the QCD transition density, i.e., 5-10 nuclear density. Such kind of energy release from the phase transition would be responsible for GRB and its possible beaming effect. The relative dense gaseous environment of GRB location and the iron line observed in the X-ray afterglow support this idea. Some predictions in this model are given. keywords: Gamma-ray Bursts, Supernova, Quark-gluon plasma

Conventional dynamic model of gamma-ray burst remnants is found to be incorrect for adiabatic blastwaves during the non-relativistic phase. A new model is derived, which is shown to be correct for both radiative and adiabatic blastwaves during both ultra-relativistic and non-relativistic phase. Our model also takes the evolution of the radiative efficiency into account. The importance of the transition from the ultra-relativistic phase to the non-relativistic phase is stressed.

Recent observations show that the temporal decay of the R-band afterglow from GRB 990123 steepened about 2.5 days after the burst. We here propose a possible explanation for such a steepening: a shock expanding in a dense medium has undergone the transition from a relativistic phase to a nonrelativistic phase. We find that this model is consistent with the observations if the medium density is about $3\times 10^6 {\rm cm}^{-3}$. By fitting our model to the observed optical and X-ray afterglow quantitatively, we further infer the electron and magnetic energy fractions of the shocked medium and find these two parameters are about 0.1 and $2\times 10^{-8}$ respectively. The former parameter is near the equipartition value while the latter is about six orders of magnitude smaller than inferred from the GRB 970508 afterglow. We also discuss possibilities that the dense medium can be produced.

The conventional generic model is deemed to explain the dynamics of $\gamma$-ray burst remnants very well, no matter whether they are adiabatic or highly radiative. However, we find that for adiabatic expansion, the model could not reproduce the Sedov solution in the non-relativistic phase, thus the model needs to be revised. In the present paper, a new differential equation is derived. The generic model based on this equation has been shown to be correct for both radiative and adiabatic fireballs, and in both ultra-relativistic and non-relativistic phase.

Millisecond pulsars with strong magnetic fields may be formed through several processes, e.g. accretion-induced collapse of magnetized white dwarfs, merger of two neutron stars. During the birth of such a pulsar, an initial fireball available for a gamma-ray burst (GRB) may occur. We here study evolution of a postburst relativistic fireball with energy injection from the pulsar through magnetic dipole radiation, and find that the magnitude of the optical afterglow from this fireball first decreases with time, subsequently flattens, and finally declines again. This may provide a natural explanation for the behavior of the lightcurve of the afterglow of GRB970228 if this burst resulted from the birth of a strongly magnetic millisecond pulsar.

We here discuss a new model of $\gamma$-ray bursts (GRBs) based on differentially rotating strange stars. Strange stars in this model and differentially rotating neutron stars in the Klu\'zniak-Ruderman model can produce extremely relativistic, variable fireballs required by GRBs and then become millisecond pulsars. The effect of such pulsars on expansion of the postburst fireballs through magnetic dipole radiation is studied. We show that these two models can explain naturally not only various features of GRBs but also light curves of afterglows.

Postburst evolution of adiabatic fireballs that produce $\gamma$-ray bursts is studied. Emphasis has been put on the transition from the highly relativistic phase to the non-relativistic phase, which, according to our calculation, should happen much earlier than previously expected. The theoretical light curves decline a little sharper at the non-relativistic stage than at early times, but still can fit the observations well. However, disagreements are obvious when $t \ge 100$ d, implying a large initial energy $E_0$, a low interstellar density $n$, or the possible existence of a persistent energy source at the center of the fireball.

A GRB afterglow has been commonly thought to be due to continuous deceleration of a postburst fireball. Many analytical models have made simplifications for deceleration dynamics of the fireball and its radiation property, although they are successful at explaining the overall features of the observed afterglows. We here propose a model for a GRB afterglow in which the evolution of a postburst fireball is in an intermediate case between the adiabatic and highly radiative expansion. In our model, the afterglow is both due to the contribution of the adiabatic electrons behind the external blastwave of the fireball and due to the contribution of the radiative electrons. In addition, this model can describe evolution of the fireball from the extremely relativistic phase to the non-relativistic phase. Our calculations show that the fireball will go to the adiabatic expansion phase after about a day if the accelerated electrons are assumed to occupy the total internal energy. In all cases considered, the fireball will go to the mildly relativistic phase about $10^4$ seconds later, and to the non-relativistic phase after several days. These results imply that the relativistic adiabatic model cannot describe the deceleration dynamics of the several-days-later fireball. The comparison of the calculated light curves with the observed results at late times may imply the presence of impulsive events or energy injection with much longer durations.

The postburst evolution of fireballs that produce $\gamma$-ray bursts is studied, assuming the expansion of fireballs to be adiabatic and relativistic. Numerical results as well as an approximate analytic solution for the evolution are presented. Due to adoption of a new relation among $t$, $R$ and $\gamma$ (see the text), our results differ markedly from the previous studies. Synchrotron radiation from the shocked interstellar medium is attentively calculated, using a convenient set of equations. The observed X-ray flux of GRB afterglows can be reproduced easily. Although the optical afterglows seem much more complicated, our results can still present a rather satisfactory approach to observations. It is also found that the expansion will no longer be highly relativistic about 4 days after the main GRB. We thus suggest that the marginally relativistic phase of the expansion should be investigated so as to check the afterglows observed a week or more later.

The afterglow of a gamma-ray burst (GRB) is commonly thought to be due to continuous deceleration of a relativistically expanding fireball in the surrounding medium. Assuming that the expansion of the fireball is adiabatic and that the density of the medium is a power-law function of shock radius, viz., $n_{ext}\propto R^{-k}$, we analytically study the effects of the first-order radiative correction and the nonuniformity of the medium on a GRB afterglow. We first derive a new relation among the observed time, the shock radius and the fireball's Lorentz factor: $t_\oplus=R/4(4-k)\gamma^2c$, and also derive a new relation among the comoving time, the shock radius and the fireball's Lorentz factor: $t_{co}=2R/(5-k)\gamma c$. We next study the evolution of the fireball by using the analytic solution of Blandford and McKee (1976). The radiation losses may not significantly influence this evolution. We further derive new scaling laws both between the X-ray flux and observed time and between the optical flux and observed time. We use these scaling laws to discuss the afterglows of GRB 970228 and GRB 970616, and find that if the spectral index of the electron distribution is $p=2.5$, implied from the spectra of GRBs, the X-ray afterglow of GRB970616 is well fitted by assuming $k=2$.

The unusal hard x-ray burster GRO J1744-28 recently discovered by the Compton Gamma-ray Observatory (GRO) can be modeled as a strange star with a dipolar magnetic field $\le 10^{11} $Gauss. When the accreted mass of the star exceeds some critical mass, its crust may break, resulting in conversion of the accreted matter into strange matter and release of energy. Subsequently, a fireball may form and expand relativistically outward. The expanding fireball may interact with the surrounding interstellar medium, causing its kinetic energy to be radiated in shock waves, producing a burst of x-ray radiation. The burst energy, duration, interval and spectrum derived from such a model are consistent with the observations of GRO J1744-28.

We propose that when neutron stars in low-mass X-ray binaries accrete sufficient mass and become millisecond pulsars, the interiors of these stars may undergo phase transitions, which excite stellar radial oscillations. We show that the radial oscillations will be mainly damped by gravitational-wave radiation instead of internal viscosity. The gravitational waves can be detected by the advanced Laser Interferometer Gravitational-Wave Observatory at a rate of about three events per year.

We consider a model in which the soft gamma-ray repeaters (SGRs) result from young, magnetized strange stars with superconducting cores. As such a strange star spins down, the quantized vortex lines move outward and drag the magnetic field tubes because of the strong coupling between them. Since the terminations of the tubes interact with the stellar crust, the dragged tubes can produce sufficient tension to crack the crust. Part of the broken platelet will be dragged into the quark core, which is only $10^4$ cm from the surface, leading to the deconfinement of crustal matter into strange quark matter and thus the release of energy. We will show that the burst energy, duration, time interval and spectrum for our model are in agreement with the observational results. The persistent X-ray emission from the SGRs can be well explained by our model.

This paper argues that conversion of neutron stars to strange stars as an origin of cosmological stars in the binaries with low-mass companions. Our model may provide an explanation why the binary millisecond pulsars seem to have same low magnetic fields.