Source author record

Xue-Feng Wu

Xue-Feng Wu 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

66works
7topics
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

66 published item(s)

preprint2026arXiv

Investigating the Anisotropy of Dispersion Measure Contribution from the Galactic Halo by Using Fast Radio Bursts

We propose a data-driven approach to reconstruct the all-sky distribution of the dispersion measure contribution from the Galactic halo ($\mathrm{DM_{halo}}$) through a spherical harmonic expansion, enabling an investigation of its possible anisotropies. Based on the NE2001 model and using 92 localized and 574 unlocalized non-repeating fast radio bursts (FRBs) at Galactic latitudes $|b|>15^\circ$, we find a significant dipole anisotropy in $\mathrm{DM_{halo}}$, pointing toward $(l=130^\circ,\, b=+5^\circ)$ with a $1σ$ uncertainty of approximately $28^\circ$. The $\mathrm{DM_{halo}}$ value in this direction is $63\pm9~\mathrm{pc~cm^{-3}}$, exceeding the all-sky mean by about $2.6σ$. This result is not significantly affected by the choice of Galactic ISM models. Furthermore, even when using a refined sample of 62 localized FRBs (excluding CHIME detections, repeaters, and unlocalized events), the dipole anisotropic structure persists, with a direction of $(l=141^\circ,\, b=+51^\circ)$ and a larger 1$σ$ uncertainty of $\sim 44^\circ$. Model comparisons using the Akaike Information Criterion and Bayesian evidence yield consistent preferences, and together they suggest that current FRB data slightly favor the existence of a dipole structure in $\mathrm{DM_{halo}}$. If this feature is not a statistical fluctuation or systematic error, its physical origin requires further investigation. Future FRB samples with larger sizes and more complete sky coverage will be essential to confirm or refute this possible anisotropic structure.

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 semi-analytical solution to the forward-reverse shock hydrodynamics of the gamma-ray burst afterglow

We extend the standard model of forward-reverse shock (FS-RS) for gamma-ray burst (GRB) afterglow to more general cases. On one hand, we derive the analytical solution to the hydrodynamics of the shocks in two limiting cases, i.e., an ultra-relativistic reverse shock case and a Newtonian reverse shock case. Based on the asymptotic solutions in these two limiting cases, we constitute a semi-analytical solution for the hydrodynamics of the shocks in the generic case, covering the mildly-relativistic reverse shock case. On the other hand, we derive the evolution of the system taking into account the condition of energy conservation which is not satisfied in the standard FS-RS model. A generic solution of semi-analytical expressions is also given. In both the extended standard FS-RS model (satisfying pressure balance condition) and the model satisfying energy conservation, we find that the results in the ultra-relativistic reverse shock case and in the early stage of the Newtonian reverse shock case are different from those in the standard FS-RS model by only a factor that close to one while the same initial conditions adopted. However, the asymptotic solutions in the limiting cases are not good approximations to those in the intermediate case. Our semi-analytical results agree well with the numerical results for a large range of model parameters, and hence can be easily employed to diagnose the physical quantities of the GRB shell and circumburst environment.

preprint2022arXiv

Bounding the Photon Mass with the Dedispersed Pulses of the Crab Pulsar and FRB 180916B

Tight limits on the photon mass have been set through analyzing the arrival time differences of photons with different frequencies originating from the same astrophysical source. However, all these constraints have relied on using the first-order Taylor expansion of the dispersion due to a nonzero photon mass. In this work, we present an analysis of the nonzero photon mass dispersion with the second-order derivative of Taylor series. If the arrival time delay corrected for all known effects (including the first-order delay time due to the plasma and photon mass effects) is assumed to be dominated by the second-order term of the nonzero photon mass dispersion, a conservative upper limit on the photon mass can be estimated. Here we show that the dedispersed pulses with the second-order time delays from the Crab pulsar and the fast radio burst FRB 180916B pose strict limits on the photon mass, i.e., $m_{γ,2} \leq5.7\times10^{-46}\;{\rm kg}\simeq3.2\times10^{-10}\; {\rm eV}/c^{2}$ and $m_{γ,2} \leq6.0\times10^{-47}\;{\rm kg}\simeq3.4\times10^{-11}\; {\rm eV}/c^{2}$, respectively. This is the first time to study the possible second-order photon mass effect.

preprint2022arXiv

Diagnosing Circumburst Environment with Multiband Gamma-Ray Burst Radio Afterglows

It has been widely recognized that gamma-ray burst (GRB) afterglows arise from interactions between GRB outflow and circumburst medium, while their evolution follows the behaviors of relativistic shock waves. Assuming the distribution of circumburst medium follows a general power-law form, that is, $n = A_{\ast} R^{-k}$, where $R$ denotes the distance from the burst, it is obvious that the value of density-distribution index $k$ can affect the behaviors of the afterglow. In this paper, we analyze the temporal and spectral behaviors of GRB radio afterglows with arbitrary $k$-values. In the radio band, a standard GRB afterglow produced by forward shock exhibits a late-time flux peak, and the relative peak fluxes as well as peak times at different frequencies show dependencies on $k$. Thus with multi-band radio peak observations, one can determine the density profile of circumburst medium by comparing the relations between peak flux/time and frequency at each observing band. Also, the effects of trans-relativistic shock waves, as well as jets in afterglows are discussed. By analyzing 31 long and 1 short GRBs with multi-band data of radio afterglows, we find that nearly half of them can be explained with uniform interstellar medium ($k=0$), $\sim 1/5$ can be constrained to exhibiting stellar wind environment ($k=2$), while less than $\sim 1/3$ samples show $0< k< 2$.

preprint2022arXiv

Direct Estimate of the Post-Newtonian Parameter and Cosmic Curvature from Galaxy-scale Strong Gravitational Lensing

Einstein's theory of general relativity (GR) has been precisely tested on solar system scales, but extragalactic tests are still poorly performed. In this work, we use a newly compiled sample of galaxy-scale strong gravitational lenses to test the validity of GR on kiloparsec scales. In order to solve the circularity problem caused by the preassumption of a specific cosmological model based on GR, we employ the distance sum rule in the Friedmann-Lema\^ıtre-Robertson-Walker metric to directly estimate the parameterized post-Newtonian (PPN) parameter $γ_{\rm PPN}$ and the cosmic curvature $Ω_k$ by combining observations of strong lensing and Type Ia supernovae. This is the first simultaneous measurement of $γ_{\rm PPN}$ and $Ω_k$ without any assumptions about the contents of the universe or the theory of gravity. Our results show that $γ_{\rm PPN}=1.11^{+0.11}_{-0.09}$ and $Ω_{k}=0.48^{+1.09}_{-0.71}$, indicating a strong degeneracy between the two quantities. The measured $γ_{\rm PPN}$, which is consistent with the prediction of 1 from GR, provides a precise extragalactic test of GR with a fractional accuracy better than 9.0\%. If a prior of the spatial flatness (i.e., $Ω_{k}=0$) is adopted, the PPN parameter constraint can be further improved to $γ_{\rm PPN}=1.07^{+0.07}_{-0.07}$, representing a precision of 6.5\%. On the other hand, in the framework of GR (i.e., $γ_{\rm PPN}=1$), our results are still marginally compatible with zero curvature ($Ω_k=-0.12^{+0.48}_{-0.36}$), supporting no significant deviation from a flat universe.

preprint2022arXiv

Gravitationally Lensed Orphan Afterglows of Gamma-Ray Bursts

The cosmological nature of gamma-ray bursts (GRBs) implies that a small portion of them could be gravitationally lensed by foreground objects during their propagation. The gravitational lensing effect on the GRB prompt emission and on-axis afterglows has been discussed, and some candidates have been found in the literature. In this work, considering the high detection rate of GRB orphan afterglows in future wide-field survey era, we investigate the gravitationally lensed orphan afterglows in view of three lens models, i.e., the point-mass model, the singular isothermal sphere model, and the Chang-Refsdal model. The structure of the GRB jet itself is also incorporated in calculating the lensed afterglow light curves. It is found that lensed optical/X-ray orphan afterglows in principle could be diagnosed through their temporal characteristics, and the optical band is the best band to observe the galaxy-lensed orphan afterglows. Moreover, the event rate for galaxy-lensed orphan afterglows is estimated to be $\lesssim$ 1.8 $\text{yr}^{-1}$ for the whole sky. If most orphan afterglows could be identified (from other transients in the survey data), the optimistic detection rates of the 2.5m Wide Field Survey Telescope of China and 8.4m Vera Rubin Observatory Legacy Survey of Space and Time for galaxy-lensed orphan afterglows in the optical band are $\lesssim$ 0.01$\sim$0.02 $\text{yr}^{-1}$ and $\lesssim$ 0.04$\sim$0.08 $\text{yr}^{-1}$, respectively.

preprint2021arXiv

Testing Fundamental Physics with Astrophysical Transients

Explosive astrophysical transients at cosmological distances can be used to place precision tests of the basic assumptions of relativity theory, such as Lorentz invariance, the photon zero-mass hypothesis, and the weak equivalence principle (WEP). Signatures of Lorentz invariance violations (LIV) include vacuum dispersion and vacuum birefringence. Sensitive searches for LIV using astrophysical sources such as gamma-ray bursts, active galactic nuclei, and pulsars are discussed. The most direct consequence of a nonzero photon rest mass is a frequency dependence in the velocity of light propagating in vacuum. A detailed representation of how to obtain a combined severe limit on the photon mass using fast radio bursts at different redshifts through the dispersion method is presented. The accuracy of the WEP has been well tested based on the Shapiro time delay of astrophysical messengers traveling through a gravitational field. Some caveats of Shapiro delay tests are discussed. In this article, we review and update the status of astrophysical tests of fundamental physics.

preprint2020arXiv

Combined Limit on the Photon Mass with Nine Localized Fast Radio Bursts

A nonzero-mass hypothesis for the photon can produces a frequency-dependent dispersion of light, which results in arrival-time differences of photons with different frequencies originating from a given transient source. Extragalactic fast radio bursts (FRBs), with their low frequency emissions, short time durations, and long propagation distances, are excellent astrophysical probes to constrain the rest mass of the photon $m_γ$. However, the derivation of a limit on $m_γ$ is complicated by the similar frequency dependences of dispersion expected from the plasma and nonzero photon mass effects. If a handful measurements of redshift for FRBs are available, the different redshift dependences of the plasma and photon mass contributions to the dispersion measure (DM) might be able to break dispersion degeneracy in testing the photon mass. For now, nine FRBs with redshift measurements have been reported, which can turn this idea into reality. Taking into account the DM contributions from both the plasma and a possible photon mass, we use the data on the nine FRBs to derive a combined limit of $m_γ\leq7.1\times10^{-51}\;{\rm kg}$, or equivalently $m_γ\leq4.0\times10^{-15}\; {\rm eV}/c^{2}$ at 68\% confidence level, which is essentially as good as or represents a factor of 7 improvement over previous limits obtained by the single FRBs. Additionally, a reasonable estimation for the DM contribution from the host galaxy, $\rm DM_{host}$, can be simultaneously achieved in our analysis. The rapid progress in localizing FRBs will further tighten the constraints on both $m_γ$ and $\rm DM_{host}$.

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

Self-organized Criticality in Multi-pulse Gamma-Ray Bursts

The variability in multi-pulse gamma-ray bursts (GRBs) may help to reveal the mechanism of underlying processes from the central engine. To investigate whether the self-organized criticality (SOC) phenomena exist in the prompt phase of GRBs, we statistically study the properties of GRBs with more than 3 pulses in each burst by fitting the distributions of several observed physical variables with a Markov Chain Monte Carlo approach, including the isotropic energy $E_{\rm iso}$, the duration time $T$ and the peak count rate $P$ of each pulse. Our sample consists of 454 pulses in 93 GRBs observed by the CGRO/BATSE satellite. The best-fitting values and uncertainties for these power-law indices of the differential frequency distributions are: $α^d_{E}=1.54 \pm 0.09$, $α^d_{T}=1.82_{-0.15}^{+0.14}$ and $α^d_{P}=2.09_{-0.19}^{+0.18}$, while the power-law indices in the cumulative frequency distributions are: $α^c_{E}=1.44_{-0.10}^{+0.08}$, $α^c_{T}=1.75_{-0.13}^{+0.11}$ and $α^c_{P}=1.99_{-0.19}^{+0.16}$. We find that these distributions are roughly consistent with the physical framework of a Fractal-Diffusive, Self-Organized Criticality (FD-SOC) system with the spatial dimension $S=3$ and the classical diffusion $β$=1. Our results support that the jet responsible for the GRBs should be magnetically dominated and magnetic instabilities (e.g., kink model, or tearing-model instability) lead the GRB emission region into the SOC state.

preprint2020arXiv

Testing the Weak Equivalence Principle and Lorentz Invariance with Multiwavelength Polarization Observations of GRB Optical Afterglows

Violations of both the weak equivalence principle (WEP) and Lorentz invariance can produce vacuum birefringence, which leads to an energy-dependent rotation of the polarization vector of linearly polarized emission from a given astrophysical source. However, the search for the birefringent effect has been hindered by our ignorance concerning the intrinsic polarization angle in different energy bands. Considering the contributions to the observed linear polarization angle from both the intrinsic polarization angle and the rotation angles induced by violations of the WEP and Lorentz invariance, and assuming the intrinsic polarization angle is an unknown constant, we simultaneously obtain robust bounds on possible deviations from the WEP and Lorentz invariance, by directly fitting the multiwavelength polarimetric data of the optical afterglows of gamma-ray burst (GRB) 020813 and GRB 021004. Here we show that at the $3σ$ confidence level, the difference of the parameterized post-Newtonian parameter $γ$ values characterizing the departure from the WEP is constrained to be $Δγ=\left(-4.5^{+10.0}_{-16.0}\right)\times10^{-24}$ and the birefringent parameter $η$ quantifying the broken degree of Lorentz invariance is limited to be $η=\left(6.5^{+15.0}_{-14.0}\right)\times10^{-7}$. These are the first simultaneous verifications of the WEP and Lorentz invariance in the photon sector. More stringent limits can be expected as the analysis presented here is applied to future multiwavelength polarization observations in the prompt gamma-ray emission of GRB

preprint2020arXiv

The Bright Reverse Shock Emission in the Optical Afterglows of Gamma-ray Bursts in a Stratified Medium

The reverse shock (RS) model is generally introduced to interpret the optical afterglows with the rapid rising and decaying, such as the early optical afterglow of GRB 990123 (which is also called optical flash). In this paper, we collected 11 gamma-ray bursts (GRBs) early optical afterglows, which have such signatures of dominant RS emission. Since the temporal slopes of the optical flashes are determined by both the medium density distribution index $k$ and the electron spectral index $p$, we apply the RS model of the thin shell case to the optical flashes and determine the ambient medium of the progenitors. We find that the $k$ value is in the range of 0 - 1.5. The $k$ value in this paper is consistent with the result in Yi et al. (2013), where the forward shock (FS) model was applied to some onset bumps. However, the method adopted in this paper is only applicable to GRB afterglows with significant sharp rising and decaying RS emission. Our results indicate that the RS model can also be applied to confirm the circumburst medium, further implying that GRBs may have diverse circumburst media.

preprint2020arXiv

The SiTian project

SiTian is an ambitious ground-based all-sky optical monitoring project, developed by the Chinese Academy of Sciences. The concept is an integrated network of dozens of 1-m-class telescopes deployed partly in China and partly at various other sites around the world. The main science goals are the detection, identification and monitoring of optical transients (such as gravitational wave events, fast radio bursts, supernovae) on the largely unknown timescales of less than 1 day; SiTian will also provide a treasure trove of data for studies of AGN, quasars, variable stars, planets, asteroids, and microlensing events. To achieve those goals, SiTian will scan at least 10,000 square deg of sky every 30 min, down to a detection limit of $V \approx 21$ mag. The scans will produce simultaneous light-curves in 3 optical bands. In addition, SiTian will include at least three 4-m telescopes specifically allocated for follow-up spectroscopy of the most interesting targets. We plan to complete the installation of 72 telescopes by 2030 and start full scientific operations in 2032.

preprint2016arXiv

Constraining the Mass of the Photon with Gamma-Ray Bursts

One of the cornerstones of modern physics is Einstein's special relativity, with its constant speed of light and zero photon mass assumptions. Constraint on the rest mass m_γ of photons is a fundamental way to test Einstein's theory, as well as other essential electromagnetic and particle theories. Since non-zero photon mass can give rise to frequency-(or energy-) dependent dispersions, measuring the time delay of photons with different frequencies emitted from explosive astrophysical events is an important and model-independent method to put such a constraint. The cosmological gamma-ray bursts (GRBs), with short time scales, high redshifts as well as broadband prompt and afterglow emissions, provide an ideal testbed for m_γ constraints. In this paper we calculate the upper limits of the photon mass with GRB early time radio afterglow observations as well as multi-band radio peaks, thus improve the results of Schaefer (1999) by nearly half an order of magnitude.

preprint2016arXiv

Constraints on the Photon Mass with Fast Radio Bursts

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

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

Fast Radio Bursts from the Inspiral of Double Neutron Stars

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

preprint2016arXiv

Limits on the Neutrino Velocity, Lorentz Invariance, and the Weak Equivalence Principle with TeV Neutrinos from Gamma-Ray Bursts

Five TeV neutrino events weakly correlated with five gamma-ray bursts (GRBs) were detected recently by IceCube. This work is an attempt to show that if the GRB identifications are verified, the observed time delays between the TeV neutrinos and gamma-ray photons from GRBs provide attractive candidates for testing fundamental physics with high accuracy. Based on the assumed associations between the TeV neutrinos and GRBs, we find that the limiting velocity of the neutrinos is equal to that of photons to an accuracy of $\sim1.9\times10^{-15}-2.5\times10^{-18}$, which is about $10^{4}-10^{7}$ times better than the constraint obtained with the neutrino possibly from a blazar flare. In addition, we set the most stringent limits up to date on the energy scale of quantum gravity for both the linear and quadratic violations of Lorentz invariance, namely $E_{\rm QG, 1}>6.3\times10^{18}-1.5\times10^{21}$ GeV and $E_{\rm QG, 2}>2.0\times10^{11}-4.2\times10^{12}$ GeV, which are essentially as good as or are an improvement of one order of magnitude over the results previously obtained by the GeV photons of GRB 090510 and the PeV neutrino from a blazar flare. Assuming that the Shapiro time delay is caused by the gravitational potential of the Laniakea supercluster of galaxies, we also place the tightest limits to date on Einstein's weak equivalence principle through the relative differential variations of the parameterized post-Newtonian parameter $γ$ values for two different species of particles (i.e., neutrinos and photons), yielding $Δγ\sim 10^{-11}-10^{-13}$. However, it should be emphasized again that these limits here obtained are at best forecast of what could be achieved if the GRB/neutrino correlations would be finally confirmed.

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

Solving the 56Ni puzzle of magnetar-powered broad-lined type Ic supernovae

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

preprint2016arXiv

Testing Einstein's Equivalence Principle with Supercluster Laniakea's Gravitational Field

Comparing the parameterized post-Newtonian parameter $γ$ values for different types of particles, or the same type of particles with different energies is an important method to test the Einstein Equivalence Principle (EEP). Assuming that the observed time delays are dominated by the gravitational potential of the Laniakea supercluster of galaxies, better results of EEP constraints can be obtained. In this paper, we apply photons from three kinds of cosmic transients, including TeV blazars, gamma-ray bursts as well as fast radio bursts to constrain EEP. With a gravitational field far more stronger than a single galaxy, we obtain 4--5 orders of magnitude more stringent than the pervious results.

preprint2016arXiv

Testing Einstein's Weak Equivalence Principle With Gravitational Waves

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

preprint2016arXiv

Testing Models for the Shallow Decay Phase of Gamma-Ray Burst Afterglows with Polarization Observations

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

preprint2016arXiv

Tests of the Einstein Equivalence Principle using TeV Blazars

The observed time delays between different energy bands from TeV blazars provide a new interesting way of testing the Einstein Equivalence Principle (EEP). If the whole time delay is assumed to be dominated by the gravitational field of the Milky Way, the conservative upper limit on the EEP can be estimated. Here we show that the strict limits on the differences of the parameterized post-Newtonian parameter $γ$ values are $γ_{\rm TeV}-γ_{\rm keV}<3.86\times10^{-3}$ for Mrk 421 and $γ_{\rm TeV}-γ_{\rm keV}<4.43\times10^{-3}$ for Mrk 501, while expanding the scope of the tested EEP energy range out to the TeV--keV range for the first time. With the small time lag from the 0.2--0.8 TeV and $>0.8$ TeV light curves of PKS 2155-304, a much more severe constraint on $γ$ differences of $\sim10^{-6}$ can be achieved, although the energy difference is of order of $\sim$ TeV. Furthermore, we can combine these limits on the energy dependence of $γ$ with the bound on the absolute $γ$ value $γ-1\sim0.3\%$ from light deflection measurements at the optical (eV) bands, and conclude that this absolute bound on $γ$ can be extended from optical to TeV energies.

preprint2016arXiv

The HII Galaxy Hubble Diagram Strongly Favors $R_{\rm h}=ct$ over $Λ$CDM

We continue to build support for the proposal to use HII galaxies (HIIGx) and giant extragalactic HII regions (GEHR) as standard candles to construct the Hubble diagram at redshifts beyond the current reach of Type Ia supernovae. Using a sample of 25 high-redshift HIIGx, 107 local HIIGx, and 24 GEHR, we confirm that the correlation between the emission-line luminosity and ionized-gas velocity dispersion is a viable luminosity indicator, and use it to test and compare the standard model $Λ$CDM and the $R_{\rm h}=ct$ Universe by optimizing the parameters in each cosmology using a maximization of the likelihood function. For the flat $Λ$CDM model, the best fit is obtained with $Ω_{\rm m}= 0.40_{-0.09}^{+0.09}$. However, statistical tools, such as the Akaike (AIC), Kullback (KIC) and Bayes (BIC) Information Criteria favor $R_{\rm h}=ct$ over the standard model with a likelihood of $\approx 94.8\%-98.8\%$ versus only $\approx 1.2\%-5.2\%$. For $w$CDM (the version of $Λ$CDM with a dark-energy equation of state $w_{\rm de}\equiv p_{\rm de}/ρ_{\rm de}$ rather than $w_{\rm de}=w_Λ=-1$), a statistically acceptable fit is realized with $Ω_{\rm m}=0.22_{-0.14}^{+0.16}$ and $w_{\rm de}= -0.51_{-0.25}^{+0.15}$ which, however, are not fully consistent with their concordance values. In this case, $w$CDM has two more free parameters than $R_{\rm h}=ct$, and is penalized more heavily by these criteria. We find that $R_{\rm h}=ct$ is strongly favored over $w$CDM with a likelihood of $\approx 92.9\%-99.6\%$ versus only $0.4\%-7.1\%$. The current HIIGx sample is already large enough for the BIC to rule out $Λ$CDM/$w$CDM in favor of $R_{\rm h}=ct$ at a confidence level approaching $3σ$.

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

Constraints on the Bulk Lorentz Factors of GRB X-Ray Flares

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

preprint2015arXiv

Cosmic Transients Test Einstein's Equivalence Principle out to GeV Energies

The Einstein Equivalence Principle (EEP) can be probed with astrophysical sources emitting simultaneously different types of neutral particles, or particles with varying energies, by testing their time of flight through the same gravitational field. Here we use the time delays between correlated photons from cosmological transients to constrain the accuracy of the EEP. We take data from two gamma-ray bursts as an example, and use, as a lower limit to the theoretical time delays between different energies, delays arising from only the gravitational field of our own galaxy. We then show that the parameterized post-Newtonian parameter $γ$ is the same for photons over energy ranges between eV and MeV and between MeV and GeV to a part in $10^{-7}$, which is at least one order of magnitude better than previous limits. Combining this bound on the wavelength dependence of $γ$ with the absolute bound $|γ-1|<0.3\%$ from light-deflection measurements at optical (eV) wavelengths, we thus extend this absolute bound on $γ$ to GeV energies.

preprint2015arXiv

GRB 080503 late afterglow re-brightening: signature of a magnetar powered merger-nova

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

preprint2015arXiv

Long GRBs as a Tool to Investigate Star Formation in Dark Matter Halos

First stars can only form in structures that are suitably dense, which can be parametrized by the minimum dark matter halo mass $M_{\rm min}$. $M_{\rm min}$ must plays an important role in star formation. The connection of long gamma-ray bursts (LGRBs) with the collapse of massive stars has provided a good opportunity for probing star formation in dark matter halos. We place some constraints on $M_{\rm min}$ using the latest $Swift$ LGRB data. We conservatively consider that LGRB rate is proportional to the cosmic star formation rate (CSFR) and an additional evolution parametrized as $(1+z)^α$, where the CSFR model as a function of $M_{\rm min}$. Using the $χ^{2}$ statistic, the contour constraints on the $M_{\rm min}$--$α$ plane show that at the $1σ$ confidence level, we have $M_{\rm min}<10^{10.5}$ $\rm M_{\odot}$ from 118 LGRBs with redshift $z<4$ and luminosity $L_{\rm iso}>1.8\times10^{51}$ erg $\rm s^{-1}$. We also find that adding 12 high-\emph{z} $(4<z<5)$ LGRBs (consisting of 104 LGRBs with $z<5$ and $L_{\rm iso}>3.1\times10^{51}$ erg $\rm s^{-1}$) could result in much tighter constraints on $M_{\rm min}$, for which, $10^{7.7}\rm M_{\odot}<M_{\rm min}<10^{11.6}\rm M_{\odot}$ ($1σ$). Through Monte Carlo simulations, we estimate that future five years of Sino-French spacebased multiband astronomical variable objects monitor (\emph{SVOM}) observations would tighten these constraints to $10^{9.7}\rm M_{\odot}<M_{\rm min}<10^{11.3}\rm M_{\odot}$. The strong constraints on $M_{\rm min}$ indicate that LGRBs are a new promising tool for investigating star formation in dark matter halos.

preprint2015arXiv

Polarization Evolution of Early Optical Afterglows of Gamma-Ray Bursts

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

preprint2015arXiv

Radio Afterglows and Host Galaxies of Gamma-Ray Bursts

Considering the contribution of the emission from the host galaxies of gamma-ray bursts (GRBs) to the radio afterglows, we investigate the effect of host galaxies on observations statistically. For the three types of events, e.g. low-luminosity, standard and high-luminosity GRBs, it is found that a tight correlation exists between the ratio of the radio flux (RRF) of host galaxy to the total radio peak emission and the observational frequency. Especially, toward lower frequencies, the contribution from the host increases significantly. The correlation can be used to get a useful estimate for the radio brightness of those host galaxies which only have very limited radio afterglow data. Using this prediction, we re-considered the theoretical radio afterglow light curves for four kinds of events, i.e. high-luminosity, low-luminosity, standard and failed GRBs, taking into account the contribution from the host galaxies and aiming at exploring the detectability of these events by the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Lying at a typical redshift of $z=1$, most of the events can be detected by FAST easily. For the less fierce low-luminosity GRBs, their radio afterglows are not strong enough to exceed the sensitivity limit of FAST at such distances. However, since a large number of low luminosity bursts actually happen very near to us, it is expected that FAST will still be able to detect many of them.

preprint2015arXiv

Testing Cosmological Models with Type Ic Super Luminous Supernovae

The use of type Ic Super Luminous Supernovae (SLSN Ic) to examine the cosmological expansion introduces a new standard ruler with which to test theoretical models. The sample suitable for this kind of work now includes 11 SLSNe Ic, which have thus far been used solely in tests involving $Λ$CDM. In this paper, we broaden the base of support for this new, important cosmic probe by using these observations to carry out a one-on-one comparison between the $R_{\rm h}=ct$ and $Λ$CDM cosmologies. We individually optimize the parameters in each cosmological model by minimizing the $χ^{2}$ statistic. We also carry out Monte Carlo simulations based on these current SLSN Ic measurements to estimate how large the sample would have to be in order to rule out either model at a $\sim 99.7\%$ confidence level. The currently available sample indicates a likelihood of $\sim$$70-80\%$ that the $R_{\rm h}=ct$ Universe is the correct cosmology versus $\sim$$20-30\%$ for the standard model. These results are suggestive, though not yet compelling, given the current limited number of SLSNe Ic. We find that if the real cosmology is $Λ$CDM, a sample of $\sim$$240$ SLSNe Ic would be sufficient to rule out $R_{\rm h}=ct$ at this level of confidence, while $\sim$$480$ SLSNe Ic would be required to rule out $Λ$CDM if the real Universe is instead $R_{\rm h}=ct$. This difference in required sample size reflects the greater number of free parameters available to fit the data with $Λ$CDM. If such SLSNe Ic are commonly detected in the future, they could be a powerful tool for constraining the dark-energy equation of state in $Λ$CDM, and differentiating between this model and the $R_{\rm h}=ct$ Universe.

preprint2015arXiv

The Age-Redshift Relationship of Old Passive Galaxies

We use 32 age measurements of passively evolving galaxies as a function of redshift to test and compare the standard model ($Λ$CDM) with the $R_{\rm h}=ct$ Universe. We show that the latter fits the data with a reduced $χ^2_{\rm dof}=0.435$ for a Hubble constant $H_{0}= 67.2_{-4.0}^{+4.5}$ km $\rm s^{-1}$ $\rm Mpc^{-1}$. By comparison, the optimal flat $Λ$CDM model, with two free parameters (including $Ω_{\rm m}=0.12_{-0.11}^{+0.54}$ and $H_{0}=94.3_{-35.8}^{+32.7}$ km $\rm s^{-1}$ $\rm Mpc^{-1}$), fits the age-\emph{z} data with a reduced $χ^2_{\rm dof}=0.428$. Based solely on their $χ^2_{\rm dof}$ values, both models appear to account for the data very well, though the optimized $Λ$CDM parameters are only marginally consistent with those of the concordance model ($Ω_{\rm m}=0.27$ and $H_{0}= 70$ km $\rm s^{-1}$ $\rm Mpc^{-1}$). Fitting the age-$z$ data with the latter results in a reduced $χ^2_{\rm dof}=0.523$. However, because of the different number of free parameters in these models, selection tools, such as the Akaike, Kullback and Bayes Information Criteria, favour $R_{\rm h}=ct$ over $Λ$CDM with a likelihood of $\sim 66.5\%-80.5\%$ versus $\sim 19.5\%-33.5\%$. These results are suggestive, though not yet compelling, given the current limited galaxy age-$z$ sample. We carry out Monte Carlo simulations based on these current age measurements to estimate how large the sample would have to be in order to rule out either model at a $\sim 99.7\%$ confidence level. We find that if the real cosmology is $Λ$CDM, a sample of $\sim 45$ galaxy ages would be sufficient to rule out $R_{\rm h}=ct$ at this level of accuracy, while $\sim 350$ galaxy ages would be required to rule out $Λ$CDM if the real Universe were instead $R_{\rm h}=ct$.

preprint2015arXiv

Utilizing the Updated Gamma-Ray Bursts and Type Ia Supernovae to Constrain the Cardassian Expansion Model and Dark Energy

We update gamma-ray burst (GRB) luminosity relations among certain spectral and light-curve features with 139 GRBs. The distance modulus of 82 GRBs at $z>1.4$ can be calibrated with the sample at $z\leq1.4$ by using the cubic spline interpolation method from the Union2.1 Type Ia supernovae (SNe Ia) set. We investigate the joint constraints on the Cardassian expansion model and dark energy with 580 Union2.1 SNe Ia sample ($z<1.4$) and 82 calibrated GRBs data ($1.4<z\leq8.2$). In $Λ$CDM, we find that adding 82 high-\emph{z} GRBs to 580 SNe Ia significantly improves the constrain on $Ω_{m}-Ω_Λ$ plane. In the Cardassian expansion model, the best fit is $Ω_{m}= 0.24_{-0.15}^{+0.15}$ and $n=0.16_{-0.52}^{+0.30}$ $(1σ)$, which is consistent with the $Λ$CDM cosmology $(n=0)$ in the $1σ$ confidence region. We also discuss two dark energy models in which the equation of state $w(z)$ is parametrized as $w(z)=w_{0}$ and $w(z)=w_{0}+w_{1}z/(1+z)$, respectively. Based on our analysis, we see that our Universe at higher redshift up to $z=8.2$ is consistent with the concordance model within $1σ$ confidence level.

preprint2014arXiv

A Comparison of Cosmological Models Using Time Delay Lenses

The use of time-delay gravitational lenses to examine the cosmological expansion introduces a new standard ruler with which to test theoretical models. The sample suitable for this kind of work now includes 12 lens systems, which have thus far been used solely for optimizing the parameters of $Λ$CDM. In this paper, we broaden the base of support for this new, important cosmic probe by using these observations to carry out a one-on-one comparison between {\it competing} models. The currently available sample indicates a likelihood of $\sim 70-80%$ that the $R_{\rm h}=ct$ Universe is the correct cosmology versus $\sim 20-30%$ for the standard model. This possibly interesting result reinforces the need to greatly expand the sample of time-delay lenses, e.g., with the successful implementation of the Dark Energy Survey, the VST ATLAS survey, and the Large Synoptic Survey Telescope. In anticipation of a greatly expanded catalog of time-delay lenses identified with these surveys, we have produced synthetic samples to estimate how large they would have to be in order to rule out either model at a $\sim 99.7%$ confidence level. We find that if the real cosmology is $Λ$CDM, a sample of $\sim 150$ time-delay lenses would be sufficient to rule out $R_{\rm h}=ct$ at this level of accuracy, while $\sim 1,000$ time-delay lenses would be required to rule out $Λ$CDM if the real Universe is instead $R_{\rm h}=ct$. This difference in required sample size reflects the greater number of free parameters available to fit the data with $Λ$CDM.

preprint2014arXiv

A double neutron star merger origin for the cosmological relativistic fading source PTF11agg?

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

preprint2014arXiv

Cosmological Tests Using GRBs, the Star Formation Rate and Possible Abundance Evolution

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

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

Time evolution of flares in GRB 130925A: jet precession in black hole accretion system

GRB 130925A, composed of three gamma-ray emission episodes and a series of orderly flares, has been detected by $Swift$, Fermi, Konus-$Wind$, and $INTEGRAL$. If the third weakest gamma-ray episode can be considered as a giant flare, we find that after the second gamma-ray episode observed by $INTEGRAL$ located at about 2000 s, there exists a positive relation between the time intervals of the adjacent flares and the time since the episode. We suggest that the second gamma-ray episode and its flares originate from the resumption of the accretion process due to the fragments from the collapsar falling back, such a relation may be related to a hyperaccretion disk around a precessed black hole (BH). We propose that the origin and time evolution of the flares, and the approximately symmetrical temporal structure and spectral evolution of the single flare, can be explained well by the jet precession model \citep{Liu2010}. In addition, the mass and spin of the BH can be constrained, which indicates a stellar-mass, fast rotating BH located in the center of GRB 130925A.

preprint2014arXiv

Variability of the giant X-ray bump in GRB 121027A and possible origin

The particular giant X-ray bump of GRB 121027A triggered by \emph{Swift} is quite different from the typical X-ray flares in gamma-ray bursts. There exhibit four parts of the observed structural variabilities in the rise and decay phase of the bump. Considering the quality of four parts of the data, we can only analyze the data from about 5300 s to about 6100 s in the bump using the stepwise filter correlation method (Gao et al. 2012), and find that the $86^{+5.9}_{-9.4}~\rm s$ periodic oscillation may exist, which is confirmed by the Lomb-Scargle method (Scargle 1982). Furthermore, a jet precession model (Liu et al. 2010) is proposed to account for such a variability.

preprint2013arXiv

A Complete Reference of the Analytical Synchrotron External Shock Models of Gamma-Ray Bursts

Gamma-ray bursts are most luminous explosions in the universe. Their ejecta are believed to move towards Earth with a relativistic speed. The interaction between this "relativistic jet" and a circum burst medium drives a pair of (forward and reverse) shocks. The electrons accelerated in these shocks radiate synchrotron emission to power the broad-band afterglow of GRBs. The external shock theory is an elegant theory, since it invokes a limit number of model parameters, and has well predicted spectral and temporal properties. On the other hand, depending on many factors (e.g. the energy content, ambient density profile, collimation of the ejecta, forward vs. reverse shock dynamics, and synchrotron spectral regimes), there is a wide variety of the models. These models have distinct predictions on the afterglow decaying indices, the spectral indices, and the relations between them (the so-called "closure relations"), which have been widely used to interpret the rich multi-wavelength afterglow observations. This review article provides a complete reference of all the analytical synchrotron external shock afterglow models by deriving the temporal and spectral indices of all the models in all spectral regimes, including some regimes that have not been published before. The review article is designated to serve as a useful tool for afterglow observers to quickly identify relevant models to interpret their data. The limitations of the analytical models are reviewed, with a list of situations summarized when numerical treatments are needed.

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 supra-massive magnetar central engine for short GRB 130603B

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

preprint2013arXiv

Bright broad-band afterglows of gravitational wave bursts from mergers of binary neutron stars

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

preprint2013arXiv

Compton Scattering of Self-Absorbed Synchrotron Emission

Synchrotron self-Compton (SSC) scattering is an important emission mechanism in many astronomical sources, such as gamma-ray bursts (GRBs) and active galactic nuclei (AGNs). We give a complete presentation of the analytical approximations for the Compton scattering of synchrotron emission with both weak and strong synchrotron self-absorption. All possible orders of the characteristic synchrotron spectral breaks ($ν_{\rm a}$, $ν_{\rm m}$, and $ν_{\rm c}$) are studied. In the weak self-absorption regime, i.e., $ν_{\rm a} < ν_c$, the electron energy distribution is not modified by the self-absorption process. The shape of the SSC component broadly resembles that of synchrotron, but with the following features: The SSC flux increases linearly with frequency up to the SSC break frequency corresponding to the self-absorption frequency $ν_{\rm a}$; and the presence of a logarithmic term in the high-frequency range of the SSC spectra makes it harder than the power-law approximation. In the strong absorption regime, i.e. $ν_{\rm a} > ν_{\rm c}$, heating of low energy electrons due to synchrotron absorption leads to pile-up of electrons, and form a thermal component besides the broken power-law component. This leads to two-component (thermal + non-thermal) spectra for both the synchrotron and SSC spectral components. For $ν_{\rm c} < ν_{\rm a} < ν_{\rm m}$, the spectrum is thermal (non-thermal) -dominated if $ν_a > \sqrt{ν_m ν_c}$ ($ν_a < \sqrt{ν_m ν_c}$). Similar to the weak-absorption regime, the SSC spectral component is broader than the simple broken power law approximation. We derive the critical condition for strong absorption (electron pile-up), and discuss a case of GRB reverse shock emission in a wind medium, which invokes $ν_{\rm a} > {\rm max} (ν_{\rm m}, ν_{\rm c})$.

preprint2013arXiv

Early Afterglows of Gamma-Ray Bursts in a Stratified Medium with a Power-Law Density Distribution

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

preprint2013arXiv

Giant X-ray Bump in GRB 121027A: Evidence for Fall-back Disk Accretion

A particularly interesting discovery in observations of GRB 121027A is that of a giant X-ray bump detected by the Swift/X-Ray Telescope. The X-ray afterglow re-brightens sharply at about 1000 s after the trigger by more than two orders of magnitude in less than 200 s. This X-ray bump lasts for more than 10 ks. It is quite different from typical X-ray flares. In this Letter we propose a fall-back accretion model to interpret this X-ray bump within the context of the collapse of a massive star for a long-duration gamma-ray burst. The required fall-back radius of about 3.5e10 cm and mass of about 0.9-2.6 solar masses imply that a significant part of the helium envelope should survive through the mass loss during the last stage of the massive progenitor of GRB 121027A.

preprint2013arXiv

Interpretation of the unprecedentedly long-lived high-energy emission of GRB 130427A

High energy photons (>100 MeV) are detected by the Fermi/LAT from GRB 130427A up to almost one day after the burst, with an extra hard spectral component being discovered in the high-energy afterglow. We show that this hard spectral component arises from afterglow synchrotron-self Compton emission. This scenario can explain the origin of >10 GeV photons detected up to ~30000s after the burst, which would be difficult to be explained by synchrotron radiation due to the limited maximum synchrotron photon energy. The lower energy multi-wavelength afterglow data can be fitted simultaneously by the afterglow synchrotron emission. The implication of detecting the SSC emission for the circumburst environment is discussed.

preprint2013arXiv

Is the late near-infrared bump in short-hard GRB 130603B due to the Li-Paczynski kilonova?

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

preprint2013arXiv

Possible High-Energy Neutrino and Photon Signals from Gravitational Wave Bursts due to Double Neutron Star Mergers

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

preprint2013arXiv

Signature of gravitational radiation in afterglow of short Gamma-Ray Bursts?

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

preprint2013arXiv

The Gamma-Ray Burst Hubble Diagram and Its Cosmological Implications

In this paper, we continue to build support for the proposal to use gamma-ray bursts (GRBs) as standard candles in constructing the Hubble Diagram at redshifts beyond the current reach of Type Ia supernova observations. We confirm that correlations among certain spectral and lightcurve features can indeed be used as luminosity indicators, and demonstrate from the most up-to-date GRB sample appropriate for this work that the $Λ$CDM model optimized with these data is characterized by parameter values consistent with those in the concordance model. Specifically, we find that $(Ω_m,Ω_Λ)\approx (0.25_{-0.06}^{+0.05}, 0.75_{-0.05}^{+0.06})$, which are consistent, to within $1σ$, with $(0.29,0.71)$ obtained from the 9-yr WMAP data. We also carry out a comparative analysis between $Λ$CDM and the $R_{\rm h}=ct$ Universe and find that the optimal $Λ$CDM model fits the GRB Hubble Diagram with a reduced $χ^2_{\rm dof}\approx 2.26$, whereas the fit using $R_{\rm h}=ct$ results in a $χ^2_{\rm dof}\approx 2.14$. In both cases, about 20% of the events lie at least $2σ$ away from the best-fit curves, suggesting that either some contamination by non-standard GRB luminosities is unavoidable, or that the errors and intrinsic scatter associated with the data are being underestimated. With these optimized fits, we use three statistical tools---the Akaike Information Criterion (AIC), the Kullback Information Criterion (KIC), and the Bayes Information Criterion (BIC)---to show that, based on the GRB Hubble Diagram, the likelihood of $R_{\rm h}=ct$ being closer to the correct model is $\sim 85-96%$, compared to $\sim 4-15%$ for $Λ$CDM.

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

Is There a Relation between Duration and E_iso in Gamma-Ray Bursts?

The system of accretion disk and black hole is usually considered as the central engine of Gamma-ray Bursts (GRBs). It is usually thought that the disk in the central engine of GRBs is the advection-dominated accretion disk, which is developed from a massive (mass M_disk) torus at radius r_disk. We find a positive correlation between the isotropic gamma-ray energy E_iso and duration (so-called T_90) for GRBs. We interpret this correlation within the advection-dominated accretion disk model, associating E_iso and T_90 with M_ disk and viscous timescale respectively.

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.

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

A Photosphere-Internal Shock Model of Gamma-Ray Bursts: Case Studies of Fermi/LAT Bursts

Radially inhomogeneous gamma-ray burst (GRB) jets release variable photospheric emission and can have internal shocks occurring above the photosphere. We generically formulate a photospheric emission model of GRBs including Compton up-scattered photospheric (UP) emission off the electrons (and positrons) in the internal shocks, and find that the photospheric emission may correspond to the traditional (Band) component at <~1 MeV and the UP emission to the high-energy emission observed by Fermi/LAT for some GRBs at >~ 10 MeV. The two components can be separate in the spectrum in some cases or can mimic a smooth broad Band spectrum in other cases. We apply our formulation to the well-studied long and short LAT GRBs, GRB 080916C, GRB 090902B, and GRB 090510, and typically find reasonable parameters for fitting the time-binned spectra, although fine tuning of several parameters is required. The observed delays of the high-energy emission with respect to the MeV emission which are large compared to the variability times are unlikely to be due to simple kinematic effects of a non-evolving jet. These delays may be attributed to the temporal evolution of the physical parameters of the jet, and thus the delay timescales could provide a potential tool for investigating the structures of GRB jets themselves and their progenitors. The difference of the delay timescales of long and short GRBs inferred from the Fermi data might be due to the differences in the progenitors of long and short GRBs. Some other properties and consequences of this model are discussed, including temporal correlations among the prompt optical, the soft X-ray, and the distinct high-energy component as well as the Band component.

preprint2011arXiv

A statistical model for the $γ$-ray variability of the Crab nebula

A statistical scenario is proposed to explain the $γ$-ray variability and flares of the Crab nebula, which were observed recently by the Fermi/LAT. In this scenario electrons are accelerated in a series of knots, whose sizes follow a power-law distribution. These knots presumably move outwards from the pulsar and have a distribution in the Doppler boost factor. The maximal electron energy is assumed to be proportional to the size of the knot. Fluctuations at the highest energy end of the overall electron distribution will result in variable $γ$-ray emission via the synchrotron process in the $\sim 100$ MeV range. Since highly boosted larger knots are rarer than smaller knots, the model predicts that the variability of the synchrotron emission increases with the photon energy. We realize such a scenario with a Monte-Carlo simulation and find that the model can reproduce both the two $γ$-ray flares over a period of $\sim$ year and the monthly scale $γ$-ray flux fluctuations as observed by the Fermi/LAT. The observed $γ$-ray spectra in both the steady and flaring states are also well reproduced.

preprint2011arXiv

A statistical model to explain the gamma-ray variability and flares of the Crab nebula

Recently the AGILE and Fermi/LAT detectors uncovered giant $γ$-ray flares from the Crab nebula. The duration of these flares is a few days. The Fermi/LAT data with monthly time binning further showed significant variability of the synchrotron tail of the emission, while the inverse Compton component was stable. The simultaneous or follow-up observations in X-ray, optical, infrared and radio bands did not find significant flux variation. Based on these observations, we propose that the $γ$-ray variability and flares are due to statistical fluctuations of knots that can accelerate electrons to $\sim$PeV energies. The maximum achievable energy of electrons is adopted to be proportional to the size of the knot, which is assumed to follow a power-law distribution. Thus the low energy electron flux will be stable due to the large number of small knots, while the high energy electron flux may experience large fluctuations. Monte Carlo realization of such a picture can reproduce the observational data quite well given proper model parameters.

preprint2011arXiv

On the High Energy Emission of the Short GRB 090510

Long-lived high-energy (>100MeV) emission, a common feature of most Fermi-LAT detected gamma-ray burst, is detected up to \sim 10^2 s in the short GRB 090510. We study the origin of this long-lived high-energy emission, using broad-band observations including X-ray and optical data. We confirm that the late > 100 MeV, X-ray and optical emission can be naturally explained via synchrotron emission from an adiabatic forward shock propagating into a homogeneous ambient medium with low number density. The Klein-Nishina effects are found to be significant, and effects due to jet spreading and magnetic field amplification in the shock appear to be required. Under the constraints from the low-energy observations, the adiabatic forward shock synchrotron emission is consistent with the later-time (t>2s) high-energy emission, but falls below the early-time (t < 2s) high energy emission. Thus we argue that an extra high energy component is needed at early times. A standard reverse shock origin is found to be inconsistent with this extra component. Therefore, we attribute the early part of the high-energy emission (t< 2s) to the prompt component, and the long-lived high energy emission (t>2s) to the adiabatic forward shock synchrotron afterglow radiation. This avoids the requirement for an extremely high initial Lorentz factor.

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.

preprint2009arXiv

Diffuse high energy neutrinos and cosmic rays from hyperflares of soft-gamma repeaters

We calculate the diffuse high energy (TeV - PeV) neutrino emission from hyperflares of Soft-Gamma Repeaters (SGRs), like the hyperflare risen from \astrobj{SGR 1806-20} on December 27 of 2004, within the framework of the fireball model. The fireball model for gamma-ray bursts (GRBs) can explain well the main features of this hyperflare and the subsequent multi-frequency terglow emission. The expected rate, $\sim 20-100$ Gpc$^{-3}$day$^{-1}$, of such hyperflares is well in excess of the GRBs rate. Our result shows that the contribution to the diffuse TeV-PeV neutrino background from such hyperflares is less than 10% of the contribution from GRBs. We also discuss the high energy cosmic rays (CRs) from these sources.

preprint2009arXiv

Early afterglows from radially structured outflows and the application to X-ray shallow decays

In the fireball model, it is more physically realistic that gamma-ray burst (GRB) ejecta have a range of bulk Lorentz factors (assuming $M\propto Γ^{-s}$). The low Lorentz factor part of the ejecta will catch up with the high Lorentz factor part when the latter is decelerated by the surrounding medium to a comparable Lorentz factor. Such a process will develop a long-lasting weak reverse shock until the whole ejecta are shocked. Meanwhile, the forward shocked materials are gradually supplied with energy from the ejecta that are catching-up, and thus the temporal decay of the forward shock emission will be slower than that without an energy supply. However, the reverse shock may be strong. Here, we extend the standard reverse-forward shock model to the case of radially nonuniform ejecta. We show that this process can be classified into two cases: the thick shell case and the thin shell case. In the thin shell case, the reverse shock is weak and the temporal scaling law of the afterglow is the same as that in Sari & Mesz (2000). However, in the thick shell case, the reverse shock is strong and thus its emission dominates the afterglow in the high energy band. Our results also show slower decaying behavior of the afterglow due to the energy supply by low Lorentz factor materials, which may help the understanding of the plateau observed in the early optical and X-ray afterglows.

preprint2009arXiv

Klein-Nishina effects on the high-energy afterglow emission of gamma-ray bursts

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