Gamma-ray bursts (GRBs) are known to be the most violent explosions in the universe, and a variety of correlations between observable GRB properties have been proposed in literature, but none of these correlations is valid for both long GRBs and short GRBs. In this paper we report the discovery of a universal correlation which is suitable for both long and short GRBs using three prompt emission properties of GRBs, i.e. the isotropic peak luminosity $L_{\rm iso}$, the peak energy of the time-integtated prompt emission spectrum $E_{\rm peak}$, and the "high signal" timescale $T_{\rm 0.45}$, $L_{\rm iso} \propto E_{\rm peak}^{1.94} T_{0.45}^{0.37}$. This universal correlation just involves properties of GRB prompt emission and does not require any information of afterglow phase, which can be used as a relatively unbiased redshift estimator. Here we use this correlation to estimate the pseudo-redshifts for short Gamma Ray Bursts and then use Lynden-Bell method to obtain a non-parametric estimate of their luminosity function and formation rate. The luminosity function is $\psi(L_0)\propto{L_0^{-0.63\pm{0.07}}}$ for dim SGRBs and $\psi(L_0)\propto{L_0^{-1.96\pm{0.28}}}$ for bright SGRBs, with the break point $6.95_{-0.76}^{+0.84}\times10^{50} erg/s$. The local formation rate of SGRBs is about 15 events $\rm Gpc^{-3}yr^{-1}$ . This universal correlation may have important implications for GRB physics, implying that the long and short GRBs should share similar radiation processes.

Luminosity of GRB 170817A is much lower than that of other sGRBs. The measurement of the superluminal movement of the radio afterglow emission confirms the presence of the relativistic jet, and the emission features can be well explained by the structured jet model. In this paper, we calculate the luminosity distribution of sGRBs and its evolution with redshift based on the structured (Gaussian) jet model, and find that the typical luminosity increase with redshift, for nearby sGRBs (such as for luminosity distance less than 200 Mpc) the typical gamma-ray luminosity is just around 10^47-10^48 erg s-1, which naturally explains the very low radiation luminosity of GRB 170817A. We derived the detection probability of sGRBs by Fermi-GBM and found that the expected detection rate of sGRBs is only about 1 yr-1 within the distance of several hundred Mpc. We explored the effect of the power-law index {\alpha} of the merger time distribution on the observed characteristics and found that it had little effect on the observed luminosity and viewing-angle distributions. However, it is very interesting that, for different values of {\alpha}, the distributions of the number of observed sGRBs are quite different, so it is possible to determine the value of {\alpha} through observed distributions of the number of sGRBs. We used the Bayesian method to make a quantitative analysis and found that the value of {\alpha} may be identified when the number of observed sGRBs with known redshifts is more than 200. Finally, we compare our results of gamma-ray luminosity distribution with sGRBs with known redshifts, and found that our results are consistent with the observation, which implies that our simulation results can reproduce the observed luminosity distribution well.

The post-merger gravitational wave radiation of the remnant formed in the double neutron star coalescence has not been directly measured, yet. We show in this work that the properties of the binary neutron stars involved in GW170817, additionally constrained by PSR J0030+0451, the lower limit on the maximum gravitational mass of non-rotating neutron star and some nuclear data, are in favor of strong post-merger gravitational wave radiation. This conclusion applies to the mergers of Galactic binary neutron star systems as well. Significant post-merger gravitational wave radiation is also preferred to improve the consistency between the maximum gravitational mass of the non-rotating neutron star inferred from GW170817/GRB170817A/AT2017gfo and the latest mass measurements of pulsars. The prominent post-merger gravitational radiation of GW170817-like events are expected to be detectable by advanced LIGO/Virgo detectors in the next decade and then shed valuable lights on the properties of the matter in the extremely high density.

For on-axis typical short gamma-ray bursts (sGRBs), the forward shock emission is usually so bright that renders the identification of kilonovae (also known as macronovae) in the early afterglow ($t<0.5$ d) phase rather challenging. This is why previously no thermal-like kilonova component has been identified at such early time except in the off-axis dim GRB 170817A associated with GW170817. Here we report the identification of an unusual optical radiation component in GRB 070809 at $t\sim 0.47$ d, thanks plausibly to the very-weak/subdominant forward shock emission. The optical emission with a very red spectrum is well in excess of the extrapolation of the X-ray emission that is distinguished by an unusually hard spectrum, which is at odds with the forward shock afterglow prediction but can be naturally interpreted as a kilonova. Our finding supports the speculation that kilonovae are ubiquitous , and demonstrates the possibility of revealing the neutron star merger origin with the early afterglow data of some typical sGRBs that take place well beyond the sensitive radius of the advanced gravitational wave detectors and hence the opportunity of organizing dedicated follow-up observations for events of interest.

Very recently the NICER collaboration has published the first-ever accurate measurement of mass and radius together for PSR J0030+0451, a nearby isolated quickly-rotating neutron star (NS). In this work we set the joint constraints on the equation of state (EoS) and some bulk properties of NSs with the data of PSR J0030+0451, GW170817 and some nuclear experiments. The piecewise polytropic expansion method and the spectral decomposition method have been adopted to parameterize the EoS. The resulting constraints are consistent with each other. Assuming the maximal gravitational mass of non-rotating NS $M_{\rm TOV}$ lies between $2.04 \rm M_{\odot}$ and $2.4 \rm M_{\odot}$, with the piecewise method the pressure at twice nuclear saturation density is measured to be $3.38^{+2.43}_{-1.50}\times 10^{34}~{\rm dyn~cm^{-2}}$ at the $90\%$ level. For a NS with canonical mass of $1.4 \rm M_\odot$, we have the moment of inertia $I_{1.4} = {1.43}^{+0.28}_{-0.13} \times 10^{38}~{\rm kg \cdot m^2}$, tidal deformability $\Lambda_{1.4} = 390_{-140}^{+320}$, radius $R_{1.4} = 12.2_{-0.9}^{+1.0}~{\rm km}$, and binding energy $BE_{1.4} = {0.16}^{+0.01}_{-0.02} \rm M_{\odot}$ at the $90\%$ level, which are improved in comparison to the constraints with the sole data of GW170817.

Assuming that the differential rotation of the massive neutron star (NS) formed in the double NS (DNS) mergers has been effectively terminated by the magnetic braking and a uniform rotation has been subsequently established (i.e., a supramassive NS is formed), we analytically derive in this work an approximated expression for the critical total gravitational mass ($M_{\rm tot,c}$) to form supramassive NS (SMNS) in the DNS mergers, benefited from some equation of state (EoS) insensitive relationships. The maximum gravitational mass of the nonrotating NSs ($M_{\rm TOV}$) as well as the dimensionless angular momentum of the remnant ($j$) play the dominant roles in modifying $M_{\rm tot,c}$, while the radius and mass differences of the premerger NSs do not. The GW170817/GRB 170817A/AT2017gfo observations have provided so far the best opportunity to quantitatively evaluate $M_{\rm TOV}$. Supposing the central engine for GRB 170817A is a black hole quickly formed in the collapse of an SMNS, we find $M_{\rm TOV}=2.13^{+0.09}_{-0.08}M_\odot$ (68.3% credibility interval, including also the uncertainties of the EoS insensitive relationships), which is consistent with the constraints set by current NS mass measurements.

The discovery of gravitational waves from compact objects coalescence opens a brand-new window to observe the universe. With more events being detected in the future, statistical examinations would be essential to better understand the underlying astrophysical processes. In this work we investigate the prospect of measuring the mass function of black holes that are merging with the neutron stars. Applying Bayesian parameter estimation for hundreds of simulated neutron star$-$black hole (NSBH) mergers, we find that the parameters for most of the injected events can be well recovered. We also take a Bayesian hierarchical model to reconstruct the population properties of the masses of black holes, in the presence of a low mass gap, both the mass gap and power-law index ($\alpha$) of black hole mass function can be well measured, thus we can reveal where the $\alpha$ is different for binary black hole (BBH) and NSBH systems. In the absence of a low mass gap, the gravitational wave data as well as the electromagnetic data can be used to pin down the nature of the merger event and then measure the mass of these very light black holes. However, as a result of the misclassification of BBH into NSBH, the measurement of $\alpha$ is more challenging and further dedicated efforts are needed.

We estimate the stellar mass for a sample of low surface brightness galaxies (LSBGs) by fitting their multiband spectral energy distributions (SEDs) to the stellar population synthesis (SPS) model. The derived stellar masses (log M*/Msun) span from 7.1 to 11.1, with a mean of log M*/Msun=8.5, which is lower than that for normal galaxies. The stellar mass-to-light ratio (M/L) in each band varies little with absolute magnitude, but increases with higher M*. This trend of M/L with M* is even stronger in bluer bands. In addition, the M/L for our LSBGs slightly declines from the r band to the longer wavelength bands. The log M/L_j (j=g,r,i,and z) have relatively tight relations with optical colors of g-r and g-i. Compared with several representative M/L-color relations (MLCRs) from literature, our MLCRs based on LSBG data are consistently among those literature MLCRs previously defined on diverse galaxy samples, and the existing minor differences between the MLCRs are more caused by the differences in the SED model ingredients including initial mass function, star formation history, and stellar population model, and the line fitting techniques, galaxy samples, and photometric zero-point as well, rather than the galaxy surface brightness itself which distinguishes LSBGs from HSBGs. Our LSBGs would be very likely to follow those representative MLCRs previously defined on diverse galaxy populations, if those main ingredients were taken into account.

While it is well known that mass transfer in binaries can pollute the surfaces of the accretors, it is still unclear whether this mechanism can reproduce the observed chemical inhomogeneities in globular clusters. We study the surface abundances of the accretors in low-mass binaries, as a first step towards understanding whether mass transfer in low-mass binaries is one of the potential origins of the aforementioned abundance anomalies in globular clusters. We use the Modules for Experiments in Stellar Astrophysics code to calculate binary evolutionary models with different initial donor masses between 0.9 and 1.9 $ \rm {M}_\odot$ for an initial metallicity of Z = 0.0034. The results show that in some low-mass binary systems, the accretors exhibit peculiar chemical patterns when they are still unevolved stars, e.g. C and O depletion; Na and N enhancement; constant Mg, Al and C+N+O. The abundance patterns of the accretors are significantly different from their initial abundances (or that of normal single stars), and can match the observed populations. These abundance patterns strongly depend not only on the initial parameters of binaries (donor mass, mass ratio, and orbital period), but also on the assumptions regarding mass-transfer efficiency and angular momentum loss. These results support the hypothesis that mass transfer in low-mass binaries is at least partly responsible for the unevolved anomalous stars in globular clusters. More work on binary evolutionary models and binary population synthesis is required to fully evaluate the contribution of this scenario.

GW190425 is the second neutron star merger event detected by the Advanced LIGO/Virgo detectors. If interpreted as a double neutron star merger, the total gravitational mass is substantially larger than that of the binary systems identified in the Galaxy. In this work we analyze the gravitational-wave data within the neutron star$-$black hole merger scenario. For the black hole, we yield a mass of $2.40^{+0.36}_{-0.32}M_\odot$ and an aligned spin of $0.141^{+0.067}_{-0.064}$. As for the neutron star we find a mass of $1.15^{+0.15}_{-0.13}M_\odot$ and the dimensionless tidal deformability of $1.4^{+3.8}_{-1.2}\times 10^{3}$. These parameter ranges are for 90\% credibility. The inferred masses of the neutron star and the black hole are not in tension with current observations and we suggest that GW190425 is a viable candidate of a neutron star$-$black hole merger event. Benefitting from the continual enhancement of the sensitivities of the advanced gravitational detectors and the increase of the number of the observatories, similar events are anticipated to be much more precisely measured in the future and the presence of black holes below the so-called mass gap will be unambiguously clarified. If confirmed, the mergers of neutron stars with (quickly rotating) low-mass black holes are likely important production sites of the heaviest r-process elements.

For some neutron stars (NSs) in the binary systems, the masses have been accurately measured. While for the isolated neutron stars (INSs), no mass measurement has been reported yet. The situation will change soon thanks to the successful performance of the Neutron Star Interior Composition Explorer (NICER), with which the radius and mass of the isolated PSR J0030+0451 can be simultaneously measured. For most INSs, no mass measurements are possible for NICER because of observational limitations. Benefiting from recent significant progress made on constraining the equation of state of NSs, in this work we propose a way to estimate the masses of the INSs with the measured gravitational redshifts. We apply our method to RX J1856.5-3754, RX J0720.4-3125, and RBS 1223, three members of "The Magnificent Seven" (M7), and estimate their masses to be $1.24_{-0.29}^{+0.29}M_{\odot}$, $1.23_{-0.05}^{+0.10}M_{\odot}$, and $1.08_{-0.11}^{+0.20}M_{\odot}$, respectively. These masses are consistent with that of binary NS systems, suggesting no evidence for experiencing significant accretion of these isolated objects.

In this work we parameterize the Equation of State of dense neutron star (NS) matter with four pressure parameters of $\{\hat{p}_1, \hat{p}_2, \hat{p}_3, \hat{p}_4\}$ and then set the combined constraints with the data of GW 170817 and the data of 6 Low Mass X-ray Binaries (LMXBs) with thermonuclear burst or alternatively the symmetry energy of the nuclear interaction. We find that the nuclear data effectively narrow down the possible range of $\hat{p}_1$, the gravitational wave data plays the leading role in bounding $\hat{p}_2$, and the LMXB data as well as the lower bound on maximal gravitational mass of non-rotating NSs govern the constraints on $\hat{p}_3$ and $\hat{p}_4$. Using posterior samples of pressure parameters and some universal relations, we further investigate how the current data sets can advance our understanding of tidal deformability ($\Lambda$), moment of inertia ($I$) and binding energy ($BE$) of NSs. For a canonical mass of $1.4M_\odot$, we have $I_{1.4} = {1.43}^{+0.30}_{-0.13} \times 10^{38}~{\rm kg \cdot m^2}$, $\Lambda_{1.4} = 390_{-210}^{+280}$ , $R_{1.4} = 11.8_{-0.7}^{+1.2}~{\rm km}$ and $BE_{1.4} = {0.16}^{+0.01}_{-0.02} M_{\odot}$ if the constraints from the nuclear data and the gravitational wave data have been jointly applied. For the joint analysis of gravitational wave data and the LMXB data, we have $I_{1.4} = {1.28}^{+0.15}_{-0.08} \times 10^{38}~{\rm kg \cdot m^2}$, $\Lambda_{1.4} = 220_{-90}^{+90}$, $R_{1.4} = 11.1_{-0.6}^{+0.7}~{\rm km}$ and $BE_{1.4} = {0.18}^{+0.01}_{-0.01} M_{\odot}$. These results suggest that the current constraints on $\Lambda$ and $R$ still suffer from significant systematic uncertainties, while $I_{1.4}$ and $BE_{1.4}$ are better constrained.

We report the detection of GeV $\gamma$-ray emission from the very-high-energy (VHE) $\gamma$-ray source VER J2227+608 associated with the "tail" region of SNR G106.3+2.7. The GeV $\gamma$-ray emission is extended and spatially coincident with molecular clouds traced by CO emission. The broadband GeV to TeV emission of VER J2227+608 can be well fitted by a single power-law function with an index of 1.90$\pm$0.04, without obvious indication of spectral cutoff toward high energies. The pure leptonic model for the $\gamma$-ray emission can be marginally ruled out by the X-ray and TeV data. In the hadronic model, the low energy content of CRs and the hard $\gamma$-ray spectrum, in combination with the center-bright source structure, suggest that VER J2227+608 may be powered by the PWN instead of shocks of the SNR. And the cutoff energy of the proton distribution needs to be higher than $\sim$ 400 TeV, which makes it an attractive PeVatron candidate. Future observations by the upcoming Large High Altitude Air Shower Observatory (LHAASO) and the Cherenkov Telescope Array in the north (CTA-North) could distinguish these models and constrain the maximum energy of cosmic rays in supernova remnants.

The precise measurement of the spectrum of protons, the most abundant component of the cosmic radiation, is necessary to understand the source and acceleration of cosmic rays in the Milky Way. This work reports the measurement of the cosmic ray proton fluxes with kinetic energies from 40 GeV to 100 TeV, with two and a half years of data recorded by the DArk Matter Particle Explorer (DAMPE). This is the first time an experiment directly measures the cosmic ray protons up to ~100 TeV with a high statistics. The measured spectrum confirms the spectral hardening found by previous experiments and reveals a softening at ~13.6 TeV, with the spectral index changing from ~2.60 to ~2.85. Our result suggests the existence of a new spectral feature of cosmic rays at energies lower than the so-called knee, and sheds new light on the origin of Galactic cosmic rays.

The DArk Matter Particle Explorer (DAMPE), a satellite-based cosmic ray and gamma-ray detector, was launched on December 17, 2015, and began its on-orbit operation on December 24, 2015. In this work we document the on-orbit calibration procedures used by DAMPE and report the calibration results of the Plastic Scintillator strip Detector (PSD), the Silicon-Tungsten tracKer-converter (STK), the BGO imaging calorimeter (BGO), and the Neutron Detector (NUD). The results are obtained using Galactic cosmic rays, bright known GeV gamma-ray sources, and charge injection into the front-end electronics of each sub-detector. The determination of the boundary of the South Atlantic Anomaly (SAA), the measurement of the live time, and the alignments of the detectors are also introduced. The calibration results demonstrate the stability of the detectors in almost two years of the on-orbit operation.

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.

The $\gamma$-ray observation of interstellar gas provides a unique way to probe the cosmic rays (CRs) outside the solar system. In this work, we use an updated version of Fermi-LAT data and recent multi-wavelength tracers of interstellar gas to re-analyze a mid-latitude region in the third Galactic quadrant and estimate the local CR proton spectrum. Two $\gamma$-ray production cross section models for $pp$ interaction, the commonly used one from Kamae et al. (2006) and the up-to-date one from Kafexhiu et al. (2014), are adopted separately in the analysis. Both of them can well fit the emissivity and the derived proton spectra roughly resemble the direct measurements from AMS-02 and Voyager 1, but rather different spectral parameters are indicated. A break at $4\pm1~{\rm GeV}\;c^{-1}$ is shown if the cross section model by Kamae et al. (2006) is adopted. The resulting spectrum is $\lesssim 20\%$ larger than the AMS-02 observation above $15~\rm GeV$ and consistent with the de-modulated spectrum within $2\%$. The proton spectrum based on the cross section model of Kafexhiu et al. (2014) is about $1.4-1.8$ times that of AMS-02 at $2-100~\rm GeV$, however the difference decreases to $20\%$ below $10~\rm GeV$ with respect to the de-modulated spectrum. A spectral break at $20\pm11~{\rm GeV}\;c^{-1}$ is required in this model. An extrapolation down to $300~\rm MeV$ is performed to compare with the observation of Voyager 1, and we find a deviation of $\lesssim 2.5\sigma$ for both the models. In general, an approximately consistent CR spectrum can be obtained using $\gamma$-ray observation nowadays, but we still need a better $\gamma$-ray production cross section model to derive the parameters accurately.

As revealed recently by the modeling of the multi-wavelength data of the emission following GW170817/GRB 170817A, there was an off-axis energetic relativistic outflow component launched by this historic double neutron star merger event. In this work we use the results of these modeling to examine the energy extraction process of the central engine. We show that the magnetic process (i.e., the Blandford-Znajek mechanism) is favored, while the neutrino process usually requires a too massive accretion disk if the duration of the central engine activity is comparable to the observed $T_{90}$ of GRB 170817A, unless the timescale of the central engine activity is less than $\sim$ 0.2s. We propose that the GRB observations are helpful to constrain the combined tidal parameter $\tilde{\Lambda}$, and by adopting the accretion disk mass distribution estimated in BZ mechanism, the $90\%$ credible interval of $\tilde{\Lambda}$ for the progenitor of GW170817 is inferred as $309-954$.

GW170817, the first neutron star merger event detected by advanced LIGO/Virgo detectors, was associated with an underluminous short duration GRB 170817A. In this work we compare the forward shock afterglow emission of GW170817/GRB 170817A to other luminous short GRBs (sGRBs) with both a known redshift and an afterglow emission lasting at least one day after the burst. In the rapid decay phase, the afterglow emission of the bright sGRBs and GW170817/GRB 170817A form a natural and continuous sequence, though separated by an observation time gap. If viewed on-axis, the forward shock afterglow emission of GW170817/GRB 170817A would be among the brightest ones detected so far. This provides a strong evidence for the GW170817-like merger origin of bright sGRBs, and suggests that the detection of the forward shock afterglow emission of most neutron star merger events are more challenging than the case of GW170817 since usually the mergers will be more distant and the viewing angles are plausibly higher.

The precise observations of Galactic cosmic ray fluxes of the secondary family, such as Li, Be, B, are expected to have significant implications on our understanding of the cosmic ray origin and propagation. Here we employ the recent very precise measurements of those species by the Alpha Magnetic Spectrometer on the International Space Station, together with their parent species (C and O), as well as the data collected by the Voyager-1 spacecraft outside the heliosphere and the Advanced Composition Explorer, we investigate the propagation of cosmic rays in the Milky Way. We find that these new data strongly favor the existence of reacceleration of cosmic rays during the propagation rather than the convective transport. We further find that for the reacceleration model, the slope of rigidity dependence of the diffusion coefficient is close to the Kolmogorov-type of turbulence, and the spectral hardenings of both the primary and secondary particles can be well described by the injection hardening rather than the propagation hardening.

HESS J1640-465 is an extended TeV $\gamma$-ray source and its $\gamma$-ray emission whether from the shell of a supernova remnant (SNR) or a pulsar wind nebula (PWN) is still under debate. We reanalyze the GeV $\gamma$-ray data in the field of HESS J1640-465 using eight years of Pass 8 data recorded by the Fermi Large Area Telescope. An extended GeV $\gamma$-ray source positionally coincident with HESS J1640-465 is found. Its photon spectrum can be described by a power-law with an index of $1.42\pm0.19$ in the energy range of 10-500 GeV, and smoothly connects with the TeV spectrum of HESS J1640-465. The broadband spectrum of HESS J1640-465 can be well fit by a leptonic model with a broken power-law spectrum of electrons with an exponential cut-off at $\sim$ 300 TeV. The spectral properties of HESS J1640-465 are broadly consistent with the characteristics of other sources identified as PWNe, such as the correlations between high-energy luminosity ratios and the physical parameters of pulsar, including spin-down luminosity $\dot{E}$ and characteristic age $\tau_c$. All these pieces of evidence support that the $\gamma$-ray emission of HESS J1640-465 may originate from the PWN powered by PSR J1640-4631 rather than the shell of the SNR G338.3-0.0.

The acceleration of cosmic ray particles and their propagation in the Milky Way and the heliosphere tangle with each other, leading to complexity and degeneracy of the modeling of Galactic cosmic rays (GCRs). The recent measurements of the GCR spectra by Voyager-1 from outside of the heliosphere gave a first direct observation of GCRs in the local interstellar (LIS) environment. Together with the high-precision data near the Earth taken by ACE and AMS-02, we derive the LIS spectra of Helium, Lithium, Beryllium, Boron, Carbon, and Oxygen nuclei from a few MeV/n to TeV/n, using a non-parameterization method. These LIS spectra are helpful in further studying the injection and propagation parameters of GCRs. The nearly 20 years of data recorded by ACE are used to determine the solar modulation parameters over the solar cycles 23 and 24, based on the force-field approximation. We find general agreements of the modulation potential with the results inferred from neutron monitors and other cosmic ray data.

So far, in total 15 double neutron star systems (DNSs) with a reliable measurement of the total gravitational mass ($M_{\rm T}$) have been detected in the Galaxy. In this work we study the distribution of $M_{\rm T}$. The data prefer the double Gaussian distribution over a single Gaussian distribution and the low and high mass populations center at $M_{\rm T}\sim 2.58M_\odot$ and $\sim 2.72M_\odot$, respectively. The progenitor stars of GW170817 have a $M_{\rm T}=2.74^{+0.04}_{-0.01}M_\odot$, falling into the high mass population. With a local neutron star merger rate of $\sim 10^{3}~{\rm Gpc^{-3}~yr^{-1}}$, supposing the $M_{\rm T}$ of those merging neutron stars also follow the double Gaussian distribution, the upcoming runs of the advanced LIGO/Virgo will soon detect some events with a $M_{\rm T}\lesssim 2.6M_\odot$ that can effectively probe the equation of state of the neutron stars and the distribution function is expected to be reliably reconstructed in the next decade.

The gravitational wave data of GW170817 favor the equation of state (EoS) models that predict compact neutron stars (NSs), consistent with the radius constraints from X-ray observations. Motivated by such a remarkable progress, we examine the fate of the remnants formed in NS mergers and focus on the roles of the angular momentum and the mass distribution of the binary NSs. In the mass shedding limit (for which the dimensionless angular momentum equals to the Keplerian value, i.e., $j=j_{\rm Kep}$), the adopted { seven EoS models, except H4 and ALF2,} yield supramassive NSs in more than half of the mergers. However, for $j\lesssim 0.7j_{\rm Kep}$, the presence or absence of a non-negligible fraction of supramassive NSs formed in the mergers depends sensitively on both the EoS and the mass distribution of the binary systems. The NS mergers with a total gravitational mass $\leq 2.6M_\odot$ are found to be able to shed valuable light on both the EoS model and the angular momentum of the remnants if supramassive NSs are still absent. We have also discussed the uncertainty on estimating the maximum gravitational mass of non-rotating NSs ($M_{\rm max}$) due to the unknown $j$ of the pre-collapse remnants. With the data of GW170817 and the assumption of the mass loss of $0.03M_\odot$, we have $M_{\rm max}<(2.19,~2.32)M_\odot$ (90\% confidence level) for $j=(1.0,~0.8)j_{\rm Kep}$, respectively.

The origin of the IceCube neutrinos is still an open question. Upper limits from diffuse gamma-ray observations suggest that the neutrino sources are either distant or hidden from gamma-ray observations. It is possible that the neutrinos are produced in jets that are formed in the core-collapsing massive stars and fail to break out, the so-called choked jets. We study neutrinos from the jets choked in the hydrogen envelopes of red supergiant stars. Fast photo-meson cooling softens the neutrino spectrum, making it difficult to explain the PeV neutrinos observed by IceCube in a one-component scenario, but a two-component model can explain the spectrum. Furthermore, we predict that a newly born jet-driven type-II supernova may be observed to be associated with a neutrino burst detected by IceCube.

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

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

The $\gamma$-ray emission from molecular clouds is widely believed to have a hadronic origin, while unequivocal evidence is still lacking. In this work, we analyze the Fermi-LAT Pass 8 publicly available data accumulated from 2008 August 4 to 2017 June 29 and report the significant detection of a characteristic $\pi^0$-decay feature in the $\gamma$-ray spectra of the molecular clouds Orion A and Orion B. This detection provides a direct evidence for the hadronic origin of their $\gamma$-ray emission.

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

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

High energy cosmic ray electrons plus positrons (CREs), which lose energy quickly during their propagation, provide an ideal probe of Galactic high-energy processes and may enable the observation of phenomena such as dark-matter particle annihilation or decay. The CRE spectrum has been directly measured up to $\sim 2$ TeV in previous balloon- or space-borne experiments, and indirectly up to $\sim 5$ TeV by ground-based Cherenkov $\gamma$-ray telescope arrays. Evidence for a spectral break in the TeV energy range has been provided by indirect measurements of H.E.S.S., although the results were qualified by sizeable systematic uncertainties. Here we report a direct measurement of CREs in the energy range $25~{\rm GeV}-4.6~{\rm TeV}$ by the DArk Matter Particle Explorer (DAMPE) with unprecedentedly high energy resolution and low background. The majority of the spectrum can be properly fitted by a smoothly broken power-law model rather than a single power-law model. The direct detection of a spectral break at $E \sim0.9$ TeV confirms the evidence found by H.E.S.S., clarifies the behavior of the CRE spectrum at energies above 1 TeV and sheds light on the physical origin of the sub-TeV CREs.

The DArk Matter Particle Explorer (DAMPE), one of the four scientific space science missions within the framework of the Strategic Pioneer Program on Space Science of the Chinese Academy of Sciences, is a general purpose high energy cosmic-ray and gamma-ray observatory, which was successfully launched on December 17th, 2015 from the Jiuquan Satellite Launch Center. The DAMPE scientific objectives include the study of galactic cosmic rays up to $\sim 10$ TeV and hundreds of TeV for electrons/gammas and nuclei respectively, and the search for dark matter signatures in their spectra. In this paper we illustrate the layout of the DAMPE instrument, and discuss the results of beam tests and calibrations performed on ground. Finally we present the expected performance in space and give an overview of the mission key scientific goals.

Short Gamma-Ray Bursts(SGRBs) are widely believed to be from mergers of binary compact objects involving at least one neutron star and hence have a broad range of spatial offsets from their host galaxies. In this work we search for possible correlations between the emission properties of 18 SGRBs and their offsets from the host galaxies. The SGRBs with and without extended emission do not show significant difference between their offset distribution, in agreement with some previous works. There are however possible correlations between the optical and X-ray afterglow emission and the offsets. The underlying physical origins are examined.

Puppis A is a very famous and extensively studied supernova remnant (SNR) that shows strong evidence of shock-cloud interaction. We re-analyze the GeV $\gamma$-ray emission of it using seven years Pass 8 data recorded by the Fermi Large Area Telescope (Fermi-LAT). The morphology of the $\gamma$-ray emission is more compatible with that of the thermal X-ray and IR emissions than the radio image, which suggests a possible correlation between the gamma-ray emitting region and dense clouds. The $\gamma$-ray spectrum in the energy range of 1-500 GeV shows a break at $7.92\pm1.91$ GeV with the photon indices of $1.81\pm0.08$ below the break and $2.53\pm0.12$ above the break, which can naturally explain the lack of TeV $\gamma$-ray emission from Puppis A. The multi-wavelength observations favor a hadronic origin for the $\gamma$-ray emission.

Intermediate mass black holes (IMBHs) with a mass between $10^{2}$ and $10^{5}$ times that of the sun, which bridges the {mass gap between the} stellar-mass black holes and the supermassive black holes, are crucial in understanding the evolution of the black holes. Although they are widely believed to exist, decisive evidence has long been absent. Motivated by the successful detection of massive stellar-mass black holes by advanced LIGO, through the gravitational wave radiation during the binary merger, in this work we investigate the prospect of detecting/identifying the lightest IMBHs (LIMBHs; the black holes $\gtrsim 100M_\odot$) with the second generation gravitational wave detectors. In general, the chance of hearing the birth of the LIMBHs is significantly higher than that to identify pre-merger IMBHs. The other formation channel of LIMBHs, where stars with huge mass/low-metallicity directly collapse, is likely "silent", so the merger-driven birth of the LIMBHs may be the only "observable" scenario in the near future. Moreover, the prospect of establishing the presence of (lightest) intermediate mass black holes in the O3 run and beyond of advanced LIGO is found quite promising, implying that such an instrument could make another breakthrough on astronomy in the near future. The joining of other detectors like advanced Virgo would only increase the detection rate.

GRB 160625B is an extremely-bright outburst with well-monitored afterglow emission. The geometry-corrected energy is high up to $\sim 5.2\times10^{52}$ erg or even $\sim 8\times 10^{52}$ erg, rendering it the most energetic GRB prompt emission recorded so far. We analyzed the time-resolved spectra of the prompt emission and found that in some intervals there were likely thermal-radiation components and the high energy emission were characterized by significant cutoff. The bulk Lorentz factors of the outflow material are estimated accordingly. We found out that the Lorentz factors derived in the thermal-radiation model are consistent with the luminosity-Lorentz factor correlation found in other bursts as well as in GRB 090902B for the time-resolved thermal-radiation components. While the spectral cutoff model yields much lower Lorentz factors that are in tension with the constraints set by the electron pair Compoton scattering process. We then suggest that these spectral cutoffs are more likely related to the particle acceleration process and that one should be careful in estimating the Lorentz factors if the spectrum cuts at a rather low energy (e.g., $\sim$ tens MeV). The nature of the central engine has also been discussed and a stellar-mass black hole is likely favored.

The modeling of Li-Paczy\'{n}ski macronova/kilonova signals gives reasonable estimate on the neutron-rich material ejected during the neutron star mergers. Usually the accretion disk is more massive than the macronova ejecta, with which the efficiencies of converting the disk mass into prompt emission of three merger-driven GRBs can hence be directly constrained. Supposing the macronovae/kilonovae associated with GRB 050709, GRB 060614 and GRB 130603B arose from radioactive decay of the r-process material, the upper limit on energy conversion efficiencies are found to be as low as $\sim 10^{-6}-10^{-4}$. Moreover, for all three events, neutrino annihilation is likely powerful enough to account for the brief gamma-ray flashes. Neutrino annihilation can also explain the "extended" emission lasting $\sim 100$ s in GRB 050709, but does not work for the one in GRB 060614. These progresses demonstrate that the macronova can serve as a novel probe of the central engine activity.

GRB 140903A, a short duration $\gamma-$ray burst (SGRB) detected by {\it Swift}, is characterized by its long-lasting radio emission among SGRBs. In addition to the $\sim 10^{6}$ s radio afterglow emission, the afterglow of GRB 140903A displays a plateau from $10^3$ s to $7 \times 10^3 $ s in the X-rays. In this work, we attribute the X-ray plateau to the energy injection into the decelerating blast wave and then model the later radio/optical/X-ray afterglow emission within the standard fireball afterglow model. The afterglow emission has been well reproduced with reasonable physical parameters, including a jet half-opening angle $\sim 0.05$. We also compare the physical properties of SGRBs with and without radio afterglow detections and find that nearby SGRBs are brighter and easier to be detected in radio. Such a phenomena has been interpreted within the afterglow model.

Recently Jin et al. (2016) reanalyzed the optical observation data of GRB 050709 and reported a line-like spectral energy distribution (SED) component observed by the Very Large Telescope at $t\sim 2.5$ days after the trigger of the burst, which had been interpreted as a broadened line signal arising from a macronova dominated by iron group. In this work we show that an optical flare origin of such a peculiar optical SED is still possible. Interestingly, even in such a model, an "unusual" origin of the late-time long-lasting Hubble Space Telescope F814W-band emission is still needed and a macronova/kilonova is the natural interpretation.

Neutron star$-$black hole (NS$-$BH) coalescences are widely believed to be promising gravitational wave sources in the era of advanced detectors of LIGO/Virgo but such binaries have never been directly detected yet. Evidence for NS$-$BH coalescences have been suggested in short and hybrid GRB observations, which are examined critically. Based on the suggested connection between the observed macronovae/kilonovae events and NS$-$BH coalescences, we get a fiducial lower limit of NS$-$BH coalescence rate density ${\cal R}_{\rm nsbh} \approx 18.8^{+12.5}_{-8.6} ~{\rm Gpc^{-3}~ yr^{-1}~ (\theta_j/0.1~{\rm rad})^{-2}}$, where $\theta_{\rm j}$ is the typical half-opening angle of the GRB ejecta. The real value of ${\cal R}_{\rm nsbh}$ is likely at least $\sim {\rm a~few}$ times larger, depending upon the equation of state of NS material and the properties of the NS$-$BH system, such as the mass and spin distribution of the black hole. If the link between macronovae/kilonovae and NS$-$BH coalescence is valid, one can expect that at design sensitivity the aLIGO/AdVirgo network will detect NS$-$BH coalescence signals at a rate of at least a dozen per year, and to consequently place constraints on certain physical properties of NS$-$BH systems.

GRB 050709 was the first short Gamma-ray Burst (sGRB) with an identified optical afterglow. In this work we report a re-analysis of the publicly available data of this event and the discovery of a Li-Paczynski macronova/kilonova that dominates the optical/IR signal at t $>$ 2.5 days. Such a signal would arise from ~0.05 M$_\odot$ r-process material launched by a compact binary merger. The implied mass ejection supports the suggestion that compact binary mergers are significant and possibly main sites of heavy r-process nucleosynthesis. We have re-analyzed all afterglow data from nearby short and hybrid GRBs. A statistical study of sGRB/macronova connection reveals that macronova may have taken place in all these GRBs though the fraction as low as 0.18 cannot be ruled out. The identification of two of the three macronova candidates in the I-band implies a more promising detection prospect for ground-based surveys.

Short-duration gamma-ray bursts (SGRBs) are widely believed to be powered by the mergers of compact binaries, such as binary neutron stars or possibly neutron star-black hole binaries. Though the prospect of detecting SGRBs with gravitational wave (GW) signals by the advanced Laser Interferometer Gravitational-Wave Observatory (LIGO)/VIRGO network is promising, no known SGRB has been found within the expected advanced LIGO/VIRGO sensitivity range for binary neutron star systems. We find, however, that the two long-short GRBs (GRB 060505 and GRB 060614) may be within the horizon of advanced GW detectors. In the upcoming era of GW astronomy, the merger origin of some long-short GRBs, as favored by the macronova signature displayed in GRB 060614, can be unambiguously tested. The model-dependent time lags between the merger and the onset of the prompt emission of the GRB are estimated. The comparison of such time lags between model predictions and the real data expected in the era of the GW astronomy would be helpful in revealing the physical processes taking place at the central engine (including the launch of the relativistic outflow, the emergence of the outflow from the dense material ejected during the merger, and the radiation of gamma rays). We also show that the speed of GWs, with or without a simultaneous test of Einstein's equivalence principle, can be directly measured to an accuracy of $\sim 3\times 10^{-8}~{\rm cm~s^{-1}}$ or even better in the advanced LIGO/VIRGO era. The Astrophysical Journal, Volume

The merger-driven Gamma-ray Bursts (GRBs) and their associated gravitational wave (GW) radiation, if both successfully detected, have some far-reaching implications, including for instance: (i) The statistical comparison of the physical properties of the short/long-short GRBs with and without GW detection can test the general origin model; (ii) Revealing the physical processes taking place at the central engine; (iii) Measuring the velocity of the Gravitational wave directly/accurately. In this work we discuss these implications in the case of possible association of GW150914/ GBM transient 150914. We compared GBM transient 150914 with other SGRBs and found that such an event {may be} a distinct outlier in some statistical diagrams, possibly due to its specific binary-black-hole merger origin. However, the presence of a "new" group of SGRBs with "unusual" physical parameters is also possible. If the outflow of GBM transient 150914 was launched by the accretion onto the nascent black hole, the magnetic activity rather than the neutrino process is likely responsible for the energy extraction and the accretion disk mass is estimated to be $\sim 10^{-5}~M_\odot$. The GW150914/GBM transient 150914 association, {if confirmed, would} provide the first opportunity to directly measure the GW velocity and its departure from the speed of the light {should be within} a factor of $\sim 10^{-17}$.

Supernova remnants (SNRs) are the most attractive candidates for the acceleration sites of Galactic cosmic rays. We report the detection of GeV $\gamma$-ray emission with the Pass 8 events recorded by Fermi Large Area Telescope (Fermi-LAT) in the vicinity of the shell type SNR CTB 37B that is likely associated with the TeV $\gamma-$ray source HESS J1713-381. The photon spectrum of CTB 37B is consistent with a power-law with an index of $1.89\pm0.08$ in the energy range of $0.5-500$ GeV, and the measured flux connects smoothly with that of HESS J1713-381 at a few hundred GeV. No significant spatial extension and time variation are detected. The multi-wavelength data can be well fitted with either a leptonic model or a hadronic one. However, parameters of both models suggest more efficient particle acceleration than typical SNRs. Meanwhile, the X-ray and $\gamma$-ray spectral properties of CTB 37B show that it is an interesting source bridging young SNRs dominated by non-thermal emission and old SNRs interacting with molecular clouds.

After the jet break at $t\sim 1.4$ days, the optical afterglow emission of the long-short burst GRB 060614 can be divided into two components. One is the power-law decaying forward shock afterglow emission. The other is an excess of flux in several multi-band photometric observations, which emerges at $\sim$4 days after the burst, significantly earlier than that observed for a supernova associated with a long-duration GRB. At $t>13.6$ days, the F814W-band flux drops faster than $t^{-3.2}$. Moreover, the spectrum of the excess component is very soft and the luminosity is extremely low. These observed signals are incompatible with those from weak supernovae, but the ejection of $\sim 0.1~M_\odot$ of $r-$process material from a black hole-neutron star merger, as recently found in some numerical simulations, can produce it. If this interpretation is correct, it represents the first time that a multi-epoch/band lightcurve of a Li-Paczynski macronova (also known as kilonova) has been obtained and black hole-neutron star mergers are sites of significant production of $r-$process elements.

The {\it Swift}-detected GRB 060614 was a unique burst that straddles an imaginary divide between long- and short-duration gamma-ray bursts (GRBs), and its physical origin has been heavily debated over the years. Recently, a distinct very-soft F814W-band excess at $t\sim 13.6$ days after the burst was identified in a joint-analysis of VLT and HST optical afterglow data of GRB~060614, which has been interpreted as evidence for an accompanying Li-Paczynski macronova (also called a kilonova). Under the assumption that the afterglow data in the time interval of $1.7-3.0$ days after the burst are due to external forward shock emission, when this assumption is extrapolated to later times it is found that there is an excess of flux in several multi-band photometric observations. This component emerges at $\sim$4 days after the burst, and it may represent the first time that a multi-epoch/band lightcurve of a macronova has been obtained. The macronova associated with GRB 060614 peaked at $t\lesssim 4$ days after the burst, which is significantly earlier than that observed for a supernova associated with a long-duration GRB. Due to the limited data, no strong evidence for a temperature evolution is found. We derive a conservative estimate of the macronova rate of $\sim 16.3^{+16.3}_{-8.2}~{\rm Gpc^{-3}}{\rm yr^{-1}}$, implying a promising prospect for detecting the gravitational wave radiation from compact object mergers by upcoming Advanced LIGO/VIRGO/KAGRA detectors (i.e., the rate is ${\cal R}_{\rm GW} \sim 0.5^{+0.5}_{-0.25}(D/200~{\rm Mpc})^{3}~{\rm yr^{-1}}$).

Long-duration ($>2$ s) $\gamma$-ray bursts that are believed to originate from the death of massive stars are expected to be accompanied by supernovae. GRB 060614, that lasted 102 s, lacks a supernova-like emission down to very stringent limits and its physical origin is still debated. Here we report the discovery of near-infrared bump that is significantly above the regular decaying afterglow. This red bump is inconsistent with even the weakest known supernova. However, it can arise from a Li-Paczy\'{n}ski macronova $-$ the radioactive decay of debris following a compact binary merger. If this interpretation is correct GRB 060614 arose from a compact binary merger rather than from the death of a massive star and it was a site of a significant production of heavy r-process elements. The significant ejected mass favors a black hole-neutron star merger but a double neutron star merger cannot be ruled out.

The short burst GRB 130912A was detected by Swift, Fermi satellites and several ground-based optical telescopes. Its X-ray light curve decayed with time normally. The optical emission, however, displayed a long term plateau, which is the longest one in current short GRB observations. In this work we examine the physical origin of the X-ray and optical emission of this peculiar event. We find that the canonical forward shock afterglow emission model can account for the X-ray and optical data self-consistently and the energy injection model that has been widely adopted to interpret the shallowly-decaying afterglow emission is not needed. We also find that the burst was born in a very-low density interstellar medium, consistent with the compact object merger model. Significant fractions of the energy of the forward shock have been given to accelerate the non-thermal electrons and amplify the magnetic fields (i.e., $\epsilon_{\rm e}\sim 0.37$ and $\epsilon_{\rm B}\sim 0.16$, respectively), which are much larger than those inferred in most short burst afterglow modeling and can explain why the long-lasting optical afterglow plateau is rare in short GRBs.

The existence of fast radio bursts (FRBs), a new type of extragalatic transients, has been established recently and quite a few models have been proposed. In this work we discuss the possible connection between the FRB sources and ultra-high energy ($>10^{18}$ eV) cosmic rays. We show that in the blitzar model and the model of merging binary neutron stars, the huge energy release of each FRB central engine together with the rather high rate of FRBs, the accelerated EeV cosmic rays may contribute significantly to the observed ones. In other FRB models including for example the merger of double white dwarfs and the energetic magnetar radio flares, no significant EeV cosmic ray is expected. We also suggest that the mergers of double neutron stars, even if they are irrelevant to FRBs, may play a non-ignorable role in producing EeV cosmic ray protons if supramassive neutron stars were formed in a good fraction of mergers and the merger rate is $\gtrsim 10^{3}~{\rm yr^{-1}~ Gpc^{-3}}$. Such a possibility will be unambiguously tested in the era of gravitational wave astronomy.

GRB 120308A, a long duration $\gamma-$ray burst detected by {\it Swift}, was distinguished by a highly-polarized early optical afterglow emission that strongly suggests an ordered magnetic field component in the emitting region. In this work we model the optical and X-ray emission in the reverse and forward shock scenario and show that the strength of the magnetic field in reverse shock region is $\sim 10$ times stronger than that in the forward shock region. Consequently the outflow powering the highly-polarized reverse shock optical emission was mildly-magnetized at a degree $\sigma \sim$ a few percent. Considering the plausible magnetic energy dissipation in both the acceleration and prompt emission phases of the Gamma-ray Burst (GRB) outflow, the afterglow data of GRB 120308A provides us the compelling evidence that at least for some GRBs a non-ignorable fraction of the energy was released in the form of Poynting-flux, confirming the finding firstly made in the reverse-forward shock emission modeling of the optical afterglow of GRB 990123 (Fan et al. 2002; Zhang et al. 2003).

We discuss the possibility of using fast radio bursts (FRBs), if cosmological, as a viable cosmic probe. We find out that the contribution of the host galaxies to the detected dispersion measures can be inapparent for the FRBs not from galaxy centers or star forming regions. The inhomogeneity of the intergalactic medium (IGM), however, causes significant deviation of the dispersion measure from that predicted in the simplified homogeneous IGM model for individual event. Fortunately, with sufficient FRBs along different sightlines but within a very narrow redshift interval (e.g., $\Delta z \sim 0.05$ or $\Delta z \sim 0.05(1+z)$), the mean from averaging observed dispersion measures does not suffer such a problem and hence may be used as a cosmic probe. We show that in the optimistic case (e.g., tens FRBs in each $\Delta z$ have been measured; the most distant FRBs were at redshift $\geq 3$; the host galaxies and the FRB sources contribute little to the detected dispersion measures) and with all the uncertainties (i.e. the inhomogeneity of the IGM, the contribution and uncertainty of host galaxies as well as the evolution and error of $f_{\rm IGM}$) considered, FRBs could help constrain the equation of state of dark energy.

GRB 131231A was detected by the Large Area Telescope onboard Fermi Space Gamma-ray Telescope. The high energy gamma-ray ($> 100$ MeV) afterglow emission spectrum is $F_\nu \propto \nu^{-0.54\pm0.15}$ in the first $\sim 1300$ s after the trigger and the most energetic photon has an energy $\sim 62$ GeV arriving at $t\sim 520$ s. With reasonable parameters of the GRB outflow as well as the density of the circum-burst medium, the synchrotron radiation of electrons or protons accelerated at an external forward shock have difficulty accounting for the data. The synchrotron self-Compton radiation of the forward shock-accelerated electrons, instead, can account for both the spectrum and temporal behavior of the GeV afterglow emission. We also show that the prospect for detecting GRB 131231A$-$like GRBs with Cherenkov Telescope Array (CTA) is promising.

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.

The principal goal of this paper is to use attempts at reconciling the Swift long gamma-ray bursts (LGRBs) with the star formation history (SFH) to compare the predictions of $\Lambda$CDM with those in the $R_{\rm h}=ct$ Universe. In the context of the former, we confirm that the latest Swift sample of GRBs reveals an increasing evolution in the GRB rate relative to the star formation rate (SFR) at high redshifts. The observed discrepancy between the GRB rate and the SFR may be eliminated by assuming a modest evolution parameterized as $(1+z)^{0.8}$---perhaps indicating a cosmic evolution in metallicity. However, we find a higher metallicity cut of $Z=0.52Z_{\odot}$ than was seen in previous studies, which suggested that LGRBs occur preferentially in metal poor environments, i.e., $Z\sim0.1-0.3Z_{\odot}$. We use a simple power-law approximation to the high-\emph{z} ($\ga 3.8$) SFH, i.e., $R_{\rm SF}\propto[(1+z)/4.8]^{\alpha}$, to examine how the high-\emph{z} SFR may be impacted by a possible abundance evolution in the Swift GRB sample. For an expansion history consistent with $\Lambda$CDM, we find that the Swift redshift and luminosity distributions can be reproduced with reasonable accuracy if $\alpha=-2.41_{-2.09}^{+1.87}$. For the $R_{\rm h}=ct$ Universe, the GRB rate is slightly different from that in $\Lambda$CDM, but also requires an extra evolutionary effect, with a metallicity cut of $Z=0.44Z_{\odot}$. Assuming that the SFR and GRB rate are related via an evolving metallicity, we find that the GRB data constrain the slope of the high-\emph{z} SFR in $R_{\rm h}=ct$ to be $\alpha=-3.60_{-2.45}^{+2.45}$. Both cosmologies fit the GRB/SFR data rather well.

GRB 130821A was detected by Fermi-GBM/LAT, Konus-Wind, SPI-ACS/INTEGRAL, RHESSI and Mars Odyssey-HEND. Although the data of GRB 130821A are very limited, we show in this work that the high energy gamma-ray emission (i.e., above 100 MeV) alone imposes tight constraint on the density profile of the circum-burst medium as well as the initial Lorentz factor of the outflow. The temporal behavior of the high energy gamma-ray emission is consistent with the forward shock synchrotron radiation model and the circum-burst medium likely has a constant-density profile. The Lorentz factor is about a few hundred, similar to other bright GRBs.

We show that the peculiar early optical and in particular X-ray afterglow emission of the short duration burst GRB 130603B can be explained by continuous energy injection into the blastwave from a supra-massive magnetar central engine. The observed energetics and temporal/spectral properties of the late infrared bump (i.e., the "kilonova") are also found consistent with emission from the ejecta launched during an NS-NS merger and powered by a magnetar central engine. The isotropic-equivalent kinetic energies of both the GRB blastwave and the kilonova are about $E_{\rm k}\sim 10^{51}$ erg, consistent with being powered by a near-isotropic magnetar wind. However, this relatively small value demands that most of the initial rotational energy of the magnetar $(\sim {\rm a~ few \times 10^{52}~ erg})$ is carried away by gravitational wave radiation. Our results suggest that (i) the progenitor of GRB 130603B would be a NS-NS binary system, whose merger product would be a supra-massive neutron star that lasted for about $\sim 1000$ seconds; (ii) the equation-of-state of nuclear matter would be stiff enough to allow survival of a long-lived supra-massive neutron star, so that it is promising to detect bright electromagnetic counterparts of gravitational wave triggers without short GRB associations in the upcoming Advanced LIGO/Virgo era.

Short Gamma-Ray Bursts (GRBs), brief intense emission of $\gamma-$rays characterized by a duration shorter than 2 seconds that are plausibly powered by the coalescence of binary neutron stars, are believed to be strong gravitational wave radiation (GWR) sources. The test of such a speculation has been thought to be impossible until the performance of the detectors like advanced LIGO. Recently there has been growing evidence for the formation of highly-magnetized neutron star (i.e., magnetar) in the double neutron star mergers. In this work we re-examine the interpretation of the X-ray plateau followed by an abrupt decline detected in some short GRB afterglows within the supramassive magnetar model and find that the maximum gravitational mass of the non-rotating neutron stars is $\sim 2.3M_{\odot}$ and the observed duration of some X-ray plateaus are significantly shorter than that expected in the magnetic dipole radiation scenario, suggesting that the collapse of the supramassive magnetars has been considerably enhanced by the energy loss via GWR. Such a result demonstrates that the signature of GWR may have already existed in current electromagnetic data of short GRBs.

Cosmological time dilation is a fundamental phenomenon in an expanding universe, which stresses that both the duration and wavelength of the emitted light from a distant object at the redshift $z$ will be dilated by a factor of $1+z$ at the observer. By using a sample of 139 \emph{Swift} long GRBs with known redshift ($z\leq8.2$), we measure the observed duration ($T_{90}$) in the observed energy range between $140/(1+z)$ keV and $350/(1+z)$ keV, corresponding to a fixed energy range of 140-350 keV in the rest frame. We obtain a significant correlation between the duration and the factor $1+z$, i.e., $T_{\rm{90}}=10.5(1+z)^{0.94\pm0.26}$, which is well consistent with that expected from cosmological time dilation effect.

Short-hard gamma-ray bursts (GRBs) are widely believed to be produced by the merger of two binary compact objects, specifically by two neutron stars or by a neutron star orbiting a black hole. According to the Li-Paczynski kilonova model, the merger would launch sub-relativistic ejecta and a near-infrared/optical transient would then occur, lasting up to days, which is powered by the radioactive decay of heavy elements synthesized in the ejecta. The detection of a late bump using the {\em Hubble Space Telescope} ({\em HST}) in the near-infrared afterglow light curve of the short-hard GRB 130603B is indeed consistent with such a model. However, as shown in this Letter, the limited {\em HST} near-infrared lightcurve behavior can also be interpreted as the synchrotron radiation of the external shock driven by a wide mildly relativistic outflow. In such a scenario, the radio emission is expected to peak with a flux of $\sim 100 \mu$Jy, which is detectable for current radio arrays. Hence, the radio afterglow data can provide complementary evidence on the nature of the bump in GRB 130603B. It is worth noting that good spectroscopy during the bump phase in short-hard bursts can test validity of either model above, analogous to spectroscopy of broad-lined Type Ic supernova in long-soft GRBs.

A nearby super-luminous burst GRB 130427A was simultaneously detected by six $\gamma$-ray space telescopes ({\it Swift}, Fermi-GBM/LAT, Konus-Wind, SPI-ACS/INTEGRAL, AGILE and RHESSI) and by three RAPTOR full-sky persistent monitors. The isotropic $\gamma-$ray energy release is of $\sim 10^{54}$ erg, rendering it the most powerful explosion among the GRBs with a redshift $z\leq 0.5$. The emission above 100 MeV lasted about one day and four photons are at energies greater than 40 GeV. We show that the count rate of 100 MeV-100 GeV emission may be mainly accounted for by the forward shock synchrotron radiation and the inverse Compton radiation likely dominates at GeV-TeV energies. In particular, an inverse Compton radiation origin is established for the $\sim (95.3,~47.3,~41.4,~38.5,~32)$ GeV photons arriving at $t\sim (243,~256.3,~610.6,~3409.8,~34366.2)$ s after the trigger of Fermi-GBM. Interestingly, the external-inverse-Compton-scattering of the prompt emission (the second episode, i.e., $t\sim 120-260$ s) by the forward-shock-accelerated electrons is expected to produce a few $\gamma-$rays at energies above 10 GeV, while five were detected in the same time interval. A possible unified model for the prompt soft $\gamma-$ray, optical and GeV emission of GRB 130427A, GRB 080319B and GRB 090902B is outlined. Implication of the null detection of $>1$ TeV neutrinos from GRB 130427A by IceCube is discussed.

In addition to the two ~1 PeV neutrinos, the IceCube Collaboration recently reported a detection of 26 neutrino candidates at energies from 30 TeV to 250 TeV, implying a confidence level of 4.3\sigma over the atmospheric background. We suggest that these TeV-PeV non-atmospheric neutrinos may originate from two groups of sources, motivated by the non-detection of neutrinos in the energy range 250 TeV- 1 PeV in current data. If intrinsic, the non-detection of 250 TeV-1 PeV neutrinos disfavors the single power-law spectrum model for the TeV-PeV non-atmospheric neutrinos at a confidence level of ~ 2\sigma. We then interpret the current neutrino data with a two-component spectrum model. One has a flat spectrum with a cutoff at the energy ~ 250 TeV and the other has a sharp peak at ~1 PeV. The former is likely via pp collision while the latter may be generated by the photomeson interaction.

Ultra-luminous infrared galaxies (ULIRGs) are the most luminous and intense starburst galaxies in the Universe. Both their star-formation rate (SFR) and gas surface mass density are very high, implying a high supernovae rate and an efficient energy conversion of energetic protons. A small fraction of these supernovae is the so-called hypernovae with a typical kinetic energy ~1e52 erg and a shock velocity >=1e9 cm/s. The strong shocks driven by hypernovae are able to accelerate cosmic ray protons up to 1e17 eV. These energetic protons lose a good fraction of their energy through proton-proton collision when ejected into very dense interstellar medium, and as a result, produce high energy neutrinos (<=5 PeV). Recent deep infrared surveys provide solid constraints on the number density of ULIRGs across a wide redshift range 0<z<2.3, allowing us to derive the flux of diffuse neutrinos from hypernovae. We find that at PeV energies, the diffuse neutrinos contributed by ULIRGs are comparable with the atmosphere neutrinos with the flux of 2e-9GeV cm^-2/s/sr, by assuming the injected cosmic ray power law spectrum with an index of -2.

GRB 120422A is a low-luminosity Gamma-ray burst (GRB) associated with a bright supernova, which distinguishes itself by its relatively short T90 ~ 5 s and an energetic X-ray tail. We analyze the Swift BAT and XRT data and discuss the physical implications. We show that the early steep decline in the X-ray light curve can be interpreted as the curvature tail of a late emission episode around 58-86 s, with a curved instantaneous spectrum at the end of the emission episode. Together with the main activity in the first ~ 20 s and the weak emission from 40 s to 60 s, the prompt emission is variable, which points towards a central engine origin, in contrast to the shock breakout origin as invoked to interpret some other nearby low-luminosity supernova GRBs. The curvature effect interpretation and interpreting the early shallow decay as the coasting external forward shock emission in a wind medium both give a constraint on the bulk Lorentz factor \Gamma to be around several. Comparing the properties of GRB 120422A and other supernova GRBs, we found that the main criterion to distinguish engine-driven GRBs from the shock breakout GRBs is the time-averaged luminosity, with a separation line of about ~ 10^48 erg s^-1.

We show that the experiential $E_{\rm p}-L$, $\Gamma-L$, $E_{\rm p}-\Gamma$ and ${\rm \bar{\eta}_\gamma}-E_{\rm p}$ correlations (where $L$ is the time-averaged luminosity of the prompt emission, $E_{\rm p}$ is the spectral peak energy, $\Gamma$ is the bulk Lorentz factor and $\bar{\eta}_\gamma$ is the emission efficiency of Gamma-ray bursts) are well consistent with the relations between the resembling parameters predicted in the photospheric radiation model of the prompt emission of Gamma-ray bursts. The time-resolved thermal radiation of GRB 090902B does follow the $E_{\rm p}-L$ and $\Gamma-L$ correlations. A reliable interpretation of the four correlations in alternative models is still lacking. These facts may point towards a photospheric origin of prompt emission of some Gamma-ray Bursts.

Gamma-ray bursts (GRBs) are brief but intense emission of soft $\gamma-$rays, mostly lasting from a few seconds to a few thousand seconds. For such kind of high energy transients, their isotropic-equivalent-energy ($E_{\rm iso}$) function may be more scientifically meaningful when compared with GRB isotropic-equivalent-luminosity function ($L_{\rm iso}$), as the traditional luminosity function refers to steady emission much longer than a few thousand seconds. In this work we for the first time construct the isotropic-equivalent-energy function for a sample of 95 bursts with measured redshifts ($z$) and find an excess of high-$z$ GRBs. Assuming that the excess is caused by a GRB luminosity function evolution in a power-law form, we find a cosmic evolution of $E_{\rm iso}\propto(1+z)^{1.80^{+0.36}_{-0.63}}$, which is comparable to that between $L_{\rm iso}$ and $z$, i.e., $L_{\rm iso}\propto(1+z)^{2.30^{+0.56}_{-0.51}}$ (both $1\sigma$). The evolution-removed isotropic-equivalent-energy function can be reasonably fitted by a broken power-law, in which the dim and bright segments are $\psi(E_{\rm iso})\propto E_{\rm iso}^{-0.27\pm0.01}$ and $\psi(E_{\rm iso})\propto E_{\rm iso}^{-0.87\pm0.07}$, respectively ($1\sigma$). For the cosmic GRB formation rate, it increases quickly in the region of $0 \leq z \lesssim 1$, and roughly keeps constant for $1\lesssim z \lesssim 4$, and finally falls with a power index of $-3.80\pm2.16$ for $z\gtrsim 4$, in good agreement with the observed cosmic star formation rate so far.

We perform a statistical analysis of the temporal and spectral properties of the latest Fermi gamma-ray bursts (GRBs) to revisit the classification of GRBs. We find that the bimodalities of duration and the energy ratio ($E_{\mathrm{peak}}$/Fluence) and the anti-correlation between spectral hardness (hardness ratio ($HR$), peak energy and spectral index) and duration ($T_{90}$) support the long/soft $-$ short/hard classification scheme for Fermi GRBs. The $HR - T_{90}$ anti-correlation strongly depends upon the spectral shape of GRBs and energy bands, and the bursts with the curved spectra in the typical BATSE energy bands show a tighter anti-correlation than those with the power-law spectra in the typical BAT energy bands. This might explain why the $HR - T_{90}$ correlation is not evident for those GRB samples detected by instruments like {\it Swift} with a narrower/softer energy bandpass. We also analyze the intrinsic energy correlation for the GRBs with measured redshifts and well defined peak energies. The current sample suggests $E_{\mathrm{p,rest}}=2455\times (E_{\mathrm{iso}}/10^{52})^{0.59}$ for short GRBs, significantly different from that for long GRBs. However, both the long and short GRBs comply with the same $E_{\mathrm{p,rest}}-L_{\mathrm{iso}}$ correlation.

In order to study the effect of dust extinction on the afterglow of gamma-ray bursts (GRBs), we carry out numerical calculations with high precision based on rigorous Mie theory and latest optical properties of interstellar dust grains, and analyze the different extinction curves produced by dust grains with different physical parameters. Our results indicate that the absolute extinction quantity is substantially determined by the medium density and metallicity. However, the shape of the extinction curve is mainly determined by the size distribution of the dust grains. If the dust grains aggregate to form larger ones, they will cause a flatter or grayer extinction curve with lower extinction quantity. On the contrary, if the dust grains are disassociated to smaller ones due to some uncertain processes, they will cause a steeper extinction curve with larger amount of extinction. These results might provide an important insight into understanding the origin of the optically dark GRBs.

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.

In this work we estimate the accretion-disk mass in the specific scenario of binary-neutron-star-merger with current observational data. Assuming that the outflows of short Gamma-ray Bursts (GRBs) are driven via neutrino-antineutrino annihilation we estimate the disk mass of about half of short bursts in the sample to be $\sim 0.01-0.1 M_\odot$, in agreement with that obtained in the numerical simulations. Massive disks ($\sim {\rm several} 0.1 M_\odot$) found in some other short GRBs may point to the more efficient magnetic process of energy extraction or the neutron star and black hole binary progenitor. Our results suggest that some short bursts may be really due to the coalescence of double neutron stars and are promising gravitational wave radiation sources. For future short GRBs with simultaneous gravitational-wave detections, the disk mass may be reliably inferred and the validity of our approach will be tested. We also propose a method to constrain the initial radius of a baryonic outflow where it is launched ($R_0$) without the need of identifying an ideal thermal spectrum. We then apply it to GRB 090510 and get that $R_0 \lesssim 6.5\times 10^{6}(\Gamma_{\rm ph}/2000)^{-4}$ cm, suggesting that the central engine is a black hole with a mass $< 22 M_\odot(\Gamma_{\rm ph}/2000)^{-4}$, where $\Gamma_{\rm ph}$ is the bulk Lorentz factor of the outflow at the photospheric radius.

Superconducting cosmic strings may be able to power energetic gamma-ray bursts, in particular at high-redshifts ($z>6$). The ultra-relativistic outflows driven by the superconducting cosmic strings, however, seem too narrow to be consistent with the afterglow data of GRB 090423 and GRB 080913.

Sw 1644+57/GRB 110328A is a remarkable cosmological X-ray outburst detected by the {\it Swift} satellite. Its early-time ($t\lesssim 0.1$ days since the trigger) X-ray emission resembles some gamma-ray bursts (GRBs), e.g., GRB 090417B. But the late-time flaring X-ray plateau lasting $> 40$ days renders it unique. We examine the possibilities that the outburst is a super-long GRB powered either by the fallback accretion onto a nascent black hole or by a millisecond pulsar, and find out that these two scenarios can address some but not all of the main observational features. We then focus on the model of tidal disruption of a (giant) star by a massive black hole. The mass of the tidal-disrupted star is estimated to be $\gtrsim$ a few solar masses. A simple/straightforward argument for a magnetic origin of the relativistic outflow is presented.

We develop two methods to estimate the bulk Lorentz factor of X-ray flare outflow. In the first method the outflow is assumed to be baryonic and is accelerated by the thermal pressure, for which the final bulk Lorentz factor is limited by the outflow luminosity as well as the initial radius of the outflow getting accelerated. Such a method may be applied to a considerable fraction of flares. The second method, based on the curvature effect interpretation of the quick decline of the flare, can give a tightly constrained estimate of the bulk Lorentz factor but can only be applied to a few giant flares. The results obtained in these two different ways are consistent with each other. The obtained bulk Lorentz factor (or just upper limit) of the X-ray flare outflows, ranging from ten to a few hundred, is generally smaller than that of the Gamma-ray Burst (GRB) outflows.

With its rapid response, {\it Swift} has revealed plenty of unexpected properties of gamma-ray bursts (GRBs). With an abundance of observations, our current understanding is only limited by complexity of early X-ray light curves. In this work, based on the public {\it Swift} data of 150 well-monitored GRBs with measured redshifts, we find some interesting global features in the rest-frame X-ray light curves. The distinct spectral behaviors between the prompt emission and the afterglow emission implies dissimilar radiation scenarios. Interestingly, an unforeseen plateau is exhibited in the prompt X-ray light curves despite the presence of complex spikes, which might indicate the presence of a steady central engine. In particular, the seemingly continuous evolution with a single power law from the prompt to the afterglow of most GRBs might place strong constraints on the theoretical models.

A possible birefringence effect that arises in quantum gravity leads to a frequency-dependent rotation of the polarization angle of linearly polarized emission from distant sources. Here we use the UV/optical polarization data of the afterglows of GRB 020813 and GRB 021004 to constrain this effect. We find an upper limit on the Gambini & Pulin birefringence parameter $| \eta | <2\times 10^{-7}$. This limit is of 3 orders better than the previous limits from observations of AGNs and of the Crab pulsar. Much stronger limits may be obtained by the future observation of polarization of the prompt $\gamma$-rays.

We present the prompt BAT and afterglow XRT data of Swift-discovered GRB081109A up to ~ 5\times 10^5 sec after the trigger, and the early ground-based optical follow-ups. The temporal and spectral indices of the X-ray afterglow emission change remarkably. We interpret this as the GRB jet first traversing the freely expanding supersonic stellar wind of the progenitor with density varying as $\rho \propto r^{-2}$. Then after approximately 300 sec the jet traverses into a region of apparent constant density similar to that expected in the stalled-wind region of a stellar wind bubble or the interstellar medium (ISM). The optical afterglow data are generally consistent with such a scenario. Our best numerical model has a wind density parameter {$A_{*} \sim 0.02$, a density of the stalled wind $n\sim 0.12 {\rm cm}^{-3}$, and a transition radius $ \sim 4.5 \times 10^{17}$ cm}. Such a transition radius is smaller than that predicted by numerical simulations of the stellar wind bubbles and may be due to a rapidly evolving wind of the progenitor close to the time of its core-collapse.

The synchrotron self-Compton (SSC) emission from Gamma-ray Burst (GRB) forward shock can extend to the very-high-energy (VHE; $E_\gamma > $100 GeV) range. Such high energy photons are rare and are attenuated by the cosmic infrared background before reaching us. In this work, we discuss the prospect to detect these VHE photons using the current ground-based Cherenkov detectors. Our calculated results are consistent with the upper limits obtained with several Cherenkov detectors for GRB 030329, GRB 050509B, and GRB 060505 during the afterglow phase. For 5 bursts in our nearby GRB sample (except for GRB 030329), current ground-based Cherenkov detectors would not be expected to detect the modeled VHE signal. Only for those very bright and nearby bursts like GRB 030329, detection of VHE photons is possible under favorable observing conditions and a delayed observation time of $\la$10 hours.

We intend to determine the type of circumburst medium and measure directly the initial Lorentz factor $\Gamma_0$ of GRB outflows. If the early X-ray afterglow lightcurve has a peak and the whole profile across the peak is consistent with the standard external shock model, the early rise profile of light curves can be used to differentiate whether the burst was born in interstellar medium (ISM) or in stellar wind. In the thin shell case, related to a sub-relativistic reverse shock, the peak time occurring after the end of the prompt emission, can be used to derive an accurate $\Gamma_0$, especially for the ISM case. The afterglow lightcurves for a flat electron spectrum $1<p<2$ have been derived analytically. In our GRB sample, we obtain $\Gamma_0 \sim 300$ for the bursts born in ISM. We did not find any good case for bursts born in stellar wind and behaving as a thin shell that can be used to constrain $\Gamma_0$ reliably.

For an unsteady baryonic gamma-ray burst (GRB) outflow, the fast and slow proton shells collide with each other and produce energetic soft gamma-ray emission. If the outflow has a significant neutron component, the ultra-relativistic neutrons initially expand freely until decaying at a larger radius. The late time proton shells ejected from the GRB central engine, after powering the regular internal shocks, will sweep these $\beta-$decay products and give rise to very bright UV/optical emission. The naked-eye optical flash from GRB 080319B, an energetic explosion in the distant universe, can be well explained in this way.

Gamma-ray Bursts (GRBs) are among the potential extragalactic sources of very-high-energy (VHE) gamma-rays. We discuss the prospects of detecting VHE gamma-rays with current ground-based Cherenkov instruments during the afterglow phase. Using the fireball model, we calculate the synchrotron self-Compton (SSC) emission from forward-shock electrons. The modeled results are compared with the observational afterglow data taken with and/or the sensitivity level of ground-based VHE instruments (e.g. STACEE, H.E.S.S., MAGIC, VERITAS, and Whipple). We find that modeled SSC emission from bright and nearby bursts such as GRB 030329 are detectable by these instruments even with a delayed observation time of ~10 hours.

In the neutron-rich internal shocks model for gamma-ray bursts (GRBs), the Lorentz factors (LFs) of ion shells are variable, and so are the LFs of accompanying neutron shells. For slow neutron shells with a typical LF of approximate tens, the typical $\beta$-decay radius is $\sim 10^{14}-10^{15}$ cm. As GRBs last long enough [$T_{90}>14(1+z)$ s], one earlier but slower ejected neutron shell will be swept successively by later ejected ion shells in the range $\sim10^{13}-10^{15}$ cm, where slow neutrons have decayed significantly. Part of the thermal energy released in the interaction will be given to the electrons. These accelerated electrons will be mainly cooled by the prompt soft $\gamma-$rays and give rise to GeV emission. This kind of GeV emission is particularly important for some very long GRBs and is detectable for the upcoming satellite {\it Gamma-Ray Large Area Space Telescope} (GLAST).

Although it is well recognized that gamma-ray burst (GRB) afterglows are obscured and reddened by dust in their host galaxies, the wavelength-dependence and quantity of dust extinction are still poorly known. Current studies on this mostly rely on fitting the afterglow spectral energy distributions (SEDs) with template extinction models. The inferred extinction (both quantity and wavelength-dependence) and dust-to-gas ratios are often in disagreement with that obtained from dust depletion and X-ray spectroscopy studies. We argue that this discrepancy could result from the prior assumption of a template extinction law. We propose an analytical formula to approximate the GRB host extinction law. With the template extinction laws self-contained, and the capability of revealing extinction laws differing from the conventional ones, it is shown that this is a powerful approach in modeling the afterglow SEDs to derive GRB host extinction.

In the synchrotron radiation model, the polarization property depends on both the configuration of the magnetic field and the geometry of the visible emitting region. Some peculiar behaviors in the X-ray afterglows of {\it Swift} gamma-ray bursts (GRBs), such as energetic flares and the plateau followed by a sharp drop, might by highly linearly-polarized because the outflows powering these behaviors may be Poynting-flux dominated. Furthermore, the broken-down of the symmetry of the visible emitting region may be hiding in current X-ray data and will give rise to interesting polarization signatures. In this work we focus on the polarization accompanying the very early sharp decline of GRB X-ray afterglows. We show that strong polarization evolution is possible in both the high latitude emission model and the dying central engine model which are used to interpret this sharp X-ray decline. It is thus not easy to efficiently probe the physical origin of the very early X-ray sharp decline with future polarimetry. Strong polarization evolution is also possible in the decline phase of X-ray flares and in the shallow decline phase of X-ray light curves characterized by chromatic X-ray VS. Optical breaks. An {\it XRT}-like detector but with polarization capability on board a {\em Swift}-like satellite would be suitable to test our predictions.

We calculate the very high energy (sub-GeV to TeV) inverse Compton emission of GRB afterglows. We argue that this emission provides a powerful test of the currently accepted afterglow model. We focus on two processes: synchrotron self-Compton (SSC) emission within the afterglow blast wave, and external inverse Compton (EIC) emission which occurs when flare photons (produced by an internal process) pass through the blast wave. We show that if our current interpretations of the Swift XRT data are correct, there should be a canonical high energy afterglow emission light curve. Our predictions can be tested with high energy observatories such as GLAST, Whipple, H.E.S.S. and MAGIC. Under favorable conditions we expect afterglow detections in all these detectors.

The dust extinction of gamma-ray bursts (GRBs) host galaxies, containing important clues to the nature of GRB progenitors and crucial for dereddening, is still poorly known. Here we propose a straightforward method to determine the extinction of GRB host galaxies by comparing the observed optical spectra to the intrinsic ones extrapolated from the X-ray spectra. The rationale for this method is from the standard fireball model: if the optical flux decay index equals to that of the X-ray flux, then there is no break frequency between the optical and X-ray bands, therefore we can derive the intrinsic optical flux from the X-ray spectra. We apply this method to three GRBs of which the optical and X-ray fluxes have the same decay indices and another one with inferred cooling break frequency, and obtain the rest-frame extinction curves of their host galaxies. The derived extinction curves are gray and do not resemble any extinction curves of local galaxies (e.g. the Milk Way, the Small/Large Magellanic Clouds, or nearby starburst galaxies). The amount of extinction is rather large (with visual extinction $A_V$ $\sim$ 1.6--3.4$\magni$). We model the derived extinction curves in terms of the silicate-graphite interstellar grain model. As expected from the ``gray'' nature of the derived extinction curve, the dust size distribution is skewed to large grains. We determine, for the first time, the local dust-to-gas ratios of GRB host galaxies using the model-derived dust parameters and the hydrogen column densities determined from X-ray absorptions.

Before 2004, nearly all GRB afterglow data could be understood in the context of the external shocks model. This situation has changed in the past two years, when it became clear that some afterglow components should be attributed to the activity of the central engine; i.e., the {\it central engine afterglow}. We review here the afterglow emission that is directly related to the GRB central engine. Such an interpretation proposed by Katz, Piran & Sari, peculiar in pre-{\it Swift} era, has become generally accepted now.

High energy emission (> tens MeV) of Gamma-Ray Bursts (GRBs) provides an important clue to understand the physical processes involved in GRBs, which may be correlated with the GRB early afterglow. A shallow decline phase has been well detected in about half {\it Swift} Gamma-ray Burst X-ray afterglows. The widely considered interpretation involves a significant energy injection and possibly time-evolving shock parameter(s). This work we calculate the synchrotron-self-Compton (SSC) radiation of such an external forward shock and show that it could explain the well-known long term high energy (i.e., tens MeV to GeV) afterglow of GRB 940217. We propose that the cooperation of Swift and GLAST will help to reveal the nature of GRBs.

In the double neutron star merger or neutron star-black hole merger model for short GRBs, the outflow launched might be mildly magnetized and neutron rich. The magnetized neutron-rich outflow will be accelerated by the magnetic and thermal pressure and may form a two component jet finally, as suggested by Vlahakis, Peng & K\"{o}nigl (2003). We show in this work that such a two component jet model could well reproduce the multi-wavelength afterglow lightcurves, in particular the X-ray flat segment, of short GRB 051221A. In this model, the central engine need not to be active much longer than the prompt $\gamma-$ray emission.

With the successful launch of the Swift Gamma-Ray Burst Explorer, people expected the prompt optical flash like GRB990123 would be easily detected. However the fact that early optical flash have not been detected for a number of GRBs indicates the reverse shock must be suppressed. Here we explore the possibility that the optical flash may arise from the internal shock. We find that, under certain circumstance, the optical flash of GRB990123 and GRB060111b can really be explained by the internal shock. For GRB041219a, the prompt optical emission was correlated with the gamma-ray emission, we explain this feature also in the internal shock scenario, the optical emission is the low energy extension of the gamma-ray emission, and we can restrict its redshift $z\sim 0.2$. As for GRB050904, we have shown in previous paper that the optical flash was produced by synchrotron radiation and the X-ray flare was produced by the synchrotron-self-Compton mechanism. Therefore we conclude that the early optical flash of GRBs can usually come from the internal shock. Meanwhile since the condition to produce the optical flash is not easily satisfied, so the optical flash like GRB990123 should not be common in GRBs. In addition, we also discussed the synchrotron-self-Compton effect in the internal shock model, and find that for different values of parameters, there would be soft gamma-ray (100 KeV), hard gamma-ray (10 MeV) and GeV flare accompanying the optical flash. For GRB like GRB990123, a GeV flare with fluence about 10^{-8} erg cm^{-2} s^{-1} is expected, which may be detected by the GLAST satellite.

Dust in the host galaxies of gamma-ray bursts (GRBs) dims and reddens their afterglow spectra. Knowledge of the nature of this dust is crucial for correcting for extinction, providing clues to the nature of GRB progenitors, and probing the interstellar medium of high-redshift galaxies as well as the nature of cosmic dust when the universe was much younger and galaxies were much less evolved. The dust and extinction properties of GRB host galaxies are still poorly known. Unlike previous work, we derive in this Letter the extinction curves for 10 GRB host galaxies without a priori assumption of any specific extinction types (such as that of the Milky Way, or the Small/Large Magellanic Clouds). It is found that there appears to exist two different types of extinction curves: one is relatively flat and gray, the other displays a steeper dependence on inverse wavelength, closely resembling that of the Milky Way but with the 2175$\Angstrom$ feature removed.

If X-ray flashes (XRFs) and X-ray rich Gamma-ray Bursts(XRRGs) have the same origin with Gamma-ray Bursts (GRBs) but are viewed from larger angles of structured jets, their early afterglows may differ from those of GRBs. When the ultra-relativistic outflow interact with the surrounding medium, there are two shocks formed, one is a forward shock, the other is a reverse shock. In this paper we calculate numerically the early afterglow powered by uniform jet, Gaussian jet and power-law jet in the forward-reverse shock scenario. A set of differential equations are used to govern the dynamical evolution and synchrotron self-Compton effect has been taken into account to calculate the emission. In uniform jets, the very early afterglows of XRRGs and XRFs are significantly lower than GRBs and the observed peak times of RS emission are longer in interstellar medium environment. The RS components in XRRGs and XRFs are difficult to be detected. But in stellar wind, the reduce of very early flux and the delay of RS peak time are not so remarkable. In nonuniform jet(Gaussian jet and power-law jet), where there are emission materials on the line of sight, the very early light curve resembles isotropic-equivalent ejecta in general although the RS flux decay index shows notable deviation if the RS is relativistic(in stellar wind).

GRB050904 is very interesting since it is by far the most distant GRB event known to date($z=6.29$). It was reported that during the prompt high energy emission phase, a very bright optical flare was detected, and it was temporal coincident with an X-ray flare. Here we use two models to explain the optical flare, One is the "late internal shock model", in which the optical flare is produced by the synchrotron radiation of the electrons accelerated by the late internal shock, and the X-ray flare is produced by the synchrotron-self-Compton mechanism. The other is the external forward-reverse shock model, in which the optical flare is from the reverse shock emission and the X-ray flare is attributed to the central engine activity. We show that with proper parameters, a bright optical flare can appear in both models. We think the "late internal shock model" is more favored since in this model the optical flash and the X-ray flare have the same origin, which provides a natural explanation of the temporal coincidence of them. In the forward-reverse shock scenario, fits to the optical flare and the late afterglow suggests that the physical parameters of the reverse shock are much different from that of forward shock, as found in modeling the optical flash of GRB 990123 previously.

Recently, X-ray emission lines have been observed in X-ray afterglows of several $\gamma$-ray bursts. It is a major breakthrough for understanding the nature of the progenitors. It is proposed that the X-ray emission lines can be well explained by the Geometry-Dominated models, but in these models the illuminating angle is much larger than that of the collimated jet of the $\gamma$-ray bursts(GRBs). For GRB 011211, we obtain the illuminating angle is about $\theta\sim45^{\circ}$, while the angle of GRB jet is only $3.6^{\circ}$, so we propose that the outflow of the GRBs with emission lines should have two distinct components. The wide component illuminates the reprocessing material, and produces the emission lines, while the narrow one produces the $\gamma$-ray bursts. The observations show that the energy for producing the emission lines is higher than that of the GRBs. In this case, when the wide component dominates the afterglows, a bump will appear in the GRBs afterglows. For GRB 011211, the emergence time of the bump is less than 0.05 days after the GRB, it is obviously too early for the observation to catch it. With the presence of the X-ray emission lines there should also be a bright emission component between the UV and the soft X-rays. These features can be tested by the $Swift$ satellite in the near future.

We explore two possible models which might give rise to bright X-ray flares in GRBs afterglows. One is an external forward-reverse shock model, in which the shock parameters of forward/reverse shocks are taken to be quite different. The other is a so called "late internal shock model", which requires a refreshed unsteady relativistic outflow generated after the prompt $\gamma-$ray emission. In the forward-reverse shock model, after the time $t_\times$ at which the RS crosses the ejecta, the flux declines more slowly than $(t_\oplus/t_\times)^{-(2+\beta)}$, where $t_\oplus$ denotes the observer's time and $\beta$ is the spectral index of the X-ray emission. In the ``late internal shock model", decaying slopes much steeper than $(t_\oplus/t_{\rm e, \oplus})^{-(2+\beta)}$ are possible if the central engine shuts down after $t_{\rm e, \oplus}$ and the observed variability timescale of the X-ray flare is much shorter than $t_{\rm e, \oplus}$. The sharp decline of the X-ray flares detected in GRB 011121, XRF 050406, GRB 050502b, and GRB 050730 rules out the external forward-reverse shock model directly and favors the "late internal shock model". These X-ray flares could thus hint that the central engine operates again and a new unsteady relativistic outflow is generated just a few minutes after the intrinsic hard burst.

Very early optical and near infrared (IR) emission was discovered accompanying the long gamma-ray burst (GRB) 041219a. We show that the optical/IR flash tracking the gamma-ray lightcurve during the prompt emission could be understood as emission from neutron-rich internal shocks, as has been suggested by Fan & Wei. The early K_s-band afterglow lightcurve after the prompt phase could be well-modeled as the superposition of a reverse shock and a forward shock component. The radio data also support the reverse shock interpretation.

For gamma-ray bursts (GRBs) born in a stellar wind, as the reverse shock crosses the ejecta, usually the shocked regions are still precipitated by the prompt MeV \gamma-ray emission. Because of the tight overlapping of the MeV photon flow with the shocked regions, the optical depth for the GeV photons produced in the shocks is very large. These high energy photons are absorbed by the MeV photon flow and generate relativistic e^\pm pairs. These pairs re-scatter the soft X-ray photons from the forward shock as well as the prompt \gamma-ray photons and power detectable high energy emission, significant part of which is in the sub-GeV energy range. Since the total energy contained in the forward shock region and the reverse shock region are comparable, the predicted sub-GeV emission is independent on whether the GRB ejecta are magnetized (in which case the reverse shock IC and synchrotron self-Compton emission is suppressed). As a result, a sub-GeV flash is a generic signature for the GRB wind model, and it should be typically detectable by the future {\em Gamma-Ray Large Area Telescope} (GLAST). Overlapping also influence neutrino emission. Besides the 10^{15} \sim 10^{17} eV neutrino emission powered by the interaction of the shock accelerated protons with the synchrotron photons in both the forward and reverse shock regions, there comes another $10^{14}$eV neutrino emission component powered by protons interacting with the MeV photon flow. This last component has a similar spectrum to the one generated in the internal shock phase, but the typical energy is slightly lower.

In this paper, within the popular internal shock scenario of GRBs, we calculate the early optical afterglow lightcurves of a neutron-fed GRB fireball for different assumed neutron fractions in the fireball and for both ISM- and wind-interaction models. The cases for both long and short GRBs are considered. We show that as long as the neutron fraction is significant (e.g. the number of neutrons is comparable to that of protons), rich afterglow signatures would show up. For a constant density (ISM) model, a neutron-rich early afterglow is characterized by a slowly rising lightcurve followed by a sharp re-brightening bump caused by collision between the leading neutron decay trail ejecta and the trailing ion ejecta. For a massive star stellar-wind model, the neutron-rich early afterglow shows an extended plateau lasting for about 100 seconds before the lightcurve starts to decay. The plateau is mainly attributed to the emission from the unshocked neutron decay trail. When the overlapping of the initial prompt $\gamma-$rays with the shocks and the trail is important, as is common for the wind model and is also possible in the ISM model under some conditions, the IC cooling effect suppresses the very early optical afterglow significantly, making the neutron-fed signature dimmer. For short GRBs powered by compact star mergers, a neutron-decay-induced step-like re-brightening is predicted, although the amplitude is not large. All these neutron-fed signatures are likely detectable by the Ultraviolet Optical Telescope (UVOT) on board the {\em Swift} observatory if GRB fireballs are indeed baryonic and neutron-rich. Close monitoring of early afterglows from 10s to 1000s of seconds, when combined with detailed theoretical modeling, could be used to potentially diagnose the existence of the neutron component in GRB fireballs.

We discuss the high enegry afterglow emission (including high energy photons, neutrinos and cosmic rays) following the 2004 December 27 Giant Flare from SGR 1806-20. If the initial outflow is relativistic with a bulk Lorentz factor \Gamma_0\sim {\rm tens}, the high-energy tail of the synchrotron emission from electrons in the forward shock region gives rise to a prominent sub-GeV emission, if the electron spectrum is hard enough and if the intial Lorentz factor is high enough. This signal could serve as a diagnosis of the initial Lorentz factor of the giant flare outflow. This component is potentially detectable by GLAST if a similar giant flare occurs in the GLAST era. With the available 10 MeV data, we constrain that \Gamma_0 < 50 if the electron distribution is a single power law. For a broken power law distribution of electrons, a higher \Gamma_0 is allowed. At energies higher than 1 GeV, the flux is lower because of a high energy cut off of the synchrotron emission component. The synchrotron self-Compton emission component and the inverse Compton scattering component off the photons in the giant flare oscillation tail are also considered, but they are found not significant given a moderate \Gamma_0 (e.g. \leq 10). The forward shock also accelerates cosmic rays to the maximum energy 10^{17}eV, and generate neutrinos with a typical energy 10^{14}eV through photomeson interaction with the X-ray tail photons. However, they are too weak to be detectable.

Recently iron emission features have been observed in several X-ray afterglows of GRBs. It is found that the energy obtained from the illuminating continuum which produces the emission lines is much higher than that of the main burst.The observation of SN-GRB association indicates a fallback disk should be formed after the supernovae explosion. The disk is optically thick and advection-dominated and dense. We suggest that the delayed injection energy after the initial main burst, much higher than energy of the main burst, causes the re-burst appearance in GRB afterglow and illuminates the region of the disk surface with $\tau\approx1$($\tau$ is the optical depth for the Thomson scatter ) and produces the iron emission line whose luminosity can be up to $10^{45}$ erg$s^{-1}$. The duration of the iron line emission can be $10^{4}-10^{5}$ s. This model can explain the appearance of re-burst and emission lines in GRB afterglow and disappearance of the iron emission lines, and also can naturally solve the problem of higher energy of the illuminating continuum than that of the main burst. This scenario is different from the models put forward to explain the emission lines before, that can be tested by SWIFT satellite.

We show that if the GRB ejecta itself is magnetized or neutron-rich, the very early afterglow of GRBs is very different from that powered by a pure ion fireball. In the Swift era, with the well monitored early afterglow data, we can potentially diagnose the ejecta composition and reveal the nature of the central GRB engine.

In the neutron-rich internal shocks model for Gamma-ray Burts (GRBs), the Lorentz factors (LFs) of ions shells are variable, so are the LFs of accompanying neutron shells. For slow neutron shells with a typical LF tens, the typical beta-decay radius reads R_{\beta,s} several 10^{14} cm, which is much larger than the typical internal shocks radius 10^{13} cm, so their impact on the internal shocks may be unimportant. However, as GRBs last long enough (T_{90}>20(1+z) s), one earlier but slower ejected neutron shell will be swept successively by later ejected ion shells in the range 10^{13}-10^{15} cm, where slow neutrons have decayed significantly. We show in this work that ion shells interacting with the beta-decay products of slow neutron shells can power a ultraviolet (UV) flash bright to 12th magnitude during the prompt gamma-ray emission phase or slightly delayed, which can be detected by the upcoming Satellite SWIFT in the near future.

In the Poynting Flux dominated outflow (the initial ratio of the electromagnetic energy flux to the particle energy flux $\sigma_0\gg1$) model for Gamma-ray bursts, nearly half of the internally dissipated magnetic energy is converted into the prompt $\gamma-$ray energy emission and the rest is converted into the kinetic energy of the outflow. Consequently, at the end of the $\gamma-$ray burst, $\sigma$ decreases significantly ($\sigma\sim 1$ or even smaller). We numerically investigate the very early reverse shock emission powered by such mildly magnetized outflows interacting with medium--uniform interstellar medium (ISM) or stellar wind (WIND). We show that for $\sigma\sim0.05-1$ and typical parameters of Gamma-ray bursts, both the ISM-ejecta interaction and the WIND-ejecta interaction can power very strong optical emission ($m_{\rm R}\sim 10-12{\rm th}$ magnitude or even brighter). Similar to the very early afterglow powered by the non-magnetized ejecta interacting with the external medium, the main difference between the ISM-ejecta interaction case and the WIND-ejecta interaction case is that, before the reverse shock crosses the ejecta, the R-band emission flux increases rapidly for the former, but for the latter it increases only slightly. (The abstract has been shortened). We suggest that the linear polarization detection of the early multi-wavelength afterglow is highly needed to see whether the outflows powering GRBs are magnetized or not.

We discuss the possible simultaneously UV/optical emission accompanying Gamma-ray bursts (GRBs). We show that as long as the intrinsic spectrum of GRB can extend to $\sim$10 GeV or higher, there is a large amount of relativistic $e^\pm$ pairs generated due to the annihilation of the soft $\gamma-$rays with the very energetic photons, which dominates over the electrons/positrons associated with the fireball, no matter the fireball is highly magnetized or not (For the highly magnetized fireball, the magnetic field is ordered, the high linear polarization of the multi-wavelength emission is expected). We find that these $e^\pm$ pairs can power an UV flash with $m\simeq 12-13{\rm th}$ magnitude, and the corresponding optical emission can be up to $m_{\rm R}\simeq15-16{\rm th}$ magnitude. Such bright UV emission can be detected by the upcoming satellite Swift, planned for launch in early 2004. The behavior of the optical-UV spectrum ($F_{\nu}\propto \nu^{5/2}$) differs significantly from that of the reverse shock emission ($F_{\nu}\propto \nu^{-\beta/2}$, $\beta \simeq 2.2$), which is a signature of the emission accompanying with GRB. The mild optical emission can be detected with the ROTSE-IIIa telescope system, if the response to the GRB alert is fast enough.

(Please note that the abstract has been significantly shorted) We investigate Gamma-ray Burst (GRB) internal shocks with moderate magnetization, with the magnetization parameter $\sigma$ ranging from 0.001 to 10. Possible magnetic dissipation in the stripped magnetized shells is also taken into account through introducing a parameter $k$ ($0<k\leq 1$), which is the ratio of the electric field strength of the downstream and the upstream. By solving the general MHD jump conditions, we show that the dynamic evolution of the shock with magnetic dissipation is different from the familiar one obtained in the ideal MHD limit. In view of the possible high degree of linear polarization of GRB 021206 and the identification of a possible highly magnetized flow in GRB 990123 and GRB 021211, we suggest that a mildly magnetized internal shock model ($0.01<\sigma<1$) with moderate magnetic dissipation is a good candidate to explain the GRB prompt emission data.

In the past two years, tremendous progresses to understand X-ray Flashes have been made. Now it is widely believed that X-ray Flashes and Gamma-ray Bursts are intrinsically the same, their much different peak energy and flux may be just due to our different viewing angles to them. Here we analytically calculate the very early afterglow of X-ray Flashes, i.e., the reverse shock emission powered by the outflows interacting with the interstellar medium. Assuming $z\sim 0.3$, we have shown that typically the R Band flux of reverse shock emission can be bright to $16\sim 17$th magnitude (the actual values are model dependent and sensitive to the initial Lorentz factor of the viewed ejecta). That emission is bright enough to be detected by the telescope on work today such as ROTSE-III or the upcoming UVOT carried by Swift Satellite, planned for launching in later 2004.

Extremely powerful emission lines are observed in the X-ray afterglow of several GRBs. The energy contained in the illuminating continuum which is responsible for the line production exceeds 10$^{51}$ erg, much higher than that of the collimated GRBs. It constrains the models which explain the production of X-ray emission lines. In this paper, We argue that this energy can come from a continuous postburst outflow. Focusing on a central engine of highly magnetized millisecond pulsar or magnetar we find that afterglow can be affected by the illuminating continuum, and therefore a distinct achromatic bump may be observed in the early afterglow lightcurves. With the luminosity of the continuous outflow which produces the line emission, we define the upper limit of the time when the bump feature appears. We argue that the reason why the achromatic bumps have not been detected so far is that the bumps should appear at the time too early to be observed.

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.

It is widely believed that outflows of gamma-ray bursts are jetted. Some people also suggest that the jets may have structures like $\epsilon(\theta)\propto\theta^{-k}$. We test the possibility of X-ray flashes coming from such jets in this paper. Both qualitative and quantitative analyses have shown that this model can reproduce most observational features for both X-ray flashes and gamma-ray bursts. Using common parameters of gamma-ray bursts, we have done both uniform and nonuniform jets' numerical calculations for their fluxes, spectra and peak energies. It seems that nonuniform jets are more appropriated to these observational properties than uniform ones. We also give a spectrum and flux fit to the most significant X-ray flash, XRF971019 by our model. We also have shown that in our model the observational ratio of gamma-ray bursts to X-ray flashes is about several.

The variability of gamma-ray burst (GRB) is thought to be correlated with its absolute peak luminosity, and this relation had been used to derive an estimate of the redshifts of GRBs. Recently Amati et al. presented the results of spectral and energetic properties of several GRBs with known redshifts. Here we analyse the properties of two group GRBs, one group with known redshift from afterglow observation, and another group with redshift derived from the luminosity - variability relation. We study the redshift dependence of various GRBs features in their cosmological rest frames, including the burst duration, the isotropic luminosity and radiated energy, and the peak energy $E_p$ of $\nu F_\nu$ spectra. We find that, for these two group GRBs, their properties are all redshift dependent, i.e. their intrinsic duration, luminosity, radiated energy and peak energy $E_p$, are all correlated with the redshift, which means that there are cosmological evolution effects on gamma-ray bursts features, and this can provide an interesting clue to the nature of GRBs. If this is true, then the results also imply that the redshift derived from the luminosity - variability relation may be reliable.

GRB021211 was first detected by HETE II and its early afterglow has been observed. There is a break in its afterglow light curve at about 12 minutes after the bursts, before the break the optical flux decays with a power-law index of about -1.6, while at late time the power-law slope is about -1 (Chornock et al. 2002). Here we will show that the afterglow light curve of GRB021211 can be explained within the framework of the standard fireball model. We show that the afterglow emission before the break time is the contribution of the emission from both the reverse shock and the forward shock, while the afterglow emission after the break time is mainly due to the forward shock emission. From the fitting we can give constraints on the parameters: the initial Lorentz factor $250\leq \gamma_0 \leq 900$, the surrounding medium density $n\geq 1.6\times 10^{-3}$ atoms ${\rm cm^{-3}}$, and $0.077\epsilon_{\rm B}^{3/16} \leq \epsilon_{\rm B}^{1/4}(\frac{\epsilon_{\rm e}}{0.1}) \leq 0.14$, where $\epsilon_{\rm B}$ and $\epsilon_{\rm e}$ are the fractions of the shock energy given to the magnetic field and electrons. We propose that since the values of $\epsilon_{\rm B}$ and $\epsilon_{\rm e}$ are somewhat smaller for GRB021211, so the peak energy of the reverse shock emission is well below the optical band, and thus it is substantially fainter than 990123 at similar epochs. Also we suggest that such a break might be a common feature in early optical afterglows.

Here we calculate the GRB afterglow light curves from a relativistic jet as seen by observers at a wide range of viewing angles from the jet axis, and the jet is uniform or non-uniform. We find that, for uniform jet the afterglow light curves for different viewing angles are somewhat different: in general, there are two breaks in the light curve, corresponding to the time $\gamma\sim (\theta_j-\theta_v)^{-1}$ and $\gamma\sim (\theta_j+\theta_v)^{-1}$ respectively. However, for non-uniform jet, the things become more complicated. For the case $\theta_v=0$, we can obtain the analytical results, for $k<8/(p+4)$ there should be two breaks in the light curve correspond to $\gamma\sim\theta_c^{-1}$ and $\gamma\sim\theta_j^{-1}$ respectively, while for $k>8/(p+4)$ there should be only one break corresponds to $\gamma\sim\theta_c^{-1}$, and this provides a possible explanation for some rapidly fading afterglows whose light curves have no breaks since the time at which $\gamma\sim\theta_c^{-1}$ is much earlier than our first observation time. For the case $\theta_v\neq 0$, our numerical results show that, the afterglow light curves are strongly affected by the values of $\theta_v$, $\theta_c$ and $k$. If the values of $\theta_v/\theta_c$ and $k$ are larger, there will be a prominent flattening in the afterglow light curve, which is quite different from the uniform jet, and after the flattening a very sharp break will be occurred at the time $\gamma\sim (\theta_v + \theta_c)^{-1}.

Whether gamma-ray bursts are highly beamed or not is a very important question, since it has been pointed out that the beaming will lead to a sharp break in the afterglow light curves during the ultra-relativistic phase, with the breaking point determined by $\Gamma \sim 1/\theta_{0}$, where $\Gamma$ is the bulk Lorentz factor and $\theta_{0}$ is the initial half opening angle of the ejecta, and such a break is claimed to be present in the light curves of some GRBs. In this paper we will examine whether all the observed breaks in GRB afterglow light curves can be explained by jet effects. Here we present a detailed calculation of the jet evolution and emission, and have obtained a simple formula of bulk Lorentz factor evolution. We show that the light curves are very smoothly steepened by jet effect, and the shape of the light curve is determined by only one parameter -- $(E/n)^{1/8} \theta_{0}^{3/4}$, where E and n are the fireball energy and surrounding medium density. We find that for GRB990123 and GRB991216, the jet model can approximately fit their light curves, and the values of $(E_{51}/n_1)^{1/8} \theta_{0}^{3/4}$ are about 0.17 and 0.22 respectively. While for GRB990510, GRB000301c, GRB000926 and GRB010222, their light curves cannot be fitted by the jet model, which suggests that the breaks may be caused by some other reasons, jet should be not the unique reason.

Sharp breaks have been observed in the afterglow light curves of several GRBs; this is generally explained by the jet model. However, there are still some uncertainties concerning this interpretation due to the unclear hydrodynamics of jet sideways expansion. Here we propose an alternative explanation to these observed breaks. If we assume that the multiwavelength spectra of GRB afterglows are not made of exact power law segments but their slope changes smoothly, i.e. $d\beta/d log\nu<0$, where $\beta$ is the spectral index, we find that this fact can very nicely explain the afterglow light curves showing breaks. Therefore we suggest that some breaks in the afterglow light curves may be caused by their curved spectra. The main feature of this interpretation is that the break time is dependent on the observed frequency, while the jet model produces achromatic breaks in the light curves. In addition, it is very important to know the position of the characteristic frequency $\nu_{c}$ in the multiwavelength spectrum at the time of the break, since it is a further discriminant between our model and the jet model. We find that although the optical light curves of seven GRB afterglows can be well fitted by the model we propose, in fact only one of them (i.e. GRB000926) can be explained in this framework, since for other ones the characteristic frequency $\nu_{c}$ is either above the optical after the break or below the optical before the break.

There has been increasing evidence that at least some GRBs are emission beamed. The beamed GRB afterglow evolution has been discussed by several authors in the ultra-relativistic case. It has been shown that the dynamics of the blast wave will be significantly modified by the sideways expansion, and there may be a sharp break in the afterglow light curves under certain circumstances. However, this is true only when the fireball is still relativistic. Here we present an analytical approach to the evolution of the beamed GRB blast wave expanding in the surrounding medium (density $n\propto r^{-s}$) in the non-relativistic case, our purpose is to explore whether the sideways expansion will strongly affect the blast wave evolution as in the relativistic case. We find that the blast wave evolution is strongly dependent on the speed of the sideways expansion. If it expands with the sound speed, then the jet angle $\theta $ increases with time as $\theta \propto ln t$, which means that the sideways expansion has little effect on the afterglow light curves, the flux $F\propto t^{-\frac{3(5\alpha -1)}{5}}$ for $s=0$ and $F\propto t^{-\frac{7\alpha +1}{3}}$ for $s=2$.It is clear that the light curve of $s=2$ is not always steeper than that of $s=0$, as in the relativistic case.We also show that if the expansion speed is a constant, then the jet angle $\theta \propto t$, and the radius $r \propto t^{0}$, in this case the sideways expansion has the most significant effect on the blast wave evolution, the flux $F\propto t^{-(5\alpha -1)}$ independent of $s$, and we expect that there should be a smooth and gradual break in the light curve.

The fireball model of gamma-ray bursts predicted that when the energetic blast wave encountered the surrounding medium, there will be afterglow emission, and the subsequent afterglow observations appeared to confirm this prediction. In this simplest fireball model, the electrons have been accelerated to a power law energy distribution in a relativistic blast wave, then they give afterglow emission through synchrotron radiation. Up to now synchrotron radiation is believed to be the main mechanism of GRB emission, however, here we will show that under some circumstances, the inverse Compton scattering (ICS) may play an important role, and can change the light curves of GRB afterglows. Here we investigate the effects of ICS in the relativistic case (the surrounding medium density $\rho \propto r^{-2}$) and in the non-relativistic case (for both $\rho =constant$ and $\rho \propto r^{-2}$), we find that in the relativistic case the effect of ICS is usually important, while in the non-relativistic case, this effect is usually unimportant, unless the surrounding medium density is very high. We show that if ICS is important, then it can flatten and steepen the light curves of GRB afterglows, and this may provide the explanation for some afterglow observations.

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.

Beaming of relativistic ejecta in GRBs has been postulated by many authors in order to reduce the total GRB energy, thus it is very important to look for the observational evidence of beaming. Rhoads (1999) has pointed out that the dynamics of the blast wave, which is formed when the beamed ejecta sweeping the external medium, will be significantly modified by the sideways expansion due to the increased swept up matter. He claimed that shortly after the bulk Lorentz factor ($\Gamma $) of the blast wave drops below the inverse of the initial opening angle ($\theta_{0}$) of the beamed ejecta, there will be a sharp break in the afterglow light curves. However, some other authors have performed numerical calculations and shown that the break of the light curve is weaker and much smoother than the one analytically predicted. In this paper we reanalyse the dynamical evolution of the jet blast wave, calculate the jet emission analytically, we find that the sharp break predicted by Rhoads will actually not exist, and for most cases the afterglow light curve will almost not be affected by sideways expansion unless the beaming angle is extremely small. We demonstrate that only when $\theta_{0}<0.1$, the afterglow light curves may be steepened by sideways expansion, and in fact there cannot be two breaks as claimed before. We have also constructed a simple numerical code to verify our conclusion.

The detection of delayed X-ray, optical and radio emission, "afterglow", associated with $\gamma$-ray bursts (GRBs) is consistent with fireball models, where the emission are produced by relativistic expanding blast wave, driven by expanding fireball at cosmogical distances. The emission mechanisms of GRB afterglow have been discussed by many authors and synchrotron radiation is believed to be the main mechanism. The observations show that the optical light curves of two observed gamma-ray bursts, GRB970228 and GRB GRB970508, can be described by a simple power law, which seems to support the synchrotron radiation explanation. However, here we shall show that under some circumstances, the inverse Compton scattering (ICS) may play an important role in emission spectrum and this may influence the temporal properties of GRB afterglow. We expect that the light curves of GRB afterglow may consist of multi-components, which depends on the fireball parameters.

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

The detection of the delayed emission in X-ray, optical and radio band, i.e. the afterglow of $\gamma$-ray bursts (GRBs), suggestes that the sources of GRBs are likely to be at cosmological distances. Here we explore the interaction of a relativistic shell with a uniform interstellar medium (ISM) and obtain the exact solution of the evolution of $\gamma$-ray burst remnants, including the radiative losses. We show that in general the evolution of bulk Lorentz factor $\gamma$ satisfies $\gamma \propto t^{-\alpha_{t}}$ when $\gamma \gg 1$, here $\alpha_{t}$ is mainly in the range $9/22 \sim 3/8$, the latter corresponds to adiabatic expansion. So it is clear that adiabatic expansion is a good approximation even when radiative loss is considered. However, in fact, $\alpha_{t}$ is slightly larger than 3/8, which may have some effects on a detailed data analysis. Synchrotron-self-absorption is also calculated and it is demonstrated that the radio emission may become optically thin during the afterglow. Our solution can also apply to the nonrelativistic case ($\gamma \sim 1$), at that time the observed flux decrease more rapidly than that in the relativistic case.