Source author record

En-Wei Liang

En-Wei Liang appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

92works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

92 published item(s)

preprint2026arXiv

Identifying the poynting-flux-dominated outflow of Fermi GRBs with non-thermal spectrum and its energy-resolved light curve fitting

The jet compositions of gamma-ray bursts (GRBs) are very important to understand the energy dissipation and radiation mechanisms, but it remains an open question in GRB physics. In this paper, we present a systematic analysis to search for 88 bright GRBs that include a total of 129 pulses observed by Fermi/GBM with redshift measured, and extract the spectra of each pulse with Band function (Band), cutoff power-law (CPL), blackbody (BB), non-dissipative photospheric (NDP), Band+BB, as well as CPL+BB. We find that 80 pulses, 35 pulses, and 14 pulses present purely non-thermal, hybrid, and thermal spectra, respectively. By focusing on those 80 pulses with purely non-thermal spectra, one can estimate the lower limits of magnetization factor ($σ$) via suppressing the pseudo-thermal component. It is found that 30 pulses in 21 GRBs are the lower limit of $σ>5$ at the photosphere by adopting $R_{0}=10^{10}$ cm. It suggests that at least the outflow of those GRB jets with high $σ$ is dominated by Poynting-flux. On the other hand, we also perform the light curve fitting with a fast-rise-exponential-decay (FRED) model for 15 bright GRBs with a high magnetization factor in our sample, and find that a correlation between pulse width ($w$) and energy of 13 GRBs really exists in their energy-resolved light curves. It is also a piece of independent evidence for those GRBs with a high value $σ$ to support the origin of the Poynting flux outflow.

preprint2025arXiv

Evidence for the transition from thermal to non-thermal emission in the prompt emission of GRB 161117A

GRB 161117A is a long-duration GRB with three main overlapping peaks. By analyzing the time-resolved spectra of its data observed with the Gamma-Ray Burst Monitor (GBM) on board the Fermi mission, we find that the spectral evolution shows a transition from thermal (single BB) to hybrid (PL$+$BB), and finally to non-thermal (Band and CPL) emissions. Such a transition suggests that the jet composition of GRB 161117A should be changed from a fireball to a Poynting-flux-dominated jet. The bulk Lorentz factor ($Γ_{\rm ph}$), radii ($R_{\rm ph}$ and $R_{0}$), magnetization factor at the central engine ($σ_0$), and dimensionless entropy ($η$) of the outflow can be inferred by invoking the observed quasi-thermal component within two models (e.g., pure fireball and hybrid). It is found that $Γ_{\rm ph}$ seems to be tracking with the light curve, and $R_{0}$ remains a constant at $\sim$ $10^{8}$ cm. The low magnetization ($1+σ_0 \sim$ 1) and high dimensionless entropy ($η\gg$ 1) during the first seven time-intervals suggest to be a pure fireball outflow. Moreover, we also estimate the lower limit of magnetization parameter at the photosphere radius ($σ_{\rm ph}\sim 1.4$ and 0.75) for late phase via the non-thermal spectra, and it indicates that the particle acceleration mechanism is dominated by internal shocks rather than magnetic dissipation processes. Finally, the $ν\barν$ annihilation mechanism of NDAF model to explain the thermal emission of GRB 161117A is also discussed.

preprint2025arXiv

Ultrahigh-Energy Gamma-ray Emission Associated with Black Hole-Jet Systems

Black holes (BH), one of the most intriguing objects in the universe, can manifest themselves through electromagnetic radiation initiated by the accretion flow. Some stellar-mass BHs drive relativistic jets when accreting matter from their companion stars, forming microquasars. Non-thermal emission from the radio to tera-electronvolt (TeV) gamma-ray band has been observed from microquasars, indicating the acceleration of relativistic particles. Here we report detection of four microquasars (SS 433, V4641 Sgr, GRS 1915+105, MAXI J1820+070) of spectrum extending to the ultrahigh-energy (UHE; photon energy $E>100$ TeV) band and one microquasar (Cygnus X-1) of spectrum approaching 100 TeV, using the Large High Altitude Air Shower Observatory (LHAASO). Notably, the total emission associated with SS 433 cannot be interpreted with a single leptonic component. In the UHE band, its emission is in spatial coincidence with a giant atomic cloud, which is consistent with a hadronic origin. An elongated source is discovered from V4641 Sgr with the spectrum continuing up to 800 TeV. The detection of UHE gamma rays demonstrates that accreting BHs and their environments can operate as extremely efficient accelerators of particles out of 1 peta-electronvolt (PeV), suggesting microquasars to be important contributors to Galactic cosmic rays especially around the `knee' region.

preprint2022arXiv

A Channel to Form Fast-spinning Black Hole--Neutron Star Binary Mergers as Multi-messenger Sources

After the successful detection of a gravitational-wave (GW) signal and its associated electromagnetic (EM) counterparts from GW170817, neutron star--black hole (NSBH) mergers have been highly expected to be the next type of multi-messenger source. However, despite the detection of several of NSBH merger candidates during the GW third observation run, no confirmed EM counterparts from these sources have been identified. The most plausible explanation is that these NSBH merger candidates were plunging events mainly because the primary BHs had near-zero projected aligned-spins based on GW observations. In view that NSs can be easily tidally disrupted by BHs with high projected aligned-spins, we study an evolution channel to form NSBH binaries with fast-spinning BHs, the properties of BH mass and spin, and their associated tidal disruption probability. We find that if the NSs are born firstly, the companion helium stars would be tidally spun up efficiently, and would thus finally form fast-spinning BHs. If BHs do not receive significant natal kicks at birth, these NSBH binaries that can merge within the Hubble time would have BHs with the projected aligned-spins $χ_{z}\gtrsim0.8$ and, hence, can certainly allow tidal disruption to happen. Even if significant BH kicks are considered for a small fraction of NSBH binaries, the projected aligned-spins of BHs are $χ_z\gtrsim0.2$. These systems can still be disrupted events unless the NSs are very massive. Thus, NS-first-born NSBH mergers would be promising multi-messenger sources. We discuss various potential EM counterparts associated with these systems and their detectability in the upcoming fourth observation run.

preprint2022arXiv

Comparison of thermodynamic behaviors of two regular-AdS black holes

Considering the negative cosmological constant of an anti-de Sitter (AdS) background as a positive thermodynamic pressure in the extended phase space, we study the thermodynamic behaviors of the regular Hayward-AdS (HAdS) black hole (BH) by investigating $P-\upsilon$ critical and cooling-heating phase transition. We compare the difference of two thermodynamic processes between the HAdS BH and the Bardeen-AdS (BAdS) BH. It is found that the phase transition of the BAdS BH tends to be more the van der Waals (vdW) phase transition. For the cooling-heating phase transition, we find the inversion curves of the HAdS BH are always higher than the BAdS BH under the same pressure and magnetic charge. We also compare the smallest existence mass, the zero-temperature remnant, and the critical magnetic charge. The results suggest that the inner horizon and the outer horizon of the Hayward BH are easier to merge, and the singularity is easier to expose.

preprint2022arXiv

Constraints on ultracompact minihalos from the extragalactic gamma-ray background observation

Ultracompact minihalo (UCMH) is a special type of dark matter halo with a very steep density profile which may form in the early universe seeded by an overdense region or a primordial black hole. Constraints on its abundance give valuable information on the power spectrum of primordial perturbation. In this work, we update the constraints on the UCMH abundance in the universe using the extragalactic gamma-ray background (EGB) observation. Comparing to previous works, we adopt the updated Fermi-LAT EGB measurement and derive constraints based on a full consideration of the astrophysical contributions. With these improvements, we place constraints on UCMH abundance 1-2 orders of magnitude better than previous results. With the background components considered, we can also attempt to search for possible additional components beyond the known astrophysical contributions.

preprint2022arXiv

Early Evolution of a Newborn Magnetar with Strong Precession Motion in GRB 180620A

The observed early X-ray plateau in the afterglow lightcurves of some gamma-ray bursts (GRBs) is attributed to the dipole radiations (DRs) of a newborn magnetar. A quasi-periodic oscillation (QPO) signal in the plateau would be strong evidence of the magnetar precession motion. By making a time-frequency domain analysis for the X-ray afterglow lightcurve of GRB 180620A, we find a QPO signal of $\sim650$ seconds in its early X-ray plateau. We fit the lightcurve with a magnetar precession model by adopting the Markov chain Monte Carlo algorithm. The observed lightcurve and the QPO signal are well represented with our model. The derived magnetic field strength of the magnetar is $B_{\rm p}= (1.02^{+0.59}_{-0.61})\times10^{15}$~G. It rapidly spins down with angular velocity evolving as $Ω_{s} \propto(1+t/τ_{\rm sd})^{-0.96}$, where $τ_{\rm sd}=9430$~s. Its precession velocity evolution is even faster than $Ω_s$, i.e. $Ω_{ p}\propto (1+t/τ_{p})^{-2.18\pm0.11}$, where $τ_{p}=2239\pm206$~s. The inferred braking index is $n=2.04$. We argue that the extra energy loss via the magnetospheric processes results in its rapid spin-down, a low braking index of the magnetar, and the strong precession motion.

preprint2022arXiv

Ehrenfest's scheme and microstructure for regular-AdS black hole in the extended phase space

The regular (Bardeen)-AdS (BAdS) black hole (BH) in the extended phase space is taken as an example for investigating the BH phase transition grade from both macroscopic and microscopic points of view. The equation of state and thermodynamic quantities of this BH are obtained. It is found that the BAdS BH phase space in the extended phase space should be a second-order phase transition near the critical point by verifying the Ehrenfest's equation, and the possibility of its first-order phase transition can be ruled out by the entropy continuity and the heat capacity mutation. The critical exponents from the microscopic structure are analytically and numerically presented with the Landau continuous phase transition theory by introducing a microscopic order-parameter.

preprint2022arXiv

Fermi Observations of GRB 220426A: a burst similar to GRB 090902B

We report on a very bright, long-duration gamma-ray burst (GRB), GRB~220426A, observed by \emph{Fermi} satellite. GRB~220426A with total duration of $T_{90}=6$~s is composed with two main pulses and some sub-peaks. The spectral analysis of this burst with Band function reveals that both the time-integrated and the time-resolved spectra are very narrow with high $α\gtrsim 0.2$ and low $β\lesssim -3.1$. It is strong reminiscent of GRB~090902B, a special GRB with identification of the photospheric emission. Then, we perform the spectral analysis of this burst based on a non-dissipated photospheric emission, which can be well modelled as the multicolor-blackbody with a cutoff power-law distribution of the thermal temperature. The spectral fittings reveal that the photospheric emission can well describe the radiation spectrum of this burst. We conclude that this burst would be a second burst in the class of GRB~090902B observed by \emph{Fermi} satellite. We also discuss the physics of photosphere and the origin of the high-energy component in GRB~220426A .

preprint2022arXiv

GeV $γ$-ray Emission of Compact Steep-Spectrum Source 4C +39.23B

Thirteen yr observation data of 4FGL J0824.9+3915 with the Large Area Telescope on board the Fermi Gamma Ray Space Telescope (Fermi/LAT) are analyzed for revisiting whether 4C +39.23B, a compact steep-spectrum (CSS) source closed to a flat-spectrum radio quasar (FSRQ) 4C +39.23A in the $γ$-ray emitting region of 4FGL J0824.9+3915, is a $γ$-ray emitter. We find that the time-integrated $γ$-ray emission of 4FGL J0824.9+3915 is overwhelmingly dominated by 4C +39.23A. It shows significant variability at a 6.7$σ$ confidence level and the average $γ$-ray flux in the 0.1--300 GeV energy band is $(1.60\pm 0.15)\times10^{-8}$ ph cm$^{-2}$ s$^{-1}$ with a power-law photon spectral index of $2.48\pm0.05$. During MJD 57500--58500, 4FGL J0824.9+3915 is in a low state with a steady $γ$-ray flux. Analyzed the Fermi/LAT observation data in this time interval, it is found that the TS values of the $γ$-ray emission from 4C +39.23A and 4C +39.23B are $\sim5$ and $\sim 31$, respectively, indicating that the $γ$-ray emission in this time interval is dominated by the CSS 4C +39.23B. The derived average flux in this time interval for 4C +39.23B is $(9.40\pm4.10)\times 10^{-9}$ ph cm $^{-2}$ s$^{-1}$ with $Γ_γ=2.45\pm0.17$. Attributing the spectral energy distribution (SED) of 4C +39.23B to the radiations from its core and extended region, we show that the SED can be represented with a two-zone leptonic model. Its $γ$-ray emission is contributed by the core region. The derived magnetic field strength and Doppler boosting factor of the core are 0.13 G and 6.5. Comparing 4C +39.23B with other $γ$-emitting CSSs and compact symmetric objects (CSOs) in the $Γ_γ-L_γ$ plane, it resembles CSSs.

preprint2022arXiv

Gravitational-wave evolution of newborn magnetars with different deformed structure

Weak and continuous gravitational-wave (GW) radiation can be produced by newborn magnetars with deformed structure and is expected to be detected by the Einstein telescope in the near future. In this work we assume that the deformed structure of a nascent magnetar is not caused by a single mechanism but by multiple time-varying quadrupole moments such as those present in magnetically induced deformation, starquake-induced ellipticity, and accretion column-induced deformation. The magnetar loses its angular momentum through accretion, magnetic dipole radiation, and GW radiation. Within this scenario, we calculate the evolution of GWs from a newborn magnetar by considering the above three deformations. We find that the GW evolution depends on the physical parameters of the magnetar (e.g., period and surface magnetic field), the adiabatic index, and the fraction of poloidal magnetic energy to the total magnetic energy. In general the GW radiation from a magnetically induced deformation is dominant if the surface magnetic field of the magnetar is large, but the GW radiation from magnetar starquakes is more efficient when there is a larger adiabatic index if all other magnetar parameters remain the same. We also find that the GW radiation is not very sensitive to different magnetar equations of state.

preprint2022arXiv

GRB 211227A as a peculiar long gamma-ray burst from compact star merger

Long-duration gamma-ray bursts (GRBs) associated with supernovae (SNe) are believed to originate from massive star core-collapse events, whereas short-duration GRBs that are related to compact star mergers are expected to be accompanied by kilonovae. GRB 211227A, which lasted about 84 s, had an initial short/hard spike followed by a series of soft gamma-ray extended emission at redshift $z=$0.228. We performed follow-up observations of the optical emission using BOOTES, LCOGT, and the Lijiang 2.4m telescope, but we detected no associated supernova signature, even down to very stringent limits at such a low redshift. We observed the host galaxy within a large error-circle and roughly estimate the physical offset of GRB 211227A as $20.47\pm14.47$ kpc from the galaxy center. These properties are similar to those of GRB 060614, and suggest that the progenitor of GRB 211227A is not favored to be associated with the death of massive stars. Hence, we propose that GRB 211227A originates from a compact star merger. Calculating pseudo-kilonova emission for this case by adopting the typical parameters, we find that any associated pseudo-kilonova is too faint to be detected. If this is the case, it explains naturally the characteristics of the prompt emission, the lack of SN and kilonova emission, and the large physical offset from the galaxy center.

preprint2022arXiv

Influence of Accretion Flow and Magnetic Charge on the Observed Shadows and Rings of the Hayward Black Hole

The feature of the observed shadows and rings of an astrophysical black hole (BH) may depend on its accretion flows and magnetic charge. We find that the shadow radii and critical impact parameters of the Hayward BH are decreased with the increase of the magnetic charge. Comparing the Schwarzschild BH with the Hayward BH using the ray-tracing method, we show that the density and deflection of lights increase with the magnetic charge, and the BH singularity does not affect the generation of the shadow. Based on three optically thin accretion flow models, the two-dimensional shadows in celestial coordinates are derived. It is found that the shadow and photon ring luminosities of a Hayward BH surrounded by infalling spherical accretion flow are dimmer than that of a static spherical accretion flow. Taking three kinds of inner radii at which the accretion flow stops radiating, we find that the observed luminosity of a Hayward BH surrounded by a thin disk accretion flow is dominated by the direct emission, and the photon ring emission has a weak influence on it. These results suggest that the size of the observed shadow is related to the space-time geometry, and the luminosities of both the shadows and rings are affected by the accretion flow property and the BH magnetic charge.

preprint2022arXiv

Modeling the Multi-band Light Curves of the Afterglows of Three Gamma-Ray Bursts and Their Associated Supernovae

There are some dozen supernovae (SNe) associated with long Gamma-ray bursts (GRBs) have been confirmed. Most of previous studies derive the physical properties of the GRB-SNe by fitting the constructed (psuedo-)bolometric light curves. However, many GRB-SNe have only a few filter data, for which the (psuedo-)bolometric light curves are very difficult to be constructed. Additionally, constructing (psuedo-)bolometric light curves rely on some assumptions. In this paper, we use the multi-band broken power-law plus $^{56}$Ni model to fit the multi-band light curves of the afterglows and the SNe (SN~2001ke, SN~2013dx, and SN~2016jca) associated with three GRBs (GRB~011121, GRB~130702A, and GRB~161219B). We find our model can account for the multi-band light curves of the three GRB-SNe (except for the late-time $z-$band light curves of two events), indicating that the model is a reliable model. The $^{56}$Ni masses we derive are higher than that in the literature. This might be due to the fact that the $^{56}$Ni masses in the literature are usually obtained by fitting the psuedo-bolometric light curves whose luminosities are usually (significantly) underestimated. We suggest that the multi-band model can not only be used to fit the multi-band light curves of GRB-SNe that have many filter observations, but also fit those having sparse data.

preprint2022arXiv

Observable characteristics of the charged black hole surrounded by thin disk accretion in Rastall gravity

The observable characteristics of the charged black hole (BH) surrounded by a thin disk accretion are investigated in the Rastall gravity. We found that the radii of the direct emission, lensing ring, and photon ring dramatically increased as the radiation field parameter increases, but they only weakly depend on the BH charge. Three positions of the radiation accretion disk relative to the BH are considered, i.e., the innermost accretion disk is closed to the radii of the innermost stable circular orbit, the photon ring of the BH, and the event horizon of the BH. The observed images in three cases respectively are obtained. It is found that the total observed flux is dominated by the direct emission, the lensing ring provides a small contribution, and the photon ring is negligible. The lensing and photon rings could not be observed in the blurred image with the EHT resolution. Our results suggest that the observable characteristics of the charged BH surrounded by the thin disk accretion in the Rastall gravity depend on both the BH space-time structure and the position of the radiating accretion disk with respect to the BH. The research of these BH images may serve as a probe for the BH-disk structure in M87$^{*}$ like nearby active galactic nuclei.

preprint2022arXiv

Probing the progenitor of high-$z$ short-duration GRB 201221D and its possible bulk acceleration in prompt emission

The growing observed evidence shows that the long- and short-duration gamma-ray bursts (GRBs) originate from massive star core-collapse and the merger of compact stars, respectively. GRB 201221D is a short-duration GRB lasting $\sim 0.1$ s without extended emission (EE) at high redshift $z=1.046$. By analyzing data observed with the Swift/BAT and Fermi/GBM, we find that a cutoff power-law model can adequately fit the spectrum with a soft $E_{\rm p}=113^{+9}_{-7}$ keV, and isotropic energy $E_{γ,iso} =1.36^{+0.17}_{-0.14}\times 10^{51}~\rm erg$. In order to reveal the possible physical origin of GRB 201221D, we adopted multi-wavelength criteria (e.g., Amati relation, $\varepsilon$-parameter, amplitude parameter, local event rate density, luminosity function, and properties of the host galaxy), and find that most of the observations of GRB 201221D favor a compact star merger origin. Moreover, we find that $\hatα$ is larger than $2+\hatβ$ in the prompt emission phase which suggests that the emission region is possibly undergoing acceleration during the prompt emission phase with a Poynting-flux-dominated jet.

preprint2022arXiv

QED and accretion flow models effect on optical appearance of Euler-Heisenberg black holes

Taking the quantum electrodynamics (QED) effect into account, we investigate the geometrical-optics appearance of the Euler-Heisenberg (EH) black hole (BH) under the different accretion flows context, which depends on the BH space-time structure and different sources of light. The more significant magnetic charge leads to the smaller shadow radius for the EH BH, while the different values of the EH parameter do not ruin it. Different features of the corresponding two-dimensional shadow images are derived for the three optically thin accretion flow models. It is shown that the total observed intensity in the static spherical accretion flow scenario leads than that of the infalling spherical accretion flow under same parameters, but the size and position of the EH BH shadows do not change in both of these accretions flows, implying that the BH shadow size depends on the geometric space-time and the shadows luminosities rely on the accretion flow morphology. Of particular interest is that a thin disk accretion model illuminated the BH, we found that the contribution of the lensing ring to the total observed flux is less than 5\%, and the photon ring is less than 2\%, indicating that the direct emission dominates the optical appearance of the EH BH. It is also believed that the optical appearance of the BH image depends on the accretion disk radiation position in this scenario, which can serve as a probe for the disk structure around the active galactic nucleus (AGN) of M87^{*} like.

preprint2022arXiv

Radioactively-Powered Gamma-Ray Transient Associated with a Kilonova from Neutron Star Merger

AAssociation of GW170817/GRB170817A/AT2017gfo provides the first direct evidence for neutron star mergers as significant sources of $r$-process nucleosynthesis. A gamma-ray transient (GRT) would be powered by the radioactive decay of the freshly-synthesized $r$-process elements. By analyzing the composition and gamma-ray opacity of the kilonova ejecta in details, we calculate the lightcurve and spectrum of the GRT for a range of spherically symmetric merger ejecta models with mass $M_{\rm ej}=0.001 \sim 0.05M_{\odot}$ and expansion velocity $v_{\rm ej}= 0.1\sim 0.4c$. It is found that the peak of the GRT lightcurve depends on $M_{\rm ej}$ and $v_{\rm ej}$ as $t_{\rm pk} \approx 0.5~{\rm days} ~ (M_{\rm ej}/0.01M_{\odot})^{1/2}(v_{\rm ej}/0.1c)^{-1}$ and $L_{\rm pk} \approx 2.0\times10^{41} ~{\rm erg~s} ^{-1} (M_{\rm ej}/0.01M_{\odot})^{1/2}(v_{\rm ej}/0.1c)$. Most radiating photons are in the $100-3000$ keV band and the spectrum peaks at about 800~keV for different nuclear physics inputs. The line features are blurred out by the Doppler broadening effect. Adopting the ejecta parameters reported in literature, we examine the detection probability of the possible GRT associated with AT2017gfo. We show that the GRT cannot be convincingly detected neither with current nor with the proposed missions in the MeV band, such as ETCC and AMEGO. The low gamma-ray flux, together with the extremely low event rate at local universe, makes a discovery of GRTs a great challenge.

preprint2022arXiv

Search for gamma-ray line signals around the black hole at the galactic center with DAMPE observation

The adiabatic growth of a black hole (BH) may enhance the dark matter (DM) density surrounding it, causing a spike in the DM density profile. The spike around the supermassive BH at the center of the Milky Way may lead to a dramatic enhancement of the gamma-ray flux of DM annihilation from the galactic center (GC). In this work, we analyze the gamma-ray data of the innermost region (i.e., the inner 1$^\circ$) of the GC to search for potential line-like signals from the BH spike. Such line-like signals could be generated in the process of DM particles annihilating into double photons. We adopt the gamma-ray data from the Dark Matter Particle Explorer (DAMPE). Although the DAMPE has a much smaller effective area than the Fermi-LAT, the gamma-ray line search can benefit from its unprecedented high energy resolution. No significant line-like signals are found in our analysis. We derive upper limits on the cross section of the annihilation based on this non-detection. We find that despite the DAMPE's small effective area for photon detection, we can still place strong constraints on the cross section ($\left<σv\right>\lesssim10^{-27}\,{\rm cm^3\,s^{-1}}$) in the spike scenario due to the very bright model-expected flux from the spike. Our results indicate that either DM does not annihilate primarily through the $γγ$ channel in the mass range we considered or no sharp density spike is present at the GC.

preprint2022arXiv

TeV and keV-MeV Excesses as Probes for Hadronic Process in BL Lacertaes

A hard TeV $γ$-ray component excess over the single-zone leptonic model prediction (TeV excess) is observed in the spectral energy distributions (SEDs) of some BL Lacs. Its origin is uncertain. We revisit this issue with four BL Lacs (1ES 0229+200, 1ES 0347--121, 1ES 1101--232, and H2356--309), in which the TeV excess is detected in their intrinsic SEDs. We represent their SEDs with a single-zone leptohadronic model, where radiations of the electrons and protons as well as the cascade electrons produced by the $γγ$ and p$γ$ interactions within their jets are considered. We show that the observed SEDs below the GeV gamma-ray band are attributed to the synchrotron radiations and self-Compton process of the primary electrons, and the TeV excess is explained with the $γ$-ray emission from the p$γ$ process via the $π^{0}$ decay. The cascade emission of the electrons produced via the $γγ$ and p$γ$ interactions results in a keV-MeV excess in the SEDs, illustrated as a bump or plateau. This extra photon field enhances the production of TeV photons from the $pγ$ process, resulting in a reduction of the proton power by about one order of magnitude. However, the derived powers are still 3--4 orders of magnitude larger than the Eddington limit, being challenged by the current black hole accretion physics. Applying our model to Mrk 421, we propose that synergic observations with current and upcoming TeV and keV-MeV telescopes for its tentative TeV and MeV excesses can give insights to the hadronic process in its jet.

preprint2022arXiv

Using the Optical--NIR Spectral Energy Distributions To Search for the Evidence of Dust Formation of 66 Supernovae

In this paper, we searched for the dust formation evidence of 66 supernovae (SNe) by using the blackbody model and the blackbody plus dust {emission} model to fit their early$-$time optical$-$near infrared (NIR) spectral energy distributions (SEDs). We find that, while the blackbody model can fit most SEDs of the SNe in our sample, the model cannot fit the SEDs of some SNe, in which the SEDs of 2 SNe (SNe~2010bq and 2012ca) show NIR excesses which can be attributed to the emission from the heated dust. We use blackbody plus dust emission model to fit the SEDs showing NIR excesses, finding that both graphite and silicate dust models can fit the SEDs, and the graphite model get reasonable temperatures or better fits. Assuming that the dust is graphite, the best-fitting temperatures (masses) of the dust of the SNe~2010bq and 2012ca are $\sim 1300-1800$ K ($\sim 0.1-3.4 \times 10^{-4}$ M$_\odot$) and $\sim 600-1000$ K ($\sim 0.6-7.5 \times 10^{-3}$ M$_\odot$), respectively. We compare the vaporization radii and the blackbody radii of the dust shells of the 2 SNe with the upper limits of the ejecta radii of the SNe at the first epochs, and demonstrate that the NIR excesses of the SEDs of the 2 SNe might be caused by the pre-existing dust.

preprint2020arXiv

Constraining the nuclear equation of state via gravitational-wave radiation of short gamma-ray burst remnants

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$σ$ 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 $α=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.

preprint2020arXiv

Diversity of gamma-ray and Radio Variabilities of Bright Blazars and Implications for gamma-ray Emission Location

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.

preprint2020arXiv

Early Optical Observations of GRB 150910A: Bright Jet Optical Afterglow and X-ray Dipole Radiation from a Magnetar Central Engine

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 ($α_{\rm R,2} \approx -1.36$), as predicted in the standard external shock model, while the X-ray emission exhibits a plateau ($α_{\rm X,1} \approx -0.36$) followed by a steep decay ($α_{\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 $η_γ \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 $η_w \geq 32\%$.

preprint2020arXiv

Evidence for gravitational-wave dominated emission in the central engine of short GRB 200219A

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 model can adequately 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 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 a millisecond magnetar, which loses its rotation energy via GW quadrupole radiation. Then, the abrupt drop decay is caused by the magnetar collapsing into a black hole before switching to EM-dominated emission. This is the first short GRB for which the X-ray emission has such an intriguing feature powered by a magnetar via GW-dominated radiation. If this is the case, one can estimate the physical parameters of a magnetar, the GW signal powered by a magnetar and the merger-nova emission are also discussed.

preprint2020arXiv

Gamma-ray Burst Spectrum with a Time-dependent Injection Rate of High-energy Electrons

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.

preprint2020arXiv

GRB 140423A: A Case of Stellar Wind to Interstellar Medium Transition in the Afterglow

We present very early ground-based optical follow-up observations of GRB~140423A, which was discovered by \emph{Swift}/BAT and by {\it Fermi}/GBM. Its broadband afterglow was monitored by {\it Swift}/XRT and ground-based optical telescopes from $T_0+$70.96~s to 4.8~d after the {\it Swift}/BAT trigger. This is one more case of prompt optical emission observation. The temporal and spectral joint fit of the multiwavelength light curves of GRB 140423A reveals that achromatic behavior is consistent with the external shock model including a transition from a stellar wind to the interstellar medium (ISM) and energy injection. In terms of the optical light curves, there is an onset bump in the early afterglow with a rising index $α_{\rm O,I} = -0.59 \pm 0.04$ (peaking at $t_{\rm peak}-T_0 \approx 206$~s). It then decays with a steep index $α_{\rm O,II} = 1.78 \pm 0.03$, and shows a steeper to flatter "transition" with $α_{\rm O,III} = 1.13 \pm 0.03$ at around $T_0 + 5000$~s. The observed X-ray afterglow reflects an achromatic behavior, as does the optical light curve. There is no obvious evolution of the spectral energy distribution between the X-ray and optical afterglow, with an average value of the photon index $Γ\approx 1.95$. This "transition" is consistent with an external shock model having the circumburst medium transition from a wind to the ISM, by introducing a long-lasting energy injection with a Lorentz factor stratification of the ejecta. The best parameters from Monte Carlo Markov Chain fitting are $E_{\rm K,iso} \approx 2.14\times10^{55}$ erg, $Γ_0 \approx 162$, $ε_e \approx 0.02$, $ε_B \approx 1.7\times10^{-6}$, $A_\ast \approx 1.0$, $R_t \approx 4.1\times10^{17}$ cm, $n \approx 11.0 \rm\ cm^{-3}$, $L_0 \approx 3.1\times10^{52} \rm\ erg\ s^{-1}$, $k \approx 1.98$, $s \approx 1.54$, and $θ_j > 0.3$ rad.

preprint2020arXiv

Is GRB 110715A the progenitor of FRB 171209?

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) σ$. 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.

preprint2020arXiv

Jet Properties of Compact Steep-Spectrum Sources and an Eddington-Ratio-Driven Unification Scheme of Jet Radiation in Active Galactic Nuclei

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

preprint2020arXiv

Leptonic or Hadronic Emission: X-Ray Radiation Mechanism of Large-scale Jet Knots in 3C 273

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.

preprint2020arXiv

Lorentz Invariance Violation Limits from the Spectral Lag Transition of GRB~190114C

The spectral lags of gamma-ray bursts (GRBs) have been viewed as the most promising probes of the possible violations of Lorentz invariance (LIV). However, these constraints usually depend on the assumption of the unknown intrinsic time lag in different energy bands and the use of a single highest-energy photon. A new approach to test the LIV effects has been proposed by directly fitting the spectral lag behavior of a GRB with a well-defined transition from positive lags to negative lags. This method simultaneously provides a reasonable formulation of the intrinsic time lag and robust lower limits on the quantum-gravity energy scales ($E_{\rm QG}$). In this work, we perform a global fitting to the spectral lag data of GRB~190114C by considering the possible LIV effects based on a Bayesian approach. We then derive limits on $E_{\rm QG}$ and the coefficients of the Standard Model Extension. The Bayes factors output in our analysis shows a very strong evidence for the spectral-lag transition in GRB~190114C. Our constraints on a variety of isotropic and anisotropic coefficients for LIV are somewhat weaker than existing bounds, but they can be viewed as comparatively robust and have the promise to complement existing LIV constraints. The observations of GRBs with higher-energy emissions and higher temporal resolutions will contribute to a better formulation of the intrinsic time lag and more rigorous LIV constraints in the dispersive photon sector.

preprint2020arXiv

Magnetar as Central Engine of Gamma-Ray Bursts: Quasi-Universal Jet, Event Rate and X-ray Luminosity Function of Dipole Radiations

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$σ$) and an opening angle $θ_{\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$σ$). The inferred local GRB rate is $ρ=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 $β_1=0.78^{+0.16}_{-0.15}$ and $β_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}$, $θ_{\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.

preprint2020arXiv

Photospheric Emission in Gamma-ray Bursts: Variability

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 $η$ (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 $η$ is genetic variable, the variability of the photospheric emission is not only limited by the jet power but also affected by $η$ 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.

preprint2020arXiv

The Energy Sources of Double-peaked superluminous supernova PS1-12cil and luminous supernova SN 2012aa

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.

preprint2020arXiv

The properties of prompt emission in short GRBs with extended emission observed by {\em Fermi}/GBM

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.

preprint2020arXiv

Very High Energy Emission and Cascade Radiation of Gamma-Ray Burst Afterglows: Homogeneous Versus Wind External Media

Recent detection of sub-TeV emission from gamma-ray bursts (GRBs) represents a breakthrough in the GRB study. The multi-wavelength data of the afterglows of GRB 190114C support the synchrotron self-Compton (SSC) origin for its sub-TeV emission. We present a comparative analysis on the SSC emission of GRB afterglows in the homogeneous and wind environment in the framework of the forward shock model. The $γγ$ absorption of very high-energy photons due to pair production within the source and the Klein-Nishina effect on the inverse-Compton scattering are considered. Generally a higher SSC flux is expected for a larger circum-burst density due to a larger Compton parameter, but meanwhile the internal $γγ$ absorption is more severer for sub-TeV emission. The flux ratio between the SSC component and the synchrotron component decreases more quickly with time in the wind medium case than that in the homogenous-density medium case. The light curves of the SSC emission are also different for the two types of media. We also calculate the cascade emission resulted from the absorbed high-energy photons. In the ISM environment with $n> 1\,\rm cm^{-3}$, the cascade synchrotron emission could be comparable to the synchrotron emission of the primary electrons in the optical band, which may flatten the optical afterglow light curve at early time ($t<1$ h). In the wind medium with $A_{\ast}> 0.1$, the cascade emission in the eV-GeV band is comparable or even larger than the emission of the primary electrons at early time.

preprint2016arXiv

Central Engine of Late-Time X-ray Flares with Internal Origin

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.

preprint2016arXiv

Comprehensive study of the X-ray flares from gamma-ray bursts observed by Swift

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.

preprint2016arXiv

Cosmic Evolution of Long Gamma-Ray Burst Luminosity

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

preprint2016arXiv

GRB 110530A: Peculiar Broad Bump and Delayed Plateau in Early Optical Afterglows

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 $ε_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.

preprint2016arXiv

Lorentz factor - Beaming Corrected Energy/Luminosity Correlations and GRB Central Engine Models

We work on a GRB sample whose initial Lorentz factors ($Γ_0$) are constrained with the afterglow onset method and the jet opening angles ($θ_{\rm j}$) are determined by the jet break time. We confirm the $Γ_0$ - $E_{γ,\rm iso}$ correlation by Liang et al. (2010), and the $Γ_0$ - $L_{γ,\rm iso}$ correlation by L{ü} et al. (2012). Furthermore, we find correlations between $Γ_0$ and the beaming corrected $γ$-ray energy ($E_γ$) and mean $γ$-ray luminosity ($L_γ$). By also including the kinetic energy of the afterglow, we find rough correlations (with larger scatter) between $Γ_0$ and the total ($γ$-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 $ν\barν$ - 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.

preprint2016arXiv

On the Late-Time Spectral Softening Found in X-ray Afterglows of Gamma-Ray Bursts

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.

preprint2016arXiv

Search for a gamma-ray line feature from a group of nearby Galaxy clusters with Fermi LAT Pass 8 data

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 $γ-$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 $γ-$rays (i.e., $\langle σv \rangle_{χχ\rightarrow γγ}$) 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 $γ-$ray data.

preprint2016arXiv

Search for gamma-ray line features from Milky Way satellites with Fermi-LAT Pass 8 data

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.

preprint2016arXiv

Similarity of Jet Radiation between Flat Spectrum Radio Quasars and GeV Narrow-Line Seyfert 1 Galaxies: A Universal $δ$-$L_{\rm c}$ Correlation

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 ($δ$) 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 $δ$-$L_{\rm c}$ correlation. This indicates that the variations of observed luminosities may be essentially due to the Doppler boosting effect. And the universal $δ$-$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 $δ$. We estimate the $δ$ 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.

preprint2016arXiv

The History of GRB Outflows: Ejection Lorentz Factor and Radiation Location of X-Ray Flares

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 ($Γ_{\rm X}$). We obtain $Γ_{\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 $Γ_{\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., $Γ\propto [t_{p}/(1+z)]^{-0.69\pm 0.06}$. These results may shed lights on the long term evolution of GRB central engines.

preprint2016arXiv

The radiative efficiency of relativistic jet and wind: A case study of GRB 070110

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 $η_γ \sim 6\%$. However, the radiative efficiency of internal plateau ($η_{\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).

preprint2016arXiv

Variabilities of Gamma-ray Bursts from Black Hole Hyper-accretion Disks

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.

preprint2016arXiv

Very Bright Prompt and Reverse Shock Emission of GRB 140512A

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ε_{\rm B, r}/ε_{\rm B, f}$ is 8187, indicating a high magnetization degree in the RS region. Measuring the relative radiation efficiency with $R_e\equivε_{\rm e, r}/ε_{\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.

preprint2015arXiv

A Correlated Study of Optical and X-ray Afterglows of GRBs

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.

preprint2015arXiv

A Tight $L_{\rm p, iso}-E'_{\rm p}-Γ_0$ Correlation of Gamma-Ray Bursts

We select a sample of 34 gamma-ray bursts (GRBs) whose $Γ_0$ values are derived with the onset peaks observed in the afterglow lightcurves (except for GRB 060218 whose $Γ_0$ is estimated with its radio data), and investigate the correlations among $Γ_0$, the isotropic peak luminosity ($L_{\rm iso}$), and the peak energy ($E_{\rm p,z}$) of the $νf_ν$ 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 $Γ_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} Γ_0^{1.32\pm 0.19}$ or $E_{\rm p,z} \propto L_{\rm iso}^{0.55\pm 0.06}Γ_0^{-0.50\pm 0.17}$. Nine other GRBs whose $Γ_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}-Γ_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σ$. We argue that the $L_{\rm iso} - E_{\rm p, z} - Γ_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.

preprint2015arXiv

Correlations among Jet, Accretion Disk, and Broad Line Region of Flat Spectrum Radio Quasars

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.

preprint2015arXiv

Correlations of Disk and Jet Emission Deviating from the Fundamental Plane

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.

preprint2015arXiv

How Bad/Good Are the External Forward Shock Afterglow Models of Gamma-Ray Bursts?

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 ($α$, $β$, $t_b$ and model parameters (energy injection index q, p, $θ_j$, $η_γ$, etc) of the GRBs in the Gold Sample, and derive constraints on the magnetization parameter $ε_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.

preprint2015arXiv

Multi-Wavelength Observations of GRB 111228A and Implications for the Fireball and its environment

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 $ε_e=(6.9\pm 0.3)\times 10^{-2}$, $ε_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 ($Γ_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}-Γ_0$ relation and bridges the typical GRBs and low luminosity GRBs in this relation.

preprint2014arXiv

Distributions of Gamma-Ray Bursts and Blazars in the $L_{\rm p}-E_{\rm p}$ Plane and Possible Implications for their Radiation Physics

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 $νf_ν$ 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.

preprint2014arXiv

Episodic jet power extracted from a spinning black hole surrounded by a neutrino-dominated accretion flow in gamma-ray bursts

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.

preprint2014arXiv

Internal Energy Dissipation of Gamma-Ray Bursts Observed with Swift: Precursors, Prompt Gamma-rays, Extended emission and Late X-ray Flares

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.

preprint2014arXiv

Jet Properties of GeV-Selected Radio-Loud Narrow-Line Seyfert 1 Galaxies and Possible Connection to Their Disk and Corona

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 (δ), 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.

preprint2014arXiv

Photosphere emission in the X-Ray Flares of Swift Gamma-Ray Bursts and Implications for the Fireball Properties

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 $Γ_{\rm ph}$ of the outflow. The derived $Γ_{\rm ph}$ range from $50$ to $150$ and a relation of $Γ_{\rm ph}$ to the thermal emission luminosity is found. It is consistent with the $Γ_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.

preprint2014arXiv

Relativistic Jet Properties of GeV-TeV Blazars and Possible Implications for the Jet Formation, Composition, and Cavity Kinematics

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.

preprint2014arXiv

The "amplitude" parameter of Gamma-Ray Bursts and its implications for GRB classification

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.

preprint2013arXiv

A Comprehensive Study of Gamma-Ray Burst Optical Emission: II. Afterglow Onset and Late Re-Brightening Components

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_{γ, 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 (Γ_0) and E_{γ, 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 Γ_0, typically several tens.

preprint2013arXiv

A Comprehensive Study of Gamma-Ray Burst Optical Emission: III. Brightness Distributions and Luminosity Functions of Optical Afterglows

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.

preprint2013arXiv

Blazar Anti-Sequence of Spectral Variation within individual Blazars: Cases for Mrk 501 and 3C 279

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.

preprint2013arXiv

GRB 121027A: long-lasting, energetic X-ray flares and clues to radiation mechanism and progenitor star

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.

preprint2013arXiv

Hyper-accreting black hole as GRB central engine. I: Baryon loading in GRB jets

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 $ν\barν$-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 $Γ-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.

preprint2013arXiv

Radiation Mechanism and Jet Composition of Gamma-Ray Bursts and GeV-TeV selected Radio Loud Active Galactic Nuclei

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 (σ) 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.

preprint2012arXiv

A Comprehensive Analysis of Fermi Gamma-ray Burst Data: II. $E_{\rm p}$-Evolution Patterns and Implications for the Observed Spectrum-Luminosity Relations

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_{γ, \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.

preprint2012arXiv

A Comprehensive Analysis of Fermi Gamma-ray Burst Data: III. Energy-Dependent T90 Distributions of GBM GRBs and Instrumental Selection Effect on Duration Classification

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.

preprint2012arXiv

A Comprehensive Study of Gamma-Ray Burst Optical Emission: I. Flares and Early Shallow Decay Component

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_{γ, \rm iso}$. The relation between isotropic luminosities of the flares and gamma-rays follows $L^{\rm F}_{\rm R, iso}\propto L_{γ, \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 $γ$-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.

preprint2012arXiv

Blazar Anti-Sequence of Spectral Variability for Individual TeV Blazars

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 ΔM in each blazar results in the observed blazar anti-sequence.

preprint2012arXiv

GRB 110721A: photosphere "death line" and the physical origin of the GRB "Band" function

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.

preprint2012arXiv

Lorentz Factor-Isotropic Luminosity/Energy Correlations of GRBS and Their Interpretation

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.

preprint2012arXiv

Observational constraints on the external shock prior emission hypothesis of GRBs

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

preprint2012arXiv

Physical origin of multi-wavelength emission of GRB 100418A and implications for its progenitor

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 (Γ= 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.

preprint2012arXiv

Radial Angular Momentum Transfer and Magnetic Barrier for Short-Type Gamma-Ray Burst Central Engine Activity

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.

preprint2012arXiv

Radiation Mechanisms and Physical Properties of GeV-TeV BL Lac Objects

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.

preprint2012arXiv

Selection Effects on the Observed Redshift Dependence of GRB Jet Opening Angles

Apparent redshift dependence of the jet opening angles ($θ_{\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 γ}$) and $θ_{\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 $θ_{\rm j}-z$ dependence. We find that in case of $L_γ\propto θ_{\rm j}^2$ good consistency between the mock and observed samples can be obtained, indicating that both $L_{\rm γ}$ and $θ_{\rm j}$ are degenerate for a flux-limited sample. The parameter set $(L_{\rm γ}, θ_{\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.

preprint2012arXiv

Statistical Properties of Multiple Optical Emission Components in Gamma-Ray Bursts and Implications

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.

preprint2011arXiv

A Comprehensive Analysis of Fermi Gamma-Ray Burst Data. I. Spectral Components and Their Possible Physical Origins of LAT/GBM GRBs

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

preprint2011arXiv

Radiation Mechanisms and Physical Properties of GeV-TeV BL Lac Objects

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 $ν_{\rm s}$ and $ν_{\rm c}$ also favors the model. It is found that the break of the electron spectrum ($γ_{\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 $γ_{\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 $ν_{\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.

preprint2010arXiv

Jet Precession Driven by Neutrino-Cooled Disc for Gamma-Ray Bursts

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 $α=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 $α=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.

preprint2010arXiv

Late-time detections of the X-ray afterglow of GRB 060729 with Chandra - the latest detections ever of an X-ray afterglow

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.

preprint2010arXiv

Simulations on High-z Long Gamma-Ray Burst Rate

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.

preprint2010arXiv

The $E_{\rm p}$-Flux Correlation in the Rising and Decaying Phases of Gamma-Ray Burst Pulses: Evidence for Viewing Angle Effect?

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 $νf_ν$ 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.

preprint2010arXiv

XRF 100316D/SN 2010bh: clue to the diverse origin of nearby supernova-associated GRBs

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.

preprint2009arXiv

A Comprehensive Analysis of Swift/XRT Data: IV. Single Power-Law Decaying Lightcurves vs. Canonical Lightcurves and Implications for a Unified Origin of X-rays

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.

preprint2009arXiv

Constraining GRB Initial Lorentz Factor with the Afterglow Onset Feature and Discovery of a Tight Gamma_0-E_iso Correlation

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

preprint2009arXiv

Towards the Properties of Long Gamma-Ray Burst Progenitors with Swift Data

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.

preprint2004arXiv

Distribution of the extensive Doppler redshift of quasars

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.

preprint2002arXiv

Statistical Shapes of the Highest Pulses in Gamma--Ray Bursts

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.