Source author record

Z. G. Dai

Z. G. Dai appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.HE astro-ph astro-ph.CO astro-ph.GA astro-ph.SR hep-ph astro-ph.EP nlin.AO physics.geo-ph

Catalog footprint

What is connected

87works

9topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2026arXiv

An Enigmatic PeVatron in an Area around HII Region G35.6$-$0.5

Identifying Galactic PeVatrons (PeV particle accelerators) from the ultra-high-energy (UHE, >100 TeV) $γ$-ray sources plays a crucial role in revealing the origin of Galactic cosmic rays. The UHE source 1LHAASO J1857+0203u is suggested to be associated with HESS J1858+020, which may be attributed to the possible PeVatron candidate supernova remnant (SNR) G35.6$-$0.4 or HII region G35.6$-$0.5. We perform detailed analysis on the very-high-energy and UHE $γ$-ray emissions towards this region with data from the Large High Altitude Air Shower Observatory (LHAASO). 1LHAASO J1857+0203u is detected with a significance of 11.6$σ$ above 100 TeV, indicating the presence of a PeVatron. It has an extension of $\sim 0.18^\circ$ with a power-law (PL) spectral index of $\sim$2.5 in 1-25 TeV and a point-like emission with a PL spectral index of $\sim$3.2 above 25 TeV. Using the archival CO and HI data, we identify some molecular and atomic clouds that may be associated with the TeV $γ$-ray emissions. Our modelling indicates that the TeV $γ$-ray emissions are unlikely to arise from the clouds illuminated by the protons that escaped from SNR G35.6$-$0.4. In the scenario that HII region G35.6$-$0.5 could accelerate particles to the UHE band, the observed GeV-TeV $γ$-ray emission could be well explained by a hadronic model with a PL spectral index of $\sim$2.0 and cutoff energy of $\sim$450 TeV. However, an evolved pulsar wind nebula origin cannot be ruled out.

preprint2026arXiv

An Ultrahigh-energy $γ$-ray Bubble Powered by a Super PeVatron

We report the detection of a $γ$-ray bubble spanning at least 100$\rm deg^2$ in ultra high energy (UHE) up to a few PeV in the direction of the star-forming region Cygnus X, implying the presence Super PeVatron(s) accelerating protons to at least 10 PeV. A log-parabola form with the photon index $Γ(E) = (2.71 \pm 0.02) + (0.11 \pm 0.02) \times \log_{10} (E/10 \ {\rm TeV})$ is found fitting the gamma-ray energy spectrum of the bubble well. UHE sources, `hot spots' correlated with very massive molecular clouds, and a quasi-spherical amorphous $γ$-ray emitter with a sharp central brightening are observed in the bubble. In the core of $\sim 0.5^{\circ}$, spatially associating with a region containing massive OB association (Cygnus OB2) and a microquasar (Cygnus X-3), as well as previously reported multi-TeV sources, an enhanced concentration of UHE $γ$-rays are observed with 2 photons at energies above 1 PeV. The general feature of the bubble, the morphology and the energy spectrum, are reasonably reproduced by the assumption of a particle accelerator in the core, continuously injecting protons into the ambient medium.

preprint2026arXiv

Constraining the Cosmic-ray Energy Based on Observations of Nearby Galaxy Clusters by LHAASO

Galaxy clusters act as reservoirs of high-energy cosmic rays (CRs). As CRs propagate through the intracluster medium, they generate diffuse $γ$-rays detectable by arrays such as LHAASO. These $γ$-rays result from proton-proton ($pp$) collisions of very high-energy cosmic rays (VHECRs) or inverse Compton (IC) scattering of positron-electron pairs created by $pγ$ interactions of ultra-high-energy cosmic rays (UHECRs). We analyzed diffuse $γ$-ray emission from the Coma, Perseus, and Virgo clusters using LHAASO data. Diffuse emission was modeled as a disk of radius $R_{500}$ for each cluster while accounting for point sources. No significant diffuse emission was detected, yielding 95\% confidence level (C.L.) upper limits on the $γ$-ray flux: for WCDA (1-25~TeV) and KM2A ($>25$~TeV), less than $(49.4, 13.7, 54.0)$ and $(1.34, 1.14, 0.40) \times 10^{-14}$~ph~cm$^{-2}$~s$^{-1}$ for Coma, Perseus, and Virgo, respectively. The $γ$-ray upper limits can be used to derive model-independent constraints on the integral energy of CRp above 10~TeV (corresponding to the LHAASO observational range $>1$~TeV under the $pp$ scenario) to be less than $(1.96, 0.59, 0.08) \times 10^{61}$~erg. The absence of detectable annuli/ring-like structures, indicative of cluster accretion or merging shocks, imposes further constraints on models in which the UHECRs are accelerated in the merging shocks of galaxy clusters.

preprint2026arXiv

Constraints on heavy decaying dark matter from 570 days of LHAASO observations

The Kilometer Square Array~(KM2A) of the Large High Altitude Air Shower Observatory (LHAASO) aims at surveying the northern gamma-ray sky at energies above 10 TeV with unprecedented sensitivity. Gamma-ray observations have long been one of the most powerful tools for dark matter searches, as e.g., high-energy gamma-rays could be produced by the decays of heavy dark matter particles. In this letter, we present the first dark matter analysis with LHAASO-KM2A, using the first 340~days of data from 1/2-KM2A and 230~days of data from 3/4-KM2A. Several regions of interest are used to search for a signal and account for the residual cosmic-ray background after gamma/hadron separation. We find no excess of dark matter signals, and thus place some of the strongest gamma-ray constraints on the lifetime of heavy dark matter particles with mass between 10^5 and 10^9~GeV. Our results with LHAASO are robust, and have important implications for dark matter interpretations of the diffuse astrophysical high-energy neutrino emission.

preprint2026arXiv

Deep view of Composite SNR CTA1 with LHAASO in $γ$-rays up to 300 TeV

The ultra-high-energy (UHE) gamma-ray source 1LHAASO J0007+7303u is positionally associated with the composite SNR CTA1 that is located at high Galactic Latitude $b\approx 10.5^\circ$. This provides a rare opportunity to spatially resolve the component of the pulsar wind nebula (PWN) and supernova remnant (SNR) at UHE. This paper conducted a dedicated data analysis of 1LHAASO J0007+7303u using the data collected from December 2019 to July 2023. This source is well detected with significances of 21$σ$ and 17$σ$ at 8$-$100 TeV and $>$100 TeV, respectively. The corresponding extensions are determined to be 0.23$^{\circ}\pm$0.03$^{\circ}$ and 0.17$^{\circ}\pm$0.03$^{\circ}$. The emission is proposed to originate from the relativistic electrons and positrons accelerated within the PWN of PSR J0007+7303. The energy spectrum is well described by a power-law with an exponential cutoff function $dN/dE = (42.4\pm4.1)(\frac{E}{20\rm\ TeV})^{-2.31\pm0.11}\exp(-\frac{E}{110\pm25\rm\ TeV})$ $\rm\ TeV^{-1}\ cm^{-2}\ s^{-1}$in the energy range from 8 TeV to 300 TeV, implying a steady-state parent electron spectrum $dN_e/dE_e\propto (\frac{E_e}{100\rm\ TeV})^{-3.13\pm0.16}\exp[(\frac{-E_e}{373\pm70\rm\ TeV})^2]$ at energies above $\approx 50 \rm\ TeV$. The cutoff energy of the electron spectrum is roughly equal to the expected current maximum energy of particles accelerated at the PWN terminal shock. Combining the X-ray and gamma-ray emission, the current space-averaged magnetic field can be limited to $\approx 4.5\rm\ μG$. To satisfy the multi-wavelength spectrum and the $γ$-ray extensions, the transport of relativistic particles within the PWN is likely dominated by the advection process under the free-expansion phase assumption.

preprint2026arXiv

Discovery of a new $γ$-ray source LHAASO J0341+5258 with emission up to 200TeV

We report the discovery of a new unidentified extended $γ$-ray source in the Galactic plane named LHAASO J0341+5258 with a pre-trial significance of 8.2 standard deviations above 25 TeV. The best fit position is R.A.$=55.34^{\circ}\pm0.11^{\circ}$ and Dec$=52.97^{\circ}\pm0.07^{\circ}$. The angular size of LHAASO J0341+5258 is $0.29^\circ \pm 0.06^\circ_{stat} \pm0.02^\circ_{sys}$. The flux above 25 TeV is about $20\%$ of the flux of Crab Nebula. Although a power-law fit of the spectrum from 10 TeV to 200 TeV with the photon index $α=2.98 \pm 0.19_{stat} \pm 0.02_{sys}$ is not excluded, the LHAASO data together with the flux upper limit at 10 GeV set by the Fermi LAT observation, indicate a noticeable steepening of an initially hard power-law spectrum %($α\leq 1.75$) spectrum with a cutoff at $\approx 50$ TeV. We briefly discuss the origin of UHE gamma-rays. The lack of an energetic pulsar and a young SNR inside or in the vicinity of LHAASO J0341+5258 challenge, but do not exclude both the leptonic and hadronic scenarios of gamma-ray production.

preprint2026arXiv

Discovery of the Ultra-high energy gamma-ray source LHAASO J2108+5157

We report the discovery of a UHE gamma-ray source, LHAASO J2108+5157, by analyzing the LHAASO-KM2A data of 308.33 live days. Significant excess of gamma-ray induced showers is observed in both energy bands of 25-100 TeV and $\gt$100 TeV with 9.5 sigma and 8.5 sigma, respectively. This source is not significantly favored as an extensive source with the angular extension smaller than the point-spread function of KM2A. The measured energy spectrum from 20 to 200 TeV can be approximately described by a power-law function with an index of -2.83$\pm$ 0.18stat. A harder spectrum is demanded at lower energies considering the flux upper limit set by Fermi-LAT observations. The position of the gamma-ray emission is correlated with a giant molecular cloud, which favors a hadronic origin. No obvious counterparts have been found, deeper multiwavelength observations will help to shed new light on this intriguing UHE source.

preprint2026arXiv

Energy calibration of LHAASO-KM2A using the cosmic ray Moon shadow

We present a precise measurement of the westward, rigidity-dependent shift of the Moon's shadow using three and a half years of cosmic-ray data collected by the Kilometer Square Array (KM2A) of the Large High Altitude Air Shower Observatory (LHAASO). These measurements enable us to calibrate the detector energy response in the range 20-260 TeV, with results showing excellent agreement with the response derived from Monte Carlo (MC) simulations of the KM2A detector. We also measure a best-fit parameter $ε= 0.015 \pm 0.08$, corresponding to a 95% confidence interval of [-14%, +17%] for the energy-scale estimation. This result establishes the exceptional accuracy of the KM2A-MC in simulating the detector's response within this energy range.

preprint2026arXiv

Energy-Dependent Shifts of Medium-Scale Anisotropies in Very-High-Energy Cosmic Rays Observed by LHAASO-KM2A

Small deviations from isotropy in the arrival directions of Galactic cosmic rays serve as a unique probe of the local magnetic environment. In this Letter, we report observations of medium-scale anisotropies (MSA) at energies above 10 TeV using the LHAASO-KM2A array. Our analysis identifies four regions of excess and four regions of deficit, each spanning angular scales of approximately ten degrees. Crucially, we detect significant energy-dependent shifts in the centroids of two excess regions: Region B and the newly identified Region $\mathrm{\widetilde{D}}$. We also characterize the energy evolution of the fractional relative intensity across both excess and deficit regions. These findings imply that the observed anisotropies are shaped by the specific realization of the local turbulent magnetic field within the cosmic ray scattering length. Such energy-dependent behaviors impose strict constraints on local turbulence models and cosmic ray propagation theories.

preprint2026arXiv

Evidence for particle acceleration approaching PeV energies in the W51 complex

The $γ$-ray emission from the W51 complex is widely acknowledged to be attributed to the interaction between the cosmic rays (CRs) accelerated by the shock of supernova remnant (SNR) W51C and the dense molecular clouds in the adjacent star-forming region, W51B. However, the maximum acceleration capability of W51C for CRs remains elusive. Based on observations conducted with the Large High Altitude Air Shower Observatory (LHAASO), we report a significant detection of $γ$ rays emanating from the W51 complex, with energies from 2 TeV to 200 TeV. The LHAASO measurements, for the first time, extend the $γ$-ray emission from the W51 complex beyond 100 TeV and reveal a significant spectrum bending at tens of TeV. By combining the ``$π^0$-decay bump" featured data from Fermi-LAT, the broadband $γ$-ray spectrum of the W51 region can be well-characterized by a simple pp-collision model. The observed spectral bending feature suggests an exponential cutoff at $\sim400$~TeV or a power-law break at $\sim200$~TeV in the CR proton spectrum, most likely providing the first evidence of SNRs serving as CR accelerators approaching the PeV regime. Additionally, two young star clusters within W51B could also be theoretically viable to produce the most energetic $γ$ rays observed by LHAASO. Our findings strongly support the presence of extreme CR accelerators within the W51 complex and provide new insights into the origin of Galactic CRs.

preprint2026arXiv

Exploring Lorentz Invariance Violation from Ultra-high-energy Gamma Rays Observed by LHAASO

Recently the LHAASO Collaboration published the detection of 12 ultra-high-energy gamma-ray sources above 100 TeV, with the highest energy photon reaching 1.4 PeV. The first detection of PeV gamma rays from astrophysical sources may provide a very sensitive probe of the effect of the Lorentz invariance violation (LIV), which results in decay of high-energy gamma rays in the superluminal scenario and hence a sharp cutoff of the energy spectrum. Two highest energy sources are studied in this work. No signature of the existence of LIV is found in their energy spectra, and the lower limits on the LIV energy scale are derived. Our results show that the first-order LIV energy scale should be higher than about 10^5 times the Planck scale M_{pl} and that the second-order LIV scale is >10^{-3}M_{pl}. Both limits improve by at least one order of magnitude the previous results.

preprint2026arXiv

Extended Very-High-Energy Gamma-Ray Emission Surrounding PSR J0622 + 3749 Observed by LHAASO-KM2A

We report the discovery of an extended very-high-energy (VHE) gamma-ray source around the location of the middle-aged (207.8 kyr) pulsar PSR J0622+3749 with the Large High Altitude Air Shower Observatory (LHAASO). The source is detected with a significance of $8.2σ$ for $E>25$~TeV assuming a Gaussian template. The best-fit location is (R.A., Dec.)$=(95^{\circ}\!.47\pm0^{\circ}\!.11,\,37^{\circ}\!.92 \pm0^{\circ}\!.09)$, and the extension is $0^{\circ}\!.40\pm0^{\circ}\!.07$. The energy spectrum can be described by a power-law spectrum with an index of ${-2.92 \pm 0.17_{\rm stat} \pm 0.02_{\rm sys} }$. No clear extended multi-wavelength counterpart of the LHAASO source has been found from the radio to sub-TeV bands. The LHAASO observations are consistent with the scenario that VHE electrons escaped from the pulsar, diffused in the interstellar medium, and scattered the interstellar radiation field. If interpreted as the pulsar halo scenario, the diffusion coefficient, inferred for electrons with median energies of $\sim160$~TeV, is consistent with those obtained from the extended halos around Geminga and Monogem and much smaller than that derived from cosmic ray secondaries. The LHAASO discovery of this source thus likely enriches the class of so-called pulsar halos and confirms that high-energy particles generally diffuse very slowly in the disturbed medium around pulsars.

preprint2026arXiv

LHAASO Detection of Ultra-High-Energy Gamma-Ray Emission toward the Giant Molecular Clouds

The $γ$-ray from Giant molecular clouds (GMCs) is regarded as the most ideal tool to perform in-situ measurement of cosmic ray (CR) density and spectra in our Galaxy. We report the first detection of $γ$-ray emissions in the very-high-energy (VHE) domain from the five nearby GMCs with a stacking analysis based on a 4.5-year $γ$-ray observation with the Large High Altitude Air Shower Observatory (LHAASO) experiment. The spectral energy distributions derived from the GMCs are consistent with the expected $γ$-ray flux produced via CR interacting with the ISM in the energy interval 1 - 100 $~\rm$ TeV. In addition, we investigate the presence of the CR spectral `knee' by introducing a spectral break in the $γ$-ray data. While no significant evidence for the CR knee is found, the current KM2A measurements from GMCs strongly favor a proton CR knee located above 0.9$~\rm$ PeV, which is consistent with the latest measurement of the CR spectrum by ground-based experiments.

preprint2026arXiv

Measurement of attenuation length of the muon content in extensive air showers from 0.3 to 30 PeV with LHAASO

The attenuation length of the muon content in extensive air showers provides important information regarding the generation and development of air showers. This information can be used not only to improve the description of such showers but also to test fundamental models of hadronic interactions. Using data from the LHAASO-KM2A experiment, the development of the muon content in high-energy air showers was studied. The attenuation length of muon content in the air showers was measured from experimental data in the energy range from 0.3 to 30 PeV using the constant intensity cut method. By comparing the attenuation length of the muon content with predictions from high-energy hadronic interaction models (QGSJET-II-04, SIBYLL 2.3d, and EPOS-LHC), it is evident that LHAASO results are significantly shorter than those predicted by the first two models (QGSJET-II-04 and SIBYLL 2.3d) but relatively close to those predicted by the third model (EPOS-LHC). Thus, the LHAASO data favor the EPOS-LHC model over the other two models. The three interaction models confirmed an increasing trend in the attenuation length as the cosmic-ray energy increases.

preprint2026arXiv

Measurement of Very-high-energy Diffuse Gamma-ray Emissions from the Galactic Plane with LHAASO-WCDA

The diffuse Galactic gamma-ray emission is a very important tool used to study the propagation and interaction of cosmic rays in the Milky Way. In this work, we report the measurements of the diffuse emission from the Galactic plane, covering Galactic longitudes from $15^{\circ}$ to $235^{\circ}$ and latitudes from $-5^{\circ}$ to $+5^{\circ}$, in an energy range of 1 TeV to 25 TeV, with the Water Cherenkov Detector Array (WCDA) of the Large High Altitude Air Shower Observatory (LHAASO). After masking the sky regions of known sources, the diffuse emission is detected with $24.6σ$ and $9.1σ$ significance in the inner Galactic plane and outer Galactic plane, respectively. The WCDA spectra in both regions can be well described by a power-law function, with spectral indices of $-2.67\pm0.05_{\rm stat}$ in the inner region and $-2.83\pm0.19_{\rm stat}$ in the outer region, respectively. Combined with the Square Kilometer Array (KM2A) measurements at higher energies, a clear softening of the spectrum is found in the inner region, with change of spectral indices by $\sim0.5$ at a break energy around $30$ TeV. The fluxes of the diffuse emission are higher by a factor of $1.5-2.7$ than the model prediction assuming local CR spectra and the gas column density, which are consistent with those measured by the KM2A. Along Galactic longitude, the spatial distribution of the diffuse emission shows deviation from that of the gas column density. The spectral shape of the diffuse emission are possibly variation in different longitude region. The WCDA measurements bridge the gap between the low-energy measurements by space detectors and the ultra-high-energy observations by LHAASO-KM2A and other experiments. These results suggest that improved modeling of the wide-band diffuse emission is required.

preprint2026arXiv

Transient Large-Scale Anisotropy in TeV Cosmic Rays due to an Interplanetary Coronal Mass Ejection

Large- or medium-scale cosmic ray anisotropy at TeV energies has not previously been confirmed to vary with time. Transient anisotropy changes have been observed below 150 GeV, especially near the passage of an interplanetary shock and coronal mass ejection containing a magnetic flux rope ejected by a solar storm, which can trigger a geomagnetic storm with practical consequences. In such events, cosmic rays provide remote sensing of the magnetic field properties. Here we report the observation of transient large-scale anisotropy in TeV cosmic ray ions using data from the Large High Altitude Air Shower Observatory (LHAASO). We analyze hourly skymaps of the transient cosmic ray intensity excess or deficit, the gradient of which indicates the direction and magnitude of transient large-scale anisotropy across the field of view. We observe enhanced anisotropy above typical hourly fluctuations with $>$5$σ$ significance during some hours of November 4, 2021, in separate data sets for four primary cosmic ray energy ranges of median energy from $E$=0.7 to 3.1 TeV. The gradient varies with energy as $E^γ$, where $γ\approx-0.5$. At a median energy $\leq$1.0 TeV, this gradient corresponds to a dipole anisotropy of at least 1\%, or possibly a weaker anisotropy of higher order. This new type of observation opens the opportunity to study interplanetary magnetic structures using air shower arrays around the world, complementing existing in situ and remote measurements of plasma properties.

preprint2025arXiv

Observation of the $γ$-ray Emission from W43 with LHAASO

In this paper, we report the detection of the very-high-energy (VHE, $ 100{\rm\ GeV} < E < 100{\rm\ TeV} $) and ultra-high-energy (UHE, $E > 100\rm\ TeV$) $γ$-ray emissions from the direction of the young star-forming region W43, observed by the Large High Altitude Air Shower Observation (LHAASO). The extended $γ$-ray source was detected with a significance of ${\sim}16\,σ$ by KM2A and ${\sim}17\,σ$ by WCDA, respectively. The angular extension of this $γ$-ray source is about 0.5 degrees, corresponding to a physical size of about 50 pc. We discuss the origin of the $γ$-ray emission and possible cosmic ray acceleration in the W43 region using multi-wavelength data. Our findings suggest that W43 is likely another young star cluster capable of accelerating cosmic rays (CRs) to at least several hundred TeV.

preprint2025arXiv

Study of Ultra-High-Energy Gamma-Ray Source 1LHAASO J0056+6346u and Its Possible Origins

We report a dedicated study of the newly discovered extended UHE $γ$-ray source 1LHAASO J0056+6346u. Analyzing 979 days of LHAASO-WCDA data and 1389 days of LHAASO-KM2A data, we observed a significant excess of $γ$-ray events with both WCDA and KM2A. Assuming a point power-law source with a fixed spectral index, the significance maps reveal excesses of ${\sim}12.65\,σ$, ${\sim}22.18\,σ$, and ${\sim}10.24\,σ$ in the energy ranges of 1--25 TeV, 25--100 TeV, and $> 100$ TeV, respectively. We use a 3D likelihood algorithm to derive the morphological and spectral parameters, and the source is detected with significances of $12.65\,σ$ by WCDA and $25.27\,σ$ by KM2A. The best-fit positions derived from WCDA and KM2A data are (R.A. = $13.96^\circ\pm0.09^\circ$, Decl. = $63.92^\circ\pm0.05^\circ$) and (R.A. = $14.00^\circ\pm0.05^\circ$, Decl. = $63.79^\circ\pm0.02^\circ$), respectively. The angular size ($r_{39}$) of 1LHAASO J0056+6346u is $0.34^\circ\pm0.04^\circ$ at 1--25 TeV and $0.24^\circ\pm0.02^\circ$ at $> 25$ TeV. The differential flux of this UHE $γ$-ray source can be described by an exponential cutoff power-law function: $(2.67\pm0.25) \times 10^{-15} (E/20\,\text{TeV})^{-1.97\pm0.10} e^{-E/(55.1\pm7.2)\,\text{TeV}} \,\text{TeV}^{-1}\,\text{cm}^{-2}\,\text{s}^{-1}$. To explore potential sources of $γ$-ray emission, we investigated the gas distribution around 1LHAASO J0056+6346u. 1LHAASO J0056+6346u is likely to be a TeV PWN powered by an unknown pulsar, which would naturally explain both its spatial and spectral properties. Another explanation is that this UHE $γ$-ray source might be associated with gas content illuminated by a nearby CR accelerator, possibly the SNR candidate G124.0+1.4.

preprint2022arXiv

Flaring-associated Complex Dynamics in Two M-dwarfs Revealed by Fast, Time-resolved Spectroscopy

Habitability of an exoplanet is believed to be profoundly affected by activities of the host stars, although the related coronal mass ejections (CMEs) are still rarely detected in solar-like and late-type stars. We here report an observational study on flares of two M-dwarfs triggered by the high-cadence survey performed by the Ground Wide-angle Camera system. In both events, the fast, time-resolved spectroscopy enables us to identify symmetric broad H$α$ emission with not only a nearly zero bulk velocity, but also a large projected maximum velocity as high as $\sim700-800\ \mathrm{km\ s^{-1}}$. This broadening could be resulted from either Stark (pressure) effect or a flaring-associated CME at stellar limb. In the context of the CME scenario, the CME mass is estimated to be $\sim4\times10^{18}$ g and $2\times10^{19}$ g. In addition, our spectral analysis reveals a temporal variation of the line center of the narrow H$α$ emission in both events. The variation amplitudes are at tens of $\mathrm{km\ s^{-1}}$, which could be ascribed to the chromospheric evaporation in one event, and to a binary scenario in the other one. With the total flaring energy determined from our photometric monitor, we show a reinforced trend in which larger the flaring energy, higher the CME mass is.

preprint2022arXiv

Quasi-periodic Eruptions from Helium Envelope of Hydrogen-deficient Stars Stripped by Supermassive Black Holes

Quasi-periodic eruptions (QPEs), which are a new kind of X-ray bursts with a recurrence time of several hours, have been detected from supermassive black holes (SMBHs) in galactic nuclei. Recently, the two QPEs discovered by the \textit{eROSITA} show asymmetric light curves with a fast rise and a slow decline. Current models cannot explain the observational characteristics of QPEs. Here we show that QPEs can be generated from the Roche lobe overflows at each periapsis passage of an evolved star orbiting an SMBH. The properties of the companion stars are constrained via analytic estimations. We find that hydrogen-deficient post-AGB stars are promising candidates. Modules for Experiments in Stellar Astrophysics (MESA) stellar evolution code is used to construct the hydrogen-deficient stars which can fulfill the requirements, as obtained through analytical estimates, to produce the properties of QPEs, including the fast-rise and slow-decay light curves, periods, energetics, and rates. Furthermore, the extreme mass ratio $\sim 10^5$ between the SMBH and the donor will lead to a phenomenon called extreme mass-ratio inspiral (EMRI), producing millihertz gravitational waves. These QPEs would be detected as EMRI sources with electromagnetic counterparts for space-based GW detectors, such as Laser Interferometer Space Antenna (LISA) and Tianqin. They would provide a new way to measure the Hubble constant and further test the Hubble constant tension.

preprint2021arXiv

Standardized long gamma-ray bursts as a cosmic distance indicator

Gamma-ray bursts (GRBs) are the most luminous explosions and can be detectable out to the edge of Universe. It has long been thought they can extend the Hubble diagram to very high redshifts. Several correlations between temporal or spectral properties and GRB luminosities have been proposed to make GRBs cosmological tools. However, those correlations cannot be properly standardized. In this paper, we select a long GRB sample with X-ray plateau phases produced by electromagnetic dipole emissions from central new-born magnetars. A tight correlation is found between the plateau luminosity and the end time of the plateau in X-ray afterglows out to the redshift $z=5.91$. We standardize these long GRBs X-ray light curves to a universal behavior by this correlation for the first time, with a luminosity dispersion of 0.5 dex. The derived distance-redshift relation of GRBs is in agreement with the standard $Λ$CDM model both at low and high redshifts. The evidence of accelerating universe from this GRB sample is $3σ$, which is the highest statistical significance from GRBs to date.

preprint2020arXiv

A Magnetar-Asteroid Impact Model for FRB 200428 Associated with an X-ray Burst from SGR 1935+2154

Very recently, an extremely bright fast radio burst (FRB) 200428 with two sub-millisecond pulses was discovered to come from the direction of the Galactic magnetar SGR 1935+2154, and an X-ray burst (XRB) counterpart was detected simultaneously. These observations favor magnetar-based interior-driven models. In this Letter, we propose a different model for FRB 200428 associated with an XRB from SGR 1935+2154, in which a magnetar with high proper velocity encounters an asteroid of mass $\sim10^{20}\,$g. This infalling asteroid in the stellar gravitational field is first possibly disrupted tidally into a great number of fragments at radius $\sim {\rm a\,\,few}$ times $10^{10}\,$cm, and then slowed around the Alfv$\acute{\rm e}$n radius by an ultra-strong magnetic field and in the meantime two major fragments of mass $\sim 10^{17}\,$g that cross magnetic field lines produce two pulses of FRB 200428. The whole asteroid is eventually accreted onto the poles along magnetic field lines, impacting the stellar surface, creating a photon-e$^\pm$ pair fireball trapped initially in the stellar magnetosphere, and further leading to an XRB. We show that this gravitationally-powered model can interpret all of the observed features self-consistently.

preprint2020arXiv

Fast radio bursts from activities of neutron stars newborn in BNS mergers: offset, birth rate and observational properties

Young neutron stars (NSs) born in core-collapse explosions are promising candidates for the central engines of fast radio bursts (FRBs), since the first localized repeating burst FRB 121102 happens in a star forming dwarf galaxy, which is similar to the host galaxies of superluminous supernovae (SLSNe) and long gamma-ray bursts (LGRBs). However, FRB 180924 and FRB 190523 are localized to massive galaxies with low rates of star formation, compared with the host of FRB 121102. Meanwhile, the offsets between the bursts and host centers are about 4 kpc and 29 kpc for FRB 180924 and FRB 190523, respectively. These properties of hosts are similar to short gamma-ray bursts \textbf{(SGRBs)}, which are produced by mergers of binary neutron star (BNS) or neutron star-black hole (NS-BH). Therefore, the NSs powering FRBs may be formed in BNS mergers. In this paper, we study the BNS merger rates, merger times, and predict their most likely merger locations for different types of host galaxies using population synthesis method. We find that the BNS merger channel is consistent with the recently reported offsets of FRB 180924 and FRB 190523. The offset distribution of short GRBs is well reproduced by population synthesis using galaxy model which is similar to GRB hosts. The event rate of FRBs (including non-repeating and repeating), is larger than those of BNS merger and short GRBs, which requires a large fraction of observed FRBs emitting several bursts. Using curvature radiation by bunches in NS magnetospheres, we also predict the observational properties of FRBs from BNS mergers, including the dispersion measure, and rotation measure. At late times ($t\geq1$yr), the contribution to dispersion measure and rotation measure from BNS merger ejecta could be neglected.

preprint2020arXiv

Galactic and cosmological fast radio bursts as scaled-up solar radio bursts

Fast radio bursts (FRBs) are bright milliseconds radio transients with large dispersion measures. Recently, FRB 200428 was detected in temporal coincidence with a hard X-ray flare from the Galactic magnetar SGR 1935+2154, which supports that at least some FRBs are from magnetar activity. Interestingly, a portion of X-ray flares from magnetar XTE J1810-197 and the Sun are also accompanied by radio bursts. Many features of Galactic FRB 200428 and cosmological FRBs resemble solar radio bursts. However, a common physical origin among FRBs, magnetar radio pulses and solar radio bursts has not yet been established. Here we report a universal correlation between X-ray luminosity and radio luminosity over twenty orders of magnitude among solar type III radio bursts, XTE J1810-197 and Galactic FRB 200428. This universal correlation reveals that the energetic electrons which produce the X-ray flares and those which cause radio emissions have a common origin, which can give stringent limits on the generation process of radio bursts. Moreover, we find similar occurrence frequency distributions of energy, duration and waiting time for solar radio bursts, SGR 1935+2154 and repeating FRB 121102, which also support the tight correlation and the X-ray flares temporally associated with radio bursts. All of these distributions can be understood by avalanche models of self-organized criticality systems. The universal correlation and statistical similarities indicate that the Galactic FRB 200428 and FRBs seen at cosmological distances can be treated as scaled-up solar radio bursts.

preprint2020arXiv

On the Distance of SGR 1935+2154 Associated with FRB 200428 and Hosted in SNR G57.2+0.8

Owing to the detection of an extremely bright fast radio burst (FRB) 200428 associated with a hard X-ray counterpart from the magnetar soft gamma-ray repeater (SGR) 1935+2154, the distance of SGR 1935+2154 potentially hosted in the supernova remnant (SNR) G57.2+0.8 can be revisited. Under the assumption that the SGR and the SNR are physically related, in this Letter, by investigating the dispersion measure (DM) of the FRB contributed by the foreground medium of our Galaxy and the local environments and combining with other observational constraints, we find that the distance of SGR 1935+2154 turns out to be $9.0\pm2.5\,$kpc and the SNR radius falls into $10$ to $18\,$pc since the local DM contribution is as low as $0-18\,$pc cm$^{-3}$. These results are basically consistent with the previous studies. In addition, an estimate for the Faraday rotation measure of the SGR and SNR is also carried out.

preprint2020arXiv

Periodic Fast Radio Bursts as a Probe of Extragalactic Asteroid Belts

The periodic activity of repeating fast radio burst (FRB) 180916.J0158+65 was recently reported by the CHIME/FRB Collaboration team. 28 bursts from this source not only show a $\sim16$-day period with an active phase of $\sim 4.0$ days but also exhibit a broken power law in differential energy distribution. In this paper, we suggest that FRB 180916.J0158+65-like periodic FRBs would provide a unique probe of extragalactic asteroid belts (EABs), based on our previously-proposed pulsar-EAB impact model, in which repeating FRBs arise from an old-aged, slowly-spinning, moderately-magnetized pulsar traveling through an EAB around another stellar-mass object. These two objects form a binary and thus the observed period is in fact the orbital period. We show that this model can be used to well interpret all of the observed data of FRB 180916.J0158+65. Furthermore, we constrain the EAB's physical properties and find that (1) the outer radius of the EAB is at least an order of magnitude smaller than that of its analogue in the solar system, (2) the differential size distribution of the EAB's asteroids at small diameters (large diameters) is shallower (steeper) than that of solar-system small objects, and (3) the two belts have a comparable mass.

preprint2020arXiv

Reply to "Rapid $^{14}$C excursion at 3372-3371 BCE not observed at two different locations"

The nuclide $^{14}$C can be produced in the atmosphere by high energy particles and $γ$-rays from high-energy phenomena. Through the carbon cycle, some of $^{14}$CO$_2$ produced in the atmosphere can be retained in annual tree rings. Four events of rapid increase of the $^{14}$C content occurred in AD 775, AD 994, BC 660 and BC 3371 were found. Recently, the data of Jull et al. (2020) was inconsistent with our records around BC 3371. We measured our sample again and found the $^{14}$C records are consistent with the value in Wang et al. (2017). Therefore, our $^{14}$C records are robust. The inconsistency may be caused by the difference of calendar ages for the wood samples, or the physical origin of the event. First, crossdating on ring width can be performed only between trees whose growth has the same environmental conditions. Because the master tree-ring for dendrochronology is lack for Chinese trees. The master tree-ring from California has to be used. Therefore, the calendar ages derived from dendrochronology may be not precise. Second, the $^{14}$C even may be not global. One evidence is the variation of $^{14}$C content around AD 1006. The $^{14}$C contents of Californian trees increase 12\textperthousand~ in two years, while Japanese trees show no $^{14}$C increase.

preprint2020arXiv

The braking index of PSR B0540-69 and the associated pulsar wind nebula emission after spin-down rate transition

In Dec. 2011 PSR B0540-69 experienced a spin-down rate transition (SRT), after which the spin-down power of the pulsar increased by ~36%. About 1000 days after the SRT, the X-ray luminosity of the associated pulsar wind nebula (PWN) was found to brighten by 32+/-8%. After the SRT, the braking index n of PSR B0540-69 changes from n=2.12 to n=0.03 and then keeps this value for about five years before rising to n=0.9 in the following years. We find that most of the current models have difficulties in explaining the measured braking index. One exceptive model of the braking index evolution is the increasing dipole magnetic field of PSR B0540-69. We suggest that the field increase may result from some instabilities within the pulsar core that enhance the poloidal component at the price of toroidal component of the magnetic field. The increasing dipole magnetic field will result in the X-ray brightening of the PWN. We fit the PWN X-ray light curve by two models: one assumes a constant magnetic field within the PWN during the brightening and the other assumes an enhanced magnetic field proportional to the energy density of the PWN. It appears that the two models fit the data well, though the later model seems to fit the data a bit better. This provides marginal observational evidence that magnetic field in the PWN is generated by the termination shock. Future high-quality and high-cadence data are required to draw a solid conclusion.

preprint2020arXiv

The observation of the Crab Nebula with LHAASO-KM2A for the performance study

As a sub-array of the Large High Altitude Air Shower Observatory (LHAASO), KM2A is mainly designed to cover a large fraction of the northern sky to hunt for gamma-ray sources at energies above 10 TeV. Even though the detector construction is still underway, a half of the KM2A array has been operating stably since the end of 2019. In this paper, we present the pipeline of KM2A data analysis and the first observation on the Crab Nebula, a standard candle in very high energy gamma-ray astronomy. We detect gamma-ray signals from the Crab Nebula in both energy ranges of 10$-$100 TeV and $>$100 TeV with high significance, by analyzing the KM2A data of 136 live days between December 2019 and May 2020. With the observations, we test the detector performance including angular resolution, pointing accuracy and cosmic ray background rejection power. The energy spectrum of the Crab Nebula in the energy range 10-250 TeV fits well with a single power-law function dN/dE =(1.13$\pm$0.05$_{stat}$$\pm$0.08$_{sys}$)$\times$10$^{-14}$$\cdot$(E/20TeV)$^{-3.09\pm0.06_{stat}\pm0.02_{sys}}$ cm$^{-2}$ s$^{-1}$ TeV$^{-1}$. It is consistent with previous measurements by other experiments. This opens a new window of gamma-ray astronomy above 0.1 PeV through which ultrahigh-energy gamma-ray new phenomena, such as cosmic PeVatrons, might be discovered.

preprint2020arXiv

Understanding FRB 200428 in the synchrotron maser shock model: consistency and possible challenge

Recently, the discovery of Galactic FRB 200428 associated with a X-ray burst (XRB) of SGR 1935+2154 has built a bridge between FRBs and magnetar activities. In this paper, we assume that the XRB occurs in the magnetar magnetosphere. We show that the observational properties of FRB 200428 and the associated XRB are consistent with the predictions of synchrotron maser emission at ultrarelativistic magnetized shocks, including radiation efficiency, similar energy occurrence frequency distributions, and simultaneous arrive times. It requires that the upstream medium is a mildly relativistic baryonic shell ejected by a previous flare. The energy injection by flares responsible for the radio bursts will produce a magnetar wind nebula, which has been used to explain the persistent radio source associated FRB 121102. We find that the radio continuum around SGR 1935+2154 can be well understood in the magnetar wind nebula model, by assuming the same energy injection rate $\dot{E} \propto t^{-1.37}$ as FRB 121102. The required baryonic mass is also estimated form the observations of FRB 121102 by GBT and FAST. By assuming the same radiation efficiency $η\sim 10^{-5}$, the total baryonic mass ejected from the central magnetar is about 0.005 solar mass. This value is much larger than the typical mass of a magnetar outer crust, but is comparable to the total mass of a magnetar crust.

preprint2019arXiv

A bright electromagnetic counterpart to extreme mass ratio inspirals

The extreme mass ratio inspiral (EMRI), defined as a stellar-mass compact object inspiraling into a supermassive black hole (SMBH), has been widely argued to be a low-frequency gravitational wave (GW) source. EMRIs providing accurate measurements of black hole mass and spin, are one of the primary interests for Laser Interferometer Space Antenna (LISA). However, it is usually believed that there are no electromagnetic (EM) counterparts to EMRIs. Here we show a new formation channel of EMRIs with tidal disruption flares as EM counterparts. In this scenario, flares can be produced from the tidal stripping of the helium (He) envelope of a massive star by an SMBH. The left compact core of the massive star will evolve into an EMRI. We find that, under certain initial eccentricity and semimajor axis, the GW frequency of the inspiral can enter LISA band within 10 $\sim$ 20 years, which makes the tidal disruption flare an EM precursor to EMRI. Although the event rate is just $2\times 10^{-4}~\rm Gpc^{-3}yr^{-1}$, this association can not only improve the localization accuracy of LISA and help to find the host galaxy of EMRI, but also serve as a new GW standard siren for cosmology.

preprint2019arXiv

Gamma-ray Emission of 60Fe and 26Al Radioactivities in our Galaxy

The isotopes $^{60}$Fe and $^{26}$Al originate from massive stars and their supernovae, reflecting ongoing nucleosynthesis in the Galaxy. We studied the gamma-ray emission from these isotopes at characteristic energies 1173, 1332, and 1809 keV with over 15 years of SPI data, finding a line flux in $^{60}$Fe combined lines of $(0.31\pm 0.06) \times 10^{-3}$ ph cm$^{-2}$ s$^{-1}$ and the $^{26}$Al line flux of $(16.8\pm 0.7) \times 10^{-4}$ ph cm$^{-2}$ s$^{-1}$ above the background and continuum emission for the whole sky. Based on the exponential-disk grid maps, we characterise the emission extent of $^{26}$Al to find scale parameters $R_0 =7.0^{+1.5}_{-1.0}$ kpc and $z_0=0.8^{+0.3}_{-0.2}$ kpc, however the $^{60}$Fe lines are too weak to spatially constrain the emission. Based on a point source model test across the Galactic plane, the $^{60}$Fe emission would not be consistent with a single strong point source in the Galactic center or somewhere else, providing a hint for a diffuse nature. We carried out comparisons of emission morphology maps using different candidate-source tracers for both $^{26}$Al and $^{60}$Fe emissions, and suggests that the $^{60}$Fe emission is more likely to be concentrated towards the Galactic plane. We determine the $^{60}$Fe/$^{26}$Al $γ$-ray flux ratio at $(18.4\pm4.2)\,\%$ , when using a parameterized spatial morphology model. Across the range of plausible morphologies, it appears possible that $^{26}$Al and $^{60}$Fe are distributed differently in the Galaxy. Using the best fitting maps for each of the elements, we constrain flux ratios in the range 0.2--0.4. We discuss its implications for massive star models and their nucleosynthesis.

preprint2016arXiv

A Triple-Energy-Source Model for Superluminous Supernova iPTF13ehe

Almost all superluminous supernovae (SLSNe) whose peak magnitudes are $\lesssim -21$ mag can be explained by the $^{56}$Ni-powered model, magnetar-powered (highly magnetized pulsar) model or ejecta-circumstellar medium (CSM) interaction model. Recently, iPTF13ehe challenges these energy-source models, because the spectral analysis shows that $\sim 2.5M_\odot$ of $^{56}$Ni have been synthesized but are inadequate to power the peak bolometric emission of iPTF13ehe, while the rebrightening of the late-time light-curve (LC) and the H$α$ emission lines indicate that the ejecta-CSM interaction must play a key role in powering the late-time LC. Here we propose a triple-energy-source model, in which a magnetar together with some amount ($\lesssim 2.5M_\odot$) of $^{56}$Ni may power the early LC of iPTF13ehe while the late-time rebrightening can be quantitatively explained by an ejecta-CSM interaction. Furthermore, we suggest that iPTF13ehe is a genuine core-collapse supernova rather than a pulsational pair-instability supernova candidate. Further studies on similar SLSNe in the future would eventually shed light on their explosion and energy-source mechanisms.

preprint2016arXiv

Could FRB 131104 Originate from the Merger of Binary Neutron Stars?

Recently, DeLaunay et al. (2016) discovered a gamma-ray transient, Swift J0644.5-5111, associated with the fast radio burst (FRB) 131104. They also reported follow-up broadband observations beginning two days after the FRB and provided upper limits on a putative afterglow of this transient. In this paper, we show that if such a transient drives a relativistic shock as in a cosmological gamma-ray burst (GRB), these upper limits are consistent with an environment of which density is much less than that of an interstellar medium but typical for the outskirts' density of a galaxy when the typical values of three microphysical parameters of the shock are taken. This appears to be inconsistent with the catastrophic event models in which the central engine of Swift J0644.5-5111 is surrounded by an interstellar medium, but together with the properties of the gamma-ray transient, favors the binary neutron star merger origin. We further constrain the physical parameters of the postmerger object by assuming that Swift J0644.5-5111 results from internal dissipation of a spinning-down pulsar wind, and we find that the postmerger object is an ultra-strongly magnetized, very rapidly rotating pulsar. This merger event should have given birth to a gravitational wave burst, an FRB, and a short GRB or an extended X-ray/gamma-ray emission if a relativistic jet of the GRB is missed. Such "triplets" would be testable in the near future with the advanced LIGO and Virgo gravitational-wave observatories.

preprint2016arXiv

eXTP -- enhanced X-ray Timing and Polarimetry Mission

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

preprint2016arXiv

Optical Transients Powered by Magnetars: Dynamics, Light Curves, and Transition to the Nebular Phase

Millisecond magnetars can be formed via several channels: core-collapse of massive stars, accretion-induced collapse of white dwarfs (WDs), double WD mergers, double neutron star (NS) mergers, and WD-NS mergers. Because the mass of ejecta from these channels could be quite different, their light curves are also expected to be diverse. We evaluate the dynamic evolution of optical transients powered by millisecond magnetars. We find that the magnetar with short spin-down timescale converts its rotational energy mostly into the kinetic energy of the transient, while the energy of a magnetar with long spin-down timescale goes into radiation of the transient. This leads us to speculate that hypernovae could be powered by magnetars with short spin-down timescales. At late times the optical transients will gradually evolve into a nebular phase because of the photospheric recession. We treat the photosphere and nebula separately because their radiation mechanisms are different. In some cases the ejecta could be light enough that the magnetar can accelerate it to a relativistic speed. It is well known that the peak luminosity of a supernova (SN) occurs when the luminosity is equal to the instantaneous energy input rate, as shown by Arnett (1979). We show that photospheric recession and relativistic motion can modify this law. The photospheric recession always leads to a delay of the peak time $t_{\mathrm{pk}}$ relative to the time $t_{\times }$ at which the SN luminosity equals the instantaneous energy input rate. Relativistic motion, however, may change this result significantly.

preprint2016arXiv

Probing the Birth of Post-merger Millisecond Magnetars by X-ray and Gamma-ray Emission

There is growing evidence that a stable magnetar could be formed from the coalescence of double neutron stars. In previous papers, we investigated the signature of formation of stable millisecond magnetars in radio and optical/ultraviolet bands by assuming that the central rapidly rotating magnetar deposits its rotational energy in the form of a relativistic leptonized wind. We found that the optical transient PTF11agg could be the first evidence for the formation of post-merger millisecond magnetars. To enhance the probability of finding more evidence for the post-merger magnetar formation, it is better to extend the observational channel to other photon energy bands. In this paper we propose to search the signature of post-merger magnetar formation in X-ray and especially gamma-ray bands. We calculate the SSC emission of the reverse shock powered by post-merger millisecond magnetars. We find that the SSC component peaks at $1\,{\rm GeV}$ in the spectral energy distribution and extends to $\gtrsim 10\,{\rm TeV}$ for typical parameters. These energy bands are quite suitable for Fermi/LAT and CTA, which, with their current observational sensitivities, can detect the SSC emission powered by post-merger magnetars up to $1\,{\rm Gpc}$. NuSTAR, sensible in X-ray bands, can detect the formation of post-merger millisecond magnetars at redshift $z\sim 1$. Future improvement in sensitivity of CTA can also probe the birth of post-merger millisecond magnetars at redshift $z\sim 1$. However, because of the $γ$-$γ$ collisions, strong high-energy emission is clearly predicted only for ejecta masses lower than $10^{-3}M_\odot$.

preprint2016arXiv

Repeating Fast Radio Bursts from Highly Magnetized Pulsars Travelling through Asteroid Belts

Very recently Spitler et al. (2016) and Scholz et al. (2016) reported their detections of sixteen additional bright bursts from the direction of the fast radio burst (FRB) 121102. This repeating FRB is inconsistent with all the catastrophic event models put forward previously for hypothetically non-repeating FRBs. Here we propose a different model, in which highly magnetized pulsars travel through asteroid belts of other stars. We show that a repeating FRB could originate from such a pulsar encountering lots of asteroids in the belt. During each pulsar-asteroid impact, an electric field induced outside the asteroid has such a large component parallel to the stellar magnetic field that electrons are torn off the asteroidal surface and accelerated to ultra-relativistic energies instantaneously. Subsequent movement of these electrons along magnetic field lines will cause coherent curvature radiation, which can account for all the properties of an FRB. In addition, this model can self-consistently explain the typical duration, luminosity, and repetitive rate of the seventeen bursts of FRB 121102. The predicted occurrence rate of repeating FRB sources may imply that our model would be testable in the next few years.

preprint2016arXiv

Reverse Shock Emission Driven By Post-Merger Millisecond Magnetar Winds: Effects of the Magnetization Parameter

The study of short-duration gamma-ray bursts provides growing evidence that a good fraction of double neutron star mergers lead to the formation of stable millisecond magnetars. The launch of Poynting flux by the millisecond magnetars could leave distinct electromagnetic signatures that reveal the energy dissipation processes in the magnetar wind. In previous studies (Wang & Dai 2013b; Wang et al. 2015), we assume that the magnetar wind becomes completely lepton-dominated so that electrons/positrons in the magnetar wind are accelerated by a diffusive shock. However, theoretical modeling of pulsar wind nebulae shows that in many cases the magnetic field energy in the pulsar wind may be strong enough to suppress diffusive shock acceleration. In this paper, we investigate the reverse shock emission as well as the forward shock emission with an arbitrary magnetization parameter $σ$ of a magnetar wind. We find that the reverse shock emission strongly depends on $σ$, and in particular, $σ\sim 0.3$ leads to the strongest reverse shock emission. Future observations would be helpful to diagnose the composition of the magnetar wind.

preprint2015arXiv

A Unified Energy-Reservoir Model Containing Contributions from $^{56}$Ni and Neutron Stars and Its Implication to Luminous Type Ic Supernovae

Most type-Ic core-collapse supernovae (CCSNe) produce $^{56}$Ni and neutron stars (NSs) or black holes (BHs). The dipole radiation of nascent NSs has usually been neglected in explaining supernovae (SNe) with peak absolute magnitude $M_{\rm peak}$ in any band are $\gtrsim -19.5$~mag, while the $^{56}$Ni can be neglected in fitting most type-Ic superluminous supernovae (SLSNe Ic) whose $M_{\rm peak}$ in any band are $\lesssim -21$~mag, since the luminosity from a magnetar (highly magnetized NS) can outshine that from a moderate amount of $^{56}$Ni. For luminous SNe Ic with $-21 \lesssim M_{\rm peak}\lesssim -19.5$~mag, however, both contributions from $^{56}$Ni and NSs cannot be neglected without serious modeling, since they are not SLSNe and the $^{56}$Ni mass could be up to $\sim 0.5 M_{\odot}$. In this paper we propose a unified model that contain contributions from both $^{56}$Ni and a nascent NS. We select three luminous SNe Ic-BL, SN~2010ay, SN~2006nx, and SN~14475, and show that, if these SNe are powered by $^{56}$Ni, the ratio of $M_{\rm Ni}$ to $M_{\rm ej}$ are unrealistic. Alternatively, we invoke the magnetar model and the hybrid ($^{56}$Ni + NS) model and find that they can fit the observations, indicating that our models are valid and necessary for luminous SNe Ic. Owing to the lack of late-time photometric data, we cannot break the parameter degeneracy and thus distinguish among the model parameters, but we can expect that future multi-epoch observations of luminous SNe can provide stringent constraints on $^{56}$Ni yields and the parameters of putative magnetars.

preprint2015arXiv

An unexpectedly low-redshift excess of Swift gamma-ray burst rate

Gamma-ray bursts (GRBs) are the most violent explosions in the Universe and can be used to explore the properties of high-redshift universe. It is believed that the long GRBs are associated with the deaths of massive stars. So it is possible to use GRBs to investigate the star formation rate (SFR). In this paper, we use Lynden-Bell's $c^-$ method to study the luminosity function and rate of \emph{Swift} long GRBs without any assumptions. We find that the luminosity of GRBs evolves with redshift as $L(z)\propto g(z)=(1+z)^k$ with $k=2.43_{-0.38}^{+0.41}$. After correcting the redshift evolution through $L_0(z)=L(z)/g(z)$, the luminosity function can be expressed as $ψ(L_0)\propto L_0^{-0.14\pm0.02}$ for dim GRBs and $ψ(L_0)\propto L_0^{-0.70\pm0.03}$ for bright GRBs, with the break point $L_{0}^{b}=1.43\times10^{51}~{\rm erg~s^{-1}}$. We also find that the formation rate of GRBs is almost constant at $z<1.0$ for the first time, which is remarkably different from the SFR. At $z>1.0$, the formation rate of GRB is consistent with the SFR. Our results are dramatically different from previous studies. Some possible reasons for this low-redshift excess are discussed. We also test the robustness of our results with Monte Carlo simulations. The distributions of mock data (i.e., luminosity-redshift distribution, luminosity function, cumulative distribution and $\log N-\log S$ distribution) are in good agreement with the observations. Besides, we also find that there are remarkable difference between the mock data and the observations if long GRB are unbiased tracers of SFR at $z<1.0$.

preprint2015arXiv

Gamma-ray Burst Cosmology

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

preprint2015arXiv

TeV-PeV Neutrino Oscillation of Low-luminosity Gamma-ray Bursts

There is a sign that long-duration gamma-ray bursts (GRBs) originate from the core collapse of massive stars. During a jet puncturing through the progenitor envelope, high energy neutrinos can be produced by the reverse shock formed at the jet head. It is suggested that low-luminosity GRBs (LL-GRBs) are possible candidates of this high energy neutrino precursor up to $\sim {\rm PeV}$. Before leaving the progenitor, these high energy neutrinos must oscillate from one flavor to another with matter effect in the envelope. Under the assumption of a power-law stellar envelope density profile $ρ\propto r^{-α}$ with an index $α$, we study the properties of ${\rm TeV-PeV}$ neutrino oscillation. We find that adiabatic conversion is violated for these neutrinos so we do certain calibration of level crossing effect. The resonance condition is reached for different energies at different radii. We notice that the effective mixing angles in matter for ${\rm PeV}$ neutrinos are close to zero so the transition probabilities from one flavor to another are almost invariant for ${\rm PeV}$ neutrinos. We plot all the transition probabilities versus energy of ${\rm TeV-PeV}$ neutrinos from the birth place to the surface of the progenitor. With an initial flavor ratio $ϕ_{ν_e}^0:ϕ_{ν_μ}^0:ϕ_{ν_τ}^0=1:2:0$, we plot how the flavor ratio evolves with energy and distance when neutrinos are still in the envelope, and further get the ratio when they reach the Earth. For ${\rm PeV}$ neutrinos, the ratio is always $ϕ_{ν_e}:ϕ_{ν_μ}:ϕ_{ν_τ}\simeq0.30:0.37:0.33$ on Earth. In addition, we discuss the dependence of the flavor ratio on energy and $α$ and get a pretty good result. This dependence may provide a promising probe of the progenitor structure.

preprint2015arXiv

The Most Luminous Supernova ASASSN-15lh: Signature of a Newborn Rapidly-Rotating Strange Quark Star

In this paper we show that the most luminous supernova discovered very recently, ASASSN-15lh, could have been powered by a newborn ultra-strongly-magnetized pulsar, which initially rotates near the Kepler limit. We find that if this pulsar is a neutron star, its rotational energy could be quickly lost as a result of gravitational-radiation-driven r-mode instability; if it is a strange quark star, however, this instability is highly suppressed due to a large bulk viscosity associated with the nonleptonic weak interaction among quarks and thus most of its rotational energy could be extracted to drive ASASSN-15lh. Therefore, we conclude that such an ultra-energetic supernova provides a possible signature for the birth of a strange quark star.

preprint2015arXiv

The near-infrared radiation background, gravitational wave background and star formation rate of Pop III and Pop II during cosmic reionization

In this paper, we obtain the NIRB and SBGWs from the early stars, which are constrained by the observation of reionization and star formation rate. We study the transition from Pop III to Pop II stars via the star formation model of different population, which takes into account the reionization and the metal enrichment evolution. We calculate the two main metal pollution channels arising from the supernova-driven protogalactic outflows and "genetic channel". We obtain the SFRs of Pop III and Pop II and their NIRB and SBGWs radiation. We predict that the upper limit of metallicity in metal-enriched IGM (the galaxies whose polluted via "genetic channel") reaches $Z_{\rm crit}=10^{-3.5}Z_{\odot}$ at $z\sim13$ ($z\sim11$), which is consistent with our star formation model. We constrain on the SFR of Pop III stars from the observation of reionization. The peak intensity of NIRB is about $0.03-0.2~nW m^{-2}{sr}^{-1}$ at $\sim 1 μm$ for $z>6$. The prediction of NIRB signal is consistent with the metallicity evolution. We also obtain the gravitational wave background from the black holes formed by these early stars. The predicted gravitational wave background has a peak amplitude of $Ω_{GW}\simeq8\times10^{-9}$ at $ν=158$ Hz for Pop II star remnants. However, the background generated by Pop III.2 stars is much less than Pop II stars, with a peak amplitude of $Ω_{GW}\simeq1.2\times10^{-11}$ at $ν=28~Hz$. The background of Pop III.1 shifts to lower frequencies, and the amplitude of $Ω_{GW}$ for Pop III.1 stars shows a minimum value at $ν\simeq 10$ Hz, due to the lack of gravitational wave signals from the stars with $140~M_{\odot}<M_\ast<260~M_{\odot}$.

preprint2014arXiv

Long GRBs are metallicity-biased tracers of star formation: evidence from host galaxies and redshift distribution

We investigate the mass distribution of long gamma-ray burst (GRB) host galaxies and the redshift distribution of long GRBs by considering that long GRBs occur in low-metallicity environments. We calculate the upper limit on the stellar mass of a galaxy which can produce long GRBs by utilizing the mass-metallicity (M-Z) relation of galaxies. After comparing with the observed GRB host galaxies masses, we find that the observed GRB host galaxy masses can fit the predicted masses well if GRBs occur in low-metallicity $12+\log\rm(O/H)_{\rm KK04}<8.7$. GRB host galaxies have low metallicity, low mass, and high star formation rate compared with galaxies of seventh data release of the Sloan Digital Sky Survey. We also study the cumulative redshift distribution of the latest \emph{Swift} long GRBs by adding dark GRBs and 10 new GRBs redshifts from TOUGH survey. The observed discrepancy between the GRB rate and the star formation history can be reconciled by considering that GRBs tend to occur in low-metallicity galaxies with $12+\log\rm(O/H)_{\rm KK04}<8.7$. We conclude that the metallicity cutoff that can produce long GRBs is about $12+\log\rm(O/H)_{\rm KK04}<8.7$ from the host mass distribution and redshift distribution.

preprint2014arXiv

Neutrino Emission in Jet Propagation Process

Relativistic jets are universal in long-duration gamma-ray burst (GRB) models. Before breaking out, they must propagate in the progenitor envelope along with a forward shock and a reverse shock forming at the jet head. Both electrons and protons will be accelerated by the shocks. High energy neutrinos could be produced by these protons interacting with stellar materials and electron-radiating photons. The jet will probably be collimated, which may have a strong effect on the final neutrino flux. Under the assumption of a power-law stellar-envelope density profile $ρ\propto r^{-α}$ with an index $α$, we calculate this neutrino emission flux by these shocks for low-luminosity GRBs (LL-GRBs) and ultra-long GRBs (UL-GRBs) in different collimation regimes, using the jet propagation framework developed by \citet{bro11}. We find that LL-GRBs and UL-GRBs are capable for detectable high energy neutrinos up to $\sim {\rm PeV}$, and obtain the final neutrino spectrum. Besides, we conclude that larger $α$ corresponds to greater neutrino flux at high energy end ($\sim {\rm PeV}$) and higher maximum neutrino energy as well. However, such differences are so small that it is not promising for us to distinguish from observations, given the energy resolution we have now.

preprint2014arXiv

Similar radiation mechanism in gamma-ray bursts and blazars: evidence from two luminosity correlations

Active galactic nuclei (AGNs) and gamma-ray bursts (GRBs) are powerful astrophysical events with relativistic jets. In this Letter the broadband spectral properties are compared between GRBs and the well-observed blazars. The distribution of GRBs are consistent with the well-known blazar sequence including the $νL_ν(5\rm GHz)-α_{\rm RX}$ and $νL_ν(5\rm GHz)-ν_{\rm peak}$ correlations, where $α_{\rm RX}$ is defined as the broadband spectral slope in radio-to-X-ray bands, and $ν_{\rm peak}$ is defined as the spectral peak frequency. Moreover, GRBs occupy the low radio luminosity end of these sequences. These two correlations suggest that GRBs could have a similar radiation process with blazars both in the prompt emission and afterglow phases, i.e., synchrotron radiation.

preprint2014arXiv

Superluminous Supernovae Powered by Magnetars: Late-time Light Curves and Hard Emission Leakage

Recently, researches performed by two groups have revealed that the magnetar spin-down energy injection model with full energy trapping can explain the early-time light curves of SN 2010gx, SN 2013dg, LSQ12dlf, SSS120810 and CSS121015, but fails to fit the late-time light curves of these Superluminous Supernovae (SLSNe). These results imply that the original magnetar-powered model is challenged in explaining these SLSNe. Our paper aims to simultaneously explain both the early- and late-time data/upper limits by considering the leakage of hard emissions. We incorporate quantitatively the leakage effect into the original magnetar-powered model and derive a new semi-analytical equation. Comparing the light curves reproduced by our revised magnetar-powered model to the observed data and/or upper limits of these five SLSNe, we found that the late-time light curves reproduced by our semi-analytical equation are in good agreement with the late-time observed data and/or upper limits of SN 2010gx, CSS121015, SN 2013dg and LSQ12dlf and the late-time excess of SSS120810, indicating that the magnetar-powered model might be responsible for these SLSNe and that the gamma ray and X-ray leakage are unavoidable when the hard photons were down-Comptonized to softer photons. To determine the details of the leakage effect and unveil the nature of SLSNe, more high quality bolometric light curves and spectra of SLSNe are required.

preprint2014arXiv

The origin of the plateau and late rebrightening in the afterglow of GRB 120326A

GRB 120326A is an unusual gamma-ray burst (GRB) which has a quite long plateau and a very late rebrightening both in X-ray and optical bands. The similar behavior of the optical and X-ray light curves suggests that they maybe have a common origin. The long plateau starts from several hundred seconds and ends at tens of thousands seconds. The peak time of the late rebrightening is about 30000 s. We analyze the energy injection model by means of numerical and analytical solutions, considering both the wind environment and ISM environment for GRB afterglows. We especially study the influence of the injection starting time, ending time, stellar wind density (or density of the circumburst environment), and injection luminosity on the shape of the afterglow light curves, respectively. We find that the light curve is largely affected by the parameters in the wind model. There is a "bump" at the late time only in the wind model too. In the wind case, it is interesting that the longer the energy injected, the more obvious the rebrightening will be. We also find the peak time of bump is determined by the stellar wind density. We use the late continuous injection model to interpret the unusual afterglow of GRB 120326A. The model can well fit the observational data, however, we find that the time scale of the injection must be larger than ten thousands seconds. This implies that the time scale of the central engine activity must be more than ten thousands seconds. This can give useful constraints on the central engine of GRBs. We consider a new born millisecond pulsar with strong magnetic field as the central engine. On the other hand, our results suggest that the circumburst environment of GRB 120326A is very likely a stellar wind.

preprint2014arXiv

Universal Behavior of X-ray Flares from Black Hole Systems

X-ray flares have been discovered in black hole systems, such as gamma-ray bursts, the tidal disruption event Swift J1644+57, the supermassive black hole Sagittarius A$^*$ at the center of our Galaxy, and some active galactic nuclei. Their occurrences are always companied by relativistic jets. However, it is still unknown whether there is a physical analogy among such X-ray flares produced in black hole systems spanning nine orders of magnitude in mass. Here we report the observed data of X-ray flares, and show that they have three statistical properties similar to solar flares, including power-law distributions of energies, durations, and waiting times, which both can be explained by a fractal-diffusive self-organized criticality model. These statistical similarities, together with the fact that solar flares are triggered by a magnetic reconnection process, suggest that all of the X-ray flares are consistent with magnetic reconnection events, implying that their concomitant relativistic jets may be magnetically dominated.

preprint2014arXiv

X-Ray Afterglow Plateaus of Long Gamma-Ray Bursts: Further Evidence for Millisecond Magnetars

Many long-duration gamma-ray bursts (GRBs) were observed by {\it Swift}/XRT to have plateaus in their X-ray afterglow light curves. This plateau phase has been argued to be evidence for long-lasting activity of magnetar (ultra-strongly magnetized neutron stars) central engines. However, the emission efficiency of such magnetars in X-rays is still unknown. Here we collect 24 long GRB X-ray afterglows showing plateaus followed by steep decays. We extend the well-known relationship between the X-ray luminosity ${L_{\mathrm{X}}}$ and spin-down luminosity $L_{\mathrm{sd}}$ of pulsars to magnetar central engines, and find that the initial rotation period $P_{0}$ ranges from 1 ms to 10 ms and that the dipole magnetic field $B$ is centered around $10^{15}$ G. These constraints not only favor the suggestion that the central engines of some long GRBs are very likely to be rapidly rotating magnetars but also indicate that the magnetar plateau emission efficiency in X-rays is close to 100%.

preprint2013arXiv

Estimating the uncorrelated dark energy evolution in the Planck era

The equation of state (EOS), $w(z)$, is the most important parameter of dark energy. We reconstruct the evolution of this EOS in a model-independent way using the latest cosmic microwave background (CMB) data from Planck and other observations, such as type Ia supernovae (SNe Ia), the baryonic acoustic oscillation measurements (SDSS, 6dF, BOSS, and WiggleZ), and the Hubble parameter value $H(z)$. The results show that the EOS is consistent with the cosmological constant at the $2σ$ confidence level, not preferring a dynamical dark energy. The uncorrelated EOS of dark energy constraints from Planck CMB data are much tighter than those from the WMAP 9-year CMB data.

preprint2013arXiv

GeV Emission during X-Ray Flares from Late Internal Shocks: Application to GRB 100728A

Recently, the GeV radiation during the X-ray flare activity in GRB 100728A was detected by Feimi/LAT. Here we study the dynamics and emission properties of collision between two homogeneous shells based on the late internal shock model. The GeV photons can be produced from X-ray flare photons up-scattered by relativistic electrons accelerated by forward-reverse shocks, where involved radiative processes include synchrotron self-Compton and crossing inverse-Compton scattering. By analytical and numerical calculations, the observed spectral properties in GRB 100728A can be well explained.

preprint2013arXiv

PTF11agg as the First Evidence for Reverse Shock Emission from a Postmerger Millisecond Magnetar

Based on the stiff equations of state of neutron stars (NS) and the discovery of high-mass NSs, a NS-NS merger will leave behind, with high probabilities, a rapidly rotating massive magnetar. The central magnetar will dissipate its rotational energy to the outflow by injecting Poynting flux, which will become lepton-dominated so that a long-lasting reverse shock (RS) is developed. We calculate the emission of the RS as well as the emission of forward shock (FS) and find that in most cases the RS emission is stronger than FS emission. It is found that the recently discovered transient, PTF11agg, can be neatly accounted for by the RS emission powered by a millisecond magnetar. Other alternative models have been considered and cannot explain the observed light curves well. We therefore suggest that PTF11agg be the first evidence for RS emission from a postmerger millisecond magnetar.

preprint2013arXiv

Self-organized criticality in X-ray flares of gamma-ray burst afterglows

X-ray flares detected in nearly half of gamma-ray burst (GRB) afterglows are one of the most intriguing phenomena in high-energy astrophysics. All the observations indicate that the central engines of bursts, after the gamma-ray emission has ended, still have long periods of activity, during which energetic explosions eject relativistic materials, leading to late-time X-ray emission. It is thus expected that X-ray flares provide important clues to the nature of the central engines of GRBs, and more importantly, unveil the physical mechanism of the flares themselves, which has so far remained mysterious. Here we report statistical results of X-ray flares of GRBs with known redshifts, and show that X-ray flares and solar flares share three statistical properties: power-law frequency distributions for energies, durations, and waiting times. All of the distributions can be well understood within the physical framework of a self-organized criticality (SOC) system. The statistical properties of X-ray flares of GRBs are similar to solar flares, and thus can both be attributed to a SOC process. Both types of flares may be driven by a magnetic reconnection process, but X-ray flares of GRBs are produced in ultra-strongly magnetized millisecond pulsars or long-term hyperaccreting disks around stellar-mass black holes.

preprint2013arXiv

Testing the local-void alternative to dark energy using galaxy pairs

The possibility that we live in a special place in the universe, close to the center of a large, radially inhomogeneous void, has attracted attention recently as an alternative to dark energy or modified gravity to explain the accelerating universe. We show that the distribution of orientations of galaxy pairs can be used to test the Copernican principle that we are not in a central or special region of Universe. The popular void models can not fit both the latest type Ia supernova, cosmic microwave background data and the distribution of orientations of galaxy pairs simultaneously. Our results rule out the void models at the $4σ$ confidence level as the origin of cosmic acceleration and favor the Copernican principle.

preprint2012arXiv

Probing Pre-galactic Metal Enrichment with High-Redshift Gamma-Ray Bursts

We explore high-redshift gamma-ray bursts (GRBs) as promising tools to probe pre-galactic metal enrichment. We utilize the bright afterglow of a Pop III GRB exploding in a primordial dwarf galaxy as a luminous background source, and calculate the strength of metal absorption lines that are imprinted by the first heavy elements in the intergalactic medium (IGM). To derive the GRB absorption line diagnostics, we use an existing highly-resolved simulation of the formation of a first galaxy which is characterized by the onset of atomic hydrogen cooling in a halo with virial temperature >10^4 K. We explore the unusual circumburst environment inside the systems that hosted Pop III stars, modeling the density evolution with the self-similar solution for a champagne flow. For minihalos close to the cooling threshold, the circumburst density is roughly proportional to (1+z) with values of about a few cm^{-3}. In more massive halos, corresponding to the first galaxies, the density may be larger, n>100 cm^{-3}. The resulting afterglow fluxes may be detectable with the JWST and VLA in the near-IR and radio wavebands, respectively, out to redshift z>20. We predict that the maximum of the afterglow emission shifts from near-IR to millimeter bands with peak fluxes from mJy to Jy at different observed times. GRBs are ideal tools for probing the metal enrichment in the early IGM, due to their high luminosities and featureless power-law spectra. The metals in the first galaxies produced by the first supernova (SN) explosions are likely to reside in low-ionization stages. We show that if the afterglow can be observed sufficiently early, analysis of the metal lines can distinguish whether the first heavy elements were produced in a PISN, or a core-collapse (Type II) SN, thus constraining the initial mass function of the first stars.

preprint2012arXiv

Spin Evolution of Millisecond Magnetars with Hyperaccreting Fallback Disks: Implications for Early Afterglows of Gamma-Ray Bursts

The shallow decay phase or plateau phase of early afterglows of gamma-ray bursts (GRBs), discovered by Swift, is currently understood as being due to energy injection to a relativistic blast wave. One natural scenario for energy injection invokes a millisecond magnetar as the central engine of GRBs, because the conventional model of a pulsar predicts a nearly constant magnetic-dipole-radiation luminosity within the spin-down timescale. However, we note that significant brightening occurs in some early afterglows, which apparently conflicts with the above scenario. Here we propose a new model to explain this significant brightening phenomena by considering a hyperaccreting fallback disk around a newborn millisecond magnetar. We show that for typical values of the model parameters, sufficient angular momentum of the accreted matter is transferred to the magnetar and spins it up. It is this spin-up that leads to a dramatic increase of the magnetic dipole radiation luminosity with time and thus significant brightening of an early afterglow. Based on this model, we carry out numerical calculations and fit well early afterglows of 12 GRBs assuming sufficiently strong fallback accretion. If the accretion is very weak, our model turns out to be the conventional energy-injection scenario of a pulsar. Therefore, our model can provide a unified explanation for the shallow decay phase, plateaus, and significant brightening of early afterglows.

preprint2011arXiv

An Evolving Stellar Initial Mass Function and the Gamma-Ray Burst Redshift Distribution

Recent studies suggest that Swift gamma-ray bursts (GRBs) may not trace an ordinary star formation history. Here we show that the GRB rate turns out to be consistent with the star formation history with an evolving stellar initial mass function (IMF). We first show that the latest Swift sample of GRBs reveals an increasing evolution in the GRB rate relative to the ordinary star formation rate at high redshifts. We then assume only massive stars with masses greater than the critical value to produce GRBs, and use an evolving stellar IMF suggested by Davé (2010) to fit the latest GRB redshift distribution. This evolving IMF would increase the relative number of massive stars, which could lead to more GRB explosions at high redshifts. We find that the evolving IMF can well reproduce the observed redshift distribution of Swift GRBs.

preprint2011arXiv

Multiband Fitting to Three Long GRBs with Fermi/LAT Data: Structured Ejecta Sweeping up a Density-Jump Medium

We present broadband (radio, optical, X-ray and GeV) fits to the afterglow light curves and spectra of three long-duration gamma-ray bursts (GRBs 080916C, 090902B, and 090926A) detected by the Gamma-Ray Burst Monitor (GBM) and Large Area Telescope (LAT) instruments on the Fermi satellite. Using the observed broadband data, we study the origin of the high energy emission, and suggest that the early-time GeV emission and the late-time radio, optical, and X-ray afterglows can be understood as being due to synchrotron emission from an external forward shock caused by structured ejecta propagating in a wind bubble jumping to a homogeneous density medium. If the ceasing time for majority of the energy injection is assumed to be close to the deceleration time of the forward shock, the structured ejecta with continuous energy injection to the forward shock can well explain the early rising feature of the GeV mission from these burst, and the density-jump medium can account for some certain plateaus or flares in the late afterglows. From our fits, we find that, on one hand, the external shock origin of the GeV photons will make the optical depth have not significant contribution to the early LAT rising part, which will loosen strong constraint of lower limits of Lorentz factor. On the other hand, these Fermi-LAT events preferentially occur in a low-density circumburst environment, in which case the Klein-Nishina cutoff will significantly suppress the Self-Synchrotron Compton (SSC) radiation. Such an environment might result from superbubbles or low-metallicity progenitor stars (which have a low mass-loss rate at late times of stellar evolution) of type Ib/c supernovae.

preprint2011arXiv

The updated luminosity correlations of gamma-ray bursts and cosmological implications

Several interesting luminosity correlations among gamma-ray burst (GRB) variables have been recently discussed extensively. In this paper, we derive the six luminosity correlations ($\tlag-L$, $V-L$, $\epkk-L$, $\epkk-E_γ$, $\trt-L$, $\epkk-E_{γ, \mathrm{iso}}$) from the light curves and spectra of the latest 116 long GRBs, including the time lag ($\tlag$) between low and high photon energy light curves, the variability ($V$) of the light curve, the peak energy of the spectrum ($\epkk$), and the minimum rise time ($\trt$) of the peaks. We find that the intrinsic scatter of the $V-L$ correlation is too large and there seems no inherent correlation between the two parameters using the latest GRB data. The other five correlations indeed exist when the sample is enlarged. The $\epkk-E_γ$ correlation has a significantly lower intrinsic scatter compared to the other correlations. We divide the full data into four redshift bins when testing possible evolution of the correlations with redshift. We find no statistically significant evidence for the redshift evolution of the luminosity correlations. To avoid the circularity problem when constraining the cosmological parameters, we simultaneously minimize $χ^2$ with respect to both correlation parameters $a$, $b$ and the cosmological parameters using the maximum likelihood method. For the flat $Λ$CDM, the best fit is $\omm=0.31^{+0.13}_{-0.10}$. We also constrain the possible evolution of the equation of state (EOS) of the dark energy using the GRBs together with the Union2 compilation of SNe Ia and the $H(z)$ data. The result is consistent with the cosmological constant at $2 σ$ confidence level and mainly due to the GRB data, the dark energy EOS shows slight deviation from -1 at $z \geq0.5$ as was persistently presented with many previous data sets.

preprint2011arXiv

Weak gravitational lensing effects on cosmological parameters and dark energy from gamma-ray bursts

Gamma-ray bursts (GRBs) are attractive objects for constraining the nature of dark energy in a way complementary to other cosmological probes, especially at high redshifts. However, the apparent magnitude of distant GRBs can be distorted by the gravitational lensing from the density fluctuations along the line of sight. We investigate the gravitational lensing effect on the cosmological parameters and dark energy equation of state from GRBs. We first calibrated the GRB luminosity relations without assuming any cosmological models. The luminosity distances of low-redshift GRBs were calibrated with the cosmography method using a latest type Ia supernova (SNe Ia) sample. We find that the gravitational lensing has non-negligible effects on the determination of cosmological parameters and dark energy. The gravitational lensing shifts the best-fit constraints on cosmological parameters and dark energy. Because high-redshift GRBs are more likely to be reduced, the most probable value of the observed matter density $Ω_M$ is slightly lower than its actual value. In the $Λ$CDM model, we find that the matter density parameter $Ω_M$ will shift from 0.30 to 0.26 after including the gravitational lensing effect. The gravitational lensing also affects the dark energy equation of state by shifting it to a more negative value. We constrain the dark energy equation of state out to redshift $z\sim 8$ using GRBs for the first time, and find that the equation of state deviates from $Λ$CDM at the $1σ$ confidence level, but agrees with $w=-1$ at the $2σ$ confidence level.

preprint2010arXiv

High-redshift star formation rate up to z~8.3 derived from gamma-ray bursts and influence of background cosmology

The high-redshift star formation rate (SFR) is difficult to measure directly even by modern approaches. Long-duration gamma-ray bursts (GRBs) can be detected to the edge of the visible universe because of their high luminorsities. The collapsar model of long gamma-ray bursts indicates that they may trace the star formation history. So long gamma-ray bursts may be a useful tool of measuring the high-redshift SFR. Observations show that long gamma-ray bursts prefer to form in a low-metallicity environment. We study the high-redshift SFR up to z~8.3 considering the Swift GRBs tracing the star formation history and the cosmic metallicity evolution in different background cosmological models including $Λ$CDM, quintessence, quintessence with a time-varying equation of state, and brane-world model. We use latest Swift GRBs including two highest-$z$ GRBs, GRB 080913 at $z=6.7$ and GRB 090423 at $z=8.3$. We find that the SFR at $z>4$ shows a steep decay with a slope of $\sim -5.0$ in $Λ$CDM. In the other three models, the high-redshift SFR is slightly different from $Λ$CDM model, and also shows a steep decay.

preprint2010arXiv

Hyperaccreting Disks around Magnetars for Gamma-Ray Bursts: Effects of Strong Magnetic Fields

(Abridged) The hyperaccreting neutron star or magnetar disks cooled via neutrino emission can be a candidate of gamma-ray burst (GRB) central engines. The strong field $\geq10^{15}-10^{16}$ G of the magnetar can play a significant role in affecting the disk properties and even lead to the funnel accretion process. We investigate the effects of strong fields on the disks around magnetars, and discuss implications of such accreting magnetar systems for GRB and GRB-like events. We discuss quantum effects of the strong fields on the disk, and use the MHD conservation equations to describe the behavior of the disk flow coupled with a large scale field, which is generated by the star-disk interaction. In general, stronger fields give higher disk densities, pressures, temperatures and neutrino luminosity, and change the electron fraction and degeneracy state significantly. A magnetized disk is always viscously stable outside the Alfvén radius, but will be thermally unstable near the Alfvén radius where the magnetic field plays a more important role in transferring the angular momentum and heating the disk than the viscous stress. The funnel accretion process will be only important for an extremely strong field, which creates a magnetosphere inside the Alfvén radius and truncates the plane disk. Because of higher temperature and more concentrated neutrino emission of the magnetar surface ring-like belt region covered by funnel accretion, the neutrino annihilation rate from the accreting magnetars can be much higher than that from accreting neutron stars without fields. Furthermore, the neutrino annihilation mechanism and the magnetically-driven pulsar wind from the magnetar surface can work together to generate and feed an ultra-relativistic jet along the stellar magnetic poles.

preprint2010arXiv

Hyperaccreting Disks around Neutrons Stars and Magnetars for GRBs: Neutrino Annihilation and Strong Magnetic Fields

Hyperaccreting disks around neutron stars or magnetars cooled via neutrino emission can be the potential central engine of GRBs. The neutron-star disk can cool more efficiently, produce much higher neutrino luminosity and neutrino annihilation luminosity than its black hole counterpart with the same accretion rate. The neutron star surface boundary layer could increase the annihilation luminosity as well. An ultra relativistic jet via neutrino annihilation can be produced along the stellar poles. Moreover, we investigate the effects of strong fields on the disks around magnetars. In general, stronger fields give higher disk densities, pressures, temperatures and neutrino luminosity; the neutrino annihilation mechanism and the magnetically-driven pulsar wind which extracts the stellar rotational energy can work together to generate and feed an even stronger ultra-relativistic jet along the stellar magnetic poles.

preprint2010arXiv

Probing the cosmographic parameters to distinguish between dark energy and modified gravity models

In this paper we investigate the deceleration, jerk and snap parameters to distinguish between the dark energy and modified gravity models by using high redshift gamma-ray bursts (GRBs) and supernovae (SNe). We first derive the expressions of deceleration, jerk and snap parameters in dark energy and modified gravity models. In order to constrain the cosmographic parameters, we calibrate the GRB luminosity relations without assuming any cosmological models using SNe Ia. Then we constrain the models (including dark energy and modified gravity models) parameters using type Ia supernovae and gamma-ray bursts. Finally we calculate the cosmographic parameters. GRBs can extend the redshift - distance relation up to high redshifts, because they can be detected to high redshifts.We find that the statefinder pair (r,s) could not be used to distinguish between some dark energy and modified gravity models, but these models could be differentiated by the snap parameter. Using the model-independent constraints on cosmographic parameters, we conclude that the \LambdaCDM model is consistent with the current data.

preprint2010arXiv

The termination shock of a magnetar wind: a possible origin of gamma-ray burst X-ray afterglow emission

Context: Swift observations suggest that the X-ray afterglow emission of some gamma-ray bursts (GRB) may have internal origins, and the conventional external shock (ES) cannot be the exclusive source of the afterglow emission. Aims: If the central compact objects of some GRBs are millisecond magentars, the magnetar winds could play an important role in the (internal) X-ray afterglow emission, which is our focus here. Methods: The dynamics and the synchrotron radiation of the termination shock (TS) of the magmnetar winds, as well as the simultaneous GRB ES, are investigated by considering the magnetization of the winds. Results: As a result of the competition between the emission of the wind TS and the GRB ES, two basic types of X-ray afterglows are predicted, i.e., the TS-dominated and the ES-dominated types. Moreover, our results also show that both of the two types of afterglows have a shallow-decay phase and a normal-decay one, as observed by the \textit{Swift} satellite. This indicates that some observed X-ray afterglows could be (internally) produced by the magnetar winds, but not necessarily GRB ESs.

preprint2009arXiv

GRB Jet Beaming Angle Statistics

Existing theory and models suggest that a Type I (merger) GRB should have a larger jet beaming angle than a Type II (collapsar) GRB, but so far no statistical evidence is available to support this suggestion. In this paper, we obtain a sample of 37 beaming angles and calculate the probability that this is true. A correction is also devised to account for the scarcity of Type I GRBs in our sample. The probability is calculated to be 83% without the correction and 71% with it.

preprint2008arXiv

Optical and gamma-ray emissions from internal forward-reverse shocks: application to GRB 080319B?

In the popular internal shock model for the prompt emission of gamma-ray bursts (GRBs), collisions between a series of relativistic shells generate lots of paired forward and reverse shocks. We show that the synchrotron emission produced by the forward and reverse shocks respectively could peak at two quite different energy bands if the Lorentz factors of these two types of shocks are significantly different with each other (e.g., one shock is relativistic and the other is Newtonian). We then investigate whether this scenario is applicable to the case of GRB 080319B and find that a bimodal distribution of the shell Lorentz factors, peaking at $\sim400$ and $\sim10^5$, is required. In addition, this scenario predicts an accompanying inverse-Compton (IC) GeV emission with a luminosity comparable to (not much higher than) that of the synchrotron MeV emission, which can be tested with future \textit{Fermi} observations.

preprint2008arXiv

Self-similar structure of magnetized ADAFs and CDAFs

(Abridged) We study the effects of a global magnetic field on viscously-rotating and vertically-integrated accretion disks around compact objects using a self-similar treatment. We extend Akizuki & Fukue's work (2006) by discussing a general magnetic field with three components ($r, ϕ, z$) in advection-dominated accretion flows (ADAFs). We also investigate the effects of a global magnetic field on flows with convection. For these purposes, we first adopt a simple form of the kinematic viscosity $ν=αc_{s}^{2}/Ω_{K}$ to study magnetized ADAFs. Then we consider a more realistic model of the kinematic viscosity $ν=αc_{s}H$, which makes the infall velocity increase but the sound speed and toroidal velocity decrease. We next use two methods to study magnetized flows with convection, i.e., we take the convective coefficient $α_{c}$ as a free parameter to discuss the effects of convection for simplicity. We establish the $α_{c}-α$ relation for magnetized flows using the mixing-length theory and compare this relation with the non-magnetized case. If $α_{c}$ is set as a free parameter, then $|v_{r}|$ and $c_{s}$ increase for a large toroidal magnetic field, while $|v_{r}|$ decreases but $|v_ϕ|$ increases (or decreases) for a strong and dominated radial (or vertical) magnetic field with increasing $α_{c}$. In addition, the magnetic field makes the $α_{c}-α$ relation be distinct from that of non-magnetized flows, and allows the $ρ\propto r^{-1}$ or $ρ\propto r^{-2}$ structure for magnetized non-accreting convection-dominated accretion flows with $α+gα_{c}< 0$ (where $g$ is the parameter to determine the condition of convective angular momentum transport).

preprint2008arXiv

X-Ray and high energy flares from late internal shocks of gamma-ray bursts

We study afterglow flares of gamma-ray bursts (GRBs) in the framework of the late internal shock (LIS) model based on a careful description for the dynamics of a pair of shocks generated by a collision between two homogeneous shells,. First, by confronting the model with some fundamental observational features of X-ray flares, we find some constraints on the properties of the pre-collision shells that are directly determined by the central engine of GRBs. Second, high energy emission associated with X-ray flares, which arises from synchrotron self-Compton (SSC) emission of LISs, is investigated in a wide parameter space. The predicted flux of high energy flares may reach as high as $\sim 10^{-8}\rm erg cm^{-2}s^{-1}$, which is likely to be detectable with the Large Area Telescope (LAT) aboard \textit{the Fermi Space Telescope}

preprint2007arXiv

Behavior of X-Ray Dust Scattering and Implications for X-Ray Afterglows of Gamma-Ray Bursts

The afterglows of gamma-ray bursts (GRBs) have commonly been assumed to be due to shocks sweeping up the circum-stellar medium. However, most GRBs have been found in dense star-forming regions where a significant fraction of the prompt X-ray emission can be scattered by dust grains. Here we revisit the behavior of dust scattering of X-rays in GRBs. We find that the features of some X-ray afterglows from minutes to days after the gamma-ray triggers are consistent with the scattering of prompt X-ray emission from GRBs off host dust grains. This implies that some of the observed X-ray afterglows (especially those without sharp rising and decaying flares) could be understood with a dust-scattering--driven emission model.

preprint2007arXiv

Echo Emission From Dust Scattering and X-Ray Afterglows of Gamma-Ray Bursts

We investigate the effect of X-ray echo emission in gamma-ray bursts (GRBs). We find that the echo emission can provide an alternative way of understanding X-ray shallow decays and jet breaks. In particular, a shallow decay followed by a "normal" decay and a further rapid decay of X-ray afterglows can be together explained as being due to the echo from prompt X-ray emission scattered by dust grains in a massive wind bubble around a GRB progenitor. We also introduce an extra temporal break in the X-ray echo emission. By fitting the afterglow light curves, we can measure the locations of the massive wind bubbles, which will bring us closer to finding the mass loss rate, wind velocity, and the age of the progenitors prior to the GRB explosions.

preprint2007arXiv

Low Luminosity Gamma-Ray Bursts as a Unique Population: Luminosity Function, Local Rate, and Beaming Factor

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

preprint2006arXiv

Tail emission from a ring-like jet: its application to shallow decays of early afterglows and to GRB 050709

Similar to the pulsar, the magnetic axis and the spin axis of the gamma-ray burst source may not lie on the same line. This may cause a ring-like jet due to collimation of the precessing magnetic axis. We analyze the tail emission from such a jet, and find that it has a shallow decay phase with temporal index equal to -1/2 if the Lorentz factor of the ejecta is not very high. This phase is consistent with the shallow decay phase of some early X-ray afterglow detected by {\it{swift}}. The ring-like jet has a tail cusp with sharp rising and very sharp decay. This effect can provide an explanation for the re-brightening and sharp decay of the X-ray afterglow of GRB 050709.

preprint2006arXiv

The shock break-out of GRB 060218/SN 2006aj

Although the link between long Gamma Ray Bursts (GRBs) and supernovae (SNe) has been established, hitherto there have been no observations of the beginning of a supernova explosion and its intimate link to a GRB. In particular, we do not know however how a GRB jet emerges from the star surface nor how a GRB progenitor explodes. Here we report on observations of the close GRB060218 and its connection to SN2006aj. In addition to the classical non-thermal emission, GRB060218 shows a thermal component in its X-ray spectrum, which cools and shifts into the optical/UV band as time passes. We interpret these features as arising from the break out of a shock driven by a mildly relativistic shell into the dense wind surrounding the progenitor. Our observations allow us for the first time to catch a SN in the act of exploding, to directly observe the shock break-out and to provide strong evidence that the GRB progenitor was a Wolf-Rayet star.

preprint2005arXiv

Fireball/Blastwave Model and Soft Gamma-ray Repeaters

Soft gamma-ray repeaters are at determined distances and their positions are known accurately. If observed, afterglows from their soft gamma-ray bursts will provide important clues to the study of the so called "classical gamma-ray bursts". On applying the popular fireball/blastwave model of classical gamma-ray bursts to soft gamma-ray repeaters, it is found that their X-ray and optical afterglows are detectable. Monitoring of the three repeaters is solicited.

preprint2004arXiv

Afterglow Light Curves of Jetted Gamma-ray Burst Ejecta in Stellar Winds

Optical and radio afterglows arising from the shocks by relativistic conical ejecta running into pre-burst massive stellar winds are revisited. Under the homogeneous thin-shell approximation and a realistic treatment for the lateral expansion of jets, our results show that a notable break exists in the optical light curve in most cases we calculated in which the physical parameters are varied within reasonable ranges. For a relatively tenuous wind which cannot decelerate the relativistic jet to cause a light curve break within days, the wind termination shock due to the ram pressure of the surrounding medium occurs at a small radius, namely, a few times 10^{17} cm. In such a structured wind environment, the jet will pass through the wind within several hours and run into the outer uniform dense medium. The resulting optical light curve flattens with a shallower drop after the jet encounters the uniform medium, and then declines deeply, triggered by runaway lateral expansion.

preprint2004arXiv

Constraining $Ω_M$ and Dark Energy with Gamma-Ray Bursts

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

preprint2003arXiv

Optical Flash of GRB 990123: constraints on the physical parameter of the reverse shock

The optical flash accompanying GRB 990123 is believed to be powered by the reverse shock of a thin shell. With the best fitted physical parameters for GRB 990123 (Panaitescu & Kumar 2001) and the assumption that the parameters in the optical flash are the same as those in the afterglow, we show that: 1) the shell is thick but not thin, and we have provided the light curve for the thick shell case which coincides with the observation; 2) the theoretical peak flux of the optical flash accounts for only 3\times10^{-4} of the observed. In order to compensate this divergency, the physical parameters the electron energy ration and the magnetic ratio $ε_e$, $ε_B$ should be 0.61, 0.39 respectively, which are much different from those in the late afterglow.

preprint2002arXiv

Beaming Effects in Gamma-Ray Bursts

Based on a refined generic dynamical model, we investigate afterglows from jetted gamma-ray burst (GRB) remnants numerically. In the relativistic phase, the light curve break could marginally be seen. However, an obvious break does exist at the transition from the relativistic phase to the non-relativistic phase, which typically occurs at time 10 to 30 days. It is very interesting that the break is affected by many parameters, especially by the electron energy fraction (xi_e), and the magnetic energy fraction (xi_B^2). Implication of orphan afterglow surveys on GRB beaming is investigated. The possible existence of a kind of cylindrical jets is also discussed.

preprint2002arXiv

Dynamical Evolution of Gamma-Ray Burst Remnants with Evolving Radiative Efficiency

In previous works, a generic dynamical model has been suggested by Huang et al., which is shown to be correct for both adiabatic and radiative blastwaves, and in both ultra-relativistic and non-relativistic phases. In deriving their equations, Huang et al. have assumed that the radiative efficiency of the fireball is constant. They then applied their model directly to realistic cases where the radiative efficiency evolves with time. In this paper, we abandon the above assumption and re-derive a more accurate dynamical equation for gamma-ray burst remnants. Numerical results show that Huang et al.'s model is accurate enough in general cases.

preprint2001arXiv

Light Curves of Jetted Gamma-Ray Burst Afterglows in Circumstellar Clouds

The afterglow emission from a spreading jet expanding in a circumstellar cloud is discussed. Prompt X-ray radiation and a strong UV flash from the reverse shock produced by the interaction of the jet with the cloud may destroy and clear the dust out to about 30 pc within the initial solid angle of the jet. As the sideways expansion of the jet becomes significant, most of the optical radiation from the high-latitude part of the jet may be absorbed by the dust outside the initial solid angle of the jet, but only the radiation from the part within the initial solid angle can be observed. We analytically show that the flux of the observational radiation decays as $\propto t^{-(p+1)}$ (where $p$ is the power-law index of the electron distribution) in the relativistic phase. This preliminary result motivates us to perform numerical calculations. Our results show that one break in the optical afterglow ligh curve extends over a factor of $\sim 3$ in time rather than one decade in time in the previous jet model. These results may provide a way to judge whether GRBs locate in dense clouds or not. Finally, we carry out a detailed modelling for the R-band afterglow of GRB 000926.

preprint2000arXiv

On the Optical Light Curves of Afterglows from Jetted Gamma-ray Burst Ejecta: Effects of Parameters

Due to some refinements in the dynamics, we can follow the overall evolution of a realistic jet numerically till its bulk velocity being as small as $βc \sim 10^{-3} c$. We find no obvious break in the optical light curve during the relativistic phase itself. However, an obvious break does exist at the transition from the relativistic phase to the non-relativistic phase, which typically occurs at time $t \sim 10^6$ --- $10^{6.5}$ s (i.e., 10 --- 30 d). The break is affected by many parameters, such as the electron energy fraction $ξ_{\rm e}$, the magnetic energy fraction $ξ_{\rm B}^2$, the initial half opening angle $θ_0$, and the medium number density $n$. Increase of any of them to a large enough value will make the break disappear. Although the break itself is parameter-dependent, afterglows from jetted GRB remnants are uniformly characterized by a quick decay during the non-relativistic phase, with power law timing index $α\geq 2.1$. This is quite different from that of isotropic fireballs, and may be of fundamental importance for determining the degree of beaming in $γ$-ray bursts observationally.

preprint1999arXiv

Overall Evolution of Realistic Gamma-ray Burst Remnant and Its Afterglow

Conventional dynamic model of gamma-ray burst remnants is found to be incorrect for adiabatic blastwaves during the non-relativistic phase. A new model is derived, which is shown to be correct for both radiative and adiabatic blastwaves during both ultra-relativistic and non-relativistic phase. Our model also takes the evolution of the radiative efficiency into account. The importance of the transition from the ultra-relativistic phase to the non-relativistic phase is stressed.

preprint1998arXiv

GRB afterglows: from ultra-relativistic to non-relativistic phase

Postburst evolution of adiabatic fireballs that produce $γ$-ray bursts is studied. Emphasis has been put on the transition from the highly relativistic phase to the non-relativistic phase, which, according to our calculation, should happen much earlier than previously expected. The theoretical light curves decline a little sharper at the non-relativistic stage than at early times, but still can fit the observations well. However, disagreements are obvious when $t \ge 100$ d, implying a large initial energy $E_0$, a low interstellar density $n$, or the possible existence of a persistent energy source at the center of the fireball.

Institution

Affiliation not imported yet

This author record came from a source that does not expose affiliation metadata. Once the author claims the profile or we enrich the record from another provider, this section will link to the concrete institution.

Topic footprint

Fields this researcher appears in

astro-ph.HE astro-ph astro-ph.CO astro-ph.GA astro-ph.SR hep-ph astro-ph.EP nlin.AO physics.geo-ph

Source provenance

Where this author record came from

arxivconfidence 95%

external id: arxiv:2408.09905:author:43:z-g-dai

Imported May 21, 2026Synced May 21, 2026

arxivconfidence 95%

external id: arxiv:2504.00601:author:43:z-g-dai

Imported May 21, 2026Synced May 21, 2026

23 works

F. Y. Wang

Researcher

F. Y. Wang contributes to research discovery and scholarly infrastructure.

Open to collaborate

19 works

D. Bastieri

Researcher

D. Bastieri contributes to research discovery and scholarly infrastructure.

Open to collaborate

19 works

F. Aharonian

Researcher

F. Aharonian contributes to research discovery and scholarly infrastructure.

Open to collaborate

19 works

J. Chang

Researcher

J. Chang contributes to research discovery and scholarly infrastructure.

Open to collaborate

Z. G. Dai

What is connected

Connect this record

See the researcher in context

Building this map preview

87 published item(s)

An Enigmatic PeVatron in an Area around HII Region G35.6$-$0.5

An Ultrahigh-energy $γ$-ray Bubble Powered by a Super PeVatron

Constraining the Cosmic-ray Energy Based on Observations of Nearby Galaxy Clusters by LHAASO

Constraints on heavy decaying dark matter from 570 days of LHAASO observations

Deep view of Composite SNR CTA1 with LHAASO in $γ$-rays up to 300 TeV

Discovery of a new $γ$-ray source LHAASO J0341+5258 with emission up to 200TeV

Discovery of the Ultra-high energy gamma-ray source LHAASO J2108+5157

Energy calibration of LHAASO-KM2A using the cosmic ray Moon shadow

Energy-Dependent Shifts of Medium-Scale Anisotropies in Very-High-Energy Cosmic Rays Observed by LHAASO-KM2A

Evidence for particle acceleration approaching PeV energies in the W51 complex

Exploring Lorentz Invariance Violation from Ultra-high-energy Gamma Rays Observed by LHAASO

Extended Very-High-Energy Gamma-Ray Emission Surrounding PSR J0622 + 3749 Observed by LHAASO-KM2A

LHAASO Detection of Ultra-High-Energy Gamma-Ray Emission toward the Giant Molecular Clouds

Measurement of attenuation length of the muon content in extensive air showers from 0.3 to 30 PeV with LHAASO

Measurement of Very-high-energy Diffuse Gamma-ray Emissions from the Galactic Plane with LHAASO-WCDA

Transient Large-Scale Anisotropy in TeV Cosmic Rays due to an Interplanetary Coronal Mass Ejection

Observation of the $γ$-ray Emission from W43 with LHAASO

Study of Ultra-High-Energy Gamma-Ray Source 1LHAASO J0056+6346u and Its Possible Origins

Flaring-associated Complex Dynamics in Two M-dwarfs Revealed by Fast, Time-resolved Spectroscopy

Quasi-periodic Eruptions from Helium Envelope of Hydrogen-deficient Stars Stripped by Supermassive Black Holes

Standardized long gamma-ray bursts as a cosmic distance indicator

A Magnetar-Asteroid Impact Model for FRB 200428 Associated with an X-ray Burst from SGR 1935+2154

Fast radio bursts from activities of neutron stars newborn in BNS mergers: offset, birth rate and observational properties

Galactic and cosmological fast radio bursts as scaled-up solar radio bursts

On the Distance of SGR 1935+2154 Associated with FRB 200428 and Hosted in SNR G57.2+0.8

Periodic Fast Radio Bursts as a Probe of Extragalactic Asteroid Belts

Reply to "Rapid $^{14}$C excursion at 3372-3371 BCE not observed at two different locations"

The braking index of PSR B0540-69 and the associated pulsar wind nebula emission after spin-down rate transition

The observation of the Crab Nebula with LHAASO-KM2A for the performance study

Understanding FRB 200428 in the synchrotron maser shock model: consistency and possible challenge

A bright electromagnetic counterpart to extreme mass ratio inspirals

Gamma-ray Emission of 60Fe and 26Al Radioactivities in our Galaxy

A Triple-Energy-Source Model for Superluminous Supernova iPTF13ehe

Could FRB 131104 Originate from the Merger of Binary Neutron Stars?

eXTP -- enhanced X-ray Timing and Polarimetry Mission

Optical Transients Powered by Magnetars: Dynamics, Light Curves, and Transition to the Nebular Phase

Probing the Birth of Post-merger Millisecond Magnetars by X-ray and Gamma-ray Emission

Repeating Fast Radio Bursts from Highly Magnetized Pulsars Travelling through Asteroid Belts

Reverse Shock Emission Driven By Post-Merger Millisecond Magnetar Winds: Effects of the Magnetization Parameter

A Unified Energy-Reservoir Model Containing Contributions from $^{56}$Ni and Neutron Stars and Its Implication to Luminous Type Ic Supernovae

An unexpectedly low-redshift excess of Swift gamma-ray burst rate

Gamma-ray Burst Cosmology

TeV-PeV Neutrino Oscillation of Low-luminosity Gamma-ray Bursts

The Most Luminous Supernova ASASSN-15lh: Signature of a Newborn Rapidly-Rotating Strange Quark Star

The near-infrared radiation background, gravitational wave background and star formation rate of Pop III and Pop II during cosmic reionization

Long GRBs are metallicity-biased tracers of star formation: evidence from host galaxies and redshift distribution

Neutrino Emission in Jet Propagation Process

Similar radiation mechanism in gamma-ray bursts and blazars: evidence from two luminosity correlations

Superluminous Supernovae Powered by Magnetars: Late-time Light Curves and Hard Emission Leakage

The origin of the plateau and late rebrightening in the afterglow of GRB 120326A

Universal Behavior of X-ray Flares from Black Hole Systems

X-Ray Afterglow Plateaus of Long Gamma-Ray Bursts: Further Evidence for Millisecond Magnetars

Estimating the uncorrelated dark energy evolution in the Planck era

GeV Emission during X-Ray Flares from Late Internal Shocks: Application to GRB 100728A

PTF11agg as the First Evidence for Reverse Shock Emission from a Postmerger Millisecond Magnetar

Self-organized criticality in X-ray flares of gamma-ray burst afterglows

Testing the local-void alternative to dark energy using galaxy pairs

Probing Pre-galactic Metal Enrichment with High-Redshift Gamma-Ray Bursts

Spin Evolution of Millisecond Magnetars with Hyperaccreting Fallback Disks: Implications for Early Afterglows of Gamma-Ray Bursts

An Evolving Stellar Initial Mass Function and the Gamma-Ray Burst Redshift Distribution

Multiband Fitting to Three Long GRBs with Fermi/LAT Data: Structured Ejecta Sweeping up a Density-Jump Medium

The updated luminosity correlations of gamma-ray bursts and cosmological implications

Weak gravitational lensing effects on cosmological parameters and dark energy from gamma-ray bursts

High-redshift star formation rate up to z~8.3 derived from gamma-ray bursts and influence of background cosmology

Hyperaccreting Disks around Magnetars for Gamma-Ray Bursts: Effects of Strong Magnetic Fields

Hyperaccreting Disks around Neutrons Stars and Magnetars for GRBs: Neutrino Annihilation and Strong Magnetic Fields

Probing the cosmographic parameters to distinguish between dark energy and modified gravity models

The termination shock of a magnetar wind: a possible origin of gamma-ray burst X-ray afterglow emission

GRB Jet Beaming Angle Statistics

Optical and gamma-ray emissions from internal forward-reverse shocks: application to GRB 080319B?

Self-similar structure of magnetized ADAFs and CDAFs

X-Ray and high energy flares from late internal shocks of gamma-ray bursts

Behavior of X-Ray Dust Scattering and Implications for X-Ray Afterglows of Gamma-Ray Bursts

Echo Emission From Dust Scattering and X-Ray Afterglows of Gamma-Ray Bursts