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 ach
romatic. 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.
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 typi
cal 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.
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 lig
htcurves 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 th
e 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_ppropto t_{p}^{-1.81+/-0.32} and L_ppropto 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 npropto 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.
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.20pm 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.
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.
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 sh
allow-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.11pm 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.15pm0.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.78pm 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.
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 signific
ant 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.
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 di
fferent 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}=2sim 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.
