Compact steep-spectrum sources (CSSs) likely represent a population of young radio-loud active galactic nuclei (AGNs) and have been identified as gamma-ray emitting sources. We present a comprehensive analysis of their gamma-ray emission observed with Fermi/LAT and establish their broadband spectral energy distributions (SEDs). We derive their jet properties by the SED fits with a two-zone leptonic model for radiations from the compact core and large-scale extended region, and explore the possible signature of a unification picture of jet radiation among subclasses of AGNs. We show that the observed gamma-rays of CSSs with significant variability are contributed by the radiation of their compact cores via the inverse Compton process of the torus photons. The derived power-law distribution index of the radiating electrons is p_1~1.5-1.8, magnetic field strength is B~0.15-0.6 G, and Doppler boosting factor is \delta~2.8-8.9. Assuming that the jet is composed of electron-positron pairs, the compact cores of CSSs are magnetized and have a high radiation efficiency, similar to that of flat spectrum radio quasars. The six CSSs on average have higher Eddington ratio and black hole mass than those non-GeV-detected CSSs, and they follow the correlation between the jet power in units of Eddington luminosity (P_jet/L_Edd) and Eddington ratio (R_Edd) with other sub-classes of AGNs, P_jet/L_Edd~R_Edd^0.52, indicating that R_Edd would be a key physical driver for the unification scheme of AGN jet radiation.

GRB 200219A is a short gamma-ray burst (GRB) with an extended emission (EE) lasting $\sim 90$s. By analyzing data observed with the {\em Swift}/BAT and {\em Fermi}/GBM, we find that a cutoff power law (CPL) model can adequately well fit the spectra of the initial short pulse with $\rm E_{p}=1387^{+232}_{-134}$ keV. More interestingly, together with the EE component and early X-ray data, it exhibits a plateau emission smoothly connected with a $\sim t^{-1}$ segment, and followed by an extremely steep decay. The short GRB composed of those three segments is unique in the {\em Swift} era, and is very difficult to explain with the standard internal/external shock model of a black hole central engine, but could be consistent with the prediction of a magnetar central engine from the merger of an NS binary. We suggest that the plateau emission followed by a $\sim t^{-1}$ decay phase is powered by the spin-down of millisecond magnetar, which loses its rotation energy via GW quadrupole radiation. Then, the abrupt drop decay is caused by magnetar collapsing into black hole before switching to the EM-dominated emission. This is the first short GRB for which the X-ray emission has such an intriguing feature powered by magnetar via GW-dominated radiation. If this is the case, one can estimate the physical parameters of a magnetar, the GW signal powered by magnetar and the merger-nova emission are also discussed.

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

A Kepler solver is an analytical method used to solve a two-body problem. In this paper, we propose a new correction method by slightly modifying the Kepler solver. The only change to the analytical solutions is that the obtainment of the eccentric anomaly relies on the true anomaly that is associated to a unit radial vector calculated by an integrator. This scheme rigorously conserves all integrals and orbital elements except the mean longitude. However, the Kepler energy, angular momentum vector and Laplace-Runge-Lenz vector for perturbed Kepler problems are slowly-varying quantities. However, their integral invariant relations give the quantities high-precision values that directly govern five slowly-varying orbital elements. These elements combined with the eccentric anomaly determine the desired numerical solutions. The newly proposed method can considerably reduce various errors for a post-Newtonian two-body problem compared with an uncorrected integrator, making it suitable for a dissipative two-body problem. Spurious secular changes of some elements or quasi-integrals in the outer solar system may be caused by short integration times of the fourth-order Runge-Kutta algorithm. However, they can be eliminated in a long integration time of 10^8 years by the proposed method, similar to Wisdom-Holman second-order symplectic integrator. The proposed method has an advantage over the symplectic algorithm in the accuracy but gives a larger slope to the phase error growth.

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

Gamma-ray burst (GRB) 150910A was detected by {\it Swift}/BAT, and then rapidly observed by {\it Swift}/XRT, {\it Swift}/UVOT, and ground-based telescopes. We report Lick Observatory spectroscopic and photometric observations of GRB~150910A, and we investigate the physical origins of both the optical and X-ray afterglows, incorporating data obtained with BAT and XRT. The light curves show that the jet emission episode lasts $\sim 360$~s with a sharp pulse from BAT to XRT (Episode I). In Episode II, the optical emission has a smooth onset bump followed by a normal decay ($\alpha_{\rm R,2} \approx -1.36$), as predicted in the standard external shock model, while the X-ray emission exhibits a plateau ($\alpha_{\rm X,1} \approx -0.36$) followed by a steep decay ($\alpha_{\rm X,2} \approx -2.12$). The light curves show obvious chromatic behavior with an excess in the X-ray flux. Our results suggest that GRB 150910A is an unusual GRB driven by a newly-born magnetar with its extremely energetic magnetic dipole (MD) wind in Episode II, which overwhelmingly dominates the observed early X-ray plateau. The radiative efficiency of the jet prompt emission is $\eta_{\gamma} \approx 11\%$. The MD wind emission was detected in both the BAT and XRT bands, making it the brightest among the current sample of MD winds seen by XRT. We infer the initial spin period ($P_0$) and the surface polar cap magnetic field strength ($B_p$) of the magnetar as $1.02 \times 10^{15}~{\rm G} \leq B_{p} \leq 1.80 \times 10^{15}~{\rm G}$ and 1~ms $\leq P_{0}v\leq 1.77$~ms, and the radiative efficiency of the wind is $\eta_w \geq 32\%$.

The physical origin of fast radio bursts (FRBs) is unknown. Young magnetars born from gamma-ray bursts (GRBs) have been suggested to be a possible central engine of FRBs. We test such a hypothesis by systematically searching for GRB-FRB spatial associations from 110 FRBs and 1440 GRBs. We find that one FRB event, FRB 171209, discovered by the Parkes telescope is spatially coincident with a historical long-duration GRB 110715A at $z=0.82$. The afterglow of GRB 110715A is consistent with being powered by a millisecond magnetar. The extragalactic dispersion measure of FRB 171209 is in excess of that contributed by the intergalactic medium, which can be interpreted as being contributed by a young supernova remnant associated with the GRB. Overall, the significance of the association is $(2.28 - 2.55) \sigma$. If the association is indeed physical, our result suggests that the magnetars associated with long GRBs can be the progenitors of at least some FRBs.

The WT scheme, a piecewise polynomial force interpolation scheme with time-step dependency, is proposed in this paper for relativistic particle-in-cell (PIC) simulations. The WT scheme removes the lowest order numerical Cherenkov instability (NCI) growth rate for arbitrary time steps allowed by the Courant condition. While NCI from higher order resonances is still present, the numerical tests show that for smaller time steps, the numerical instability grows much slower than using the optimal time step found in previous studies. The WT scheme is efficient for improving the quality and flexibility of relativistic PIC simulations.

Transition Radiation (TR) plays an important role in particle identification in high-energy physics and its characteristics provide a feasible method of energy calibration in the energy range up to 10 TeV, which is of interest for dark matter searches in cosmic rays. In a Transition Radiation Detector (TRD), the TR signal is superimposed onto the ionization energy loss signal induced by incident charged particles. In order to make the TR signal stand out from the background of ionization energy loss in a significant way, we optimized both the radiators and the detector. We have designed a new prototype of regular radiator optimized for a maximal TR photon yield, combined with the Side-On TRD which is supposed to improve the detection efficiency of TR. We started a test beam experiment with the Side-On TRD at Conseil Europ\'{e}en pour la Recherche Nucl\'{e}aire (CERN), and found that the experimental data is consistent with the simulation results.

Although the physical origin of prompt emission in gamma-ray bursts (GRBs) remains inconclusive, previous studies have considered the synchrotron radiation of relativistic electrons as a promising mechanism. These works usually adopted a invariable injection rate of electrons ($Q$) which may be discordant with that in a Poynting-flux dominated jet. In a Poynting-flux dominated jet (e.g., ICMART model, Zhang & Yan 2011), the number of magnetic reconnections occurred simultaneously may grow rapidly with time and results in an increase of $Q$ with time. This paper is dedicated to study the synchrotron radiation spectrum in this scenario. It is found that the radiation spectrum would obviously get harder if an increasing $Q$ is adopted and a Band-like radiation spectrum can be obtained if the increase of $Q$ is fast enough. The latter is related to the fact that a bump-shape rather than a power-law spectrum appears in the low-energy regime of the obtained electron spectrum. This effect can strongly harden the low-energy radiation spectrum. It indicates that an increasing $Q$ can help to alleviate the "fast-cooling problem" of synchrotron radiation for GRBs. Our studies also reveal that a Poynting-flux dominated jet with a large emission radius, a small length of the magnetic reconnection region, or a low-minimum energy of injected electron would prefer to form a Band-like radiation spectrum. We suggest that the Band spectrum found in GRBs may be the synchrotron emission of the electrons with a bump-shape distribution in its low-energy regime.

Early shallow-decaying X-ray afterglows of gamma-ray bursts (GRBs) may be attributed to the dipole radiations of newly-born magnetars. Assuming that the GRB jets powered by magnetars are quasi-universal, we find that the jet structure can be parameterized as a uniform jet with a luminosity of $\log L_{\rm j}/{\rm erg\ s^{-1}}=52.68^{+0.76}_{-0.33}$ (1$\sigma$) and an opening angle $\theta_{\rm j}=2.10_{-1.28}^{+1.90}$ (50\% confidence level), surrounding by a power-law decay component with an index of ${-4.00^{+0.27}_{-0.37}}$ (1$\sigma$). The inferred local GRB rate is $\rho=9.6$ Gpc$^{-3}$ yr$^{-1}$ by including both the typical GRBs and LL-GRBs as the same population. The typical viewing angle is $3.3^{o}$, and may be $20^{o}\sim30^{o}$ for LL-GRBs. The X-ray luminosity function of the dipole radiation wind can be empirically described by a broken power-law function with indices $\beta_1=0.78^{+0.16}_{-0.15}$ and $\beta_2>1.6$ broken at $\log L_{b, w}/{\rm erg\ s^{-1}}=48.51^{+0.53}_{-0.65}$. In case of that the wind outflow is collimated and co-axial with the GRB jet, we find that the wind structure is similar to the GRB jet, i.e., $\log L_{\rm c, w}/{\rm erg\ s^{-1}}=48.38^{+0.30}_{-0.48}$, $\theta_{\rm c, w}={2.65^{o}}_{-1.19^{o}}^{+0.1.73^{o}}$, and $k_{\rm w}=4.57^{+1.21}_{-0.75}$. The observed correlation between the prompt gamma-ray luminosity and X-ray luminosity of the wind may be resulted from the viewing angle effect in such a jet-wind system. Discussion on survey with the X-ray instruments on board the {\em Einstein\ Probe} mission in the soft X-ray band for the jet and wind emission is also presented.

Violent multi-wavelength variabilities are observed in gamma-ray-selected blazars. We present an analysis of long-term light curves for eight bright blazars to explore the co-variation pattern in the gamma-ray and radio bands. We extract their gamma-ray light curves and spectra with data observed by the Fermi/LAT since 2008. We find diverse co-variation patterns between the gamma-ray and radio (at 43 GHz) fluxes in these sources. The gamma-ray and radio fluxes of 3C 454.3 and PKS 1633+382 are correlated without any time-lag, suggesting that they are from the same radiation region. Similar correlation is also observed in 3C 273 and PKS 1222+216, but the radio flux is lag behind the gamma-ray flux approximately ~160 days and ~290 days, respectively. This likely suggests that their gamma-ray emission regions are located at the upstream of their radio cores at 43 GHz. The gamma-ray and radio fluxes of the other four blazars are not correlated, implying that the gamma-ray and radio emission may be from different regions in their jets. The gamma-ray light curves of the eight blazars can be decomposed into some long timescale variability components and fast spike flares. We propose that they may be attributed to the central engine activity and the magnetic reconnection process or turbulence in the local emission region, respectively.

A comprehensively theoretical analysis on the broadband spectral energy distributions (SEDs) of large-scale jet knots in 3C 273 is presented for revealing their X-ray radiation mechanism. We show that these SEDs cannot be explained with a single electron population model when the Doppler boosting effect is either considered or not. By adding a more energetic electron (the leptonic model) or proton (the hadronic model) population, the SEDs of all knots are well represented. In the leptonic model, the electron population that contributes the X-ray emission is more energetic than the one responsible for the radio-optical emission by almost two orders of magnitude; the derived equipartition magnetic field strengths (B_eq) are ~0.1 mG. In the hadronic model, the protons with energy of ~20 PeV are required to interpret the observed X-rays; the B_eq values are several mG, larger than that in the leptonic model. Based on the fact that no resolved substructures are observed in these knots and the fast cooling-time of the high-energy electrons is difficult to explain the observed X-ray morphologies, we argue that two distinct electron populations accelerated in these knots are unreasonable and their X-ray emission would be attributed to the proton synchrotron radiation accelerated in these knots. In case of these knots have relativistic motion towards the observer, the super-Eddington issue of the hadronic model could be avoided. Multiwavelength polarimetry and the gamma-ray observations with high resolution may be helpful to discriminate these models.

The observed internal plateau of X-ray emission in some short GRBs suggests the formation of a remnant supra-massive magnetar following a double neutron star (NS) merger. In this paper, we assume that the rotational energy is lost mainly via gravitational wave (GW) radiation instead of magnetic dipole (MD) radiation, and present further constraints on the NS nuclear equation of state (EoS) via mass quadrupole deformation and r-mode fluid oscillations of the magnetar. We present two short GRBs with measured redshifts, 101219A and 160821B, whose X-ray light curves exhibit an internal plateau. This suggests that a supra-massive NS may survive as the central engine. By considering twelve NS EoSs, within the mass quadrupole deformation scenario we find that the GM1, DD2, and DDME2 models give an $M_{\rm p}$ band falling within the 2$\sigma$ region of the proto-magnetar mass distribution for $\varepsilon=0.01$. This is consistent with the constraints from the MD radiation dominated model of rotational energy loss. However, for an r-mode fluid oscillation model with $\alpha=0.1$ the data suggest that the NS EOS is close to the Shen and APR models, which is obviously different from the MD radiation dominated and mass quadrupole deformation cases.

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

Short GRBs with extended emission (EE) that are composed initially of a short-hard spike and followed by a long-lasting EE, are thought to be classified as a subsection of short GRBs. The narrow energy band available during the {\em Swift} era combined with a lack of spectral information prevented discovery of the intrinsic properties of those events. In this paper, we performed a systematic search of short GRBs with EE by using all available {\em Fermi}/GBM data. The search identified 26 GBM-detected short GRBs with EE that are similar to GRB 060614 observed by {\em Swift}/BAT. We focus on investigating the spectral and temporal properties for both the hard spike and the EE components of all 26 GRBs, and explore differences and possible correlations between them. We find that while the peak energy ($E_{\rm p}$) of the hard spikes is a little bit harder than that of the EE, but their fluences are comparable. The harder $E_{\rm p}$ seems to correspond to a larger fluence and peak flux with a large scatter for both the hard spike and EE components. Moreover, the $E_{\rm p}$ of both the hard spikes and EE are compared to other short GRBs. Finally, we also compare the properties of GRB 170817A with those short GRBs with EE and find no significant statistical differences between them. We find that GRB 170817A has the lowest $E_{\rm p}$, likely because it was off-axis.

We propose a self-similar (SS) solution to hydrodynamic non-relativistic flow behind a spherical strong blast wave (BW) passing through a homogeneous plasma with efficient relativistic particle acceleration at the shock front. The flow is described by an ideal two-fluid model with a relativistic component so that the post-shock gas has an effective SS adiabatic index $ \gamma $ varying from $ 5/3 $ to $ 4/3 $. This solution is calculated numerically and compared with the standard Sedov solution. We find that the BW center in our solution is dominated by the relativistic component with $ \gamma =4/3 $ for the divergence of expansion there, and the relativistic component dominates the interior for a moderate acceleration efficiency at the shock front. The overall efficiency of relativistic particle acceleration can be enhanced by a factor of $ 2 $ due to the slower adiabatic energy loss rate of the relativistic component during expansion. Tendency of the dominance by the relativistic component may be common in expanding astrophysical two-fluid systems such as supernova remnants, lobes of radio galaxies.

The first detected gravitational wave GW170817 from a binary neutron star merger is associated with an important optical transient AT 2017gfo, which is a direct observation of kilonova. Recent observations suggest that the remnant compact object of the binary neutron star merger associated with GW170817/GRB 170817A may be a stable long-lived magnetized neutron star. In this situation, there would be a pulsar wind nebula (PWN) embedded inside the dynamic ejecta. The PWN emission may be absorbed by the ejecta or leak out of the system. We study the effect of the PWN emission on the observed light curves and radiation spectra. Different from previous works, the absorption and leakage of the PWN emission are all involved in our model, where the absorption of the PWN emission heats up the ejecta and alters its radiation. It is found that the characteristic emission of the embedded PWN quickly evolves. For the multiband and long-term observations of AT 2017gfo, we find that the dynamic ejecta with a PWN emission can fit the observational data very well, especially for the light curves at $t\sim 5$ days and those in the late phase. In addition, our model can naturally generate the thermal to nonthermal spectrum evolution of AT 2017gfo. Our fitting result suggests that a PWN is embedded in the AT 2017gfo.

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

The second observational campaign of gravitational waves organized by the LIGO/Virgo Collaborations has led to several breakthroughs such as the detection of gravitational wave signals from merger systems involving black holes or neutrons stars. During O2,14 gravitational wave alerts were sent to the astronomical community with sky regions covering mostly over hundreds of square degrees. Among them, 6 have been finally confirmed as real astrophysical events. Since 2013, a new set of ground-based robotic telescopes called GWAC and its pathfinder mini-GWAC have been developed to contribute to the various challenges of themulti-messenger and time domain astronomy. The GWAC system is built up in the framework of the ground-segment system of the SVOM mission that will be devoted to the study of the multi-wavelength transient sky in the next decade. During O2, only the mini-GWAC telescopenetwork was fully operational. Due to the wide field of view and fast automatic follow-up capabilities of the mini-GWAC telescopes, they were well adapted to efficiently cover the sky localization areas of the gravitational wave event candidates. In this paper, we present the mini-GWAC pipeline we have set up to respond to the GW alerts and we report our optical follow-up observations of 8 GW alerts detected during the O2 run. Our observations provided the largest coverage of the GW localization areas in a short latency made by any optical facility. We found tens of optical transient candidates in our images, but none of those could be securely associated with any confirmed black hole-black hole merger event. Based on this first experience and the near future technical improvements of our network system, we will be more competitive to detect the optical counterparts from some gravitational wave events that will be detected during the upcoming O3 run, especially those emerging from binary neutron star mergers.

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

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

The Gamma-Ray Integrated Detectors (GRID) is a space mission concept dedicated to monitoring the transient gamma-ray sky in the energy range from 10 keV to 2 MeV using scintillation detectors onboard CubeSats in low Earth orbits. The primary targets of GRID are the gamma-ray bursts (GRBs) in the local universe. The scientific goal of GRID is, in synergy with ground-based gravitational wave (GW) detectors such as LIGO and VIRGO, to accumulate a sample of GRBs associated with the merger of two compact stars and study jets and related physics of those objects. It also involves observing and studying other gamma-ray transients such as long GRBs, soft gamma-ray repeaters, terrestrial gamma-ray flashes, and solar flares. With multiple CubeSats in various orbits, GRID is unaffected by the Earth occultation and serves as a full-time and all-sky monitor. Assuming a horizon of 200 Mpc for ground-based GW detectors, we expect to see a few associated GW-GRB events per year. With about 10 CubeSats in operation, GRID is capable of localizing a faint GRB like 170817A with a 90% error radius of about 10 degrees, through triangulation and flux modulation. GRID is proposed and developed by students, with considerable contribution from undergraduate students, and will remain operated as a student project in the future. The current GRID collaboration involves more than 20 institutes and keeps growing. On August 29th, the first GRID detector onboard a CubeSat was launched into a Sun-synchronous orbit and is currently under test.

We perform a detailed analysis and simulations on the transient detection capability in the near-ultraviolet (NUV) band by focusing on some major local transient events. These events include the tidal disruption event due to a supermassive blackhole, the shock breakout of a core-collapse supernova and the flare of a late-type star. Our simulations show that a set of small wide-angle NUV cameras can allow us to detect and study numerous galactic and extra-galactic transient events. Based on the analysis and simulations, here we propose a space-based NUV sky patrol mission by updating the proposal that was originally submitted to the Chinese Space Station mission in 2011. The mission proposed here is composed of a set of eight small wide-field NUV cameras each with a diameter of 20cm. The total sky area simultaneously covered by the NUV cameras is as large as 3000$\mathrm{deg^2}$. The survey cadence ranges from 30 to 300s. The transient events are required to be detected by a dedicated on-board software in real time.

The structured jet is involved to explain the afterglows and even the prompt emission of GRB 170817A. In this paper, we stress that for a precessing jet, the jet structure in the prompt emission phase and that in the afterglow phase may be different. The jet structure in the afterglow phase can be non-uniform even if a narrow-uniform jet is presented in the prompt emission phase. We estimate the jet structure in the afterglow phase under the situation that a narrow-uniform-precessing jet is launched from the central engine of gamma-ray burst. With different precession angles, it is found that the structured jet can be roughly described as follows: a narrow uniform core with power-law wings and sharp cut-off edges, a Gaussian profile, a ring shape, or other complex profile in energy per solid angle. Correspondingly, the afterglows for our obtained structured jets are also estimated. We find that the estimates of the intrinsic kinetic energy, the electron index, and the jet opening angle based on the afterglows formed in a precessing system may be incorrect. Our obtained structured jet is likely to be revealed by future observations for a fraction of gravitational wave detected merging compact binary systems (e.g., black hole-neutron star mergers).

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 event GW170817/GRB 170817A, discovered via the successful joint observation of its gravitational wave radiation and its multi-wavelength electromagnetic counterparts, was the first definite "smoking-gun" from the merger of two neutron stars (NSs). However, the remnant of the merger remains unknown. Piro et al. recently claimed that a low-significance X-ray variability in GRB 170817A. By systematically comparing the properties of variability in the afterglow of GRB 170817A and X-ray flares in GRB afterglows, we find that this X-ray variability seems to share similar statistical correlations with X-ray flares in GRB afterglows. We further investigate several possible merger product scenarios to see whether they can produce the observed X-ray variability in GRB 170817A. The first scenario invokes a stable magnetar as the central engine producing the later X-ray variability via differential rotation or fall-back accretion onto the NS. The second scenario invokes a black hole as the central engine with a fall-back accretion process. The final scenario is a central engine with a long-lived supra-massive NS. We find that the first two scenarios have difficulty producing the later X-ray variability, which requires either an impractical NS magnetic field or an extraordinarily large stellar envelope and an extremely long accretion timescale. However, the third scenario seems to be consistent with observations, and the later X-ray variability can be produced by the magnetosphere which is expelled following the collapse of the NS with a $B_p\in(3.6, 13.5)\times10^{13}$ G.

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

The merging process of binary neutron stars is the source of the quadrupole gravitational-wave (GW) radiation. Its remnant can be either a supra-massive or a stable NS. Such evidence of magnetar signature has been supported indirectly by observed X-ray plateau of some gamma-ray bursts (GRBs) afterglow. Recently, Xue et al. (2019) discovered an X-ray transient CDF-S XT2 that is claimed to be powered by a stable magnetar from the merger of double NS. In this paper, we revisit the X-ray emission of CDF-S XT2 and find that it is more consistent with a supra-massive magnetar central engine, surviving thousands of seconds to collapse black hole. We present the comparisons of the X-ray plateau luminosity, break time, and the parameters of magnetar between CDF-S XT2 and other short GRBs with internal plateau samples. By adapting the collapse time to constrain the equation of state (EOS), we find that three EOSs (GM1, DD2, and DDME2) are consistent with the observational data. On the other hand, if the most released rotation energy of magnetar is dominated by GW radiation, we also constrain the upper limit of ellipticity of NS for given EOS, and it is in a range in $[0.89-1.8]\times 10^{-3}$. Its GW signal cannot be detected by aLIGO or even for more sensitive Einstein Telescope in the future.

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

Gamma-ray burst (GRB) central engines and jet production mechanisms are still open questions. Assuming that the shallow decay segments of canonical X-ray afterglow lightcurves of {\em Swift} GRBs are attributed to the magnetic dipole (MD) radiations of newly-born magnetars, we derive the parameters of the magnetars and explore their possible relations to jet and MD wind emission. We show that the magnetar initial spin period ($P_0$) are tightly correlated with the jet energy ($E_{\rm jet}$), which is almost proportional to the wind energy ($E_{\rm wind}$). Our least square fits yield $P_0\propto E^{-0.36\pm 0.03}_{\rm jet}$ and $E_{\rm wind}\propto E^{0.91\pm 0.07}_{\rm jet}$. These relations may imply that a magnetar with faster rotating speed can power a more energetic GRB, and energy partition between the jet and wind may be quasi-universal. Although the $P_0-E_{\rm jet}$ relation is driven by a few sub-energetic GRBs in our sample, our Monte Carlo simulation analysis shows that sample selection biases from instrumental flux limits and contaminations of the bright jet afterglows cannot make this correlation. Within this jet-wind paradigm, we propose that GRB 101225A-like ultra-long GRBs, whose prompt gamma-ray/X-ray lightcurves are featured as a long-lasting plateau with a sharp drop, may be the orphan MD wind emission being due to misalignment of their jet axis to the light of sight. Brief discussion on the orphan MD wind emission and its association with the gravitational wave radiation of newly-born magnetars is presented.

The long-lasting X-ray plateau emission in long gamma-ray bursts (GRBs) shows observational evidence for ongoing energy injection, which may be from magnetar spindown due to energy released via either magnetic dipole (MD) or gravitational wave (GW) radiation. In this paper, by systematically analyzing the {\em Swift}/XRT light curves detected before 2018 July, we find 45 light curves with a measured redshift that monotonically decay as a smooth broken power law. By assuming that the central engines of these GRBs are newly born magnetars, we measure the braking index $n$ of putative millisecond magnetars, due to MD and GW radiations. The inferred braking indices are not close to 3 or 5, but range between them with a normal distribution ($n_{\rm c}=4.02\pm 0.11$). We define a dimensionless parameter $\Re$, which is the ratio between the MD and GW components, and find that the energy released via magnetar spindown in most GRBs of our sample is dominated by GW radiation for $P_0=3$ ms and $\epsilon=0.005$ and 0.01. On the other hand, we find that $\Re$ and the braking index $n$ seem to be anticorrelated within a large systematic error at $t=0$, but depend on the values of the parameters $P_0$ and $\epsilon$. These results suggest that the contribution of GW radiation cannot be ignored, and that a larger braking index leads to GWs dominating the energy released during magnetar spindown if indeed magnetars are operating in some long GRBs.

We report first results from the CHinese Exoplanet Searching Program from Antarctica (CHESPA)---a wide-field high-resolution photometric survey for transiting exoplanets carried out using telescopes of the AST3 (Antarctic Survey Telescopes times 3) project. There are now three telescopes (AST3-I, AST3-II, and CSTAR-II) operating at Dome A---the highest point on the Antarctic Plateau---in a fully automatic and remote mode to exploit the superb observing conditions of the site, and its long and uninterrupted polar nights. The search for transiting exoplanets is one of the key projects for AST3. During the Austral winters of 2016 and 2017 we used the AST3-II telescope to survey a set of target fields near the southern ecliptic pole, falling within the continuous viewing zone of the TESS mission \citep{Ricker10}. The first data release of the 2016 data, including images, catalogs and lightcurves of 26578 bright stars ($7.5\le i \le15$) was presented in \citet{Zhang18}. The best precision, as measured by the RMS of the lightcurves at the optimum magnitude of the survey ($i=10$), is around 2\,mmag. We detect 222 objects with plausible transit signals from these data, 116 of which are plausible transiting exoplanet candidates according to their stellar properties as given by the TESS Input Catalog \citep{Stassun17}, Gaia DR2 \citep{Gaia18} and TESS-HERMES spectroscopy \citep{Sharma18}. With the first data release from TESS expected in late 2018, this candidate list will be a timely for improving the rejection of potential false positives.

Located at Dome A, the highest point of the Antarctic plateau, the Chinese Kunlun station is considered to be one of the best ground-based photometric sites because of its extremely cold, dry, and stable atmosphere(Saunders et al. 2009). A target can be monitored from there for over 40 days without diurnal interruption during a polar winter. This makes Kunlun station a perfect site to search for short-period transiting exoplanets. Since 2008, an observatory has been built at Kunlun station and three telescopes are working there. Using these telescopes, the AST3 project has been carried out over the last six years with a search for transiting exoplanets as one of its key programs (CHESPA). In the austral winters of 2016 and 2017, a set of target fields in the Southern CVZ of TESS (Ricker et al. 2009) were monitored by the AST3-II telescope. In this paper, we introduce the CHESPA and present the first data release containing photometry of 26,578 bright stars (m_i < 15). The best photometric precision at the optimum magnitude for the survey is around 2 mmag. To demonstrate the data quality, we also present a catalog of 221 variables with a brightness variation greater than 5 mmag from the 2016 data. Among these variables, 179 are newly identified periodic variables not listed in the AAVSO databasea), and 67 are listed in the Candidate Target List(Stassun et al. 2017). These variables will require careful attention to avoid false-positive signals when searching for transiting exoplanets. Dozens of new transiting exoplanet candidates will be also released in a subsequent paper(Zhang et al. 2018b).

In this paper we propose the model that the coalescence of primordial black holes (PBHs) binaries with equal mass $M \sim 10^{28}$g can emit luminous gigahertz (GHz) radio transient, which may be candidate sources for the observed fast radio bursts (FRBs), if at least one black hole holds appropriate amount of net electric charge $Q$. Using a dimensionless quantity for the charge $q = Q/\sqrt{G}M$, our analyses infer that $q\sim O(10^{-4.5})$ can explain the FRBs with released energy of order $O(10^{40}) {\rm ergs}$. With the current sample of FRBs and assuming a distribution of charge $\phi(q)$ for all PBHs, we can deduce that its form is proportional to $q^{-3.0\pm0.1}$ for $q\geq 7.2\times10^{-5}$ if PBHs are sources of the observed FRBs. Furthermore, with the proposed hypothetical scenario and by estimating the local event rate of FRBs $\sim 2.6 \times 10^3 {\rm Gpc}^{-3} {\rm yr}^{-1}$, one derives a lower bound for the fraction of PBHs (at the mass of $10^{28}$g) against that of matter $f_{\rm PBH}(10^{28}{\rm g})$ $\gtrsim 10^{-5}$ needed to explain the rate. With this inspiring estimate, we expect that future observations of FRBs can help to falsify their physical origins from the PBH binaries coalescences. In the future, the gravitational waves produced by mergers of small black holes can be detected by high frequency gravitational wave detectors. We believe that this work would be a useful addition to the current literature on multimessenger astronomy and cosmology.

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

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

Planetary systems can survive the stellar evolution, as evidenced by the atmospheric metal pollution and dusty disks of single white dwarfs. Recent observations show that 1 to 4 percent of single white dwarfs are accompanied by dusty disks, while the occurrence rate of metal pollution is about 25 to 50 percent . The dusty disks and metal pollution have been associated with accretion of remnant planetary systems around white dwarfs, yet the relation between these two phenomena is still unclear. Here we suggest an evolutionary scenario to link the two observational phenomena. By analyzing a sample of metal polluted white dwarfs, we find that the mass accretion rate onto the white dwarf generally follows a broken power law decay, which matches well with the theoretical prediction, if assuming dust accretion is primarily driven by Poynting-Robertson drag and the dust source is primarily delivered via dynamically falling asteroids perturbed by a Jovian planet. The presence of disks is mainly at the early stage (about 0.1 to 0.7 billion years) of the whole process of metal pollution, which is detectable until about 8 billion years, naturally explaining the fraction (about 2-16 percent) of metal-polluted white dwarfs having dusty disks. The success of this scenario also implies that the configuration of an asteroid belt with an outer gas giant might be common around stars of several solar masses.

Using atomistic simulations, we characterize the adsorption process of organic molecules on carbon nanoparticles, both of which have been reported to be abundant in the interstellar medium (ISM). It is found that the aromatic organics are adsorbed more readily than the aliphatic ones. This selectivity would favor the formation of polycyclic aromatic hydrocarbons (PAHs) or fullerene-like structures in the ISM due to structural similarity. It is also observed in our simulations that the molecules form a monolayer over the nanoparticle surface before stacking up in aggregates. This suggests a possible layer-by-layer formation process of onion-like nanostructures in the ISM. These findings reveal the possible role of carbon nanoparticles as selective catalysts that could provide reaction substrates for the formation of interstellar PAHs, high-fullerenes and soots from gas-phase molecules.

We report the first confirmation of a hot Jupiter discovered by the Transiting Exoplanet Survey Satellite (TESS) mission: HD 202772A b. The transit signal was detected in the data from TESS Sector 1, and was confirmed to be of planetary origin through radial-velocity measurements. HD 202772A b is orbiting a mildly evolved star with a period of 3.3 days. With an apparent magnitude of V = 8.3, the star is among the brightest known to host a hot Jupiter. Based on the 27days of TESS photometry, and radial velocity data from the CHIRON and HARPS spectrographs, the planet has a mass of 1.008+/-0.074 M_J and radius of 1.562+/-0.053 R_J , making it an inflated gas giant. HD 202772A b is a rare example of a transiting hot Jupiter around a quickly evolving star. It is also one of the most strongly irradiated hot Jupiters currently known.

4C +49.22 is a gamma-ray flat spectrum radio quasar with a bright and knotty jet. We investigate the properties of the core and large-scale knots by using their spectral energy distributions (SEDs). Analyzing its Fermi/LAT data in the past 8 years, a long-term steady gamma-ray emission component is found besides bright outbursts. For the core region, the gamma-ray emission together with the simultaneous emission in the low-energy bands at different epochs is explained with the single-zone leptonic model. The derived magnetization parameters and radiation efficiencies of the radio-core jet decrease as gamma-ray flux decays, likely indicating that a large part of the magnetic energy is converted to the kinetic energy of particles in pc-scale. For the large-scale knots, their radio-optical-X-ray SEDs can be reproduced with the leptonic model by considering the inverse Compton scattering of cosmic microwave background photons. The sum of the predicted gamma-ray fluxes of these knots is comparable to that observed with LAT at 10^{24} Hz of the steady gamma-ray component, indicating that the steady gamma-ray emission may be partially contributed by these large-scale knots. This may conceal the flux variations of the low-level gamma-ray emission from the radio-core. The derived bulk Lorentz factors of the knots decrease along the distance to the core, illustrating as deceleration of jet in large-scale. The powers of the core and knots are roughly in the same order, but the jet changes from highly magnetized at the core region into particle-dominated at the large-scale knots.

The radiation mechanism of the prompt emission of gamma-ray bursts (GRBs) remains an open question. Although their spectra are usually well fitted with the empirical Band function, which is widely believed to be fully non-thermal and interpreted as an optically thin synchrotron emission, accumulating evidence shows that a thermal component actually exists. In this paper, a multi-color blackbody (mBB) model is proposed for the time-integrated spectrum of GRB 081221 by assuming a power-law distribution of the thermal luminosities with temperature, which manifests photospheric emissions from a different radius and/or angle. The effects of the minimum temperature $kT_{\min }$, the maximum temperature $kT_{\max }$ and the power law index $m$ of the luminosity distribution of an mBB are discussed. The fitting to the time-integrated spectrum during the bright phase (from 20s to 30s since the trigger) of GRB 081221 by the mBB model yields $kT_{\min }$ = $4.4\pm 0.3$ keV, $kT_{\max }$ = 57.0 $_{-1.4}^{+1.6}$ keV, and $m=-0.46_{-0.06}^{+0.05}$. When the time bin is small enough, the time-resolved spectra of GRB 081221 are well fitted with a series of single-temperature blackbodies. Our results imply the prompt emission of GRB 081221 is dominated by the photosphere emission and its time-integrated spectrum is a superposition of pure blackbody components at different times, indicating that some empirical Band spectra may be interpreted as mBB if the temperature is widely distributed.

A newly-born magnetar is thought to be central engine of some long gamma-ray bursts (GRBs). We investigate the evolution of the electromagnetic (EM) emission from the magnetic dipole (MD) radiation wind injected by spin-down of a newly-born magnetar via both quadrupole gravitational-wave (GW) and MD radiations. We show that the EM luminosity evolves as $L_{\rm em}\propto (1+t/\tau_c)^{\alpha}$, and $\alpha$ is $-1$ and $-2$ in the GW and MD radiation dominated scenarios, respectively. Transition from the GW to MD radiation dominated epoch may show up as a smooth break with slope changing from $-1$ to $-2$. If the magnetar collapses to a black hole before $\tau_c$, the MD radiation should be shut down, then the EM light curve should be a plateau followed by a sharp drop. The expected generic light curve in this paradigm is consistent with the canonical X-ray light curve of {\em Swift} long GRBs. The X-ray emission of several long GRBs are identified and interpreted as magnetar spin-down via GW or MD, as well as constrain the physical parameters of magnetar. The combination of MD emission and GRB afterglows may make the diversity of the observed X-ray light curves. This may interpret the observed chromatic behaviors of the X-ray and optical afterglow light curves and the extremely low detection rate of a jet-like break in the X-ray afterglow light curves of long GRBs.

The spectral evolution and spectral lag behavior of 92 bright pulses from 84 gamma-ray bursts (GRBs) observed by the Fermi GBM telescope are studied. These pulses can be classified into hard-to-soft pulses (H2S, 64/92), H2S-dominated-tracking pulses (21/92), and other tracking pulses (7/92). We focus on the relationship between spectral evolution and spectral lags of H2S and H2S-dominated-tracking pulses. %in hard-to-soft pulses (H2S, 64/92) and H2S-dominating-tracking (21/92) pulses. The main trend of spectral evolution (lag behavior) is estimated with $\log E_p\propto k_E\log(t+t_0)$ (${\hat{\tau}} \propto k_{\hat{\tau}}\log E$), where $E_p$ is the peak photon energy in the radiation spectrum, $t+t_0$ is the observer time relative to the beginning of pulse $-t_0$, and ${\hat{\tau}}$ is the spectral lag of photons with energy $E$ with respect to the energy band $8$-$25$ keV. For H2S and H2S-dominated-tracking pulses, a weak correlation between $k_{{\hat{\tau}}}/W$ and $k_E$ is found, where $W$ is the pulse width. We also study the spectral lag behavior with peak time $t_{\rm p_E}$ of pulses for 30 well-shaped pulses and estimate the main trend of the spectral lag behavior with $\log t_{\rm p_E}\propto k_{t_p}\log E$. It is found that $k_{t_p}$ is correlated with $k_E$. We perform simulations under a phenomenological model of spectral evolution, and find that these correlations are reproduced. We then conclude that spectral lags are closely related to spectral evolution within the pulse. The most natural explanation of these observations is that the emission is from the electrons in the same fluid unit at an emission site moving away from the central engine, as expected in the models invoking magnetic dissipation in a moderately-high-$\sigma$ outflow.

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

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

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

The occurrence rates of bright X-ray flares in z<1 gamma-ray bursts (GRBs) with or without observed supernovae (SNe) association were compared. Our Sample I: the z<1 long GRBs (LGRBs) with SNe association (SN-GRBs) and with early Swift/X-Ray Telescope (XRT) observations, consists of 18 GRBs, among which only two GRBs have bright X-ray flares. Our Sample II: for comparison, all the z<1 LGRBs without observed SNe association and with early Swift/XRT observations, consists of 45 GRBs, among which 16 GRBs present bright X-ray flares. Thus, the study indicates a lower occurrence rate of bright X-ray flares in Sample I (11.1%) than in Sample II (35.6%). In addition, if dim X-ray fluctuations are included as flares, then 16.7% of Sample I and 55.6% of Sample II are found to have flares, again showing the discrepancy between these two samples. We examined the physical origin of these bright X-ray flares and found that most of them are probably related to the central engine reactivity. To understand the discrepancy, we propose that such a lower occurrence rate of flares in the SN-GRB sample may hint at an energy partition among the GRB, SNe, and X-ray flares under a saturated energy budget of massive star explosion.

We present our optical photometric and spectroscopical observations of GRB 140629A. A redshift of $z=2.275\pm0.043$ is measured through the metal absorption lines in our spectroscopic data. Using our photometric data and multiple observational data observed with other telescopes, we show that its optical light curve is well interpreted with the standard forward shock models in the thin shell case. Its optical-X-ray afterglow spectrum is jointly fitted with a single power-law function, yielding a photon index of $-1.90\pm 0.05$. The optical extinction and neutral hydrogen absorption of the GRB host galaxy are negligible. The fit to the light curve with the standard models shows that the ambient density is $60\pm 9$ cm$^{-3}$ and the GRB radiating efficiency is as low as $\sim 0.24\%$, likely indicating a baryonic-dominated ejecta of this GRB. This burst agrees well with the $L_{\rm p, iso}-E_p^{'}-\Gamma_0$ relation, but confidently violates those empirical relations involving geometric corrections (or jet break time). This gives rise to an issue of possible selection effect on these relations since the jet opening angle of this GRB is extremely narrow (0.04 rad).

Gamma-ray Burst (GRB) collimation has been inferred with the observations of achromatic steepening in GRB light curves, known as jet breaks. Identifying a jet break from a GRB afterglow lightcurve allows a measurement of the jet opening angle and true energetics of GRBs. In this paper, we reinvestigate this problem using a large sample of GRBs that have an optical jet break which is consistent with being achromatic in the X-ray band. Our sample includes 99 GRBs from February 1997 to March 2015 that have optical and, for Swift GRBs, X-ray lightcurves that are consistent with the jet break interpretation. Out of 99 GRBs we have studied, 55 GRBs are found to have temporal and spectral behaviors both before and after the break consistent with the theoretical predictions of the jet break models, respectively. These include 53 long/soft (Type II) and 2 short/hard (Type I) GRBs. Only 1 GRB is classified as the candidate of a jet break with energy injection. Another 41 and 3 GRBs are classified as the candidates with the lower and upper limits of the jet break time, respectively. The typical beaming correction factor $f_b^{-1} \sim 1000$ for Type II GRBs, suggesting an even higher total GRB event rate density in the universe. Both isotropic and jet-corrected energies have a wide span in their distributions. We also investigate several empirical correlations (Amati, Frail, Ghirlanda and Liang-Zhang) previously discussed in the literature. We find that in general most of these relations are less tight than before. The existence of early jet breaks and hence small opening angle jets, which were detected in the {\em Swfit era}, is most likely the source of scatter. If one limits the sample to jet breaks later than $10^4$ s, the Liang-Zhang relation remains tight and the Ghirlanda relation still exists. These relations are derived from Type II GRBs, and Type I GRBs usually deviate from them.

Gamma-ray burst (GRB) 120729A was detected by Swift/BAT and Fermi/GBM, and then rapidly observed by Swift/XRT, Swift/UVOT, and ground-based telescopes. It had a single long and smooth \gamma-ray emission pulse, which extends continuously to the X-rays. We report Lick/KAIT observations of the source, and make temporal and spectral joint fits of the multiwavelength light curves of GRB 120729A. It exhibits achromatic light-curve behavior, consistent with the predictions of the external shock model. The light curves are decomposed into four typical phases: onset bump (Phase I), normal decay (Phase II), shallow decay (Phase III), and post-jet break (Phase IV). The spectral energy distribution (SED) evolves from prompt \gamma-ray emission to the afterglow with photon index from $\Gamma_{\rm \gamma}=1.36$ to $\Gamma \approx 1.75$. There is no obvious evolution of the SED during the afterglow. The multiwavelength light curves from \gamma-ray to optical can be well modeled with an external shock by considering energy injection, and a time-dependent microphysics model with $\epsilon_B\propto t^{\alpha_B}$ for the emission at early times, $T < T_{\rm 0} + 157$~s. Therefore, we conclude that both the prompt \gamma-ray emission and afterglow of GRB 120729A have the same external shock physical origin. Our model indicates that the $\epsilon_B$ evolution can be described as a broken power-law function with $\alpha_{\rm B,1} = 0.18 \pm 0.04$ and $\alpha_{\rm B,2} = 0.84 \pm 0.04$. We also systematically investigate single-pulse GRBs in the Swift era, finding that only a small fraction of GRBs (GRBs 120729A, 051111, and 070318) are likely to originate from an external shock for both the prompt \gamma-ray emission and afterglow.

3C 207 is a lobe-dominant radio galaxy with one sided jet and the bright knots in kpc-Mpc scale were resolved in the radio, optical and X-ray bands. It was confirmed as a gamma-ray emitter with Fermi/LAT, but it is uncertain whether the gamma-ray emission region is the core or knots due to the low spatial resolution of Fermi/LAT. We present an analysis of its Fermi/LAT data in the past 9 years. Different from the radio and optical emission from the core, it is found that the gamma-ray emission is steady without detection of flux variation over 2 sigma confidence level. This likely implies that the gamma-ray emission is from its knots. We collect the radio, optical, and X-ray data of knot-A, the closest knot from the core at 1 arcsec, and compile its spectral energy distribution (SED). Although the single-zone synchrotron+SSC+IC/CMB model by assuming knot-A at rest can reproduce the SED in the radio-optical-X-ray band, the predicted gamma-ray flux is lower than the LAT observations and the derived magnetic field strength deviates the equipartition condition with 3 orders of magnitude. Assuming that knot-A is relativistically moving, its SED from radio to gamma-ray bands would be well represented with the single-zone synchrotron+SSC+IC/CMB model under the equipartition condition. These results likely suggest that the gamma-ray emission may be from knot-A via the IC/CMB process and the knot should have relativistical motion. The jet power derived from our model parameters is also roughly consistent with the kinetic power estimated with the radio data.

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

We study the first electromagnetic pulse after the gravitational wave chirp signal, focusing on the profile and duration. It is found that the light curve, especially the steep decay (SD) phase, can be very different by adopting different viewing angle $\theta_{\rm view}$ on the jet shell. For an on-axis jet with a power-law radiation spectrum, the observed flux in the SD is proportional to $t_{\rm{obs}}^{-2-\beta}$ with $\beta$ being the spectral index and $t_{\rm{obs}}$ being the observer time. Here, $t_{\rm{obs}}=0$ is set at the observed time of the jet ejected from the central engine. The SD may become steep by increasing $\theta_{\rm view}$. We also study the bolometric luminosity $L$ from a jet shell with a non-power-law radiation spectrum. For an on-axis jet, $L{\propto}t_{\rm{obs}}^{-3}$ is found in the SD. However, the SD is steeper than $L{\propto}t_{\rm{obs}}^{-3}$ for the radiation from an off-axis jet. The higher value of $\theta_{\rm view}$ is, the steeper of SD would be. Then, we suggest that the SD phase can be used to discriminate an off-axis jet from an on-axis jet. The reason for above behaviors is discussed. In addition, we find that the duration of first electromagnetic pulse is close to its peak time, especially for $\theta_{\rm{view}}\sim20^\circ$. This result is consistent with that found in GW~170817/GRB~170817A. Thus, the jet corresponding to the prompt emission of GRB~170817A should be immediately ejected after the merger. Our results also reveal that the duration of the first electromagnetic pulse can provide the information of the time to search gravitational waves.

Double neutron star (DNS) merger events are promosing candidates of short Gamma-ray Burst (sGRB) progenitors as well as high-frequecy gravitational wave (GW) emitters. On August 17, 2017, such a coinciding event was detected by both the LIGO-Virgo gravitational wave detector network as GW170817 and Gamma-Ray Monitor on board NASA's {\it Fermi} Space Telescope as GRB 170817A. Here we show that the fluence and spectral peak energy of this sGRB fall into the lower portion of the distributions of known sGRBs. Its peak isotropic luminosity is abnormally low. The estimated event rate density above this luminosity is at least $190^{+440}_{-160} {\rm Gpc^{-3} \ yr^{-1}}$, which is close to but still below the DNS merger event rate density. This event likely originates from a structured jet viewed from a large viewing angle. There are similar faint soft GRBs in the {\it Fermi} archival data, a small fraction of which might belong to this new population of nearby, low-luminosity sGRBs.

We analyze the long-term lightcurve of 3C 454.3 observed with Fermi/LAT and investigate its relation to the flux in the radio, optical, and X-ray bands. By fitting the 1-day binned GeV lightcurve with multiple Gaussian function (MGF), we propose that the typical variability timescale in the GeV band is 1--10 days. The GeV flux variation is accompanied by the spectral variation characterized as flux-tracking, i.e., "harder when brighter". The GeV flux is correlated with the optical and X-ray fluxes, and a weak correlation between gamma-ray flux and radio flux is also observed. The gamma-ray flux is not correlated with the optical linear polarization degree for the global lightcurves, but they show a correlation for the lightcurves before MJD 56000. The power density spectrum of the global lightcurve shows an obvious turnover at ~7.7 days, which may indicate a typical variability timescale of 3C 454.3 in the gamma-ray band. This is also consistent with the derived timescales by fitting the global lightcurve with MGF. The spectral evolution and an increase of the optical linear polarization degree along with the increase of the gamma-ray flux may indicate that the radiation particles are accelerated and the magnetic field is ordered by the shock processes during the outbursts. In addition, the nature of 3C 454.3 may be consistent with the self-organized criticality system, similar to Sagittarius A, and thus the outbursts are from plasmoid ejections driven by magnetic reconnection. This may further support the idea that the jet radiation regions are magnetized.

Central engine of gamma-ray bursts (GRBs) may be intermittent and launch several episodes of ejecta separated by a long quiescent interval. In this scenario, an external shock is formed due to the propagation of the first launched ejecta into the circum-burst medium and the later launched ejecta may interact with the external shock at later period. Owing to the internal dissipation, the later launched ejecta may be observed at a later time ($t_{\rm{jet}}$). In this paper, we study the relation of $t_{\rm{b}}$ and $t_{\rm{jet}}$, where $t_{\rm{b}}$ is the collision time of the later launched ejecta with the formed external shock. It is found that the relation of $t_{\rm{b}}$ and $t_{\rm{jet}}$ depends on the bulk Lorentz factor ($\Gamma_{\rm{jet}}$) of the later launched ejecta and the density ($\rho$) of the circum-burst medium. If the value of $\Gamma_{\rm{jet}}$ or $\rho$ is low, the $t_{\rm{b}}$ would be significantly larger than $t_{\rm{jet}}$. However, the $t_{\rm{b}}\sim t_{\rm{jet}}$ can be found if the value of $\Gamma_{\rm{jet}}$ or $\rho$ is significantly large. Our results can explain the large lag of the optical emission relative to the $\gamma$-ray/X-ray emission in GRBs, e.g., GRB~111209A. For GRBs with a precursor, our results suggest that the energy injection into the external shock and thus more than one external-reverse shock may appear in the main prompt emission phase. According to our model, we estimate the Lorentz factor of the second launched ejecta in GRB~160625B.

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

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

The discovery of the first electromagnetic counterpart to a gravitational wave signal has generated follow-up observations by over 50 facilities world-wide, ushering in the new era of multi-messenger astronomy. In this paper, we present follow-up observations of the gravitational wave event GW170817 and its electromagnetic counterpart SSS17a/DLT17ck (IAU label AT2017gfo) by 14 Australian telescopes and partner observatories as part of Australian-based and Australian-led research programs. We report early- to late-time multi-wavelength observations, including optical imaging and spectroscopy, mid-infrared imaging, radio imaging, and searches for fast radio bursts. Our optical spectra reveal that the transient source afterglow cooled from approximately 6400K to 2100K over a 7-day period and produced no significant optical emission lines. The spectral profiles, cooling rate, and photometric light curves are consistent with the expected outburst and subsequent processes of a binary neutron star merger. Star formation in the host galaxy probably ceased at least a Gyr ago, although there is evidence for a galaxy merger. Binary pulsars with short (100 Myr) decay times are therefore unlikely progenitors, but pulsars like PSR B1534+12 with its 2.7 Gyr coalescence time could produce such a merger. The displacement (about 2.2 kpc) of the binary star system from the centre of the main galaxy is not unusual for stars in the host galaxy or stars originating in the merging galaxy, and therefore any constraints on the kick velocity imparted to the progenitor are poor.

A hyperaccreting stellar-mass black hole has been proposed as the candidate central engine of gamma-ray bursts (GRBs). The rich observations of GRBs by \textit{Fermi} and \textit{Swift} make it possible to constrain the central engine model by comparing the model predications against data. This paper is dedicated to studying the temporal evolution of central engine parameters for both prompt emission and afterglow phases. We consider two jet launching mechanisms, i.e., $\nu\bar{\nu}$ annihilations and the Blandford-Znajek (BZ) processe, and obtain analytical solutions to these two models. We then investigate the black hole central engine parameters, such as the jet power, the dimensionless entropy $\eta$, and the central engine parameter $\mu_0=\eta (1+\sigma_0)$ (where $\sigma_0$ is the initial magnetization of the engine) at the base of the jet. The black hole may be spun up by accretion, or spun down by the BZ process, leaving imprints in GRB lightcurves. Usually, a BZ jet is more powerful and is likely responsible for the late time central engine activities. However, an initially non-spinning black hole central engine may first launch a thermal "fireball" via neutrino annihilations, and then launch a Poynting-flux-dominated jet via the BZ process. Multiple flares, giant bumps, and plateaus in GRB afterglows can be well produced as the result of late time accretion onto the black hole.

GRB 160625B is an extremely bright GRB with three distinct emission episodes. By analyzing its data observed with the GBM and LAT on board the {\em Fermi} mission, we find that a multi-color black body (mBB) model can be used to fit the spectra of initial short episode (Episode I) very well within the hypothesis of photosphere emission of a fireball model. The time-resolved spectra of its main episode (Episode II), which was detected with both GBM and LAT after a long quiet stage ($\sim 180$ seconds) of the initial episode, can be fitted with a model composing of an mBB component plus a cutoff power-law (CPL) component. This GRB was detected again in the GBM and LAT bands with a long extended emission (Episode III) after a quiet period of $\sim 300$ seconds. The spectrum of Episode III is adequately fitted with a CPL plus a single power-law models, and no mBB component is required. These features may imply that the emission of three episodes are dominated by distinct physics process, i.e., Episode I is possible from cocoon emission surrounding the relativistic jet, Episode II may be from photosphere emission and internal shock of relativistic jet, and Episode III is contributed by internal and external shocks of relativistic jet. On the other hand, both X-ray and optical afterglows are consistent with standard external shocks model.

The X-ray plateau followed by a steep decay ("internal plateau") has been observed in both long and short gamma-ray burst (GRBs), implying a millisecond magnetar operating in some GRBs. The sharp decay at the end of plateau, marking the abrupt cessation of the magnetar central engine, has been considered as the collapse of a supra-massive magnetar to a black hole (BH) when it spins down. If "internal plateau" is indeed the evidence of a magnetar central engine, a natural expectation is a signature from the new-born BH in some candidates. In this work, we find that GRB 070110 is a particular case, which shows a small X-ray bump following its "internal plateau". We interpret the plateau with a spin-down supra-massive magnetar and the X-ray bump with a fall-back BH accretion. This indicates that the new-born BH is likely active in some GRBs. Therefore, GRB 070110-like events may provide a further support to the magnetar central engine model and enable us to investigate the properties of the magnetar as well as the new-born BH.

We present a comprehensive analysis of multiple wavelength observational data of the first GeV-selected narrow-line Seyfert 1 galaxy PMN J0948+0022. We derive its lightcurves in the \gamma-ray and X-ray bands from the data observed with Fermi/LAT and Swift/XRT, and make the optical and radio lightcurves by collecting the data from the literature. These lightcurves show significant flux variations. With the LAT data we show that this source is analogue to typical flat spectrum radio quasars in the $L_\gamma-\Gamma_\gamma$ plane, where $L_\gamma$ and $\Gamma_\gamma$ are the luminosity and spectral index in the LAT energy band. The gamma-ray flux is correlated with the V-band flux with a lag of ~44 days, and a moderate quasi-periodic oscillation (QPO) with a periodicity of ~490 days observed in the LAT lightcurve. A similar QPO signature is also found in the V-band lightcurve. The gamma-ray flux is not correlated with the radio flux in 15 GHz, and no similar QPO signature is found in a confidence level of 95%. Possible mechanisms of the QPO are discussed. We propose that gravitational wave observations in the future may clarify the current plausible models for the QPO.

The jet compositions, central engines, and progenitors of gamma-ray bursts (GRBs) remain open questions in GRB physics. Applying broadband observations, including GRB prompt emission and afterglow properties derived from {\em Fermi} and {\em Swift} data, as well as from Keck host-galaxy observations, we address these questions for the peculiar, bright GRB 110731A. By using the pair-opacity method, we derive $\Gamma_{0}>190$ during the prompt emission phase. Alternatively, we derive $\Gamma_{0} \approx 580$ and $\Gamma_{0} \approx 154$ by invoking the early-afterglow phase within the homogeneous density and wind cases, respectively. On the other hand, nondetection of a thermal component in the spectra suggests that the prompt emission is likely powered by dissipation of a Poynting-flux-dominated jet leading to synchrotron radiation in an optically thin region. The nondetection of a jet break in the X-ray and optical bands allows us to place a lower limit on the jet opening angle $\theta_{\rm j}>5.5^{\circ}$. Within a millisecond magnetar central engine scenario, we derive the period $P_{0}$ and polar magnetic field strength $B_{\rm p}$, which have extreme (but still allowed) values. The moderately short observed duration (7.3\,s) and relatively large redshift ($z=2.83$) places the burst as a "rest-frame short" GRB, so the progenitor of the burst is subject to debate. Its relatively large $f_{{\rm eff}, z}$ parameter (ratio of the 1\,s peak flux of a pseudo-GRB and the background flux) and a large physical offset from a potential host galaxy suggest that the progenitor of GRB 110731A may be a compact-star merger.

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.

We compile a sample of spectral energy distribution (SED) of 12 GeV radio galaxies (RGs), including eight FR I RGs and four FR II RGs. These SEDs can be represented with the one-zone leptonic model. No significant unification as expected in the unification model is found for the derived jet parameters between FR I RGs and BL Lacertae objects (BL Lacs) and between FR II RGs and flat spectrum radio quasars (FSRQs). However, on average FR I RGs have the larger gamma_b (break Lorentz factor of electrons) and lower B (magnetic field strength) than FR II RGs, analogous to the differences between BL Lacs and FSRQs. The derived Doppler factors (delta) of RGs are on average smaller than that of balzars, which is consistent with the unification model that RGs are the misaligned parent populations of blazars with smaller delta. On the basis of jet parameters from SED fits, we calculate their jet powers and the powers carried by each component, and compare their jet compositions and radiation efficiencies with blazars. Most of the RG jets may be dominated by particles, like BL Lacs, not FSRQs. However, the jets of RGs with higher radiation efficiencies tend to have higher jet magnetization. A strong anticorrelation between synchrotron peak frequency and jet power is observed for the GeV RGs and blazars in both the observer and co-moving frames, indicating that the "sequence" behavior among blazars, together with the GeV RGs, may be dominated by the jet power intrinsically.

We derive a simple analytical formula to describe the evolution of spectral index $\beta$ in the steep decay phase shaped by the curvature effect with assumption that the spectral parameters and Lorentz factor of jet shell is the same for different latitude. Here, the value of $\beta$ is estimated in 0.3$-$10keV energy band. For a spherical thin shell with a cutoff power law (CPL) intrinsic radiation spectrum, the spectral evolution can be read as a linear function of observer time. For the situation with Band function intrinsic radiation spectrum, the spectral evolution may be complex. If the observed break energy of radiation spectrum is larger than 10keV, the spectral evolution is the same as that shaped by jet shells with a CPL spectrum. If the observed break energy is less than 0.3keV, the value of $\beta$ would be a constant. Others, the spectral evolution can be approximated as a logarithmal function of the observer time in generally.

The curvature effect may be responsible for the steep decay phase observed in gamma-ray bursts. For testing the curvature effect with observations, the zero time point $t_0$ adopted to plot observer time and flux on a logarithmic scale should be appropriately selected. In practice, however, the true $t_0$ cannot be directly constrained from the data. Then, we move $t_0$ to a certain time in the steep decay phase, which can be easily identified. In this situation, we derive an analytical formula to describe the flux evolution of the steep decay phase. The analytical formula is read as $F_\nu\propto (1+\tilde t_{\rm obs}/{\tilde t_c})^{-\alpha}$ with $\alpha(\tilde{t}_{\rm obs})=2+{\int_{0}^{\log (1+\tilde{t}_{\rm obs}/{\tilde{t}_c})} {\beta(\tau)d[\log(1+\tau/{\tilde{t}_c})]}}/{\log (1 + {\tilde t}_{\rm obs}/{{\tilde t}_c})}$, where $F_\nu$ is the flux observed at frequency $\nu$, $\tilde t_{\rm obs}$ is the observer time by setting zero time point $t_0$ at a certain time in the steep decay phase, $\beta$ is the spectral index estimated around $\nu$, and ${\tilde t}_c$ is the decay timescale of the phase with $\tilde{t}_{\rm obs}{\geqslant}0$. We test the analytical formula with the data from numerical calculations. It is found that the analytical formula presents a well estimation about the evolution of flux shaped by the curvature effect. Our analytical formula can be used to confront the curvature effect with observations and estimate the decay timescale of the steep decay phase.

We compare Particle-in-Cell simulation results of relativistic electron-ion shear flows with different bulk Lorentz factors, and discuss their implications for spine-sheath models of blazar versus gamma-ray burst (GRB) jets. Specifically, we find that most properties of the shear boundary layer scale with the bulk Lorentz factor: the lower the Lorentz factor, the thinner the boundary layer, and the weaker the self-generated fields. Similarly, the energized electron spectrum peaks at an energy near the ion drift energy, which increases with bulk Lorentz factor, and the beaming of the accelerated electrons gets narrower with increasing Lorentz factor. This predicts a strong correlation between emitted photon energy, angular beaming and temporal variability with the bulk Lorentz factor. Observationally, we expect systematic differences between the high-energy emissions of blazars and GRB jets.

We aim to test if a blazar candidate of uncertain-type (BCU) in the third Fermi active galactic nuclei catalog (3LAC) can be potentially classified as a BL Lac object or a flat spectrum radio quasar (FSRQ) by performing a statistical analysis of its broadband spectral properties. We find that 34% of the radio-selected BCUs (583 BCUs) are BL Lac-like and 20% of them are FSRQ-like, which maybe within 90% level of confidence. Similarly, 77.3% of the X-ray selected BCUs (176 BCUs) are evaluated as BL Lac-like and 6.8% of them may be FSRQ-like sources. And 88.7% of the BL Lac-like BCUs that have synchrotron peak frequencies available are high synchrotron peaked BL Lacs in the X-ray selected BCUs. The percentages are accordingly 62% and 7.3% in the sample of 124 optical-selected BCUs. The high ratio of source numbers of the BL Lac-like to the FSRQ-like BCUs in the X-ray and optically selected BCU samples is due to the selection effect. Examining the consistency between our evaluation and spectroscopic identification case by case with a sample of 78 radio-selected BCUs, it is found that the statistical analysis and its resulting classifications agree with the results of the optical follow-up spectroscopic observations. Our observation campaign for high-|\rho_{s}| BCUs , i.e., |\rho_{s}|>0.8, selected with our method is ongoing.

We present time-resolved spectral analysis of the steep decay segments of 29 bright X-ray flares of gamma-ray bursts (GRBs) observed with the Swift/X-ray telescope, and model their lightcurves and spectral index evolution behaviors with the curvature effect model. Our results show that the observed rapid flux decay and strong spectral index evolution with time can be well fit with this model, and the derived characteristic timescales ($t_c$) are in the range of $33\sim 264$ seconds. Using an empirical relation between the peak luminosity and the Lorentz factor derived from the prompt gamma-rays, we estimate the Lorentz factors of the flares ($\Gamma_{\rm X}$). We obtain $\Gamma_{\rm X}=17\sim 87$ with a median value of $52$, which is smaller than the initial Lorentz factors of prompt gamma-ray fireballs. With the derived $t_c$ and $\Gamma_{\rm X}$, we constrain the radiating regions of 13 X-ray flares, yielding $R_{\rm X}=(0.2\sim 1.1)\times 10^{16}$ cm, which are smaller than the radii of the afterglow fireballs at the peak times of the flares. A long evolution feature from prompt gamma-ray phase to the X-ray epoch is found by incorporating our results with a sample of GRBs whose initial Lorentz factors are available in literatures, i.e., $\Gamma\propto [t_{p}/(1+z)]^{-0.69\pm 0.06}$. These results may shed lights on the long term evolution of GRB central engines.

We present Particle-in-Cell simulation results of relativistic shear boundary layers between electron-ion and electron-positron plasmas and discuss their potential applications to astrophysics. Specifically, we find in the case of a fast electron-positron spine surrounded by a slow-moving or stationary electron-ion sheath, lepton acceleration proceeds in a highly anisotropic manner due to electromagnetic fields created at the shear interface. While the highest-energy leptons still produce a beaming pattern (as seen in the quasi-stationary frame of the sheath) of order 1/{\Gamma}, where {\Gamma} is the bulk Lorentz factor of the spine, for lower-energy particles, the beaming is much less pronounced. This is in stark contrast to the case of pure electron-ion shear layers, in which anisotropic particle acceleration leads to significantly narrower beaming patterns than 1/{\Gamma} for the highest-energy particles. In either case, shear-layer acceleration is expected to produce strongly angle-dependent lepton (and hence, emanating radiation) spectra, with a significantly harder spectrum in the forward direction than viewed from larger off-axis angles, much beyond the regular Doppler boosting effect from a co-moving isotropic lepton distribution. This may solve the problem of the need for high (and apparently arbitrarily chosen) minimum Lorentz factors of radiating electrons, often plaguing current blazar and GRB jet modeling efforts.

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

GRB 160821B is a short gamma-ray burst (GRB) at redshift $z=0.16$, with a duration less than 1 second and without detection of any "extended emission" up to more than 100 seconds in both {\em Swift}/BAT and {\em Fermi}/GBM bands. An X-ray plateau with a sharp drop 180 seconds after the BAT trigger was observed with {\em Swift}/XRT. No supernova or kilo-nova signature was detected. Assuming the central engine of this SGRB is a recently born supra-massive magnetar, we can explain the SGRB as jet radiation and its X-ray plateau as the internal energy dissipation of the pulsar wind as it spins down. We constrain its surface magnetic field as $B_{\rm p}<3.12\times 10^{16}$ G and initial spin period as $P_0< 8.5\times 10^{-3}$ seconds. Its equation of state is consistent with the GM1 model with $M_{\rm TOV} \sim 2.37 M_\odot$ and ellipticity $\epsilon<0.07$. Its gravitational wave (GW) radiation may be detectable with the future Einstein Telescope, but is much weaker than the current detectability limit of advanced-LIGO. The GW radiation of such an event would be detectable by advanced-LIGO if it occurred at a distance of 100 Mpc ($z=0.023$).

We summarize basic observational properties of gamma-ray bursts (GRBs), including prompt emission properties, afterglow properties, and classification schemes. We also briefly comment on the current physical understanding of these properties.

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

This work focuses on a sample of seven extremely late-time X-ray flares with peak time $t_{\rm p} > 10^4 {\rm s}$, among which two flares can be confirmed as the late-time activity of central engine. The main purpose is to investigate the mechanism of such late-time flares based on the internal origin assumption. In the hyper-accreting black hole (BH) scenario, we study the possibility of two well-known mechanisms as the central engine to power such X-ray flares, i.e., the neutrino-antineutrino annihilation and the Blandford-Znajek (BZ) process. Our results show that the annihilation luminosity is far below the observational data. Thus, the annihilation mechanism cannot account for such late-time flares. For the BZ process, if the role of outflows is taken into consideration, the inflow mass rate near the horizon will be quite low such that the magnetic field will probably be too weak to power the observed X-ray flares. We therefore argue that, for the late-time flares with internal origin, the central engine is unlikely to be associated with BHs. On the contrary, a fast rotating neutron star with strong bipolar magnetic fields may be responsible for such flares.

We report our observations of very bright prompt optical and reverse shock (RS) optical emission of GRB 140512A and analyze its multi-wavelength data observed with the {\em Swift} and {\em Fermi} missions. It is found that the joint optical-X-ray-gamma-ray spectrum with our first optical detection (R=13.09 mag) at $T_0+136$ seconds during the second episode of the prompt gamma-rays can be fit by a single power-law with index $-1.32\pm 0.01$. Our empirical fit to the afterglow lightcurves indicates that the observed bright optical afterglow with R=13.00 mag at the peak time is consistent with predictions of the RS and forward shock (FS) emission of external shock models. Joint optical-X-ray afterglow spectrum is well fit with an absorbed single power-law, with an index evolving with time from $-1.86\pm 0.01$ at the peak time to $-1.57\pm 0.01$ at late epoch, which could be due to the evolution of the ratio of the RS to FS emission fluxes. We fit the lightcurves with standard external models, and derive the physical properties of the outflow. It is found that the ratio $R_B\equiv\epsilon_{\rm B, r}/\epsilon_{\rm B, f}$ is 8187, indicating a high magnetization degree in the RS region. Measuring the relative radiation efficiency with $R_e\equiv\epsilon_{\rm e, r}/\epsilon_{\rm e, f}$, we have $R_e= 0.02$, implying the radiation efficiency of the RS is much lower than that in FS. We also show that the $R_B$ of GRBs 990123, 090102, and 130427A are similar to that of GRB 140512A and their apparent difference may be mainly attributed to the difference of the jet kinetic energy, initial Lorentz factor, and medium density among them.

The Chinese Small Telescope ARray (CSTAR) is the first Chinese astronomical instrument placed in Antarctica. It is a group of four identical, fully automatic $14.5\,\rm{cm}$ telescopes, with an field of view (FOV) of $20\,\rm{deg^2}$ centered on the South Celestial Pole. Placed at Antarctic Dome A, CSTAR is designed to provide high-cadence photometry for site monitoring and variable sources detection. During the 2008 observing season, CSTAR has taken high-precision photometric data for 18,145 stars around the South Celestial Pole. At $i\,=\,7.5$ and $12$, the photometric precision reaches $\sim 8$ mmag and $\sim 30$ mmag with a cadence of 20s or 30s, respectively. Using robust detection method, we have found 15 stellar flares on 13 sources, including two classified variables. We have also found a linear relation between the decay times and the total durations of the stellar flares. The details of all detected flares along with their stellar properties are presented in this work.

With 91 months of the publicly available Fermi-LAT Pass 8 data, we analyze the gamma-ray emission from the Milky Way satellites to search for potential line signals due to the annihilation of dark matter particles into double photons. The searched targets include a sample of dwarf spheroidal galaxies, the Large Magellanic Cloud (LMC) and Small Magellanic Cloud (SMC). No significant line emission has been found neither in the stacked dwarf galaxy sample nor in the direction of LMC/SMC. The corresponding upper limits on the cross section of DM annihilation into two photons are derived. Compared with results of previous gamma-ray line searches with the Pass 7 data, the current constraints on the line emission from dwarf spheroidal galaxies has been significantly improved in a wide energy range. With the rapid increase of the sample of dwarf spheroidal galaxies (candidates), we expect that the sensitivity of gamma ray line searches will be significantly improved in the near future.

The emission from black hole binaries (BHBs) and active galactic nuclei (AGNs) displays significant aperiodic variabilities. The most promising explanation for these variabilities is the propagating fluctuations in the accretion flow. It is natural to expect that the mechanism driving variabilities in BHBs and AGNs may operate in a black hole hyper-accretion disk, which is believed to power gamma-ray bursts (GRBs). We study the variabilities of jet power in GRBs based on the model of propagating fluctuations. It is found that the variabilities of jet power and the temporal profile of erratic spikes in this scenario are similar to those in observed light curves of prompt gamma-ray emission of GRBs. Our results show that the mechanism driving X-ray variabilities in BHBs and AGNs may operate in the central engine to drive the variabilities of GRBs.

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

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

We report our very early optical observations of GRB 110530A and investigate its jet properties together with its X-ray afterglow data. A peculiar broad onset bump followed by a plateau is observed in its early R band afterglow lightcurve. The optical data in the other bands and the X-ray data are well consistent with the temporal feature of the R band lightcurve. Our joint spectral fits of the optical and X-ray data show that they are in the same regime, with a photon index of $\sim 1.70$. The optical and X-ray afterglow lightcurves are well fitted with the standard external shock model by considering a delayed energy injection component. Based on our modeling results, we find that the radiative efficiency of the GRB jet is $\sim 1\%$ and the magnetization parameter of the afterglow jet is $<0.04$ with the derived extremely low $\epsilon_B$ (the fraction of shock energy to magnetic field) of $(1.64\pm 0.25)\times 10^{-6}$. These results indicate that the jet may be matter dominated. Discussion on delayed energy injection from accretion of late fall-back material of its pre-supernova star is also presented.

By modeling the broadband spectral energy distributions (SEDs) of a typical flat spectrum radio quasar (FSRQ, 3C 279) and two GeV narrow-line Seyfert 1 galaxies (NLS1s, PMN J0948+0022 and 1H 0323+342) in different flux stages with the one-zone leptonic models, we find a universal correlation between their Doppler factors ($\delta$) and peak luminosities ($L_{\rm c}$) of external Compton scattering bumps. Compiling a combined sample of FSRQs and GeV NLS1s, it is found that both FSRQs and GeV NLS1s in different stages and in different sources well follow the same $\delta$-$L_{\rm c}$ correlation. This indicates that the variations of observed luminosities may be essentially due to the Doppler boosting effect. And the universal $\delta$-$L_{\rm c}$ relation between FSRQs and GeV NLS1s in different stages may be further evidence that the particle acceleration and radiation mechanisms for the two kinds of sources are similar. In addition, by replacing $L_{\rm c}$ with the observed luminosity in the Fermi/LAT band ($L_{\rm LAT}$), this correlation holds, and it may serve as an empirical indicator of $\delta$. We estimate the $\delta$ values with $L_{\rm LAT}$ for 484 FSRQs in the Fermi/LAT Catalog and they range from 3 to 41, with a median of 16, which are statistically consistent with the values derived by other methods.

Galaxy clusters are the largest gravitationally bound objects in the universe and may be suitable targets for indirect dark matter searches. With 85 months of Fermi-LAT Pass 8 publicly available data, we analyze the gamma-ray emission in the directions of 16 nearby Galaxy Clusters with an unbinned likelihood analysis. No globally statistically-significant $\gamma-$ray line feature is identified and a tentative line signal may be present at $\sim 43$ GeV. The 95\% confidence level upper limits on the velocity-averaged cross section of dark matter particles annihilating into double $\gamma-$rays (i.e., $\langle \sigma v \rangle_{\chi\chi\rightarrow \gamma\gamma}$) are derived. Unless very optimistic boost factors of dark matter annihilation in these Galaxy Clusters have been assumed, such constraints are much weaker than the bounds set by the Galactic $\gamma-$ray data.

X-ray flares are generally supposed to be produced by the later central engine activities, and may share the similar physical origin with prompt emission of gamma-ray bursts (GRBs). In this paper, we have analyzed all significant X-ray flares from the GRBs observed by {\em Swift} from April 2005 to March 2015. The catalog contains 468 bright X-ray flares, including 200 flares with redshifts. We obtain the fitting results of X-ray flares, such as start time, peak time, duration, peak flux, fluence, peak luminosity, and mean luminosity. The peak luminosity decreases with peak time, following a power-law behavior $L_p \propto T_{peak,z}^{-1.27}$. The flare duration increases with peak time. The 0.3-10 keV isotropic energy of X-ray flares distribution is a lognormal peaked at $10^{51.2}$ erg. We also study the frequency distributions of flare parameters, including energies, durations, peak fluxes, rise times, decay times and waiting times. Power-law distributions of energies, durations, peak fluxes, and waiting times are found in GRB X-ray flares and solar flares. These distributions could be well explained by a fractal-diffusive, self-organized criticality model. Some theoretical models basing on magnetic reconnection have been proposed to explain X-ray flares. Our result shows that the relativistic jets of GRBs may be Poynting-flux dominated.

The cosmic evolution of gamma-ray burst (GRB) luminosity is essential for revealing the GRB physics and for using GRBs as cosmological probes. We investigate the luminosity evolution of long GRBs with a large sample of 258 {\em Swift}/BAT GRBs. Parameterized the peak luminosity of individual GRBs evolves as $L_{\rm p}\propto{\rm }(1+z)^{k}$, we get $k=1.49\pm0.19$ using the non-parametric $\tau$ statistics method without considering observational biases of GRB trigger and redshift measurement. By modeling these biases with the observed peak flux and characterizing the peak luminosity function of long GRBs as a smoothly broken power-law with a break that evolves as $L_{\rm b}\propto (1+z)^{k_{\rm b}}$, we obtain $k_{\rm b}=1.14^{+0.99}_{-0.47}$ through simulations based on assumption that the long GRB rate follows the star formation rate (SFR) incorporating with cosmic metallicity history. The derived $k$ and $k_b$ values are systematically smaller than that reported in previous papers. By removing the observational biases of the GRB trigger and redshift measurement based on our simulation analysis, we generate mock {\em complete} samples of 258 and 1000 GRBs to examine how these biases affects on the $\tau$ statistics method. We get $k=0.94\pm 0.14$ and $k=0.80\pm 0.09$ for the two samples, indicating that these observational biases may lead to overestimate the $k$ value. With the large uncertain of $k_b$ derived from our simulation analysis, one even cannot convincingly argue a robust evolution feature of the GRB luminosity.

Strong spectral softening has been revealed in the late X-ray afterglows of some gamma-ray bursts (GRBs). The scenario of X-ray scattering around circum-burst dusty medium has been supported by previous works due to its overall successful prediction of both the temporal and spectral evolution of some X-ray afterglows. To further investigate the observed feature of spectral softening, we now systematically search the X-ray afterglows detected by X-Ray Telescope (XRT) of Swift and collect twelve GRBs with significant late-time spectral softening. We find that dust scattering could be the dominant radiative mechanism for these X-ray afterglows regarding their temporal and spectral features. For some well observed bursts with high-quality data, their time-resolved spectra could be well produced within the scattering scenario by taking into account the X-ray absorption from circum-burst medium. We also find that during spectral softening the power-law index in the high energy end of the spectra does not vary much. The spectral softening is mainly manifested by the spectral peak energy continually moving to the soft end.

Observations of very early multi-wavelength afterglows are critical to reveal the properties of the radiating fireball and its environment as well as the central engine of gamma-ray bursts (GRBs). We report our optical observations of GRB 111228A from 95 sec to about 50 hours after the burst trigger and investigate its properties of the prompt gamma-rays and the ambient medium using our data and the data observed with {\em Swift} and {\em Fermi} missions. Our joint optical and X-ray spectral fits to the afterglow data show that the ambient medium features as low dust-to-gas ratio. Incorporating the energy injection effect, our best fit to the afterglow lightcurves with the standard afterglow model via the Markov Chain Monte Carlo (MCMC) technique shows that $\epsilon_e=(6.9\pm 0.3)\times 10^{-2}$, $\epsilon_B=(7.73\pm 0.62)\times 10^{-6}$, $E_{\rm K}=(6.32\pm 0.86)\times 10^{53}\rm erg$, $n=0.100\pm 0.014$ cm$^{-3}$. The low medium density likely implies that the afterglow jet may be in a halo or in a hot ISM. Achromatic shallow decay segment observed in the optical and X-ray bands is well explained with the long-lasting energy injection from the central engine, which would be a magnetar with a period of about 1.92 ms inferred from the data. The $E_p$ of its time-integrated prompt gamma-ray spectrum is $\sim 26$ KeV. Using the initial Lorentz factor ($\Gamma_0=476^{+225}_{-237}$) derived from our afterglow model fit, it is found that GRB 111228A satisfies the $L_{\rm iso}-E_{\rm p,z}-\Gamma_0$ relation and bridges the typical GRBs and low luminosity GRBs in this relation.

We propose a new way of generating magnetized supersonic jets using a ring laser to irradiate a flat surface target. Using 2D FLASH code simulations which include the Biermann Battery term, we demonstrate that strong toroidal fields can be generated and sustained downstream in the collimated jet outflow far from the target surface. The field strength can be controlled by varying the ring laser separation, thereby providing a versatile laboratory platform for studying the effects of magnetic field in a variety of astrophysical settings.

We select a sample of 34 gamma-ray bursts (GRBs) whose $\Gamma_0$ values are derived with the onset peaks observed in the afterglow lightcurves (except for GRB 060218 whose $\Gamma_0$ is estimated with its radio data), and investigate the correlations among $\Gamma_0$, the isotropic peak luminosity ($L_{\rm iso}$), and the peak energy ($E_{\rm p,z}$) of the $\nu f_\nu$ spectrum in the cosmological rest frame. An analysis of pair correlations among these observables well confirms the results reported by the previous papers. More interestingly, a tight correlation among $L_{\rm iso}$, $E_{\rm p,z}$, and $\Gamma_0$ is found from a multiple regression analysis, which takes the form of $L_{\rm iso} \propto E_{\rm p,z}^{1.34\pm 0.14} \Gamma_0^{1.32\pm 0.19}$ or $E_{\rm p,z} \propto L_{\rm iso}^{0.55\pm 0.06}\Gamma_0^{-0.50\pm 0.17}$. Nine other GRBs whose $\Gamma_0$ are derived via the pair production opacity constraint also follow such a correlation. Excluding GRB 060218, the $L_{\rm iso}-E_{\rm p,z}-\Gamma_0$ correlation is valid, and it even holds in the jet co-moving frame. However, GRB 060218 deviates the $L^{'}_{\rm iso}-E^{'}_{\rm p}$ relation of typical GRBs in the jet co-moving frame with $3\sigma$. We argue that the $L_{\rm iso} - E_{\rm p, z} - \Gamma_0$ correlation may be more physical than the $L_{\rm iso} - E_{\rm p,z}$ correlation, since physically the relationship between the observed $L_{\rm iso}$ and $E_{\rm p,z}$ not only depends on radiation physics, but also depends on the bulk motion of the jet. We explore the possible origins of this correlation and discuss its physical implications for understanding GRB jet composition and radiation mechanism.

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

The external forward shock (EFS) models have been the standard paradigm to interpret the broad-band afterglow data of gamma-ray bursts (GRBs). One prediction of the models is that some afterglow temporal breaks at different energy bands should be achromatic. Observations in the Swift era have revealed chromatic afterglow behaviors at least in some GRBs, casting doubts on the EFS origin of GRB afterglows. In this paper, we perform a systematic study to address the question: how bad/good are the external forward shock models? Our sample includes 85 GRBs well-monitored X-ray and optical lightcurves. Based on how well the data abide by the EFS models, we categorize them as: Gold sample: (Grade I and II) include 45/85 GRBs. They show evidence of, or are consistent with having, an achromatic break. The temporal/spectral behaviors in each afterglow segment are consistent with the predictions (closure relations) of the EFS models. Silver sample: (Grade III and IV) include 37/85 GRBs. They are also consistent with having an achromatic break, even though one or more afterglow segments do not comply with the closure relations. Bad sample: (Grade V), 3/85 shows direct evidence of chromatic behaviors, suggesting that the EFS models are inconsistent with the data. These are included in the Bad sample. We further perform statistical analyses of various observational properties ($\alpha$, $\beta$, $t_b$ and model parameters (energy injection index q, p, $\theta_j$, $\eta_\gamma$, etc) of the GRBs in the Gold Sample, and derive constraints on the magnetization parameter $\epsilon_B$ in the EFS. Overall, we conclude that the simplest EFS models can account for the multi-wavelength afterglow data of at least half of the GRBs. When more advanced modeling (e.g., long-lasting reverse shock, structured jets) is invoked, up to $>90 \%$ of the afterglows may be interpreted within the framework of the ESF models.

The variability of accretion rate, which is believed to induce the aperiodic variability of X-ray emission from disk, may affect the energy injection into the jet. In this spirit, a correlation between disk emission and jet emission can be formed even if the mean luminosity of disk emission keeps constant. In this work, these correlations are found in the situation that the luminosity of disk emission is variable and kept with a constant mean value. The obtained correlations may be shallower than that of the fundamental plane of black hole activity. In addition, the slope of correlation may increase with increasing observed frequency of jet emission. For the luminosities spacing with three days, the slope of correlation decreases with increasing black hole mass. The deviation of our found correlations from that of the fundamental plane is related to the suppression of variability in the jet emission in comparison with that in the disk emission. This mechanism may work in some sources in which shallower correlations have been reported. Moreover, it implies that luminosities used to estimate the relation of fundamental plane should cover an appropriate timescale, in which the variability of jet emission is not significantly suppressed.

The SEDs of 18 GeV FSRQs are collected and compiled from literature, in which both the jet emission and the accretion disk radiation can be observed, in order to investigate the correlations among their jet power (P_jet), accretion disk luminosity (L_disk), and luminosity of broad line region (BLR, L_BLR). On the basis of the SED fits with the jet radiation and accretion disk radiation models, we calculate P_jet and L_disk. No correlation between P_jet with either L_disk or L_BLR is found. With a sub-sample of L_BLR for 13 GeV FSRQs, it is observed that L_BLR is strongly correlated with their L_disk. We also study the BLR covering factors of the GeV FSRQs in our sample, averagely which are smaller than that of the large samples of radio-loud and radio-quiet quasars. P_jet of some GeV FSRQs is higher than L_disk, but P_jet of all the GeV FSRQs is lower than the accretion power of black hole (BH), which is estimated by \dot{M}c^2=L_disk/0.1, indicating that the total accretion power of BH is sufficient to drive the jets in these sources; however the uncorrelation between L_disk and P_jet of the GeV FSRQs may suggest that their jets are launched by the BZ process via extracting the rotational energy of BH. Using the L_BLR-L_disk relation of the GeV FSRQs, we estimate L_disk of a BL Lac sample with their L_BLR. A comparison of L_BLR and Eddington ratio (L_disk/L_Edd) among BL Lacs, very radio-loud NLS1 galaxies, and FSRQs is also presented. It is found that along with the BL Lac-NLS1-FSRQ sequence L_BLR and L_disk/L_Edd increase, which may correspond to the change of the accretion disk structure and the transformation of the dominant mechanism for jet launching. This is also consistent with the division of their parent populations, i.e., low-excitation radio galaxies and high-excitation radio galaxies.

The Chinese Small Telescope ARray (CSTAR) has observed an area around the Celestial South Pole at Dome A since 2008. About $20,000$ light curves in the i band were obtained lasting from March to July, 2008. The photometric precision achieves about 4 mmag at i = 7.5 and 20 mmag at i = 12 within a 30 s exposure time. These light curves are analyzed using Lomb--Scargle, Phase Dispersion Minimization, and Box Least Squares methods to search for periodic signals. False positives may appear as a variable signature caused by contaminating stars and the observation mode of CSTAR. Therefore the period and position of each variable candidate are checked to eliminate false positives. Eclipsing binaries are removed by visual inspection, frequency spectrum analysis and locally linear embedding technique. We identify 53 eclipsing binaries in the field of view of CSTAR, containing 24 detached binaries, 8 semi-detached binaries, 18 contact binaries, and 3 ellipsoidal variables. To derive the parameters of these binaries, we use the Eclipsing Binaries via Artificial Intelligence (EBAI) method. The primary and the secondary eclipse timing variations (ETVs) for semi-detached and contact systems are analyzed. Correlated primary and secondary ETVs confirmed by false alarm tests may indicate an unseen perturbing companion. Through ETV analysis, we identify two triple systems (CSTAR J084612.64-883342.9 and CSTAR J220502.55-895206.7). The orbital parameters of the third body in CSTAR J220502.55-895206.7 are derived using a simple dynamical model.

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

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

The observed spectral energy distributions of five GeV-selected narrow-line Seyfert 1 (NLS1) galaxies are fitted with a model including the radiation ingredients from the relativistic jet, the accretion disk, and the corona. We compare the properties of these GeV NLS1 galaxies with flat spectrum radio quasars (FSRQs), BL Lacertae objects (BL Lacs) and radio-quiet (RQ) Seyfert galaxies, and explore possible hints for jet-disk/corona connection. Our results show that the radiation physics and the jet properties of the GeV NLS1 galaxies resemble that of FSRQs. The luminosity variations of PMN J0948+0022 and 1H 0323+342 at GeV band is tightly correlated with the beaming factor ({\delta}), similar to that observed in FSRQ 3C 279. The accretion disk luminosities and the jet powers of the GeV NLS1 galaxies cover both the ranges of FSRQs and BL Lacs. With detection of bright corona emission in 1H 0323+342, we show that the ratio of the corona luminosity (L_corona) to the accretion disk luminosity (L_d) is marginally within the high end of this ratio distribution for a RQ Seyfert galaxy sample, and the variation of jet luminosity may connect with Lcorona. However, it is still unclear whether a system with a high L_corona/L_d ratio prefers to power a jet.

We carried out two-dimensional high-resolution simulations to study the effect of dust feedback on the evolution of vortices induced by massive planets in protoplanetary disks. Various initial dust to gas disk surface density ratios ($0.001$ -- $0.01$) and dust particle sizes (Stokes number $4\times10^{-4}$ -- $0.16$) are considered. We found that while dust particles migrate inwards, vortices are very effective in collecting them. When dust density becomes comparable to gas density within the vortex, a dynamical instability is excited and it alters the coherent vorticity pattern and destroys the vortex. This dust feedback effect is stronger with higher initial dust/gas density ratio and larger dust grain. Consequently, we found that the disk vortex lifetime can be reduced up to a factor of 10. We discuss the implications of our findings on the survivability of vortices in protoplanetary disks and planet formation.

X-ray flares of gamma-ray bursts (GRBs) are usually observed in the soft X-ray range and the spectral coverage is limited. In this paper, we present an analysis of 32 GRB X-ray flares that are simultaneously observed by both BAT and XRT on board the Swift mission, so a joint spectral analysis with a wider spectral coverage is possible. Our results show that the joint spectra of 19 flares are fitted with the absorbed single power-law or the Band function models. More interestingly, the joint spectra of the other 13 X-ray flares are fitted with the absorbed single power-law model plus a black body (BB) component. Phenomenally, the observed spectra of these 13 flares are analogous to several GRBs with a thermal component, but only with a much lower temperature of $kT=1\sim 3$ keV. Assuming that the thermal emission is the photosphere emission of the GRB fireball, we derive the fireball properties of the 13 flares that have redshift measurements, such as the bulk Lorentz factor $\Gamma_{\rm ph}$ of the outflow. The derived $\Gamma_{\rm ph}$ range from $50$ to $150$ and a relation of $\Gamma_{\rm ph}$ to the thermal emission luminosity is found. It is consistent with the $\Gamma_0-L_{\rm iso}$ relation that are derived for the prompt gamma-ray emission. We discuss the physical implications of these results within the content of jet composition and radiation mechanism of GRBs and X-ray flares.

We present the optical luminosity function (LF) of gamma-ray bursts (GRBs) estimated from a uniform sample of 58 GRBs from observations with the Robotic Optical Transient Search Experiment III (ROTSE-III). Our GRB sample is divided into two sub-samples: detected afterglows (18 GRBs), and those with upper limits (40 GRBs). The $R$ band fluxes 100s after the onset of the burst for these two sub-samples are derived. The optical LFs at 100s are fitted by assuming that the co-moving GRB rate traces the star-formation rate. The detection function of ROTSE-III is taken into account during the fitting of the optical LFs by using Monte Carlo simulations. We find that the cumulative distribution of optical emission at 100s is well-described with an exponential rise and power-law decay (ERPLD), broken power-law (BPL), and Schechter LFs. A single power-law (SPL) LF, on the other hand, is ruled out with high confidence.

We present a spectral analysis for a sample of redshift known GRBs observed with {\em Fermi}/GBM. Together with the results derived from our systematical spectral energy distribution modeling with the leptonic models for a {\em Fermi}/LAT blazar sample, we compare the distributions of the GRBs and the blazars by plotting the synchrotron peak luminosity ($L_{\rm s}$) and the corresponding peak photon energy $E_{\rm s}$ of blazars in the $L_{\rm p}-E_{\rm p}$ plane of GRBs, where $L_{\rm p}$ and $E_{\rm p}$ are the peak luminosity and peak photon energy of the GRB time-integrated $\nu f_\nu$ spectrum, respectively. The GRBs are in the high-$L_{\rm p}$, high-$E_{\rm p}$ corner of the plane and a tight $L_{\rm p}-E_{\rm p}$ relation is found, i.e., $L_{\rm p}\propto E_{\rm p}^{2.13^{+0.54}_{-0.46}}$. Both FSRQs and LBLs are clustered in the low-$E_{\rm p}$, low-$L_{\rm p}$ corner. IBLs and HBLs have $E_{\rm s}\sim 2\times 10^{-3} - 10^{2}$ keV and $L_{\rm s} \sim 10^{44} - 10^{47}$ erg s$^{-1}$, but no dependence of $L_{\rm s}$ on $E_{\rm s}$ is found. We show that the tight $L_p-E_p$ relation of GRBs is potentially explained with the synchrotron radiation of fast-cooling electrons in a highly magnetized ejecta, and the weak anti-correlation of $L_{\rm s}-E_{\rm s}$ for FSRQs and LBLs may be attributed to synchrotron radiation of slow-cooling electrons in a moderately magnetized ejecta. The distributions of IBLs and HBLs in the $L_{\rm p}-E_{\rm p}$ plane may be interpreted with synchrotron radiation of fast-cooling electrons in a matter-dominated ejecta. These results may present a unified picture for the radiation physics of relativistic jets in GRBs and blazars within the framework of the leptonic synchrotron radiation models.

It was suggested that the relativistic jets in gamma-ray bursts (GRBs) are powered via the Blandford-Znajek (BZ) mechanism or the annihilation of neutrinos and anti-neutrinos from a neutrino cooling-dominated accretion flow (NDAF). The advection and diffusion of the large-scale magnetic field of a NDAF is calculated, and the external magnetic field is found to be dragged inward efficiently by the accretion flow for a typical magnetic Prandtl number P_m=1. The maximal BZ jet power can be ~10^53-10^54 erg/sec for an extreme Kerr black hole, if an external magnetic field with 10^14 Gauss is advected by the NDAF. This is roughly consistent with the field strength of the disk formed after a tidal disrupted magnetar. The accretion flow near the black hole horizon is arrested by the magnetic field if the accretion rate is below than a critical value for a given external field. The arrested accretion flow fails to drag the field inward and the field strength decays, and then the accretion re-starts, which leads to oscillating accretion. The typical timescale of such episodic accretion is in an order of one second. This can qualitatively explain the observed oscillation in the soft extend emission of short-type GRBs.

We jointly analyze the gamma-ray burst (GRB) data observed with BAT and XRT on board the Swift mission to present a global view on the internal energy dissipation processes in GRBs, including precursors, prompt gamma-ray emission, extended soft gamma-ray emission, and late X-ray flares. The Bayesian block method is utilized to analyze the BAT lightcurves to identify various emission episodes. Our results suggest that these emission components likely share a same physical origin, which is repeated activation of the GRB central engine. What we observe in the gamma-ray band may be the tip-of-iceberg of more extended underlying activities. The precursor emission, which is detected in about 10% of {\em Swift} GRBs, is preferably detected in those GRBs that have a massive star core-collapse origin. The soft extended emission (EE) tail, on the other hand, is preferably detected in those GRBs that have a compact star merger origin. Bright X-ray emission is detected during the BAT quiescent phases prior to subsequent gamma-ray peaks, implying that X-ray emission may be detectable prior the BAT trigger time. Future GRB alert instruments with soft X-ray capability would be essential to reveal the early stage of GRB central engine activities, sheding light into jet composition and jet launching mechanism in GRBs.

Traditionally gamma-ray bursts (GRBs) are classified in the $T_{90}$-hardness ratio two-dimensional plane into long/soft and short/hard GRBs. In this paper, we suggest to add the "amplitude" of GRB prompt emission as the third dimension as a complementary criterion to classify GRBs, especially those of short durations. We define three new parameters $f$, $f_{\rm eff}$ and $f_{\rm eff,z}$ as ratios between the measured/simulated peak flux of a GRB/pseudo-GRB and the flux background, and discuss the applications of these parameters to GRB classification. We systematically derive these parameters to find that most short GRBs are likely not "tip-of-iceberg" of long GRBs. However, one needs to be cautious if a short GRB has a relatively small $f$ (e.g. $f<1.5$), since the chance for an intrinsically long GRB to appear as a "disguised" short GRB is higher. Based on avaialble data, we quantify the probability of a disguised short GRB below a certain $f$ value is as $P (<f)\sim 0.78^{+0.71}_{-0.4} f^{-4.33\pm 1.84}$. By progressively "moving" a long GRB to higher redshifts through simulations, we also find that most long GRBs would show up as rest-frame short GRBs above a certain redshift.

We fit the spectral energy distributions (SEDs) of a GeV-TeV FSRQ sample with the leptonic model. Their gamma_min of the relativistic electron distributions, which significantly affect the estimates of the jet properties, are constrained, with a typical value of 48. Their jet power, magnetized parameter, radiation efficiency, and jet production/radiation rates per central black hole (BH) mass are derived and compared to that of BL Lacs. We show that he FSRQ jets may be dominated by the Poynting flux and have a high radiation efficiency, whereas the BL Lac jets are likely dominated by particles and have a lower radiation efficiency than FSRQs. Being different from BL Lacs, the jet powers of FSRQs are proportional to their central BH masses. The jet production and radiation rates of the FSRQs distribute in narrow ranges and are correlated with each other, whereas no similar feature is found for the BL Lacs. We also show that the jet power is correlated with the cavity kinetic power, and the magnetic field energy in the jets may provide the cavity kinetic energy of FSRQs and the kinetic energy of cold protons in the jets may be crucial for cavity kinetic energy of BL Lacs. We suggest that the dominating formation mechanism of FSRQ jets may be the BZ process, but BL Lac jets may be produced via the BP and/or BZ processes, depending on the structures and accretion rates of accretion disks.

We present Particle-in-Cell simulation results of relativistic shear flows for hybrid positron-electron-ion plasmas and compare to those for pure e+e- and pure e-ion plasmas. Among the three types of relativistic shear flows, we find that only hybrid shear flow is able to energize the electrons to form a high-energy spectral peak plus a hard power-law tail. Such electron spectra are needed to model the observational properties of gamma-ray bursts.

Gamma-ray bursts (GRBs) and GeV-TeV selected radio loud Active Galactic Nuclei (AGNs) are compared based on our systematic modeling of the observed spectral energy distributions of a sample of AGNs with a single-zone leptonic model. We show that the correlation between the jet power (P_{jet}) and the prompt gamma-ray luminosity (L_{jet}) of GRBs is consistent, within the uncertainties, with the correlation between jet power and the synchrotron peak luminosity (L_{s, jet}) of flat spectrum radio quasars (FSRQs). Their radiation efficiencies (\varepsilon) are also comparable (>10% for most sources), which increase with the bolometric jet luminosity (L_{bol,jet}) for FSRQs and with the L_{jet} for GRBs with similar power-law indices. BL Lacs do not follow the P_{jet}-L_{s, jet} relation of FSRQs. They have lower \varepsilon and L_{bol, jet} values than FSRQs, and a tentative L_{bol, jet}-\varepsilon relation is also found, with a power-law index being different from that of the FSRQs. The magnetization parameters (\sigma) of FSRQs are averagely larger than that of BL Lacs. They are anti-correlated with $\varepsilon$ for the FSRQs, but positive correlated with \varepsilon for the BL Lacs. GeV Narrow-line Seyfert 1 galaxies potentially share similar properties with FSRQs. Based on the analogy between GRBs and FSRQs, we suggest that the prompt gamma-ray emission of GRBs is likely produced by synchrotron process in a magnetized jet with high radiation efficiency, similar to FSRQs. The jets of BL Lacs, on the other hand, are less efficient and are likely more matter dominated.

We continue our systematic statistical study on optical afterglow data of gamma-ray bursts (GRBs). We present the apparent magnitude distributions of early optical afterglows at different epochs (t= 10^2 s, t = 10^3 s, and 1 hour) for the optical lightcurves of a sample of 93 GRBs (the global sample), and for sub-samples with an afterglow onset bump or a shallow decay segment. For the onset sample and shallow decay sample we also present the brightness distribution at the peak time t_{p} and break time t_{b}, respectively. All the distributions can be fit with Gaussian functions. We further perform Monte Carlo simulations to infer the luminosity function of GRB optical emission at the rest-frame time 10^3 seconds, t_{p}, and t_{b}, respectively. Our results show that a single power-law luminosity function is adequate to model the data, with indices -1.40+/-0.10, -1.06+/- 0.16, and -1.54\+/- 0.22, respectively. Based on the derived rest-frame 10^3 s luminosity function, we generate the intrinsic distribution of the R-band apparent magnitude M_{R} at the observed time 10^{3} seconds post trigger, which peaks at M_{R}=22.5 mag. The fraction of GRBs whose R-band magnitude is fainter than 22 mag, and 25 mag and at the observer time 10^3 seconds are ~63% and ~25%, respectively. The detection probabilities of the optical afterglows with ground-based robotic telescopes and UVOT onboard {Swift} are roughly consistent with that inferred from this intrinsic M_{R} distribution, indicating that the variations of the dark GRB fraction among the samples with different telescopes may be due to the observational selection effect, although the existence of an intrinsically dark GRB population cannot be ruled out.

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

Supersonic plasma outflows driven by multi-beam, high-energy lasers, such as Omega and NIF, have been and will be used as platforms for a variety of laboratory astrophysics experiments. Here we propose a new way of launching high density and high velocity, plasma jets using multiple intense laser beams in a hollow ring formation. We show that such jets provide a more flexible and versatile platform for future laboratory astrophysics experiments. Using high resolution hydrodynamic simulations, we demonstrate that the collimated jets can achieve much higher density, temperature and velocity when multiple laser beams are focused to form a hollow ring pattern at the target, instead of focused onto a single spot. We carried out simulations with different ring radii and studied their effects on the jet properties. Implications for laboratory collisionless shock experiments are discussed.

A hyper-accreting stellar-mass black hole has been long speculated as the best candidate of central engine of gamma-ray bursts (GRBs). Recent rich observations of GRBs by space missions such as Swift and Fermi pose new constraints on GRB central engine models. In this paper, we study the baryon loading processes of a GRB jet launched from a black hole central engine. We consider a relativistic jet powered by $\nu \bar\nu$-annihilation or by the Blandford-Znajek (BZ) mechanism. We consider baryon loading from a neutrino-driven wind from a neutrino-cooling-dominated accretion flow. For a magnetically dominated BZ jet, we consider neutron-drifting from the magnetic wall surrounding the jet and subsequent positron capture and proton-neutron inelastic collisions. The minumim baryon loads in both types of jet are calculated. We find that in both cases, a more luminous jet tends to be more baryon poor. A neutrino-driven "fireball" is typically "dirtier" than a magnetically dominated jet, while a magnetically dominated jet can be much cleaner. Both models have the right scaling to interpret the empirical $\Gamma-L_{\rm iso}$ relation discovered recently. Since some neutrino-driven jets have too much baryon loading as compared with the data, we suggest that at least a good fraction of GRBs should have a magnetically dominated central engine.

Four well-sampled spectral energy distributions (SEDs) of PMN J0948+0022 are fitted with the syn+SSC+EC model to derive the physical parameters of its jets and to investigate the spectral variations of its SEDs. A tentative correlation between the peak luminosity L_c and peak frequency nu_c of its inverse Compton (IC) bump is found in both the observer and co-moving frames, indicating that the variations of luminosity are accompanied with the spectral shift. A correlation between L_c and delta is found, and thus the magnification of the external photon field by the bulk motion of the radiation regions is an essential reason for the spectral variation since the IC bump of PMN J0948+0022 is dominated by the EC process.

GRB 121027A is un-usual with its extremely long-lasting, energetic X-ray flares. The total energy release in X-ray flares is about one order of magnitude higher than prompt gamma-rays, making it special from most long GRBs. We show that while the prompt gamma-ray emission satisfies the empirical E_{iso}-E_{p} relation of typical long GRBs, the X-ray flares, whose spectra can be fit with a cutoff-power-law model with well-constrained E_p, significantly deviate from such a relation. Nonetheless, a time-resolved spectral analysis of X-ray flares suggest that the X-ray emission is consistent with the L_{iso}-E_{p} relation of long GRBs. We constrain the minimum Lorentz factor of the X-ray flares to be ~14, which is consistent with the Gamma-L_{iso} relation. Our results imply that prompt gamma-ray emission and late X-ray flares share the similar radiation mechanism, but originate from the outflows with different Lorentz factors. We search for similar GRBs from the Swift GRB archives, and find that the z=6.29 GRB 050904 is a carbon copy of GRB 121027A. The long-lasting, energetic X-ray flares in these GRBs demand significant accretion at late times, which point towards a large-radius progenitor star.

The jet properties of Mrk 501 and 3C 279 are derived by fitting the broadband spectral energy distributions with the lepton models. The derived gamma_b (the break Lorenz factor of the electron distribution) are 10^4-10^6 for Mrk 501 and 200~600 for 3C 279 and the magnetic field strength (B) of Mrk 501 is usually one order of magnitude lower than that of 3C 279, but their Doppler factors (delta) are comparable. A spectral variation feature that the peak luminosity is correlated with the peak frequency, which is opposite to the blazar sequence, is observed in the two sources. We find that (1) the peak luminosities of the two bumps in SEDs for Mrk 501 depend on gamma_b in both the observer and co-moving frames, but they are not correlated with B and delta; (2)the luminosity variation of 3C 279 is dominated by the external Compton (EC) peak and its peak luminosity is correlated with gamma_b and delta, but anti-correlated with B. These results suggest that gamma_b may govern the spectral variation of Mrk 501 and delta and B would be responsible for the spectral variation of 3C 279. The narrow distribution of gamma_b and the correlation of gamma_b and B in 3C 279 would be due to the cooling by the EC process and the strong magnetic field. Based on our brief discussion, we propose that this spectral variation feature may be originated from the instability of the corona but not from the variation of the accretion rate as that for the blazar sequence.

We present data for relativistic hot electron production by the Texas Petawatt Laser irradiating solid Au targets with thickness between 1 and 4 mm. The experiment was performed at the short focus target chamber TC1 in July 2011, with laser energies around 50 J. We measured hot electron spectra out to 50 MeV which show a narrow peak around 10 - 20 MeV plus high energy exponential tail. The hot electron spectral shape differs from those reported for other PW lasers. We did not observe direct evidence of positron production above background.

Using 2.5-dimensional Particle-in-Cell simulations, we study the kinetic physics of relativistic shear flow boundary in collisionless electron-positron (e+e-) plasmas. We find efficient magnetic field generation and particle energization at the shear boundary, driven by streaming instabilities across the shear interface and sustained by the shear flow. Nonthermal, anisotropic high-energy particles are accelerated across field lines to produce a power-law tail, truncated at energies below the shear Lorentz factor. These results have important implications for the dissipation and radiation of jets in blazars, gamma-ray bursts and other relativistic outflows.

The durations (T90) of 315 GRBs detected with Fermi/GBM (8-1000 keV) by 2011 September are calculated using the Bayesian Block method. We compare the T90 distributions between this sample and those derived from previous/current GRB missions. We show that the T90 distribution of this GRB sample is bimodal, with a statistical significance level being comparable to those derived from the BeppoSAX/GRBM sample and the Swift/BAT sample, but lower than that derived from the CGRO/BATSE sample. The short-to-long GRB number ratio is also much lower than that derived from the BATSE sample, i.e., 1:6.5 vs 1:3. We measure T90 in several bands, i.e., 8-15, 15-25, 25-50, 50-100, 100-350, and 350-1000 keV, to investigate the energy-dependence effect of the bimodal T90 distribution. It is found that the bimodal feature is well observed in the 50-100 and 100-350 keV bands, but is only marginally acceptable in the 25-50 keV and 350-1000 keV bands. The hypothesis of the bimodality is confidently rejected in the 8-15 and 15-25 keV bands. The T90 distributions in these bands are roughly consistent with those observed by missions with similar energy bands. The parameter T90 as a function of energy follows \bar T90 \propto E^{-0.20\pm 0.02} for long GRBs. Considering the erratic X-ray and optical flares, the duration of a burst would be even much longer for most GRBs. Our results, together with the observed extended emission of some short GRBs, indicate that the central engine activity time scale would be much longer than T90} for both long and short GRBs and the observed bimodal T90 distribution may be due to an instrumental selection effect.

The broadband SEDs of four gamma-ray NLS1s are compiled and explained with the leptonic model. It is found that their characteristics and fitting parameters of the observed SEDs are more like FSRQs than BL Lacs.

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

The prompt emission spectra of gamma-ray bursts (GRBs) usually have a dominant component that is well described by a phenomenological "Band" function. The physical origin of this spectral component is debated. Although the traditional interpretation is synchrotron radiation of non-thermal electrons accelerated in internal shocks or magnetic dissipation regions, a growing trend in the community is to interpret this component as modified thermal emission from a dissipative photosphere of a GRB fireball. We analyze the time dependent spectrum of GRB 110721A detected by {\em Fermi} GBM and LAT, and pay special attention to the rapid evolution of the peak energy $E_p$. We define a "death line" of thermally-dominated dissipative photospheric emission in the $E_p - L$ plane, and show that $E_p$ of GRB 110721A at the earliest epoch has a very high $E_p \sim 15$ MeV that is beyond the "death line". Together with the finding that an additional "shoulder" component exists in this burst that is consistent with a photospheric origin, we suggest that at least for some bursts, the "Band" component is not from a dissipative photosphere, but must invoke a non-thermal origin (e.g. synchrotron or inverse Compton) in the optically thin region of a GRB outflow. We also suggest that the rapid "hard-to-soft" spectral evolution is consistent with the quick discharge of magnetic energy in a magnetically-dominated outflow in the optically thin region.

Well-sampled optical lightcurves of 146 gamma-ray bursts (GRBs) are compiled from the literature. By empirical fitting we identify eight possible emission components and summarize the results in a "synthetic" lightcurve. Both optical flare and early shallow-decay components are likely related to long-term central engine activities. We focus on their statistical properties in this paper. Twenty-four optical flares are obtained from 19 GRBs. The isotropic R-band energy is smaller than 1% of $E_{\gamma, \rm iso}$. The relation between isotropic luminosities of the flares and gamma-rays follows $L^{\rm F}_{\rm R, iso}\propto L_{{\gamma}, \rm iso}^{1.11\pm 0.27}$. Later flares tend to be wider and dimmer, i.e., $w^{\rm F}\sim t^{\rm F}_{\rm p}/2$ and $L^{\rm F}_{\rm R, iso}\propto [t^{\rm F}_{\rm p}/(1+z)]^{-1.15\pm0.15}$. The detection probability of the optical flares is much smaller than that of X-ray flares. An optical shallow decay segment is observed in 39 GRBs. The relation between the break time and break luminosity is a power-law, with an index of $-0.78\pm 0.08$, similar to that derived from X-ray flares. The X-ray and optical breaks are usually chromatic, but a tentative correlation is found. We suggest that similar to the prompt optical emission that tracks $\gamma$-rays, the optical flares are also related to the erratic behavior of the central engine. The shallow decay component is likely related to a long-lasting spinning-down central engine or piling up of flare materials onto the blastwave. Mixing of different emission components may be the reason of the diverse chromatic afterglow behaviors.

We present a time-resolved spectral analysis of 51 long and 11 short bright GRBs observed with the {\em Femri}/GBM, paying special attention to $E_{\rm p}$ evolution within a same burst. Among 8 single-pulse long GRBs, 5 show hard-to-soft evolution, while 3 show intensity-tracking. The multi-pulse long GRBs have more complicated patterns. Among the GRBs whose time-resolved spectrum is available for the first pulse, almost half (15/32 GRBs) show clear hard-to-soft evolution, and the other half (17/32 GRBs) show clear intensity-tracking. Later pulses typically show the tracking behavior, although a hard-to-soft evolution pattern was identified in the 2nd pulse of 2 GRBs whose pulses are well separated. Statistically, the hard-to-soft evolution pulses tend to be more asymmetric than the intensity-tracking ones, with a steeper rising wing than the falling wing. Short GRBs have $E_{\rm p}$ tracking intensity exclusively with the 16ms time resolution analysis. We performed a simulation analysis, and suggest that at least for some bursts, the late intensity-tracking pulses could be a consequence of overlapping hard-to-soft pulses. However, the fact that the intensity-tracking pattern exists in the first pulse of multi-pulse long GRBs and some single-pulse GRBs suggest that intensity tracking is an independent component, which may operate in some late pulses as well. For the GRBs with measured redshifts, we present a time-resolved $E_{\rm p}-L_{\gamma, \rm iso}$ correlation analysis and show that the scatter of the correlation is comparable to that of the global Amati/Yonetoku relation. We discuss the predictions of various radiation models regarding $E_{\rm p}$ evolution, as well as the possibility of a precession jet in GRBs. It seems that the data pose great challenge to all these models, and hold the key to unveil the physics of GRB prompt emission.

Using the gamma-ray bursts simultaneously detected by Swift/BAT and Fermi/GBM we performed a joint spectral and temporal analysis of the prompt emission data and confirm the rough correlation between the BAT-band photon index Gamma_BAT and the peak spectral energy Epeak. With the redshift known sub-sample, we derived the isotropic gamma-ray energy E_gamma,iso and also confirm the E_gamma,iso - Epeak,rest relation, with a larger scatter than the Amati sample but consistent with GBM team analyses. We also compare the T_90 values derived in the GBM band with those derived in the BAT band and find that for long GRBs the BAT T_90 is usually longer than the GBM T_90, while for short GRBs the trend reverses. This is consistent with the soft/hard nature of long/short GRBs and suggests the importance of an energy-dependent temporal analysis of GRBs.

[abridged] We present results of modeling the SED and multiwavelength variability of the bright FSRQ PKS1510-089 with our time-dependent multizone Monte Carlo/Fokker-Planck code (Chen et al. 2001). As primary source of seed photons for inverse Compton scattering, we consider radiation from the broad line region (BLR), from the molecular torus, and the local synchrotron radiation (SSC). Different scenarios are assessed by comparing simulated light curves and SEDs with one of the best flares by PKS1510-089, in March 2009. The time-dependence of our code and its correct handling of light travel time effects allow us to fully take into account the effect of the finite size of the active region, and in turn to fully exploit the information carried by time resolved observed SEDs, increasingly available since the launch of Fermi. We confirm that the spectrum adopted for the external radiation has an important impact on the modeling of the SED, in particular for the lower energy end of the Compton component, observed in the X-ray band, which in turn is one of the most critical bands to assess the differences between EC and SSC emission. In the context of the scenario presented here, where the flaring is caused by the increase of the number of relativistic electrons ascribed to the effect of the interaction of a portion of the jet (blob) with a shock, we can not firmly discriminate the three main scenarios for gamma-ray emission. However, results show clearly the differences produced by a more realistic treatment of the emitting source in the shape of SEDs and their time variability over relevant, observable time-scales, and demonstrate the crucial importance of time-dependent multi-zone models to advance our understanding of the physics of these sources, by taking full advantage of the wealth of information offered by the high quality data of current multiwavelength campaigns.

BL Lac objects are the best candidates to study the jet properties since their spectral energy distributions (SEDs) are less contaminated by the emission from the accretion disk and external Compton processes. We compile the broadband SEDs observed with Fermi/LAT and other instruments from literature for 24 TeV BL Lac objects. Two SEDs, which are identified as a low or high state according to its flux density at 1 TeV, are available for each of ten objects. These SEDs can be explained well with the synchrotron+synchrotron-self-Compton model. We constrain the magnetic filed strength B and the Doppler factor delta of the radiation region by incorporating the chi^{2}-minimization technique and the gamma-ray transparency condition. Twenty-four out of the 34 SEDs can constrain B and delta in 1 sigma significance level and others only present the limits for both B and delta. The typical values of $B$ and delta are 0.1 ~ 0.6 G and 10 ~ 35, respectively. The derived values of gamma_b are significantly different among sources and even among the low and high states of a given source. Prominent flux variations with a clear spectral shift are observed and a tentative correlation between the ratio of the flux density at 1 TeV and the ratio of gamma_b in the low and high states is presented, likely indicating that the relativistic shocks for the electron acceleration may be responsible for the flux variations and the spectral shift. A weak anti-correlation between the jet power and the mass of the central black hole is observed, i.e., P_{jet} propto M^{-1}_{BH}, which disfavors the scenario of a pure accretion-driven jet. Implications for the blazar sequence and the intergalactic magnetic field from our results are also briefly discussed.

We compile from literature the broadband SEDs of twelve TeV blazars observed simultaneously or quasi-simultaneously with Fermi/LAT and other instruments. Two SEDs are available for each of the objects and the state is identified as a low or high state according to its flux density at GeV/TeV band. The observed SEDs of BL Lac objects (BL Lacs) are fitted well with the synchrotron + synchrotron-self-Compton (syn+SSC) model, whereas the SEDs of the two flat spectrum radio quasars (FSRQs) need to include the contributions of external Compton scattering. In this scenario, it is found that the Doppler factor delta of FSRQs is smaller than that of BL Lacs, but the magnetic field strength B of FSRQs is larger than that of BL Lacs. The increase of the peak frequency of the SEDs is accompanied with the increase of the flux for the individual sources, which seems opposite to the observational phenomena of the blazar sequence. We refer this phenomenon to blazar anti-sequence of spectral variability for individual TeV blazars. However, both the blazar sequence from FSRQs to BL Lacs and blazar anti-sequence of the spectral variability from low state to high state are accompanied by an increase of the break Lorentz factor of the electron's spectrum gamma_b and a decrease of B. We propose a model in which the mass accretion rate is the driving force behind both the blazar sequence for ensembles of blazars and the blazar anti-sequence for individual blazars. Specifically we suggest that the differences in mass accretion rate of different blazars produce the observed blazar sequence, but \Delta M in each blazar results in the observed blazar anti-sequence.

GRB 100418A is a long burst at z=0.624 without detection of any associated supernova (SN). Its lightcurves in both the prompt and afterglow phases are similar to GRB 060614, a nearby long GRB without an associated SN. We analyze the observational data of this event and discuss the possible origins of its multi-wavelength emission. We show that its joint lightcurve at 1 keV derived from Swift BAT and XRT observations is composed of two distinguished components. The first component, whose spectrum is extremely soft (\Gamma = 4.32), ends with a steep decay segment, indicating the internal origin of this component. The second component is a slowly-rising, broad bump which peaks at ~10^5 seconds post the BAT trigger. Assuming that the late bump is due to onset of the afterglow, we derive the initial Lorentz factor (Gamma_0) of the GRB fireball and find that it significantly deviates from the relation between the Gamma_0 and Eiso of typical GRBs. We also check whether it follows the same anti-correlation between X-ray luminosity and the break time observed in the shallow decay phase of many typical GRBs, which is usually regarded as a signal of late energy injection from the GRB central engine. However, we find that it does not obey this correlation. We propose that the late bump could be contributed by a two-component jet. We fit the second component with an off-axis jet model for a constant medium density and find the late bump can be represented by the model. The derived jet half-opening angle is 0.30 rad and the viewing angle is 0.315 rad. The medium density is 0.05 cm^-3, possibly suggesting that it may be from a merger of compact stars. The similarity between GRBs 060614 and 100418A may indicate that the two GRBs are from the same population and the late bump observed in the two GRBs may be a signal of a two-component jet powered by the GRB central engine.

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

Well-sampled optical lightcurves of 146 gamma-ray bursts (GRBs) are complied from the literature. Multiple optical emission components are extracted with power-law function fits to these lightcurves. We present a systematical analysis for statistical properties and their relations to prompt gamma-ray emission and X-ray afterglow for each component. We show that peak luminosity in the prompt and late flares are correlated and the evolution of the peak luminosity may signal the evolution of the accretion rate. No tight correlation between the shallow decay phase/plateau and prompt gamma-ray emission is found. Assuming that they are due to a long-lasting wind injected by a compact object, we show that the injected behavior favors the scenarios of a long-lasting wind after the main burst episode. The peak luminosity of the afterglow onset is tightly correlated with Eiso, and it is dimmer as peaking later. Assuming that the onset bump is due to the fireball deceleration by the external medium, we examine the Gamma_0-Eiso relation and find that it is confirmed with the current sample. Optical re-brightening is observed in 30 GRBs in our sample. It shares the same relation between the width and the peak time as found in the onset bump, but no clear correlation between the peak luminosity and Eiso as observed in the onset bumps is found. Although its peak luminosity also decays with time, the slope is much shallower than that of the onset peak. We get L t^{-1}_{p}$, being consistent with off-axis observations to an expanding external fireball in a wind-like circum medium. The late re-brightening may signal another jet component. Mixing of different emission components may be the reason for the observed chromatic breaks in different energy bands.

An intriguing hypothesis, that there exists a decaying X-ray emission component before the GRB trigger, has been suggested in order to explain the shallow decay phase of the X-ray afterglow detected in many Swift GRBs. If this "prior emission" is from an external shock, one would expect a corresponding optical emission component during the GRB prompt emission phase. In this paper we apply the available prompt optical emission data (both detections and upper limits) to constrain such a scenario. We fit the shallow and normal decay segments of the XRT light curves in our sample with a T_\Delta-shifted single power law, and extrapolate the X-ray flux back to the time of the early optical observations. We then use the synchrotron spectrum predicted by the standard external shock model to extrapolate from the X-ray flux to the optical band, and obtain the possible range of the predicted optical flux. Finally, we compare the predictions with the observations. In the cases where later optical data are available, we also compare the shapes of the optical light curves to the predicted optical light curves from the external shock prior emission model. We find that for a good fraction of GRBs (4 out of 8; up to December 2006), the available data already impose severe constraints on the hypothesis. In particular, the expected optical flux from the prior external shock model is higher than what the data allow. We conclude that if the shallow-decay X-ray component were from a prior emission component, it would have to be of an internal origin with optical flux suppressed.

Apparent redshift dependence of the jet opening angles ($\theta_{\rm j}$) of gamma-ray bursts (GRBs) is observed from current GRB sample. We investigate whether this dependence can be explained with instrumental selection effects and observational biases by a bootstrapping method. Assuming that (1) the GRB rate follows the star formation history and the cosmic metallicity history and (2) the intrinsic distributions of the jet-corrected luminosity ($L_{\rm \gamma}$) and $\theta_{\rm j}$ are a Gaussian or a power-law function, we generate a mock {\em Swift}/BAT sample by considering various instrumental selection effects, including the flux threshold and the trigger probability of BAT, the probabilities of a GRB jet pointing to the instrument solid angle and the probability of redshift measurement. Our results well reproduce the observed $\theta_{\rm j}-z$ dependence. We find that in case of $L_{\gamma}\propto \theta_{\rm j}^2$ good consistency between the mock and observed samples can be obtained, indicating that both $L_{\rm \gamma}$ and $\theta_{\rm j}$ are degenerate for a flux-limited sample. The parameter set $(L_{\rm \gamma}, \theta_{\rm j})=(4.9\times 10^{49} \rm {erg\ s}^{-1},\ 0.054 {rad})$ gives the best consistency for the current {\em Swift} GRB sample. Considering the beaming effect, the derived intrinsic local GRB rate accordingly is $2.85\times 10^2$ Gpc$^{-3}$ yr$^{-1}$, inferring that $\sim 0.59%$ of Type Ib/c SNe may be accompanied by a GRB.

Spectral fits to M87 core data from radio to hard x-ray are generated via a specially selected software suite, comprised of the HARM GRMHD accretion disk model and a 2D Monte Carlo radiation transport code. By determining appropriate parameter changes necessary to fit x-ray quiescent and flaring behavior of M87's core, we assess the reasonableness of various flaring mechanisms. This shows that an accretion disk model of M87's core out to 28 GM/c^2 can describe the inner emissions. High spin rates show GRMHD-driven polar outflow generation, without citing an external jet model. Our results favor accretion rate changes as the dominant mechanism of x-ray flux and index changes, with variations in density of approximately 20% necessary to scale between the average x-ray spectrum and flaring or quiescent spectra. The best fit parameters are black hole spin a/M > 0.8 and maximum accretion flow density n <= 3x10^7 cm^-3, equivalent to horizon accretion rates between m_dot = M_dot/M_dot_Edd ~ 2x10^-6 and 1x10^-5 (with M_dot_Edd defined assuming a radiative efficiency eta = 0.1). These results demonstrate that the immediate surroundings of M87's core are appropriate to explain observed x-ray variability.

We report simulation results of pair production by ultra-intense lasers irradiating a gold target using the GEANT4 Monte-Carlo code. Certain experimental features of the positron and electron energy spectra are reproduced, as well as trends with regard to target thickness and hot electron temperature Te. For Te in the range 5-10 MeV, the optimal target thickness for pair production is found to be about 3 mm. Further Monte-Carlo simulations may aid in the optimization of laser-driven positron sources.

Broadband spectral energy distributions (SEDs) are complied from literature for 24 TeV BL Lac objects. Two SEDs, which are identified as low and high states with their flux densities at $1\ {\rm TeV}$ ($F_{1 TeV}$), are available for each of 11 objects. We fit the SEDs with the single-zone synchrotron + synchrotron-self-Compton (SSC) model and determine the physical parameters of the jets. The model well represents the SEDs, and the observed relation between two peak frequencies $\nu_{\rm s}$ and $\nu_{\rm c}$ also favors the model. It is found that the break of the electron spectrum ($\gamma_{\rm b}$) is different among sources, but magnetic field strength distributes narrowly in $0.1-0.6 $ G, indicating that the shocks in the jets are significantly different among sources and the magnetic field may be independent of the shocks. The ratio of $F_{1 TeV}$ is correlated with the ratio of the $\gamma_{\rm b}$ in the low and high states, indicating that the relativistic shocks in the jets may be responsible for the observed flux variations and the spectral shift. The ratio of the two peak luminosities $L_{\rm c}/L_{\rm s}$ is anti-correlated with $\nu_{\rm s}$ in the co-moving frame, which is possibly due to the Klein-Nishina effect. No excess in the GeV band due to the interaction between the TeV photons and the extragalactic background light is observed, implying that the strength of the intergalactic magnetic field would be much larger than $ 10^{-16}$ G. An anti-correlation between jet power and the mass of the central black hole is observed, i.e., $P_{\rm jet}\propto M^{-1}_{\rm BH}$, disfavoring the scenario of a pure accretion-driven jet. We suggest that the spin energy extraction may be significant for powering jets in these sources, implying that the black holes with smaller mass may spin more rapidly.

(abridged) We present a new time-dependent multi-zone radiative transfer code and its application to study the SSC emission of Mrk 421. The code couples Fokker-Planck and Monte Carlo methods, in a 2D geometry. For the first time all the light travel time effects (LCTE) are fully considered, along with a proper treatment of Compton cooling, which depends on them. We study a set of simple scenarios where the variability is produced by injection of relativistic electrons as a `shock front' crosses the emission region. We consider emission from two components, with the second one either being pre-existing and co-spatial and participating in the evolution of the active region, or spatially separated and independent, only diluting the observed variability. Temporal and spectral results of the simulation are compared to the multiwavelength observations of Mrk 421 in March 2001. We find parameters that can adequately fit the observed SEDs and multiwavelength light curves and correlations. There remain however a few open issues, most notably: i) systematic soft intra-band X-ray lags. ii) The quadratic correlation between the TeV and X-ray flux during the flare decay has not been reproduced. These features are among those affected by the spatial extent and geometry of the source. The difficulty of producing hard X-ray lags is exacerbated by a bias towards soft lags caused by the combination of energy dependent radiative cooling time-scales and LCTE. About the second emission component, our results strongly favor the scenario where it is co-spatial and it participates in the flare evolution, suggesting that different phases of activity may occur in the same region. The cases presented in this paper represent only an initial study, and despite their limited scope they make a strong case for the need of true time-dependent and multi-zone modeling.

We use Monte Carlo/Fokker-Planck simulations to study the X-ray time lags. Our results show that soft lags will be observed as long as the decay of the flare is dominated by radiative cooling, even when acceleration and cooling timescales are similar. Hard lags can be produced in presence of a competitive achromatic particle energy loss mechanism if the acceleration process operates on a timescale such that particles are slowly moved towards higher energy while the flare evolves. In this type of scenario, the {\gamma} -ray/X-ray quadratic relation is also reproduced.

We present large scale 3D particle-in-cell (PIC) simulations to examine particle energization in magnetic reconnection of relativistic electron-positron (pair) plasmas. The initial configuration is set up as a relativistic Harris equilibrium without a guide field. These simulations are large enough to accommodate a sufficient number of tearing and kink modes. Contrary to the non-relativistic limit, the linear tearing instability is faster than the linear kink instability, at least in our specific parameters. We find that the magnetic energy dissipation is first facilitated by the tearing instability and followed by the secondary kink instability. Particles are mostly energized inside the magnetic islands during the tearing stage due to the spatially varying electric fields produced by the outflows from reconnection. Secondary kink instability leads to additional particle acceleration. Accelerated particles are, however, observed to be thermalized quickly. The large amplitude of the vertical magnetic field resulting from the tearing modes by the secondary kink modes further help thermalizing the non-thermal particles generated from the secondary kink instability. Implications of these results for astrophysics are briefly discussed.

We present a systematic analysis of the spectral and temporal properties of 17 GRBs co-detected by GBM and LAT on board the Fermi satellite by May 2010. We performed a time-resolved spectral analysis of all the bursts with the finest temporal resolution allowed by statistics, in order to avoid temporal smearing of different spectral components. We found that the time-resolved spectra of 14 out of 17 GRBs are best modeled with the Band function over the entire Fermi spectral range, which may suggest a common origin for emissions detected by LAT and GBM. GRB 090902B and GRB 090510 require the superposition between an MeV component and an extra power law component, with the former having a sharp cutoff above E_p. For GRB 090902B, this MeV component becomes progressively narrower as the time bin gets smaller, and can be fit with a Planck function as the time bin becomes small enough. In general, we speculate that phenomenologically there may be three elemental spectral components : (I) a Band-function component (e.g. in GRB 080916C) that extends in a wide energy range and does not narrow with reducing time bins, which may be of the non-thermal origin; (II) a quasi-thermal component (e.g. in GRB 090902B) with the spectra progressively narrowing with reducing time bins; and (III) another non-thermal power law component extending to high energies. The spectra of different bursts may be decomposed into one or more of these elemental components. We compare this sample with the BATSE sample and investigate some correlations among spectral parameters. We discuss the physical implications of the data analysis results for GRB prompt emission, including jet compositions (matter-dominated vs. Poynting-flux-dominated outflow), emission sites (internal shock, external shock or photosphere), as well as radiation mechanisms (synchrotron, synchrotron self-Compton, or thermal Compton upscattering).

We present results of models of the physical space and parameters of the accretion disk of Sagittarius A*, as well as simulations of its emergent spectrum. This begins with HARM, a 2D general relativistic magneto-hydrodynamic (GRMHD) model, specifically set up to evolve the space around a black hole. Data from HARM are then fed into a 2D Monte-Carlo (MC) code which generates and tracks emitted photons, allowing for absorption and scattering before they escape the volume.

X-ray Flash (XRF) 100316D, a nearby super-long under-luminous burst with a peak energy E_p \sim 20 keV, was detected by Swift and was found to be associated with an energetic supernova SN 2010bh. Both the spectral and the temporal behavior of this burst are rather similar to that of XRF 060218, except that the latter was associated with a "less energetic" SN 2006aj and had a prominent soft thermal emission component in the spectrum. We analyze the spectral and temporal properties of this burst, and interpret the prompt gamma-ray emission and the early X-ray plateau emission as synchrotron emission from a dissipating Poynting-flux-dominated outflow, probably powered by a magnetar with a spin period of $P \sim 10$ ms and the polar cap magnetic field $B_{\rm p} \sim 3\times 10^{15}$ G. The energetic supernova SN 2010bh associated with this burst is, however, difficult to interpret within the slow magnetar model, which implies that the nascent magnetar may spin much faster with an initial rotation period $\sim 1$ ms, and thus suggests a delay between the core collapse and the emergence of the relativistic magnetar wind from the star. The diverse behaviors of low-luminosity GRBs and their associated SNe may be understood within a unified picture that invokes different initial powers of the central engine and different delay times between the core collapse and the emergence of the relativistic jet from the star.

Recent Swift observations suggest that the traditional long vs. short GRB classification scheme does not always associate GRBs to the two physically motivated model types, i.e. Type II (massive star origin) vs. Type I (compact star origin). We propose a new phenomenological classification method of GRBs by introducing a new parameter epsilon=E_{gamma, iso,52}/E^{5/3}_{p,z,2}, where E_{\gamma,iso} is the isotropic gamma-ray energy (in units of 10^{52} erg), and E_{p,z} is the cosmic rest frame spectral peak energy (in units of 100 keV). For those short GRBs with "extended emission", both quantities are defined for the short/hard spike only. With the current complete sample of GRBs with redshift and E_p measurements, the epsilon parameter shows a clear bimodal distribution with a separation at epsilon ~ 0.03. The high-epsilon region encloses the typical long GRBs with high-luminosity, some high-z "rest-frame-short" GRBs (such as GRB 090423 and GRB 080913), as well as some high-z short GRBs (such as GRB 090426). All these GRBs have been claimed to be of the Type II origin based on other observational properties in the literature. All the GRBs that are argued to be of the Type I origin are found to be clustered in the low-epsilon region. They can be separated from some nearby low-luminosity long GRBs (in 3sigma) by an additional T_{90} criterion, i.e. T_{90,z}<~ 5 s in the Swift/BAT band. We suggest that this new classification scheme can better match the physically-motivated Type II/I classification scheme.

GRB 090426 is a short duration burst detected by Swift ($T_{90}\sim 1.28$ s in the observer frame, and $T_{90}\sim 0.33$ s in the burst frame at $z=2.609$). Its host galaxy properties and some $\gamma$-ray related correlations are analogous to those seen in long duration GRBs, which are believed to be of a massive-star origin (so-called Type II GRBs). We present the results of its early optical observations with the 0.8-m TNT telescope at Xinglong observatory, and the 1-m LOAO telescope at Mt. Lemmon Optical Astronomy Observatory in Arizona. Our well-sampled optical afterglow lightcurve covers from $\sim 90$ seconds to $\sim 10^4$ seconds post the GRB trigger. It shows two shallow decay episodes that are likely due to energy injection, which end at $\sim 230$ seconds and $\sim 7100$ seconds, respectively. The decay slopes post the injection phases are consistent with each other ($\alpha\simeq 1.22$). The X-ray afterglow lightcurve appears to trace the optical, although the second energy injection phase was missed due to visibility constraints introduced by the {\em Swift} orbit. The X-ray spectral index is $\beta_X\sim 1.0$ without temporal evolution. Its decay slope is consistent with the prediction of the forward shock model. Both X-ray and optical emission is consistent with being in the same spectral regime above the cooling frequency ($\nu_c$). The fact that $\nu_c$ is below the optical band from the very early epoch of the observation provides a constraint on the burst environment, which is similar to that seen in classical long duration GRBs. We therefore suggest that death of a massive star is the possible progenitor of this short burst.

A time-resolved spectral analysis for a sample of 22 intense, broad GRB pulses from the BATSE GRB sample is presented. We fit the spectra with the Band function and investigate the correlation between the observed flux (F) and the peak energy (E_p) of the $\nu f_\nu$ spectrum in the rising and decaying phases of these pulses. Two kinds of E_p evolution trends, i.e., hard-to-soft (the two-third pulses in our sample) and $E_{\rm p}$-tracing-$F$ (the one-third pulses in our sample) are observed in pulses from different GRBs and even from different pulses of the same burst. No dependence of spectral evolution feature on the pulse shape is found. A tight $F-E_{\rm p}$ positive correlation is observed in the decaying phases, with a power-law index $\sim 2.2$, which is much shallower than that expectation of the curvature effect. In the rising phase, the observed $F$ is either correlated or anti-correlated with $E_{\rm p}$, depending on the spectral evolution feature, and the power-law index of the correlation is dramatically different among pulses. More than $80%$ of the low energy photon indices in the time-resolved spectra whose $E_{\rm p}$ is anti-correlated with $F$ during the rising phase violate the death line of the synchrotron radiation, disfavoring the synchrotron radiation model for these gamma-rays. The $F-E_{\rm p}$ correlation, especially for those GRBs with $E_{\rm p}$-tracking-$F$ spectral evolution, may be due to the viewing angle and jet structure effects. In this cenario, the observed $F-E_{\rm p}$ correlation in the rising phase may be due to the line of sight from off-beam to on-beam toward a structured jet (or jitter), and the decaying phase is contributed by both the on-beam emission and the decayed photons from high latitude of the GRB fireball,resulting in a shallower slope of the observed $F-E_{\rm p}$ correlation than that predicted by the pure curvature effect.

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

Since the launch of Swift satellite, the detections of high-z (z>4) long gamma-ray bursts (LGRBs) have been rapidly growing, even approaching the very early Universe (the record holder currently is z=8.3). The observed high-z LGRB rate shows significant excess over that estimated from the star formation history. We investigate what may be responsible for this high productivity of GRBs at high-z through Monte Carlo simulations, with effective Swif/BAT trigger and redshift detection probabilities based on current Swift/BAT sample and CGRO/BATSE LGRB sample. We compare our simulations to the Swift observations via log N-log P, peak luminosity (L) and redshift distributions. In the case that LGRB rate is purely proportional to the star formation rate (SFR), our simulations poorly reproduce the LGRB rate at z>4, although the simulated log N-log P distribution is in good agreement with the observed one. Assuming that the excess of high-z GRB rate is due to the cosmic metallicity evolution or unknown LGRB rate increase parameterized as (1+z)^delta, we find that although the two scenarios alone can improve the consistency between our simulations and observations, incorporation of them gives much better consistency. We get 0.2<epsilon<0.6 and delta<0.6, where epsilon is the metallicity threshold for the production of LGRBs. The best consistency is obtained with a parameter set (epsilon, delta)=(~0.4, ~0.4), and BAT might trigger a few LGRBs at z~14. With increasing detections of GRBs at z>4 (~15% of GRBs in current Swift LGRB sample based on our simulations), a window for very early Universe is opening by Swift and up-coming SVOM missions.

We report on 5 Chandra observations of the X-ray afterglow of the Gamma-Ray Burst GRB 060729 performed between 2007 March and 2008 May. In all five observations the afterglow is clearly detected. The last Chandra pointing was performed on 2008-May-04, 642 days after the burst - the latest detection of a GRB X-ray afterglow ever. A reanalysis of the Swift XRT light curve together with the three detections by Chandra in 2007 reveals a break at about 1.0 Ms after the burst with a slight steepening of the decay slope from alpha = 1.32 to 1.61. This break coincides with a significant hardening of the X-ray spectrum, consistent with a cooling break in the wind medium scenario, in which the cooling frequency of the afterglow crosses the X-ray band. The last two Chandra observations in 2007 December and 2008 May provide evidence for another break at about one year after the burst. If interpreted as a jet break, this late-time break implies a jet half opening angle of about 14 degrees for a wind medium. Alternatively, this final break may have a spectral origin, in which case no jet break has been observed and the half-opening angle of the jet of GRB 060729 must be larger than about 15 degrees for a wind medium. We compare the X-ray afterglow of GRB 060729 in a wind environment with other bright X-ray afterglows, in particular GRBs 061121 and 080319B, and discuss why the X-ray afterglow of GRB 060729 is such an exceptionally long-lasting event.

The onset of GRB afterglow is characterized by a smooth bump in the early afterglow lightcurve. We make an extensive search for such a feature. Twenty optically selected GRBs and 12 X-ray selected GRBs are found, among which 17 optically selected GRBs and 2 X-ray-selected GRBs have redshift measurements. We fit the lightcurves with a smooth broken power-law and measure the temporal characteristic timescales of the bumps at FWHM. Strong mutual correlations among these timescales are found, and a dimmer and broader bump tends to peak at a later peak time. The ratio of rising to decaying timescales is almost universal among bursts, but the ratio of the rising time to the peak time varies from 0.3~1. The E_iso is tightly correlated with the peak luminosity and the peak time of the bump in the burst frame. Assuming that the bumps signal the deceleration of the GRB fireballs in a constant density medium, we calculate the initial Lorentz factor (Gamma_0) and the deceleration radius (R_dec) of the GRBs in the optical-selected sample. It is found that Gamma_0 are typically a few hundreds, and the typical deceleration radius is R_dec~10^{17} cm. More intriguingly, a tight correlation between the Gamma_0 and E_iso is found, namely Gamma_0 ~ 195 E_iso, 52}^{0.27} (satisfied for both the optical and X-ray z-known samples). It is helpful to understand GRB physics, and may serve as an indicator of Gamma_0. We find that the early bright X-rays are usually dominated by a different component from the external shock emission, but occasionally (for one case) an achromatic deceleration feature is observed. Components in X-rays would contribute to the diversity of the observed X-ray lightcurves (abridge).

The observed radio-optical-X-ray spectral energy distributions (SEDs) of 22 hot spots and 45 knots in the jets of 35 active galactic nuclei are complied from literature and modeled with single-zone lepton models. It is found that the observed luminosities at 5 GHz (L_5GHz) and at 1 keV (L_1keV) are tightly correlated, and the two kinds of sources can be roughly separated with a division of L_1keV=L_5GHz. Our SED fits show that the mechanisms of the X-rays are diverse. Considering the sources at rest, the synchrotron-self-Compton (SSC) scattering would dominate the IC process. This model can interpret the X-rays of some hot spots with a magnetic field strength (B_ssc^delta=1) being consistent with the equipartition magnetic field (B_eq^delta=1) in one order of magnitude, but an unreasonably low B_ssc^delta=1 is required to model the X-rays for all knots. The ratio R_B=B_eq^delta=1/B_ssc^delta=1 is greater than 1 and it is tightly anti-correlated with ratio R_L= L_1keV/L_5GHz for both the knots and the hot spots. We propose that the deviation may be due to the neglect of the relativistic bulk motion for these sources. Considering this effect, we show that the IC/CMB model well explains the X-ray emission for most sources. Both B_eq' and delta are tentatively correlated with R_L. Corrected by the beaming effect, the L'_5GHz-L'_1keV relations for the two kinds of sources are even tighter than the observed ones. These facts suggest that, under the equipartition condition, the observational differences of the X-rays from the knots and hot spots may be mainly due the differences on the Doppler boosting effect and the co-moving magnetic field of the two kinds of sources. Our IC scattering models predict a prominent GeV-TeV component in the SEDs for some sources, which are detectable with H.E.S.S. and Fermi/LAT.

Both analytical and numerical works show that magnetic reconnection must occur in hot accretion flows. This process will effectively heat and accelerate electrons. In this paper we use the numerical hybrid simulation of magnetic reconnection plus test-electron method to investigate the electron acceleration and heating due to magnetic reconnection in hot accretion flows. We consider fiducial values of density, temperature, and magnetic parameter $\beta_e$ (defined as the ratio of the electron pressure to the magnetic pressure) of the accretion flow as $n_{0} \sim 10^{6} {\rm cm^{-3}}$, $T_{e}^0\sim 2\times 10^9 {\rm K}$, and $\beta_e=1$. We find that electrons are heated to a higher temperature $T_{e}=5\times 10^9$K, and a fraction $\eta\sim 8%$ of electrons are accelerated into a broken power-law distribution, $dN(\gamma)\propto \gamma^{-p}$, with $p\approx 1.5$ and 4 below and above $\sim 1$ MeV, respectively. We also investigate the effect of varying $\beta$ and $n_0$. We find that when $\beta_e$ is smaller or $n_0$ is larger, i.e, the magnetic field is stronger, $T_e$, $\eta$, and $p$ all become larger.

By analyzing the Swift/XRT lightcurves detected before 2009 July, we find 19 cases that monotonously decay as a single power law (SPL) with an index of 1 ~ 1.7 from tens (or hundreds) to ~ 10^5 seconds post the GRB trigger, apparently different from the canonical lightcurves characterized by a shallow-to-normal decay transition. No statistical difference is found in their prompt gamma-rays, and the X-ray properties of the two samples are also similar, although the SPL sample tend to have a slightly lower NH value of the host galaxies and larger energy release compared with the canonical sample. The SPL XRT lightcurves in the burst frame gradually merge into a conflux.The normal decay segment for the canonical sample has the same feature. Similar to the normal decay segment, the SPL lightcurves satisfy the closure relations of external shock models. If the X-rays are the afterglow of the GRB fireball, our results indicate that the shallow decay would be due to energy injection and the total energy budget after injection for both samples of GRBs is comparable. More intriguing, we find that a prior X-ray emission model proposed by Yamazaki is more straightforward to interpret the data. We show that the zero times (T0) of the X-rays atisfy a log-normal distribution, and the negligible T0's of the SPL sample are consistent with being the tail of T0 distributions at low end. Referenced to T0, the canonical XRT lightcurves well trace the SPL lightcurves. The T0's of the canonical lightcurves in our analysis are usually much larger than the offsets of the known precursors from the main GRBs, indicating the X-rays would be an emission component from external shocks prior to GRB trigger. The lack of detection of a jet-like break in most XRT lightcurves implies that the opening angle of its jet would be usually large.

We investigate the properties of both the prompt and X-ray afterglows of gamma-ray bursts (GRBs) in the burst frame with a sample of 33 Swift GRBs. Assuming that the steep decay segment in the canonical X-ray afterglow lightcurves is due to the curvature effect, we fit the lightcurves with a broken power-law to derive the zero time of the last emission epoch of the prompt emission (t1) and the beginning as well as the end time of the shallow decay segment (t2 and t3).We show that both the isotropic peak gamma-ray luminosity and gamma-ray energy are correlated with the isotropic X-ray energy of the shallow decay phase and the isotropic X-ray luminosity at t2. We infer the properties of the progenitor stars based on a model proposed by Kumar et al. who suggested that both the prompt gamma-rays and the X-ray afterglows are due to the accretions of different layers of materials of the GRB progenitor star by a central black hole (BH). We find that most of the derived masses of the core layers are 0.1-5 solar mass with a radius of 10^8-10^10 cm. The rotation parameter is correlated with the burst duration, being consistent with the expectation of collapsar models. The estimated radii and the masses of the fall-back materials for the envelope layers are 10^10-10^12 cm and 10^-3~1 solar mass, respectively. The average accretion rates in the shallow decay phase are correlated with those in the prompt gamma-ray phase, but they are much lower. The derived radii of the envelope are smaller than the photospheric radii of Wolf-Rayet (WR) stars. It is interesting that the assembled mass density profile for the bursts in our sample is also well consistent with the simulation for a pre-supernova star with 25 solar mass.

We present results of simulations of the spectrum of the accretion flow onto the supermassive black hole in our Galactic Centre, Sagittarius A*, generated with a coupling of Monte-Carlo (MC) radiation and general relativistic magnetohydrodynamic (GRMHD) codes. In our modeling, we use the 2D HARM GRMHD code to first model the physical parameters of the disk, then feed its results into our 2D MC photon transport code. We will discuss results obtained which fit radio, IR, and Chandra-obtained flaring or quiescent x-ray data points, as well as the validity of the amount of scaling of input parameters (density, temperature, and magnetic field) required to fit these points. HARM output will be used to suggest whether the scaling is within reasonable limits.

(Abridged) The two high-redshift gamma-ray bursts, GRB 080913 at z=6.7 and GRB 090423 at z=8.3, recently detected by Swift appear as intrinsically short, hard GRBs. They could have been recognized by BATSE as short/hard GRBs should they have occurred at z <= 1. We perform a more thorough investigation on two physically distinct types (Type I/II) of cosmological GRBs and their observational characteristics. We reiterate the definitions of Type I/II GRBs and review the observational criteria and their physical motivations. Contrary to the traditional approach of assigning the physical category based on the gamma-ray properties (duration, hardness, and spectral lag), we take an alternative approach to define the Type I and Type II Gold Samples using several criteria that are more directly related to the GRB progenitors, and study the properties of the two Gold Samples and compare them with the traditional long/soft and short/hard samples. We find that the Type II Gold Sample reasonably tracks the long/soft population, although it includes several intrinsically short (shorter than 1s in the rest frame) GRBs. The Type I Gold Sample only has 5 GRBs, 4 of which are not strictly short but have extended emission. Other short/hard GRBs detected in the Swift era represent the BATSE short/hard sample well, but it is unclear whether all of them belong to Type I. We suggest that some (probably even most) high-luminosity short/hard GRBs instead belong to Type II. We suggest that GRB 080913 and GRB 090423 are more likely Type II events. We re-emphasize the importance of invoking multiple observational criteria, and cautiously propose an operational procedure to infer the physical origin of a given GRB with available multiple observational criteria, with various caveats laid out.

We derive analytic formulas for the radiation power output when electrons are accelerated by a relativistic comoving kinetic Poynting flux, and validate these analytic results with Particle-In-Cell simulations. We also derive analytically the critical frequency of the radiation spectrum. Potential astrophysical applications of these results are discussed. A quantitative model of gamma-ray bursts based on the breakout of kinetic Poynting flux is presented.

We present a systematic temporal and spectral study of all Swift-XRT observations of GRB afterglows discovered between 2005 January and 2007 December. After constructing and fitting all light curves and spectra to power-law models, we classify the components of each afterglow in terms of the canonical X-ray afterglow and test them against the closure relations of the forward shock models for a variety of parameter combinations. The closure relations are used to identify potential jet breaks with characteristics including the uniform jet model with and without lateral spreading and energy injection, and a power-law structured jet model, all with a range of parameters. With this technique, we survey the X-ray afterglows with strong evidence for jet breaks (~12% of our sample), and reveal cases of potential jet breaks that do not appear plainly from the light curve alone (another ~30%), leading to insight into the missing jet break problem. Those X-ray light curves that do not show breaks or have breaks that are not consistent with one of the jet models are explored to place limits on the times of unseen jet breaks. The distribution of jet break times ranges from a few hours to a few weeks with a median of ~1 day, similar to what was found pre-Swift. On average Swift GRBs have lower isotropic equivalent gamma-ray energies, which in turn results in lower collimation corrected gamma-ray energies than those of pre-Swift GRBs. Finally, we explore the implications for GRB jet geometry and energetics.

We present results from 2.5-dimensional Particle-in-Cell simulations of the interaction of nonlinear Alfven waves with thin current sheets in relativistic plasmas. We find that the Alfven waves cause the current sheet to bend and kink and increase its dissipation. The electrons are eventually heated to form a double Maxwellian, with the hotter Maxwellian caused by the current sheet dissipation and cooler Maxwellian caused by the Alfven turbulence cascade. These results may have important implications for the kinetic dissipation of MHD turbulence in which both nonlinear Alfven waves and current sheets are present, such as turbulence in accretion flows driven by the saturated magnetorotational instability (MRI).

We propose a star-quake model to understand X-ray flares of both long and short Gamma-ray bursts (GRBs) in a solid quark star regime. Two kinds of central engines for GRBs are available if pulsar-like stars are actually (solid) quark stars, i.e., the SNE-type GRBs and the SGR-type GRBs. It is found that a quark star could be solidified about 10^3 to 10^6 s later after its birth if the critical temperature of phase transition is a few MeV, and then a new source of free energy (i.e., elastic and gravitational ones, rather than rotational or magnetic energy) could be possible to power GRB X-ray flares.

The apparent spectral evolution observed in the steep decay phase of many GRB early afterglows raises a great concern of the high-latitude "curvature effect" interpretation of this phase. However, previous curvature effect models only invoked a simple power law spectrum upon the cessation of the prompt internal emission. We investigate a model that invokes the "curvature effect" of a more general non-power-law spectrum and test this model with the Swift/XRT data of some GRBs. We show that one can reproduce both the observed lightcurve and the apparent spectral evolution of several GRBs using a model invoking a power-law spectrum with an exponential cut off. GRB 050814 is presented as an example.

We report on the detection by Swift of GRB 080913, and subsequent optical/near-infrared follow-up observations by GROND which led to the discovery of its optical/NIR afterglow and the recognition of its high-z nature via the detection of a spectral break between the i' and z' bands. Spectroscopy obtained at the ESO-VLT revealed a continuum extending down to lambda = 9400 A, and zero flux for 7500 A < lambda<9400 A, which we interpret as the onset of a Gunn-Peterson trough at z=6.695+-0.025 (95.5% conf. level), making GRB 080913 the highest redshift GRB to date, and more distant than the highest-redshift QSO. We note that many redshift indicators which are based on promptly available burst or afterglow properties have failed for GRB 080913. We report on our follow-up campaign and compare the properties of GRB 080913 with bursts at lower redshift. In particular, since the afterglow of this burst is fainter than typical for GRBs, we show that 2m-class telescopes can identify most high-redshift GRBs.

The intriguing observations of Swift/BAT X-ray flash XRF 060218 and the BATSE-BeppoSAX gamma-ray burst GRB 980425, both with much lower luminosity and redshift compared to other observed bursts, naturally lead to the question of how these low-luminosity (LL) bursts are related to high-luminosity (HL) bursts. Incorporating the constraints from both the flux-limited samples observed with CGRO/BATSE and Swift/BAT and the redshift-known GRB sample, we investigate the luminosity function for both LL- and HL-GRBs through simulations. Our multiple criteria, including the log N - log P distributions from the flux-limited GRB sample, the redshift and luminosity distributions of the redshift-known sample, and the detection ratio of HL- and LL- GRBs with Swift/BAT, provide a set of stringent constraints to the luminosity function. Assuming that the GRB rate follows the star formation rate, our simulations show that a simple power law or a broken power law model of luminosity function fail to reproduce the observations, and a new component is required. This component can be modeled with a broken power, which is characterized by a sharp increase of the burst number at around L < 10^47 erg s^-1}. The lack of detection of moderate-luminosity GRBs at redshift ~0.3 indicates that this feature is not due to observational biases. The inferred local rate, rho_0, of LL-GRBs from our model is ~ 200 Gpc^-3 yr^-1 at ~ 10^47 erg s^-1, much larger than that of HL-GRBs. These results imply that LL-GRBs could be a separate GRB population from HL-GRBs. The recent discovery of a local X-ray transient 080109/SN 2008D would strengthen our conclusion, if the observed non-thermal emission has a similar origin as the prompt emission of most GRBs and XRFs.

Preliminary results of our analysis on the extended emission of short/medium duration GRBs observed with Swift/BAT are presented. The Bayesian blocks algorithm is used to analyze the burst durations and the temporal structure of the lightcurves in different energy bands. We show here the results of three bursts (GRBs 050724, 061006 and 070714B) that have a prominent soft extended emission component in our sample. The extended emission of these bursts is a continuous, flickering-liked component, lasting $\sim 100$ seconds post the GRB trigger at 15-25 keV bands. Without considering this component, the three bursts are classified as short GRBs, with $T_{90}=2\sim 3$ seconds. GRB 060614 has an emission component similar to the extended emission, but this component has pulse-liked structure, possibly indicating that this emission component is different from that observed in GRBs 050724, 061006, and 070714B. Further analysis on the spectral evolution behavior of the extended emission component is on going.

We use electromagnetic Particle-in-Cell codes to investigate the two-stream instability and Weibel instability when two beams of plasma collide head-on relativistically. We compare the dissipation of kinetic energies in the two cases. We also study the evolution of the resultant electron energy distribution function produced by such dissipation.

We present a systematic survey of the temporal and spectral properties of all GRB X-ray afterglows observed by Swift-XRT between January 2005 and July 2007. We have constructed a catalog of all light curves and spectra and investigate the physical origin of each afterglow segment in the framework of the forward shock models by comparing the data with the closure relations. We search for possible jet-like breaks in the lightcurves and try to explain some of the "missing" X-ray jet breaks in the lightcurves.

We derive analytic formulas for the power output and critical frequency of radiation by electrons accelerated by relativistic kinetic Poynting flux, and validate these results with Particle-In-Cell plasma simulations. We find that the in-situ radiation power output and critical frequency are much below those predicted by the classical synchrotron formulae. We discuss potential astrophysical applications of these results.

The Swift/XRT data of 179 GRBs (from 050124 to 070129) and the optical afterglow data of 57 pre- and post-Swift GRBs are analyzed, in order to systematically investigate the jet-like breaks in the X-ray and optical afterglow lightcurves. We find that not a single burst can be included in the ``Platinum'' sample, in which the data satisfy all the criteria of a jet break. By releasing one or more requirements to define a jet break, some candidates of various degrees could be identified. In the X-ray band, 42 out of the 103 well-sampled X-ray lightcurves have a decay slope of the post-break segment >1.5 (``Bronze'' sample), and 27 of them also satisfy the closure relations of the forward models (``Silver'' sample). The numbers of the ``Bronze'' and ``Silver'' candidates in the optical lightcurves are 27 and 23, respectively. Thirteen bursts have well-sampled optical and X-ray lightcurves, but only seven cases are consistent with an achromatic break, but even in these cases only one band satisfies the closure relations (``Gold'' sample). The observed break time in the XRT lightcurves is systematically earlier than that in the optical bands. All these raise great concerns in interpreting the jet-like breaks as jet breaks and further inferring GRB energetics from these breaks. By assuming that these breaks are jet breaks, we perform a similar analysis as previous work to calculate the jet opening angle (theta_j) and energetics (E_k) with the ``Silver'' and ``Gold'' jet break candidates. The derived E_K distribution reveals a much larger scatter than the pre-Swift sample. A tentative anti-correlation between theta_j and E_{K,iso} is found for both the pre-Swift and Swift GRBs, indicating that the E_K could still be quasi-universal, if the breaks in discussion are indeed jet breaks(abridge).

An early steep decay component following the prompt GRBs is commonly observed in {\em Swift} XRT light curves, which is regarded as the tail emission of the prompt gamma-rays. Prompted by the observed strong spectral evolution in the tails of GRBs 060218 and 060614, we present a systematic time-resolved spectral analysis for the {\em Swift} GRB tails detected between 2005 February and 2007 January. We select a sample of 44 tails that are bright enough to perform time-resolved spectral analyses. Among them 11 tails are smooth and without superimposing significant flares, and their spectra have no significant temporal evolution. We suggest that these tails are dominated by the curvature effect of the prompt gamma-rays due to delay of propagation of photons from large angles with respect to the line of sight . More interestingly, 33 tails show clear hard-to-soft spectral evolution, with 16 of them being smooth tails directly following the prompt GRBs,while the others being superimposed with large flares. We focus on the 16 clean, smooth tails and consider three toy models to interpret the spectral evolution. The curvature effect of a structured jet and a model invoking superposition of the curvature effect tail and a putative underlying soft emission component cannot explain all the data. The third model, which invokes an evolving exponential spectrum, seems to reproduce both the lightcurve and the spectral evolution of all the bursts, including GRBs 060218 and 060614. More detailed physical models are called for to understand the apparent evolution effect.

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

We analyze the properties of the shallow decay segment with a sample of 53 long Swift GRBs. We show that the distributions of its characteristics are log-normal or normal, and its isotropic X-ray energy (E_{iso,X}) is linearly correlated with the prompt gamma-ray energy, but with a steeper photon spectrum except for some X-ray flashes. No significant spectral evolution is observed from this phase to the follow-up phase, and the follow-up phase is usually consistent with the external shock models, implying that this segment likely due to a refreshed external shock. Within this scenario, the data are generally consistent with a roughly constant injection luminosity up to the end of this phase $t_b$. A positive correlation between E_{iso, X} and t_b also favors this scenario. Among the 13 bursts that have well-sampled optical light curves, 6 have an optical break around t_b and the breaks are consistent with being achromatic. However, the other 7 bursts either do not show an optical break or have a break at a different epoch than t_b. This raises a concern to the energy injection scenario, suggesting that the optical and X-ray emissions may not be the same component at least for a fraction of bursts. There are 4 significant outliers in the sample, and this shallow decau phase in these bursts is immediately followed by a very steep decay after t_b, which is inconsistent with any external shock model. The optical data of these bursts evolve independently from the X-ray data. These X-ray plateaus likely have an internal origin and demand continuous operation of a long-term GRB central engine. We conclude that the observed shallow decay phase likely has diverse physical origins(Abridged).

We discuss our radio (Australia Telescope Compact Array and Australian Long Baseline Array) and X-ray (XMM-Newton) monitoring observations of the unusual ultraluminous supernova SN 1978K in NGC 1313 at ~25 years after the explosion. SN 1978K is a rare example of a Type IIn supernova that has remained bright enough to have long-term X-ray and radio observations. The observations probe the dense medium that was ejected by the progenitor star prior to its explosion; the star might have been a Luminous Blue Variable. The radio imaging shows that the source remains compact, but it may be marginally resolved. The radio monitoring shows deviations from a smooth decay suggesting that inhomogeneities are present in the radio emitting region. It appears that a major change occurred in the mass-loss rate of the progenitor star shortly before the supernova event. The X-ray emission between 2000 and 2006 is consistent with the radiation coming from two shocks. All the X-ray data can be fit using the same model (with no systematic evolution or short-term variability) but this has a surprising requirement: the X-ray emitting regions have a very large abundance of helium. This would be consistent with the X-ray emitting shocks being located in a helium-rich layer that was ejected by the progenitor star, or helium-rich material was ejected in the supernova explosion. The unusual properties of the supernova motivated a search for an associated gamma-ray burst (GRB). We show that SN 1978K was inside the ~4 sigma error box of GRB 771029. If this association is correct, the GRB was exceptionally underluminous. However, the quality of the gamma-ray burst locations at that time was poor, and this is likely just a chance alignment.

We present a detailed analysis of Swift multi-wavelength observations of GRB 070110 and its remarkable afterglow. The early X-ray light curve, interpreted as the tail of the prompt emission, displays a spectral evolution already seen in other gamma-ray bursts. The optical afterglow shows a shallow decay up to ~2 d after the burst, which is not consistent with standard afterglow models. The most intriguing feature is a very steep decay in the X-ray flux at ~20 ks after the burst, ending an apparent plateau. The abrupt drop of the X-ray light curve rules out an external shock as the origin of the plateau in this burst and implies long-lasting activity of the central engine. The temporal and spectral properties of the plateau phase point towards a continuous central engine emission rather than the episodic emission of X-ray flares. We suggest that the observed X-ray plateau is powered by a spinning down central engine, possibly a millisecond pulsar, which dissipates energy at an internal radius before depositing energy into the external shock.

The absence of a supernova accompanying the nearby long GRB 060614 poses a great puzzle about the progenitor of this event and challenges the current GRB classification scheme. This burst displays a short-hard emission episode followed by extended soft emission with strong spectral evolution. Noticing that this burst has an isotropic gamma-ray energy only ~8 times that of GRB 050724, a good candidate of merger-type short GRBs, we generate a ``pseudo'' burst that is ~8 times less energetic than GRB 060614 based on the spectral properties of GRB 060614 and the Ep ~ Eiso^{1/2} (Amati) relation. We find that this pseudo-burst would have been detected by BATSE as a marginal short-duration GRB, and would have properties in the Swift BAT and XRT bands similar to GRB 050724. This suggests that GRB 060614 is likely a more intense event in the traditional short-hard GRB category as would be detected by BATSE. Events like GRB 060614 that seem to defy the traditional short vs. long classification of GRBs may require modification of our classification terminology for GRBs. By analogy with supernova classifications, we suggest that GRBs be classified into Type I (typically short and associated with old populations) and Type II (typically long and associated with young populations). We propose that GRB 060614 belongs to Type I, and predict that similar events will be detected in elliptical galaxies.

We review recent PIC simulation results which show that double-sided irradiation of a thin over-dense plasma slab with ultra-intense laser pulses from both sides can lead to sustained comoving Poynting flux acceleration of electrons to energies much higher than the conventional ponderomotive limit. The result is a robust power-law electron momentum spectrum similar to astrophysical sources. We discuss future ultra-intense laser experiments that may be used to simulate astrophysical particle acceleration.

The robotic ROTSE-III telescope network detected prompt optical emission contemporaneous with the gamma-ray emission of Swift events GRB051109A and GRB051111. Both datasets have continuous coverage at high signal-to-noise levels from the prompt phase onwards, thus the early observations are readily compared to the Swift XRT and BAT high energy detections. In both cases, the optical afterglow is established, declining steadily during the prompt emission. For GRB051111, there is evidence of an excess optical component during the prompt emission. The component is consistent with the flux spectrally extrapolated from the gamma-rays, using the gamma-ray spectral index. A compilation of spectral information from previous prompt detections shows that such a component is unusual. The existence of two prompt optical components - one connected to the high-energy emission, the other to separate afterglow flux, as indicated in GRB051111 - is not compatible with a simple ``external-external'' shock model for the GRB and its afterglow.

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

We summarize latest PIC simulation results on the radiation from Poynting jets and strongly magnetized collisionless shocks when a Poynting jet runs into cold ambient medium. We find that in all cases the radiative power output is much below that predicted by synchrotron radiation and the critical frequency is also much lower than the synchrotron critical frequency. We discuss the implications for the interpretation of GRB prompt emission data.

We systematically analyze the prompt emission and the early afterglow data of a sample of 31 GRBs detected by {\em Swift} before September 2005, and estimate the GRB radiative efficiency. BAT's narrow band inhibits a precise determination of the GRB spectral parameters, and we have developed a method to estimate these parameters with the hardness ratio information. The shallow decay component commonly existing in early X-ray afterglows, if interpreted as continuous energy injection in the external shock, suggests that the GRB efficiency previously derived from the late-time X-ray data were not reliable. We calculate two radiative efficiencies using the afterglow kinetic energy E_K derived at the putative deceleration time t_{dec}) and at the break time (t_b) when the energy injection phase ends, respectively. At t_b XRFs appear to be less efficient than normal GRBs. However, when we analyze the data at t_{dec} XRFs are found to be as efficient as GRBs. Short GRBs have similar radiative efficiencies to long GRBs despite of their different progenitors. Twenty-two bursts in the sample are identified to have the afterglow cooling frequency below the X-ray band. Assuming \epsilon_e = 0.1, we find \eta_\gamma(t_b) usually <10% and \eta_\gamma (t_{dec}) varying from a few percents to > 90%. Nine GRBs in the sample have the afterglow cooling frequency above the X-ray band for a very long time. This suggests a very small \epsilon_B and/or a very low ambient density n.

Very early observations with the Swift satellite of gamma-ray burst (GRB) afterglows reveal that the optical component is not detected in a large number of cases. This is in contrast to the bright optical flashes previously discovered in some GRBs (e.g. GRB 990123 and GRB 021211). Comparisons of the X-ray afterglow flux to the optical afterglow flux and prompt gamma-ray fluence is used to quantify the seemingly deficient optical, and in some cases X-ray, light at these early epochs. This comparison reveals that some of these bursts appear to have higher than normal gamma-ray efficiencies. We discuss possible mechanisms and their feasibility for explaining the apparent lack of early optical emission. The mechanisms considered include: foreground extinction, circumburst absorption, Ly-alpha blanketing and absorption due to high redshift, low density environments, rapid temporal decay, and intrinsic weakness of the reverse shock. Of these, foreground extinction, circumburst absorption, and high redshift provide the best explanations for most of the non-detections in our sample. There is tentative evidence of suppression of the strong reverse shock emission. This could be because of a Poynting-flux-dominated flow or a pure non-relativistic hydrodynamical reverse shock.

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

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

We review recent developments in particle acceleration by Poynting flux using plasma kinetic simulations, and discuss their potential applications to gamma-ray burst phenomenology

We report the Compton scattering emission from the Poynting flux acceleration of electron- positron plasma simulated by the 2-1/2 dimensional particle-in-cell(PIC) code. We show these and other remarkable properties of Poynting flux acceleration and Compton spectral output, and discuss the agreement with the observed spectra of GRBs and XRFs.

The X-ray light curves of many gamma-ray bursts (GRBs) observed by the Swift X-Ray Telescope (XRT) have a very steep-decay component (tail) following the prompt gamma-rays in the early phase and have some erratic flares occurring at a time from ~10^2 up to ~10^5 seconds. Based on the assumption that these tails and flares are of internal shock origin and that their decline behaviors are dominated by the curvature effect of the fireball, we present a self-consistency test for this scenario with a sample of 36 prompt-emission-tails/flare-tails in 22 GRB XRT light curves. Our results show that the t0 of the prompt emission tails and the tails of well-separated flares are usually at the rising segment of the last pulse of the prompt emission or the corresponding X-ray flare, being self-consistent with the expectation of the internal dissipation models for the prompt emission and X-ray flares. Our results indicate that each X-ray flare forms a distinct new episode of central engine activity and the GRB central engine remains active after the prompt emission is over, sometimes up to ~1 day after the GRB trigger (e.g. GRB 050502B & GRB 050724). This challenges the conventional central engine models and calls for new ideas to re-start the central engine. We further show that the on-set time of the late central engine activity does not depend on the GRB duration. We also identify a minority group of GRBs whose combined BAT-XRT light curves are smoothly connected without an abrupt transition between the prompt emission and the afterglow. These GRBs may have an external origin for both the prompt emission and the afterglow.

Several gamma-ray burst (GRB) luminosity indicators (LIs)have been proposed. Unlike SNe Ia, calibration of GRB LIs using a low-redshift sample is difficult. Based on the Bayesian theory, here we propose an approach to calibrate these LIs without introducing a low-redshift GRB sample. The essential points of our approach include, (1) calibrate the power-law indices in the LIs with a sample of GRBs in a narrow redshift range (Delta z); and (2) marginalize the coefficient of the LIs over a reasonable range of cosmological parameters. We take our newly discovered multivariable GRB LIs as an example and test the validity of our approach through simulations. We show that while the coefficient strongly depends on the cosmological parameters, the power-law indices do not as long as Delta z is small enough. The selection of Delta z for a particular GRB sample could be judged according to the size and the observational uncertainty of the sample. There is no preferable redshift to perform the calibration of the indices, while a lower redshift is preferable for the marginalization of the coefficient. The best strategy would be to collect GRBs within a narrow redshift bin around a fiducial intermediate redshift (e.g. z~1 or z~2), since the observed GRB redshift distribution is found to peak around z=(1 - 2.5). Our simulation suggests that with the current observational precision of measuring GRB isotropic energy, spectral break energy, and the optical temporal break time, 25 GRBs within a redshift bin of Delta z ~ 0.30 would give fine calibration to our luminosity indicator(abridged).

We review PIC simulation results of GRB emissions for both the Poynting flux and shock scenarios

We review recent PIC simulation results which show that double-sided irradiaton of a thin overdense plasma slab by ultra-intense laser pulses from both sides can lead to sustained comoving acceleration of surface electrons to energies much higher than the conventional ponderomotive limit. The acceleration stops only when the electrons drift transversely out of the laser beam. We show results of parameter studies based on this concept and discuss future laser experiments that can be used to test these computer results.

We present a detail analysis on the spectral lags of the short gamma-ray bursts (GRBs) and compare them with that of the long GRBs by using the CGRO/BATSE GRB Catalog. Our results are as follows. (1)The spectral lag distribution of the short GRBs is significantly different from that of the long GRBs. Excluding the statistical fluctuation effect, a proportion of ~17% of the short GRBs have a negative spectral lag, i.e., the hard photons being lag behind the soft photons. We do not find any peculiar features from their light curves to distinguish these bursts from those with a positive spectral lag. We argue that a more physical mechanism dominated the hard lag may be hid behind the morphological features of the light curves. This should be a great challenge to the current GRB models. We notice that this proportion is consistent with the proportion of short GRBs correlated with nearby galaxies newly discovered by Tanvir et al., although it is unclear if these short GRBs are indeed associated with the sources originated at low redshift. (2)While the spectral lags of the long GRBs are strongly correlated with the pulse durations, they are not for the short GRBs. However, the ratios of the spectral lag to the pulse duration for the short and long GRBs normally distributed at 0.023 and 0.046, respectively, with the sample width, indicating that the curvature effect alone could not explain the difference of the spectral lags between the two types of GRBs. The hydrodynamic timescales of the outflows and the radiative processes at work in GRBs might also play an important role as suggested by Daigne and Mochkovitch.

We present the results of a systematical analysis of the intrinsic optical afterglow light curves for a complete sample of gamma-ray bursts (GRBs) observed in the period from Feb. 1997 to Aug. 2005. These light curves are generally well-sampled, with at least four detections in the $R$ band. The redshifts of all the bursts in the sample are available. We derive the intrinsic $R$ band afterglow lightcurves (luminosity versus time within the cosmic proper rest frame) for these GRBs, and discover a fact that they essentially follow two universal tracks after 2 hours since the GRB triggers. The optical luminosities at 1 day show a clear bimodal distribution, peaking at 1.4*10^{46} ergs~s^{-1} for the luminous group and 5.3*10^{44} ergs~s^{-1} for the dim group. About 75% of the GRBs are in the luminous group, and the other 25% belong to the dim group. While the luminous group has a wide range of redshift distribution, the bursts in the dim group all appear at a redshift lower than 1.1.

We estimate the luminosity evolution and formation rate for over 900 GRBs by using redshift and luminosity data calculated by Band, Norris, $&$ Bonnell (2004) via the lag-luminosity correlation. By applying maximum likelihood techniques, we are able to infer the true distribution of the parent GRB population's luminosity function and density distributions in a way that accounts for detector selection effects. We find that after accounting for data truncation, there still exists a significant correlation between the average luminosity and redshift, indicating that distant GRBs are on average more luminous than nearby counterparts. This is consistent with previous studies showing strong source evolution and also recent observations of under luminous nearby GRBs. We find no evidence for beaming angle evolution in the current sample of GRBs with known redshift, suggesting that this increase in luminosity can not be due to an evolution of the collimation of gamma-ray emission. The resulting luminosity function is well fit with a single power law of index $L'^{-1.5}$, which is intermediate between the values predicted by the power-law and Gaussian structured jet models. We also find that the GRB comoving rate density rises steeply with a broad peak between $1<z<2$ followed by a steady decline above $z> 3$. This rate density qualitatively matches the current estimates of the cosmic star formation rate, favoring a short lived massive star progenitor model, or a binary model with a short delay between the formation of the compact object and the eventual merger.

We present a multi-variable regression analysis to several observable quantities for a sample of 15 gamma-ray bursts (GRBs). The observables used in the analysis includes the isotropic gamma-ray energy (E_{gamma,iso}), the peak energy of the nu F_nu spectrum in the rest frame (E_p'), and the rest frame break time of the optical afterglow light curves(t_b'). A strong dependence of E_{gamma, iso on E_p' and t_b' is derived. We also extend the analysis to the isotropic afterglow energies in the X-ray and the optical bands, respectively, and find that they are essentially not correlated with E_p' and t_b'. Regarding the E_{gamma, iso(E_p', t_b') relationship as a luminosity indicator, we explore the possible constraints on the cosmological parameters using the GRB sample. Our results indicate 0.13<omega_M<0.49 and 0.50<Omega_Lambda<0.85 for a flat Universe. The decelerating factor and its cosmological evolution, and the transition redshift are also investigated. Through Monte Carlo simulations, we find that the GRB sample satisfying our relationship observationally tends to be a soft and bright one, and that the contraints on the cosmological parameters can be much improved either by the enlargement of the sample size or by the increase of the observational precision. Furthermore, the detections of a few high redshift GRBs satisfying the correlation could greatly tighten the constraints (abridged).

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

Plasma outflows from gamma-ray bursts (GRB), pulsar winds, relativistic jets, and ultra-intense laser targets radiate high energy photons. However, radiation damping is ignored in conventional PIC simulations. In this letter, we study the radiation damping effect on particle acceleration via Poynting fluxes in two-and-half-dimensional particle-in-cell (PIC) plasma simulation of electron-positron plasmas. Radiation damping force is self-consistently calculated for each particle and reduces the net acceleration force. The emitted radiation is peaked within a few degrees from the direction of Poynting flux and strongly linear-polarized.

We study the radiation damping effect on the relativistic acceleration of electron-positron plasmas with two-and-half-dimensional particle-in-cell (PIC) simulation. Particles are accelerated by Poynting flux via the diamagnetic relativistic pulse accelerator (DRPA), and decelerated by the self-consistently solved radiation damping force. With $\Omega_{ce}/\omega_{pe}\geq 10$, the Lorentz factor of the highest energy particles reaches gamma>100, and the acceleration still continues. The emitted radiation is peaked within few degrees from the direction of Poynting flux and strongly linearly polarized, which may be detectable in gamma-ray burst(GRB) observations. We also show that the DRPA is insensitive to the initial supporting currents.

The radiation damping effect on the diamagnetic relativistic pulse accelerator (DRPA) is studied in two-and-half dimensional Particle-in-Cell (PIC) simulation with magnetized electron-positron plasmas. Self-consistently solved radiation damping force converts particle energy to radiation energy. The DRPA is still robust with radiation, and the Lorentz factor of the most high energy particles reach more than two thousand before they decouple from the electromagnetic pulse. Resulted emitted power from the pulse front is lower in the radiative case than the estimation from the non-radiative case due to the radiation damping. The emitted radiation is strongly linearly polarized and peaked within few degrees from the direction of Poynting flux.

Assuming that the observed gamma-ray burst (GRB) rate as a function of redshift is proportional to a corrected star formation rate, we derive the empirical distribution of the viewing angles of long BATSE GRBs, $P^{\rm em}(\theta)$, and the distribution of these bursts in the plane of $\theta$ against redshift, $P^{\rm em}(\theta, z)$, by using a tight correlation between $E_{\gamma}$) and $E_{\rm p}^{'}$). Our results show that $P^{\rm em}(\theta)$ is well fitted by a log-normal distribution centering at $\log \theta/{\rm rad}=-0.76$ with a width of $\sigma_{\log \theta}=0.57$. We test different universal structured jet models by comparing model predictions with our empirical results. To make the comparisons reasonable, an "effective" threshold, which corresponds to the sample selection criteria of the long GRB sample, is used. We find that the predictions of a two-Gaussian jet model are roughly consistent with our empirical results. A brief discussion shows that cosmological effect on the $E_{\gamma}-E_{\rm p}^{'}$ relation does not significantly affect our results, but sample selection effects on this relationship might significantly influence our results.

The narrowness of the distribution of the peak energy of $\nu F_{\nu}$ spectrum of gamma-ray bursts (GRBs) and the unification of GRB population are great puzzles yet to be solved. We investigate the two puzzles based on the global spectral behaviors of different GRB population in the $HR-E_{\rm{p}}$ plane (HR the spectral hardness ratio) with BATSE and HETE-2 observations. It is found that long GRBs and XRFs observed by HETE-2 seem to follow the same sequence in the $HR-E_{\rm{p}}$ plane, with the XRFs at the low end of this sequence. The long and short GRBs observed by BATSE follow significantly different sequences in the $HR-E_{\rm p}$ plane, with most of the short GRBs having a larger hardness ratio than the long GRBs at a given $E_{\rm{p}}$. These results indicate that the global spectral behaviors of the long GRB sample and the XRF sample are similar, while that of short GRBs is different. The short GRBs seem to be a unique subclass of GRBs, and they are not the higher energy extension of the long GRBs (abridged).

We present two-dimensional inviscid hydrodynamic simulations of a protoplanetary disk with an embedded planet, emphasizing the evolution of potential vorticity (the ratio of vorticity to density) and its dependence on numerical resolutions. By analyzing the structure of spiral shocks made by the planet, we show that progressive changes of the potential vorticity caused by spiral shocks ultimately lead to the excitation of a secondary instability. We also demonstrate that very high numerical resolution is required to both follow the potential vorticity changes and identify the location where the secondary instability is first excited. Low-resolution results are shown to give the wrong location. We establish the robustness of a secondary instability and its impact on the torque onto the planet. After the saturation of the instability, the disk shows large-scale non-axisymmetry, causing the torque on the planet to oscillate with large amplitude. The impact of the oscillating torque on the protoplanet's migration remains to be investigated.

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

We make an analysis of the distribution of the extensive Doppler redshift defined as $\widetilde{z}_{Dopp}=(z_{abs}-z_{em})/(1+z_{em})$ with a sample of 1317 absorption redshifts available from 401 quasars. The analysis reveals a bi-peak structure in the distribution, with one component located at $\widetilde{z}_{Dopp}\simeq 0.00$ and the other at $\widetilde{z}% _{Dopp}\simeq -0.01$. A study of some possible causes suggests that the structure can be well interpreted: while the absorbers inside the same cluster as the quasar concerned could contribute to the first component and the less populated space between clusters could explain the gap between the two peaks, the typical distance between clusters could account for the second component. If the bi-peak structure is true, which needs to be confirmed by complete samples, one would be able to obtain some useful cosmological information.

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

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

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

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

The diverse and complex light curves of gamma-ray bursts (GRBs) remain an outstanding astrophysical mystery. Here we report the results of 2-1/2-dimensional particle-in-cell (PIC) simulations of the relativistic expansion of magnetized electron-positron plasmas. When the simulation is carried to >150 light-crossing time of the initial plasma, the plasma pulse reproduces many of the GRB features. Remarkably, the plasma pulse bifurcates repeatedly, leading to a complex, multi-peak structure at late times, resembling GRB profiles.

Two-and-a-half-dimensional particle-in-cell plasma simulations are used to study the particle energization in expanding magnetized electron-positron plasmas with slab geometry. When the magnetized relativistic plasma with high temperature (initial electron and positron temperature are $k_{B}T_{e}=k_{B}T_{p}=5MeV$) is expanding into a vacuum, the electromagnetic (EM) pulse with large amplitude is formed and the surface plasma particles are efficiently accelerated in the forward direction owing to the energy conversion from the EM field to the plasma particles. We find that the behavior of the DRPA (Diamagnetic Relativistic Pulse Accelerator) depends strongly on the ratio of the electron plasma frequency to the cyclotron frequency $\omega_{pe}/\Omega_{e}$ and the initial plasma thickness. In the high $\omega_{pe}/\Omega_{e}$ case, the EM pulse is rapidly damped and the plasma diffuses uniformly without forming density peaks because the initial thermal energy of the plasma is much larger than the field energy. On the contrary, in the low $\omega_{pe}/\Omega_{e}$ case, the field energy becomes large enough to energize all the plasma particles, which are confined in the EM pulse and efficiently accelerated to ultrarelativistic energies. We also find that a thicker initial plasma increases the maximum energy of the accelerated particles.

Charge separation effects in the expansion of magnetized relativistic electron-ion plasmas into a vacuum are examined using 2-1/2-dimensional particle-in-cell plasma simulations. The electrostatic field at the plasma surface decelerates electrons and accelerates ions. A fraction of the surface electrons are trapped and accelerated by the pondermotive force of the propagating electromagnetic pulse, a mechanism we call the DRPA (diamagnetic relativistic pulse accelerator). This charge separation is enhanced as the initial plasma temperature is decreased. The overall energy gain of the plasma particles through the expansion strongly depends on the initial plasma temperature. Moreover, the electrons become relatively less energized and the ions more energized as the plasma temperature decreases.

We analyze the time profiles of individual gamma-ray burst (GRB) pulses, that are longer than 2 s, by modelling them with analytical functions that are based empirical descriptions of GRB spectral evolution. These analytical profiles are independent of the emission mechanism and can be used to model both the rise and decay profiles Using this method, we have studied a sample of 77 individual GRB pulses, allowing us to examine the fluence, pulse width, asymmetry, and rise and decay power-law distributions. We find that the rise phase is best modelled with a power law of average index $r = 1.31 \pm 0.11$ and that the average decay phase has an index o.f $d = 2.39 \pm 0.12$. We also find that the ratio between the rise and decay times (the pulse asymmetry) exhibited by the GRB pulse shape has an average value of 0.47 which varies little from pulse to pulse and is independent of pulse duration or intensity. We compare these parameters with those predicted to occur if individual pulse shapes are created purely by relativistic curvature effects in the context of the fireball model, a process that makes specific predictions about the shape of GRB pulses. The decay index distribution obtained from our sample shows that the average GRB pulse fades faster than the value predicted by curvature effects, with only 39% of our sample being consistent with the curvature model. We discuss several refinements of the relativistic curvature scenario that could naturally account for these observed deviations.

Using a 2-1/2-dimensional particle-in-cell (PIC) code to simulate the relativistic expansion of a magnetized collisionless plasma into a vacuum, we report a new mechanism in which the magnetic energy is efficiently converted into the directed kinetic energy of a small fraction of surface particles. We study this mechanism for both electron-positron and electron-ion (mi/me=100, me is the electron rest mass) plasmas. For the electron-positron case the pairs can be accelerated to ultra-relativistic energies. For electron-ion plasmas most of the energy gain goes to the ions.

The statistical shapes of the highest pulse have been studied by aligned method. A wavelet package analysis technique and a developed pulse--finding algorithm have been applied to select the highest pulse from burst profiles observed by BATSE on board CGRO from 1991 April 21 to 1999 January 26. The results of this work show that the statistical shapes of the highest pulses are related to energy: the higher the energy, the narrower the pulse. However, the characteristic structures of the pulses have nothing to do with energy, which strongly supports the previous conclusion that the temporal profiles in different channels are self--similar. The characteristic structures of the pulses can be well described by a model proposed by Norris et al. (1996). The fitting parameters are: $t_r$=0.12, $t_d$=0.16, $\upsilon =1.09$, the ratio of $t_r$ to $t_d$ for the pulse is 0.75. The result leads to our conjecture that the mechanisms of bursts in different gamma-ray bands might be the same. The shock, either an internal or an external one, producing the pulse, might emit photons over the four energy channels in the same way.

The observed delay in the arrival times between high and low energy photons in gamma-ray bursts (GRBs) has been shown by Norris et al. to be correlated to the absolute luminosity of a GRB. Despite the apparent importance of this spectral lag, there has yet to be a full explanation of its origin. We put forth that the lag is directly due to the evolution of the GRB spectra. In particular, as the energy at which the GRB's $\nu F_{\nu}$ spectra is a maximum ($E_{pk}$) decays through the four BATSE channels, the photon flux peak in each individual channel will inevitably be offset producing what we measure as lag. We test this hypothesis by measuring the rate of $E_{pk}$ decay ($\Phi_{o}$) for a sample of clean single peaked bursts with measured lag. We find a direct correlation between the decay timescale and the spectral lag, demonstrating the relationship between time delay of the low energy photons and the decay of $E_{pk}$. This implies that the luminosity of a GRB is directly related to the burst's rate of spectral evolution, which we believe begins to reveal the underlying physics behind the lag-luminosity correlation. We discuss several possible mechanisms that could cause the observed evolution and its connection to the luminosity of the burst.

Many classes of gamma-ray sources, such as gamma-ray bursts, blazars, Seyfert galaxies, and galactic black hole sources are surrounded by large amounts of gas and dust. X-rays and gamma-rays that traverse this material will be attenuated by Compton scattering and photoelectric absorption. One signature of an intervening scattering cloud is radiation-hardening by electrons that have been scattered and heated by the incident radiation, as illustrated by a Monte Carlo calculation. Compton scattering provides backscattered photons that will attenuate subsequent gamma rays through \gamma\gamma pair-production processes. We calculate the pair efficiency for a cloud illuminated by gamma-ray burst radiation. An analytic calculation of the flux of X-rays and gamma rays Thomson scattered by an intervening cloud is presented. Illuminated clouds near GRBs will form relativistic plasmas containing large numbers of electron-positron pairs that can be detected within ~1-2 days of the explosion before expanding and dissipating. Localized regions of pair annihilation radiation in the Galaxy could reveal gamma-ray sources embedded in dense clouds, or sites of past GRB explosions.

We discuss the space-and-time-dependent Monte Carlo code we have developed to simulate the relativistic radiation output from compact astrophysical objects, coupled to a Fokker-Planck code to determine the self-consistent lepton populations. We have applied this code to model the emission from a magnetized neutron star accretion shell near the Alfven radius, reprocessing the radiation from the neutron sar surface. We explore the parameter space defined by the accretion rate, stellar surface field and the level of wave turbulence in the shell. Our results are relevant to the emission from atoll sources, soft-X-ray transient X-ray binaries containing weakly magnetized neutron stars, and to recently suggested models of accretion-powered emission from anomalous X-ray pulsars.

Through studying the hard lags between the soft (3.3 - 5.8 keV) and hard (13.0 - 41.0 keV) photons of the 0.5 - 10 Hz QPOs in GRS 1915+105, we have classified them into three types: 0.5 - 2.0 Hz QPOs, 2.0 - 4.5 Hz QPOs, and 4.5 - 10 Hz QPOs. They are closely related to different temporal and spectral states. The first type of QPOs (0.5 - 2 Hz) have positive hard lags at both the QPO fundamental and first harmonic frequencies. These QPOs were observed in the quiescent soft state. The second type of QPOs (2 - 4.5 Hz), which were also detected in the quiescent soft state, have opposite signs of hard lags at the QPO fundamental and first harmonic frequencies. The third type (4.5 - 10 Hz), which showed up in medium soft quiescent/out-burst states, do not have significant higher harmonic peaks. There is a smooth transition between these three types of QPO behaviors. We did not detect 0.5 - 10 Hz QPOs in the very soft state. We discuss some of the implications of these results.

We study correlations between various hardness ratios of gamma-ray bursts (GRBs) and investigate if there are any differences between the two classes of the objects in the distributions of the ratios. The results suggest that, statistically, the slope of the higher part of the spectrum of the long duration bursts has nothing to do with that of the lower part; emissions at higher energy bands from the bursts of both short and long duration classes must be significantly different for different sources, while radiations at lower energy bands of the objects are relatively similar; the spectrum of the short duration bursts must be harder than that of the long duration bursts, confirming what the well-known hardness-duration correlation reveals; the profiles of the spectra between the long duration bursts must be more similar than that between the short duration bursts. The long duration bursts would share more common properties than the short duration bursts. A possible interpretation is proposed with the concept of the Doppler boosting in the relativistic beaming model in AGNs.

In this letter, the old issue of whether redshifts of quasars are of cosmological origin is investigated. We make a plot of absorption redshifts versus emission redshifts for quasars with large amounts of data. Our study shows that, almost all absorption redshifts are smaller than the corresponding emission redshifts. The relation between the absorption and emission redshifts predicted by current cosmological models is well obeyed. The result confirms that redshifts of quasars are indeed distance indicators. It might be the most obvious evidence found so far to be consistent with the cosmological interpretation of quasar redshifts.

We present a theoretical model for the alternating phase lags associated with QPO fundamental and harmonic frequencies observed in some Galactic black-hole candidates. We assume that the accretion flow exhibits a transition from an outer cool, optically thick accretion disk to a hot, inner advection-dominated accretion flow (ADAF), and that the QPOs are related to small-scale oscillations of the accretion rate and the transition radius. We present an analytical estimate of the expected phase lags at the fundamental and first harmonic frequencies of the QPOs and perform detailed time-dependent Monte-Carlo simualtions of the radiation transport in the oscillating ADAF / cool disk system. We find that this model is well suited to reproduce alternating phase lags between the fundamental and the first harmonic. It also naturally explains the trend observed in GRS 1915+105 that, as the soft X-ray luminosity increases, the QPO frequency increases and the phase lag associated with the QPO fundamental frequency changes sign from positive to negative. The relation between the disk temperature and the QPO frequency observed in GRS 1915+105 is consistent with a secular instability modulating the disk evaporation at the transition radius.

We report the results of our multi-wavelength observations of GRS 1758-258 made in August 1997. The energy bands include radio, millimeter, X-ray, and gamma-ray. The observations enable us to obtain a complete spectrum of the source over an energy range of 2 - 500 keV. The spectrum shows that GRS 1758-258 was in its hard state. It is well fitted by the Sunyaev-Titarchuk (ST) Compton scattering model. The spectrum is also fit by a power law with an exponential cutoff (PLE) plus a soft black-body component. The temperature of the soft component is about 1.2 keV, and the energy flux is less than 1.5% of the total X- and gamma-ray flux. The deduced hydrogen column density is in the range of (0.93 - 2.0) 10^{22} cm^{-2}. No significant iron lines are detected. The radio emission has a flat energy spectrum. The daily radio, X-ray and gamma-ray light curves show that GRS 1758-258 was stable during the observation period, but was highly variable on smaller time scales in X- and gamma-rays. The power density spectra are typical for the low-state, but we find the photon flux for the 5 to 10 keV band to be more variable than that in the other two energy bands (2 - 5 keV and 10 - 40 keV). Harmonically spaced quasi-periodic oscillations (QPOs) are observed in the power spectra. The phase lags between the hard photons and the soft photons have a flat distribution over a wide range of frequencies. A high coherence of about 1.0 (0.01 - 1 Hz) between the hard photons and the soft photons is also obtained in our observations. We compare these results with two variation models. Our millimeter observations did not reveal any conclusive signatures of an interaction between the jet from GRS 1758-258 and the molecular cloud that lies in the direction of GRS 1758-258.

Using the data from the Rossi X-ray Timing Explorer (RXTE), we report the different energy dependence of the variability of the four persistent hard X-ray sources in the low-hard state: Cygnus X-1, GX 339-4, GRS 1758-258 and 1E 1740.7-2942. Cygnus X-1 is found to have a flatter power density spectrum (PDS) shape at higher energies. The other three sources have energy independent PDS shapes. The energy dependence of the overall variability (the integrated rms amplitude) varies from source to source and from observation to observation. 1E~1740.7-2942, for example, has a variability generally increasing with energy while GX 339-4 has a decreasing variability. A general trend is found in the four sources that the integrated rms amplitude anti-correlates with the X-ray flux. We compare these distinct energy dependent behaviors with several emission models. None of the models can fully explain all the features that we have found.

The recent simultaneous detection of optical, X-ray and gamma-ray photons from GRB990123 during the burst provides the first broadband multi-wavelength characterization of the burst spectrum and evolution. Here we show that a direct correlation exists between the time-varying gamma-ray spectral shape and the prompt optical emission. This combined with the unique signatures of the time-resolved spectra of GRB990123 convincingly supports earlier predictions of the saturated Comptonization model. Contrary to other suggestions, we find that the entire continuum from optical to gamma-rays can be generated from a single source of leptons (electrons and pairs). The optical flux only appears to lag the gamma-ray flux due to the high initial Thomson depth of the plasma. Once the plasma has completely thinned out, the late time afterglow behavior of our model is the same as in standard models based on the Blandford-McKee (1976) solution.

GX 339-4 is an unusual black hole candidate in that it is a persistent source, being detected by X-ray telescopes most of the time, but it also has nova-like flaring states. In 1996 we performed radio, optical, X-ray, and gamma-ray observations of GX 339-4 when it was in a hard state (= soft X-ray low state). Here we present a brief summary of some of the results.

As part of our multiwavelength campaign of GX 339-4 observations in 1996 we present our radio, X-ray, and gamma-ray observations made in July, when the source was in a hard state (= soft X-ray low state). The radio observations were made at the time when there was a possible radio jet. We show that the radio spectrum was flat and significantly variable, and that the radio spectral shape and amplitude at this time were not anomalous for this source. Daily light curves from our pointed observation July 9-23 using OSSE, from BATSE, and from the ASM on RXTE also show that there was no significant change in the X- and gamma-ray flux or hardness during the time the possible radio jet-like feature was seen. The higher energy portion of our pointed RXTE observation made July 26 can be equally well fit using simple power law times exponential (PLE) and Sunyaev-Titarchuk (ST) functions. An additional soft component is required, as well as a broad emission feature centered on 6.4 keV. This may be an iron line that is broadened by orbital Doppler motions and/or scattering off a hot medium. Its equivalent width is 600 eV. Our simplistic continuum fitting does not require an extra reflection component. Both a PLE and a ST model also fit our OSSE spectrum on its own. Although the observations are not quite simultaneous, combining the RXTE and CGRO spectra we find that the PLE model easily fits the joint spectrum. However, the ST model drops off too rapidly with increasing energies to give an acceptable joint fit.

As part of our multiwavelength campaign of GX 339-4 observations in 1996 we present the rapid X-ray variability observed July 26 using the RXTE when the source was in a hard state (= soft X-ray low state). We found that the source was extremely variable, with many bright flares. The flares have relatively symmetric time profiles. There are a few time intervals where the flux rises steadily and then drops suddenly, sometimes to a level lower than the average before the increase. Hardness ratios showed that the source was slightly softer when the flux was brighter. The power density spectra (PDS) were also complicated and we found that broken power laws do not provide adequate fits to any of them. Instead a pair of zero-centered Lorentzians gives a good general description of the shape of the PDS. We found several quasi-periodic oscillations (QPO), including some that are harmonically spaced with the most stable frequency at 0.35 Hz. While the overall rms variability of the source was close to being constant throughout the observation (29% integrating between 0.01 and 50 Hz), there is a small but significant change in the PDS shape with time. More importantly, we show that the soft 2-5 keV band is more variable than the harder 5-10 and 10-40 keV bands, which is unusual for this source and for other black hole candidates. Cross correlation functions (CCF) between these bands show that the light curve for the 10-40 keV band lags that of the 2-5 keV band by 5 msec.

The time-dependent Comptonized output of a cool soft X-ray source drifting inward through an inhomogeneous hot inner disk or corona is numerically simulated. We propose that this scenario can explain from first principles the observed trends in the hard time lags and power spectra of the rapid aperiodic variability of the X-ray emission of Galactic black-hole candidates.

We have analyzed the evolution of the spectral hardness parameter Epk as a function of fluence in gamma-ray bursts. We fit 41 pulses within 26 bursts with the trend reported by Liang & Kargatis (1996) which found that Epk decays exponentially with respect to photon fluence. We also fit these pulses with a slight modification of this trend, where Epk decays linearly with energy fluence. In both cases, we found the set of 41 pulses to be consistent with the trend. For the latter trend, which we believe to be more physical, the distribution of the decay constant is roughly log-normal, with a mean of 1.75 +/- 0.07 and a FWHM of 1.0 +/- 0.1. Regarding an earlier reported invariance in the decay constant among different pulses in a single burst, we found probabilities of 0.49 to 0.84 (depending on the test used) that such invariance would occur by coincidence, most likely due to the narrow distribution of decay constant values among pulses.

If gamma-ray bursts are associated with dense star-forming regions in young galaxies, photoelectric absorption by the dense circumburster material (CBM) will occur. As the burst evolves, the surrounding material is photoionized, leading to fluorescence line emission and reduced photoelectric absorption opacity. We have analyzed this process in detail, accounting for the time-dependent photoelectric absorption, photoionization and fluorescence line emission from the CBM. We find that even if GRBs are hosted in dense star-forming regions, photoionization of the GRB environment leads to a constant, but very weak level of delayed fluorescence line emission on timescales of weeks to years after the burst. A temporally evolving iron~K edge absorption feature can serve as diagnostic tool to reveal the density structure of the CBM and may provide an opportunity for redshift measurements. We also investigated whether photoelectric absorption could be responsible for the spectral evolution of the low-energy slopes of some bright BATSE $\gamma$-ray bursts displaying extremely hard spectra below the peak energy, inconsistent with the optically-thin synchrotron shock model. We find that a very strong metal enrichment (~ 100 times solar-system abundances) in the gamma-ray burst environment and a rather peculiar spatial distribution of the CBM would be necessary in order to account for the observed hard spectra below a few 100 keV and their temporal evolution.

We present a detailed study of the observable effects of photoionization and Comptonization of line and continuum radiation from a cold accretion disk with a thin, warm, photoionized transition layer in the framework of self-consistent accretion-disk corona models for Galactic black-hole candidates. We use an iterative method to treat the non-linear radiation feedback between the transition layer and the hot corona numerically using a Monte-Carlo Comptonization code in combination with the photoionization and line transfer code XSTAR. The subset of the parameter space allowed in self-consistent accretion-disk corona systems on energetic grounds is explored, checking for the emergence of emission lines and/or absorption edges, with special emphasis on the spectral range $E \lesssim 1$ keV where such features might become observable with the advent of the AXAF satellite. Comparing our model calculations to the broadband spectrum of GX~339-4, we find good agreement with the observed spectral features and discuss how the future detection of the predicted features at lower energies can be used to constrain parameters.

We investigate the effect of inverse-Compton scattering of flares of soft radiation in different geometries of a hot, Comptonizing region and a colder accretion disk around a solar-mass black hole. The photon-energy dependent light curves, their Fourier transforms, power spectra and Fourier-period dependent time lags of hard photons with respect to softer photons are discussed. On the basis of a comparison with existing data we find arguments against Comptonization of external soft radiation as well as Comptonization in a homogeneous medium as dominant mechanisms for the rapid aperiodic variability in Galactic black-hole candidates. Possible further observational tests for the influence of Comptonization on the rapid aperiodic variability of Galactic black-hole candidates are suggested.

The spectral properties of a hybrid steady thermal-nonthermal isotropic plasma are studied. The involved emission mechanisms are thermal bremsstrahlung, nonthermal bremsstrahlung, cyclotron and nonthermal synchrotron. Formulas are derived to calculate the emissivities and absorption coefficients of all these processes. Calculations show that the nonthermal synchrotron and the thermal bremsstrahlung dominate over most of the frequency range and that the cyclotron and the nonthermal bremsstrahlung are important only at the two ends of the spectrum. The combination of the two dominant processes could generate the synchrotron self-absorbed spectrum, the Rayleigh-Jeans spectrum and an absorbed flatter spectrum which is seen in some solar flare radio spectra. Absorbed spectra steeper than the three types of spectra are rare. The conditions for these three possible spectral types are given. Cyclotron harmonic features are seen in the emergent spectra from the hybrid plasmas. These harmonic features, useful in determining the magnetic field, exist only under certain conditions which are explored in this paper.

The time-resolved spectra of gamma-ray bursts (GRBs) remain in conflict with many proposed models for these events. After proving that most of the bursts in our sample show evidence for spectral "shape-shifting", we discuss what restrictions that BATSE time-resolved burst spectra place on current models. We find that the synchrotron shock model does not allow for the steep low-energy spectral slope observed in many bursts, including GRB 970111. We also determine that saturated Comptonization with only Thomson thinning fails to explain the observed rise and fall of the low-energy spectral slope seen in GRB 970111 and other bursts. This implies that saturated Comptonization models must include some mechanism which can cause the Thomson depth to increase intially in pulses.

Recent studies have found that the spectral evolution of pulses within gamma-ray bursts (GRBs) is consistent with simple radiative cooling. Perhaps more interesting was a report that some bursts may have a single cooling rate for the multiple pulses that occur within it. We determine the probability that the observed "cooling rate invariance" is purely coincidental by sampling values from the observed distribution of cooling rates. We find a 0.1-26% probability that we would randomly observe a similar degree of invariance based on a variety of pulse selection methods and pulse comparison statistics. This probability is sufficiently high to warrant skepticism of any intrinsic invariance in the cooling rate.

In a hybrid thermal-nonthermal plasma, we find that the dominant emission and absorption mechanisms are synchrotron by nonthermal electrons and bremsstrahlung by thermal electrons. These two processes significantly change the spectrum from inverse Compton scatterings at low energies. We also find that Coulomb collisions are effective in cooling down the lower energy electrons but do not significantly alter the emission pattern. Compton cooling is more effective in changing emission and absorption coefficients when the photon energy density is high.

Using Monte Carlo methods to simulate the inverse Compton scattering of soft photons, we model the spectrum of the Galactic black hole candidate Cyg X-1, which shows evidence of a nonthermal tail extending beyond a few hundred keV. We assume an ad hoc sphere of leptons, whose energy distribution consists of a Maxwellian plus a high energy power-law tail, and inject 0.5 keV blackbody photons. The spectral data is used to constrain the nonthermal plasma fraction and the power-law index assuming a reasonable Maxwellian temperature and Thomson depth. A small but non-negligible fraction of nonthermal leptons is needed to explain the power-law tail.

We report evidence that the asymptotic low-energy power law slope alpha (below the spectral break) of BATSE gamma-ray burst photon spectra evolves with time rather than remaining constant. We find a high degree of positive correlation exists between the time-resolved spectral break energy E_pk and alpha. In samples of 18 "hard-to-soft" and 12 "tracking" pulses, evolution of alpha was found to correlate with that of the spectral break energy E_pk at the 99.7% and 98% confidence levels respectively. We also find that in the flux rise phase of "hard-to-soft" pulses, the mean value of alpha is often positive and in some bursts the maximum value of alpha is consistent with a value > +1. BATSE burst 3B 910927, for example, has a alpha_max equal to 1.6 +/- 0.3. These findings challenge GRB spectral models in which alpha must be negative of remain constant.

We consider stochastic particle acceleration in plasmas around stellar mass black holes to explain the emissions above 1 MeV from Galactic black hole candidates. We show that for certain parameter regimes, electrons can overcome Coulomb losses and be accelerated beyond the thermal distribution to form a new population, whose distribution is broad and usually not a power law; the peak energy of the distribution is determined by the balance between acceleration and cooling, with particles piling up around it. Radiation by inverse Compton scattering off the thermal (from background) and non-thermal (produced by acceleration) particles can in principle explain the hard X-ray to gamma-ray emissions from black hole candidates. We present model fits of Cyg X-1 and GRO J0422 in 50 keV -- 5 MeV region observed with OSSE and COMPTEL.

We study quantitatively the formation and radiation acceleration of electron-positron pair plasmoids produced by photon-photon collisions near Galactic black holes (GBHs). The terminal ejecta velocity is found to be completely determined by the total disk luminosity, proton loading factor and disk size, with no dependence on the initial velocity. We discuss potential applications to the recently discovered Galactic superluminal sources GRS1915, GROJ1655 and possibly other GBHs.